&pgstar

pgstar_interval = 1
pgstar_show_age_in_years = .true.
pgstar_show_age_in_seconds = .false.
pgstar_sleep = 0.0


! white_on_black flags -- true means white foreground color on black background
file_white_on_black_flag = .true.
file_device = 'png'        ! options 'png' and 'vcps' for png and postscript respectively
file_extension = 'png'     ! common names are 'png' and 'ps'


!file_white_on_black_flag = .false.
!file_device = 'vcps'      ! options 'png' and 'vcps' for png and postscript respectively
!file_extension = 'ps'     ! common names are 'png' and 'ps'


!------------------------------------------------------------------------------------

! mass fractions profile plot

Abundance_win_flag = .false.

! window properties
Abundance_win_width = 10
Abundance_win_aspect_ratio = 0.75
         
Abundance_xleft = 0.15
Abundance_xright = 0.85
Abundance_ybot = 0.15
Abundance_ytop = 0.85
Abundance_txt_scale = 0.8
Abundance_title = ''

! isotopes to plot

Abundance_num_isos_to_show = 20

Abundance_which_isos_to_show(1)  = 'h1'
Abundance_which_isos_to_show(2)  = 'h2'
Abundance_which_isos_to_show(3)  = 'he3'
Abundance_which_isos_to_show(4)  = 'he4'
Abundance_which_isos_to_show(5)  = 'li7'
Abundance_which_isos_to_show(6)  = 'be7'
Abundance_which_isos_to_show(7)  = 'be9'
Abundance_which_isos_to_show(8)  = 'be10'
Abundance_which_isos_to_show(9)  = 'b8'
Abundance_which_isos_to_show(10) = 'c12'
Abundance_which_isos_to_show(11) = 'c12'
Abundance_which_isos_to_show(12) = 'c13'
Abundance_which_isos_to_show(13) = 'n13'
Abundance_which_isos_to_show(14) = 'n14'
Abundance_which_isos_to_show(15) = 'n15'
Abundance_which_isos_to_show(16) = 'o14'
Abundance_which_isos_to_show(17) = 'o15'
Abundance_which_isos_to_show(18) = 'o16'
Abundance_which_isos_to_show(19) = 'o17'
Abundance_which_isos_to_show(20) = 'o18'


! number and size of isotope labels along curves
num_abundance_line_labels = 4
Abundance_line_txt_scale_factor = 0.8


! number and size of isotopes on legend         
Abundance_legend_max_cnt = 0
Abundance_legend_txt_scale_factor = 0.6

! xaxis name and orientation
Abundance_xaxis_name = 'mass' 
Abundance_xaxis_reversed = .false.
         
! xaxis limits 
Abundance_xmin = 0.0 
Abundance_xmax = 15.0 
         
! yaxis limits 
Abundance_log_mass_frac_min = -3.0 
Abundance_log_mass_frac_max =  0.3
         
! file output
Abundance_file_flag = .false.
Abundance_file_dir = 'Abundance'
Abundance_file_prefix = 'abund_'
Abundance_file_width = -1        ! (inches) negative means use same value as for window
Abundance_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------


! Network window
! ==============

      ! Shows N vs Z for isotopes in current net

      ! ::

    Network_win_flag = .false.

      ! ::

    Network_win_width = 10
    Network_win_aspect_ratio = 0.75

      ! ::

    Network_xleft = 0.15
    Network_xright = 0.85
    Network_ybot = 0.15
    Network_ytop = 0.85
    Network_txt_scale = 1.0
    Network_title = 'Network'

      ! axis limits -- to override system default selections

      ! ::

    Network_nmin = -101d0
    Network_nmax = -101d0
    Network_zmin = -101d0
    Network_zmax = -101d0

      ! Show mass fraction as a coloured square

      ! ::

    Network_show_mass_fraction = .true.

      ! Limits on mass fractions to show

      ! ::

    Network_log_mass_frac_min = -5.0d0
    Network_log_mass_frac_max = 0.0d0

      ! Label Y axis with element name

      ! ::

    Network_show_element_names = .true.

      ! Enables calling a subroutine to add extra information to a plot
      ! see ``$MESA_DIR/star/other/pgstar_decorator.f90``

      ! ::

    Network_use_decorator = .true.

      ! Add colorbar

      ! ::

    Network_show_colorbar = .false.

      ! file output

      ! ::

    Network_file_flag = .false.
    Network_file_dir = 'png'
    Network_file_prefix = 'Network_'
    Network_file_interval = 5
    Network_file_width = -1
    Network_file_aspect_ratio = -1

!------------------------------------------------------------------------------------


! Production window
! =================

      ! Plots the abundance for each stable isotope divided by the initial stable composition

      ! ::

    Production_win_flag = .false.

      ! ::

    Production_win_width = 10
    Production_win_aspect_ratio = 0.75

      ! ::

    Production_xleft = 0.15
    Production_xright = 0.85
    Production_ybot = 0.15
    Production_ytop = 0.85
    Production_txt_scale = 1.0
    Production_title = 'Production'

      ! axis limits -- to override system default selections

      ! ::

    Production_amin = -101d0
    Production_amax = -101d0
    Production_ymin = -101d0
    Production_ymax = -101d0

      ! Mass range of star to show

      ! ::

    Production_min_mass = -101d0
    Production_max_mass = -101d0

      ! Min log mass fraction an isotope must have to be shown

      ! ::

    Production_min_mass_frac = -5.0d0

      ! Add labels with element names

      ! ::

    Production_show_element_names = .true.

      ! Enables calling a subroutine to add extra information to a plot
      ! see ``$MESA_DIR/star/other/pgstar_decorator.f90``

      ! ::

    Production_use_decorator = .true.

      ! file output

      ! ::

    Production_file_flag = .false.
    Production_file_dir = 'png'
    Production_file_prefix = 'Production_'
    Production_file_interval = 5
    Production_file_width = -1
    Production_file_aspect_ratio = -1

!------------------------------------------------------------------------------------

! power profile plot

Power_win_flag = .false.

Power_win_width = 10
Power_win_aspect_ratio = 0.75
Power_title = ''      

Power_xleft = 0.15
Power_xright = 0.85
Power_ybot = 0.15
Power_ytop = 0.85
Power_txt_scale = 1.0
Power_title = ' '

Power_xaxis_name = 'mass'
Power_xaxis_reversed = .false.

Power_legend_max_cnt = 16
Power_legend_txt_scale_factor = 1.0 ! relative to other text

! power xaxis limits -- to override system default selections
Power_xmin = 0.0 
Power_xmax = 1.0

! power yaxis limits -- to override system default selections
Power_ymin = -5.0 !   -101d0 ! only used if /= -101d0
Power_ymax = 10.0 ! -101d0 ! only used if /= -101d0

! file output
Power_file_flag = .false.
Power_file_dir = 'png'
Power_file_prefix = 'power_'
Power_file_interval = 5 ! output when mod(model_number,Power_file_interval)==0
Power_file_width = -1 ! (inches) negative means use same value as for window
Power_file_aspect_ratio = -1 ! negative means use same value as for window




!------------------------------------------------------------------------------------

! Kippenhahn diagram history plot

Kipp_win_flag = .false.

! window properties
Kipp_win_width = 10
Kipp_win_aspect_ratio = 0.75
Kipp_txt_scale = 0.9
Kipp_title = ''      

! y axis limits
Kipp_mass_max = 15.0
Kipp_mass_min = 0 
Kipp_show_mass_boundaries = .true.

! x axis limits
Kipp_xmax = -101              ! maximum step number.  negative means use default.
Kipp_xmin = 0         ! minimum step number.  negative means use default.
Kipp_max_width = -1         ! only used if > 0.  causes xmin to move with xmax.

Kipp_show_mixing = .true.
Kipp_show_burn = .true.
Kipp_show_luminosities = .true.

! file output
Kipp_file_flag = .false.
Kipp_file_dir = 'Kipp'
Kipp_file_prefix = 'kipp_'
Kipp_file_interval = 5           ! output when mod(model_number,Kipp_file_interval)==0
Kipp_file_width = -1        ! (inches) negative means use same value as for window
Kipp_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------

!  HR diagram history plot

HR_win_flag = .false.

HR_win_width = 10
HR_win_aspect_ratio = 1.0 

HR_xleft = 0.15
HR_xright = 0.85
HR_ybot = 0.15
HR_ytop = 0.85
HR_txt_scale = 0.4
HR_title = ''

! axis limits 
HR_logT_min = 3.3
HR_logT_max = 4.0  
HR_logL_min = 1.0  
HR_logL_max = 4.0  


HR_logL_margin = 0.1
HR_logT_margin = 0.1
HR_dlogT_min = -1
HR_dlogL_min = -1

HR_step_min = 0      ! only plot models with model number >= this
HR_step_max = 999999 ! only plot models with model number <= this

show_HR_classical_instability_strip = .false.
show_HR_Mira_instability_region = .false.
show_HR_WD_instabilities = .false.

show_HR_target_box = .false.
HR_target_n_sigma = -3 ! -n means show sig 1..n
HR_target_logL = 0
HR_target_logL_sigma = 0
HR_target_logT = 0
HR_target_logT_sigma = 0

show_HR_annotation1 = .false.
show_HR_annotation2 = .false.
show_HR_annotation3 = .false.

HR_fname = '' ! file name for extra HR data

! Enables calling a subroutine to add extra information to a plot
! see `$MESA_DIR/star/other/pgstar_decorator.f90`
HR_use_decorator = .false.

! file output
HR_file_flag = .false.
HR_file_dir = 'HR'
HR_file_prefix = 'hr_'
HR_file_interval = 5      ! output when `mod(model_number,HR_file_interval)==0`
HR_file_width = -1        ! (inches) negative means use same value as for window
HR_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------

! central temperature and central density history plot

TRho_win_flag = .false.

TRho_win_width = 10
TRho_win_aspect_ratio = 0.75 
         
TRho_xleft = 0.15
TRho_xright = 0.85
TRho_ybot = 0.15
TRho_ytop = 0.85
TRho_txt_scale = 0.7
TRho_title = ' '      

! axis limits 
TRho_logT_min = 7.0 ! only used if /= -101d0
TRho_logT_max = 8.5 ! only used if /= -101d0
TRho_logRho_min = 1.0 ! only used if /= -101d0
TRho_logRho_max = 7.0 ! only used if /= -101d0


! file output
TRho_file_flag = .false.
TRho_file_dir = 'trho'
TRho_file_prefix = 'trho_'
TRho_file_interval = 5      ! output when `mod(model_number,TRho_file_interval)==0`
TRho_file_width = -1        ! (inches) negative means use same value as for window
TRho_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------

! temperature and density profile plot

TRho_Profile_win_flag = .false.

TRho_Profile_win_width = 10
TRho_Profile_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
         
TRho_Profile_xleft = 0.15
TRho_Profile_xright = 0.85
TRho_Profile_ybot = 0.15
TRho_Profile_ytop = 0.85
TRho_Profile_txt_scale = 0.7
TRho_Profile_title = ' '      
         
TRho_switch_to_Column_Depth = .false. 
TRho_switch_to_mass = .false.
            
show_TRho_Profile_legend = .true.
! TRho_Profile_legend_coord = 0.55
! TRho_Profile_legend_fjust = 0.0
! TRho_Profile_legend_disp1 = -20.5
! TRho_Profile_legend_del_disp = -1.3

 TRho_Profile_legend_coord = 0.07
 TRho_Profile_legend_fjust = 0.0
 TRho_Profile_legend_disp1 = -2.0
 TRho_Profile_legend_del_disp = -1.3


show_TRho_Profile_text_info = .false.
 TRho_Profile_text_info_xfac = 0.77 ! controls x location
 TRho_Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text
 TRho_Profile_text_info_yfac = 0.6 ! controls y location of 1st line
 TRho_Profile_text_info_dyfac = -0.04 ! controls line spacing
         
show_TRho_Profile_mass_locs = .false.
show_TRho_accretion_mesh_borders = .false.
show_TRho_Profile_kap_regions = .false.
show_TRho_Profile_eos_regions = .false.
show_TRho_Profile_degeneracy_line = .true.
show_TRho_Profile_Pgas_Prad_line = .true.
show_TRho_Profile_burn_lines = .true.
show_TRho_Profile_burn_labels = .true.
      
! axis limits
TRho_Profile_xmin = -12.0
TRho_Profile_xmax = 7.0
TRho_Profile_ymin = 3.5
TRho_Profile_ymax = 9.0        
         
! file output
TRho_Profile_file_flag = .false.
TRho_Profile_file_dir = 'TRho'
TRho_Profile_file_prefix = 'trho_'
TRho_Profile_file_width = -1 ! (inches) negative means use same value as for window
TRho_Profile_file_aspect_ratio = -1 ! negative means use same value as for window


! these are shown if show_TRho_Profile_mass_locs = .true.         
! set all the entries
profile_mass_point_q = -1
profile_mass_point_color_index = 1
profile_mass_point_symbol = -6
profile_mass_point_symbol_scale = 1.7
profile_mass_point_str = ''
profile_mass_point_str_clr = 1
profile_mass_point_str_scale = 1.0
         
! set defaults      
num_profile_mass_points = 3 ! max is defined in star_def (max_num_profile_mass_points)

profile_mass_point_q(1) = 0.5
profile_mass_point_color_index(1) = 1
profile_mass_point_symbol(1) = -6
profile_mass_point_str(1) = '  0.5 M\d\(0844)\u'
profile_mass_point_str_clr(1) = 1
         
profile_mass_point_q(2) = 0.95
profile_mass_point_color_index(2) = 1
profile_mass_point_symbol(2) = -6
profile_mass_point_str(2) = '  0.95 M\d\(0844)\u'
profile_mass_point_str_clr(3) = 1
         
profile_mass_point_q(3) = 0.999
profile_mass_point_color_index(3) = 1
profile_mass_point_symbol(3) = -6
profile_mass_point_str(3) = '  0.999 M\d\(0844)\u'
profile_mass_point_str_clr(3) = 1


!------------------------------------------------------------------------------------

! profiles of interesting quantities

Profile_Panels1_win_flag = .false.

Profile_Panels1_num_panels = 3

! window properties
Profile_Panels1_win_width = 10
Profile_Panels1_win_aspect_ratio = 1.0
Profile_Panels1_txt_scale = 0.6
Profile_Panels1_title = '' 

! show mass points and mixing 
Profile_Panels1_show_grid = .true.
Profile_Panels1_show_mix_regions_on_xaxis = .true.

! x-axis limits and properties
!Profile_Panels1_xaxis_name = 'mass'
!Profile_Panels1_xmin = 0.6
!Profile_Panels1_xmax = 1.8
!Profile_Panels1_xaxis_reversed = .false.

Profile_Panels1_xaxis_name = 'mass'
Profile_Panels1_xmin = 0.0
Profile_Panels1_xmax = 15.0 
Profile_Panels1_xaxis_reversed = .false.

! contents of panel 1
Profile_Panels1_yaxis_name(1) = 'logT'   
Profile_Panels1_ymin(1) = -101d0 ! only used if /= -101d0
Profile_Panels1_ymax(1) = -101d0 ! only used if /= -101d0

Profile_Panels1_other_yaxis_name(1) = 'opacity'         
Profile_Panels1_other_ymin(1) = 0.0d0 ! only used if /= -101d0
Profile_Panels1_other_ymax(1) = -101d0 ! only used if /= -101d0


! contents of panel 2
Profile_Panels1_yaxis_name(2) = 'logRho'        
Profile_Panels1_ymin(2) = -101d0 ! only used if /= -101d0
Profile_Panels1_ymax(2) = -101d0 ! only used if /= -101d0

Profile_Panels1_other_yaxis_name(2) = 'luminosity'    
Profile_Panels1_other_ymin(2) = -101d0 ! only used if /= -101d0
Profile_Panels1_other_ymax(2) = -101d0 ! only used if /= -101d0


! contents of panel 3
Profile_Panels1_yaxis_name(3) = 'log_thermal_time_to_surface'             
Profile_Panels1_ymin(3) = -101d0 ! only used if /= -101d0
Profile_Panels1_ymax(3) = -101d0 ! only used if /= -101d0

Profile_Panels1_other_yaxis_name(3) = 'log_t_thermal'
Profile_Panels1_other_ymin(3) = -101d0    ! only used if /= -101d0
Profile_Panels1_other_ymax(3) = -101d0 ! only used if /= -101d0

!Profile_Panels1_other_yaxis_name(3) = 'logtau'
!Profile_Panels1_other_ymin(3) = -101d0    ! only used if /= -101d0
!Profile_Panels1_other_ymax(3) = -101d0 ! only used if /= -101d0



! file output
Profile_Panels1_file_flag = .false.
Profile_Panels1_file_interval = 1     ! output when mod(model_number,Grid1_file_interval)==0
Profile_Panels1_file_dir = 'Thermo'
Profile_Panels1_file_prefix = 'thermo_'
Profile_Panels1_file_width = -1        ! (inches) negative means use same value as for window
Profile_Panels1_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------

! more profiles of interesting quantities

Profile_Panels2_win_flag = .false.

Profile_Panels2_num_panels = 3

Profile_Panels2_show_grid = .true.
Profile_Panels2_show_mix_regions_on_xaxis = .true.

! window properties
Profile_Panels2_win_width = 10
Profile_Panels2_win_aspect_ratio = 1.0
Profile_Panels2_txt_scale = 0.6
Profile_Panels2_title = ''      

! x-axis limits and properties
!Profile_Panels2_xaxis_name = 'mass'
!Profile_Panels2_xmin = 0.6  ! 0.5704600
!Profile_Panels2_xmax = 1.52  ! 0.5705005  ! -101d0
!Profile_Panels2_xaxis_reversed = .false.

Profile_Panels2_xaxis_name = 'logP'
Profile_Panels2_xmin = 1.0
Profile_Panels2_xmax = 9.0 
Profile_Panels2_xaxis_reversed = .true.

! contents of panel 1
Profile_Panels2_yaxis_name(1) = 'entropy'
Profile_Panels2_ymin(1) = 0.0   ! only used if /= -101d0
Profile_Panels2_ymax(1) = 100.0 ! only used if /= -101d0

!Profile_Panels2_other_yaxis_name(1) = ''
Profile_Panels2_other_yaxis_name(1) = 'vel_km_per_s'
Profile_Panels2_other_ymin(1) = -101d0 ! only used if /= -101d0
Profile_Panels2_other_ymax(1) = -101d0 ! only used if /= -101d0

! contents of panel 2
Profile_Panels2_yaxis_name(2) = 'log_j_rot'   
Profile_Panels2_ymin(2) = -101d0 ! only used if /= -101d0
Profile_Panels2_ymax(2) = -101d0 ! only used if /= -101d0

Profile_Panels2_other_yaxis_name(2) = ''   
Profile_Panels2_other_ymin(2) = -101d0 ! only used if /= -101d0
Profile_Panels2_other_ymax(2) = -101d0 ! only used if /= -101d0

! contents of panel 3
Profile_Panels2_yaxis_name(3) = 'gamma1' 
Profile_Panels2_ymin(3) = 1.0 ! only used if /= -101d0
Profile_Panels2_ymax(3) = 1.7 ! only used if /= -101d0

Profile_Panels2_other_yaxis_name(3) = ''
Profile_Panels2_other_ymin(3) = 1.0  ! only used if /= -101d0
Profile_Panels2_other_ymax(3) = 1.7   ! 1.45 ! only used if /= -101d0


! file output
Profile_Panels2_file_flag = .false.
Profile_Panels2_file_interval = 1     ! output when mod(model_number,Grid1_file_interval)==0
Profile_Panels2_file_dir = 'Profile2'
Profile_Panels2_file_prefix = 'ye_'
Profile_Panels2_file_width = -1        ! (inches) negative means use same value as for window
Profile_Panels2_file_aspect_ratio = -1 ! negative means use same value as for window


!------------------------------------------------------------------------------------

! Text_Summary windows

Text_Summary1_win_flag = .false.
Text_Summary1_win_width = 10
Text_Summary1_win_aspect_ratio = 0.15

Text_Summary1_xleft = 0.01
Text_Summary1_xright = 0.99
Text_Summary1_ybot = 0.0
Text_Summary1_ytop = 1.0
Text_Summary1_txt_scale = 0.95
Text_Summary1_title = ''

Text_Summary1_num_rows = 5 ! <= 20
Text_Summary1_num_cols = 6 ! <= 20
Text_Summary1_name(:,:) = ''

Text_Summary1_name(1,1) = 'time_step'
Text_Summary1_name(1,2) = 'num_zones'
Text_Summary1_name(1,3) = 'species'
Text_Summary1_name(1,4) = 'dynamic_timescale'
Text_Summary1_name(1,5) = 'kh_timescale'
Text_Summary1_name(1,6) = 'nuc_timescale'

Text_Summary1_name(2,1) = 'star_mass'
Text_Summary1_name(2,2) = 'star_mdot'
Text_Summary1_name(2,3) = 'center_T'
Text_Summary1_name(2,4) = 'center_Rho'
Text_Summary1_name(2,5) = 'max_T'
Text_Summary1_name(2,6) = 'mass_loc_of_max_eps_nuc'

Text_Summary1_name(3,1) = 'total_angular_momentum'
Text_Summary1_name(3,2) = 'effective_T'
Text_Summary1_name(3,3) = 'photosphere_L'
Text_Summary1_name(3,4) = 'photosphere_r'
Text_Summary1_name(3,5) = 'log_g'
Text_Summary1_name(3,6) = 'center_h1'

Text_Summary1_name(4,1) = 'surf_avg_v_rot'   ! km/sec rotational velocity at equator
Text_Summary1_name(4,2) = 'surface h1'
Text_Summary1_name(4,3) = 'surface he4'
Text_Summary1_name(4,4) = 'surface c12'
Text_Summary1_name(4,5) = 'center_he4' 
Text_Summary1_name(4,6) = 'center_n14' 

Text_Summary1_name(5,1) = 'he_core_mass'
Text_Summary1_name(5,2) = 'co_core_mass'
Text_Summary1_name(5,3) = 'one_core_mass'
Text_Summary1_name(5,4) = 'center_c12' 
Text_Summary1_name(5,5) = 'center_o16'
Text_Summary1_name(5,6) = 'center_ne20'


!------------------------------------------------------------------------------------

! global grid plot settings
pgstar_title_scale = 1.3
pgstar_title_disp = 1.1
pgstar_title_coord = 0.5
pgstar_title_fjust = 0.5

pgstar_age_scale = 0.8
pgstar_age_disp = 3.0
pgstar_age_coord = 0.0
pgstar_age_fjust = 0.0

pgstar_model_scale = 0.8
pgstar_model_disp = 3.0
pgstar_model_coord = 1.0
pgstar_model_fjust = 1.0

pgstar_xaxis_label_disp = 2.5
pgstar_left_yaxis_label_disp = 2.0
pgstar_right_yaxis_label_disp = 2.5



!------------------------------------------------------------------------------------

! grid of all the above plots

Grid1_win_flag = .true.
Grid1_win_width = 16
Grid1_win_aspect_ratio = 0.6

! reset the defaults

Grid1_plot_name(:) = ''
Grid1_plot_row(:) = 1           ! number from 1 at top
Grid1_plot_rowspan(:) = 1       ! plot spans this number of rows
Grid1_plot_col(:) =  1          ! number from 1 at left
Grid1_plot_colspan(:) = 1       ! plot spans this number of columns
Grid1_plot_pad_left(:) = 0.0    ! fraction of full window width for padding on left
Grid1_plot_pad_right(:) = 0.0   ! fraction of full window width for padding on right
Grid1_plot_pad_top(:) = 0.0     ! fraction of full window height for padding at top
Grid1_plot_pad_bot(:) = 0.0     ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(:) = 0.7 ! multiply txt_scale for subplot by this


Grid1_title = ''

Grid1_num_cols = 3  ! divide plotting region into this many equal width cols
Grid1_num_rows = 2  ! divide plotting region into this many equal height rows
Grid1_num_plots = 6 ! must be <= 10


Grid1_plot_name(1) = 'Text_Summary1'
Grid1_plot_row(1) = 1           ! number from 1 at top
Grid1_plot_rowspan(1) = 1       ! plot spans this number of rows
Grid1_plot_col(1) =  1          ! number from 1 at left
Grid1_plot_colspan(1) = 6       ! plot spans this number of columns

Grid1_plot_pad_left(1) = 0.0    ! fraction of full window width for padding on left
Grid1_plot_pad_right(1) = 0.0   ! fraction of full window width for padding on right
Grid1_plot_pad_top(1) = -0.04   ! fraction of full window height for padding at top
Grid1_plot_pad_bot(1) = 0.39    ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(1) = 0.8 ! multiply txt_scale for subplot by this

Grid1_plot_name(2) = 'Abundance' 
Grid1_plot_row(2) = 1           ! number from 1 at top
Grid1_plot_rowspan(2) = 1       ! plot spans this number of rows
Grid1_plot_col(2) =  1          ! number from 1 at left
Grid1_plot_colspan(2) = 1       ! plot spans this number of columns

Grid1_plot_pad_left(2) = -0.085 ! fraction of full window width for padding on left
Grid1_plot_pad_right(2) = 0.06   ! fraction of full window width for padding on right
Grid1_plot_pad_top(2) = 0.05    ! fraction of full window height for padding at top
Grid1_plot_pad_bot(2) = 0.03    ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(2) = 0.7 ! multiply txt_scale for subplot by this


Grid1_plot_name(3) = 'Production' 
Grid1_plot_row(3) = 1           ! number from 1 at top
Grid1_plot_rowspan(3) = 1       ! plot spans this number of rows
Grid1_plot_col(3) =  2          ! number from 1 at left
Grid1_plot_colspan(3) = 1       ! plot spans this number of columns

Grid1_plot_pad_left(3) = 0.01   ! fraction of full window width for padding on left
Grid1_plot_pad_right(3) = 0.0  ! fraction of full window width for padding on right
Grid1_plot_pad_top(3) = 0.05     ! fraction of full window height for padding at top
Grid1_plot_pad_bot(3) = 0.03     ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(3) = 0.5  ! multiply txt_scale for subplot by this


Grid1_plot_name(4) = 'Kipp' 
Grid1_plot_row(4) = 2           ! number from 1 at top
Grid1_plot_rowspan(4) = 1       ! plot spans this number of rows
Grid1_plot_col(4) =  1          ! number from 1 at left
Grid1_plot_colspan(4) = 1       ! plot spans this number of columns

Grid1_plot_pad_left(4) = -0.08   ! fraction of full window width for padding on left
Grid1_plot_pad_right(4) = 0.06   ! fraction of full window width for padding on right
Grid1_plot_pad_top(4) = 0.05     ! fraction of full window height for padding at top
Grid1_plot_pad_bot(4) = 0.01    ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(4) = 0.6  ! multiply txt_scale for subplot by this


Grid1_plot_name(5) = 'Network' 
Grid1_plot_row(5) = 2           ! number from 1 at top
Grid1_plot_rowspan(5) = 1       ! plot spans this number of rows
Grid1_plot_col(5) =  2          ! number from 1 at left
Grid1_plot_colspan(5) = 1       ! plot spans this number of columns

Grid1_plot_pad_left(5) = 0.01   ! fraction of full window width for padding on left
Grid1_plot_pad_right(5) = 0.0   ! fraction of full window width for padding on right
Grid1_plot_pad_top(5) = 0.05     ! fraction of full window height for padding at top
Grid1_plot_pad_bot(5) = 0.01    ! fraction of full window height for padding at bottom
!Grid1_txt_scale_factor(5) = 0.8  ! multiply txt_scale for subplot by this


Grid1_plot_name(6) = 'Profile_Panels1'
Grid1_plot_row(6) = 1           ! number from 1 at top
Grid1_plot_rowspan(6) = 2       ! plot spans this number of rows
Grid1_plot_col(6) =  3          ! number from 1 at left
Grid1_plot_colspan(6) = 1       ! plot spans this number of columns

Grid1_plot_pad_left(6) = 0.05    ! fraction of full window width for padding on left
Grid1_plot_pad_right(6) = 0.0  ! fraction of full window width for padding on right
Grid1_plot_pad_top(6) = 0.05     ! fraction of full window height for padding at top
Grid1_plot_pad_bot(6) = 0.0      ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(6) = 0.7  ! multiply txt_scale for subplot by this



! file output
Grid1_file_flag = .false.
Grid1_file_dir = 'grid_png'
Grid1_file_prefix = 'grid_png'
Grid1_file_interval = 100     ! output when mod(model_number,Grid1_file_interval)==0
Grid1_file_width = -1        ! (inches) negative means use same value as for window
!Grid1_file_width = 21       ! (inches) negative means use same value as for window
Grid1_file_aspect_ratio = -1 ! negative means use same value as for window

!------------------------------------------------------------------------------------


/ ! end of pgstar namelist

