[Mesa-users] Entropy Profiles and relaxation

Bill Paxton paxton at kitp.ucsb.edu
Sat Nov 28 14:07:24 EST 2020


> On Nov 28, 2020, at 11:00 AM, Francis Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
>> Also is there a way to get multiple entropy profiles on one plot with pgplot?
> 

I sometimes like to see 2 things together in a single plot such as gradT and grada, or avg_charge_He and avg_charge_H.

Profile_Panels are my way to do this.    Here’s an example.

Profile_Panels1_num_panels = 5

Profile_Panels1_yaxis_name(1) = 'logRho'
Profile_Panels1_other_yaxis_name(1) = 'logT'
Profile_Panels1_other_ymin(1) = -101d0
Profile_Panels1_other_ymax(1) = -101d0

Profile_Panels1_yaxis_name(2) = 'avg_charge_He'
Profile_Panels1_ymin(2) = -0.1
Profile_Panels1_ymax(2) = 2.2
Profile_Panels1_other_yaxis_name(2) = 'avg_charge_H'
Profile_Panels1_other_ymin(2) = -0.1
Profile_Panels1_other_ymax(2) = 1.1

Profile_Panels1_yaxis_name(3) = 'opacity'
Profile_Panels1_other_yaxis_name(3) = 'gradr'
Profile_Panels1_other_ymin(3) = -101d0
Profile_Panels1_other_ymax(3) = -101d0

Profile_Panels1_yaxis_name(4) = 'gradT'
Profile_Panels1_ymin(4) = -0.07
Profile_Panels1_ymax(4) = 0.6
Profile_Panels1_other_yaxis_name(4) = 'grada'
Profile_Panels1_other_ymin(4) = -0.07
Profile_Panels1_other_ymax(4) = 0.6

Profile_Panels1_yaxis_name(5) = 'entropy'
Profile_Panels1_ymin(5) = -101d0
Profile_Panels1_ymax(5) = -101d0
Profile_Panels1_other_yaxis_name(5) = 'prad_div_pgas'
Profile_Panels1_other_ymin(5) = -101d0
Profile_Panels1_other_ymax(5) = -101d0

! x-axis limits and properties
Profile_Panels1_xaxis_name = 'logtau'
Profile_Panels1_xmin = -101d0
Profile_Panels1_xmax = 8
Profile_Panels1_xaxis_reversed = .true.

!Profile_Panels1_xaxis_name = 'zone'
!Profile_Panels1_xmin = 15
!Profile_Panels1_xmax = 270
!Profile_Panels1_xaxis_reversed = .true.

!Profile_Panels1_show_grid = .true.
Profile_Panels1_show_mix_regions_on_xaxis = .true.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20201128/5cb8824a/attachment.htm>


More information about the Mesa-users mailing list