[mesa-users] show two profiles

[Bill Paxton]

Hi Jing,

Have you tried using these controls for the Profile window?

         ! axis choices -- use names from mesa/data/star_data/profile_columns.list
         ! you can also use element names or reaction category names
         ! here is a partial list of the available choices:
            ! mass, logT, logRho, logP, logR, entropy, csound, log_opacity, 
            ! eta, mu, grada, gradr, tau, chiRho, abar, ye, luminosity, ....
            ! h1, he3, c12, ...
            ! pp, cno, tri_alfa, .... (see rates/public/rates_def.f)
         Profile_xaxis_name = 'mass'
         Profile_xmin = -101 ! only used if > -100
         Profile_xmax = -101 ! only used if > -100
         
         Profile_yaxis_name = 'logT'
         Profile_yaxis_reversed = .false.
         Profile_ymin = -101 ! only used if > -100
         Profile_ymax = -101 ! only used if > -100        
         Profile_dymin = -1 
         
         Profile_other_yaxis_name = '' 
         Profile_other_yaxis_reversed = .false.
         Profile_other_ymin = -101 ! only used if > -100
         Profile_other_ymax = -101 ! only used if > -100
         Profile_other_dymin = -1
         
         ! extra axis choices
            ! in your run_star_extras routine extras_finish_step,
            ! store the values to be plotted for each cell in the model.
            !
            !    	do k=1,s% nz ! up to 9 extra values for plotting
            !  		s% profile_extra(k,1) = value1(k)
            !  		s% profile_extra(k,2) = value2(k)
            !  		s% profile_extra(k,3) = value3(k)
            !  	end do
            !  	s% profile_extra_name(1) = 'label for value1 axis'
            !  	s% profile_extra_name(2) = 'label for value2 axis'
            !  	s% profile_extra_name(3) = 'label for value3 axis'
            !
            ! in your inlist you can set Profile_xaxis_name,
            ! Profile_yaxis_name, or Profile_other_yaxis_name
            ! to any of 'profile_extra1', 'profile_extra2', 'profile_extra3', ...


-Bill





On Aug 20, 2011, at 10:28 AM, Jing Luan wrote:

> Dear Mesa-users,
> 
> There are three profiles showing on the main pgplot window under the T-Rho profile. How to change the y-axis of these profiles please? I currently only know how to change the x-axis in inlist_project;
> 
> For example, I would like to plot two profile-figures showing on the main window, each figure containing two profiles:
> 
> 1,  M(>r) (i.e. mass outside radius r) and Mh(>r) (i.e. hydrogen mass outside r) as functions of logR
> 2, logRho and \epsilon_nuc as functions of logR;
> 
> The special thing is that there is only M(<r) and s% xa(j,k) directly readable from MESA's data pointer s, I can calculate M(>r) and Mh(>r) in run_star_extras.f, but how could I let pgplot read in my own data to make the first profile-figure please?
> 
> Thank you very much :-)
> 
> -- 
> Sincerely
> Jing
> 
> Ph.D candidate at physics.caltech
> email: jingluan at caltech.edu
> address: MC350-17,Caltech,1200 E.California Blvd
>         Pasadena, CA 91125
> 
>