[mesa-users] GYRE output fails?

[Pieter Degroote]

Hi all,

I tried writting GYRE output with MESA version 5232 but it fails for 
some reason. The easiest to reproduce my problem is by adding the 
following two lines to a test_suite example inlist (e.g. 7M_prems_to_AGB):

write_pulse_info_with_profile = .true.
pulse_info_format = 'GYRE'

If I however do

write_pulse_info_with_profile = .true.
pulse_info_format = 'FGONG'

Everything works fine (I attach the two versions of the inlist). I 
recall that this worked fine in previous versions of MESA. Am I 
forgetting something obvious here? The error that I get is:

  failed to open LOGS/profile3.data.GYRE
  failed in do_write_gyre
  failed to open LOGS/profile3.data.GYRE
                                         s% model_number         100
                                 prev_num_iounits_in_use           3
                              current_num_iounits_in_use           4
                                                       i           1
STOP finish_step
DATE: 2013-07-25
TIME: 09:34:29


It doesn't fail when writing profile1 and profile2, though those files 
are empty.

Thanks and cheers,
Pieter
-------------- next part --------------
! inlist_7M_prems_to_AGB

&star_job
      
      create_pre_main_sequence_model = .true.

      change_net = .true.      
      new_net_name = 'o18_and_ne22.net'


/ ! end of star_job namelist

&controls

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 50
         max_number_retries = 100
        
      write_pulse_info_with_profile = .true.
      pulse_info_format = 'GYRE'


      initial_mass = 7.00
      initial_z = 0.001
      
      mixing_length_alpha = 1.73

      log_L_upper_limit = 4.3

      max_model_number = 4000
      history_interval = 1
      terminal_cnt = 10
      write_header_frequency = 10

      overshoot_f_above_nonburn = 0.014
      overshoot_f_below_nonburn = 0.014
      overshoot_f_above_burn_h = 0.014
      overshoot_f_below_burn_h = 0.014
      overshoot_f_above_burn_he = 0.014
      overshoot_f_below_burn_he = 0.014
      overshoot_f_above_burn_z = 0.014
      overshoot_f_below_burn_z = 0.014

      which_atm_option = 'simple_photosphere'

! mass loss       
      RGB_wind_scheme = 'Reimers'
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.0003d0  

/ ! end of controls namelist

-------------- next part --------------
! inlist_7M_prems_to_AGB

&star_job
      
      create_pre_main_sequence_model = .true.

      change_net = .true.      
      new_net_name = 'o18_and_ne22.net'


/ ! end of star_job namelist

&controls

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 50
         max_number_retries = 100
        
      write_pulse_info_with_profile = .true.
      pulse_info_format = 'FGONG'


      initial_mass = 7.00
      initial_z = 0.001
      
      mixing_length_alpha = 1.73

      log_L_upper_limit = 4.3

      max_model_number = 4000
      history_interval = 1
      terminal_cnt = 10
      write_header_frequency = 10

      overshoot_f_above_nonburn = 0.014
      overshoot_f_below_nonburn = 0.014
      overshoot_f_above_burn_h = 0.014
      overshoot_f_below_burn_h = 0.014
      overshoot_f_above_burn_he = 0.014
      overshoot_f_below_burn_he = 0.014
      overshoot_f_above_burn_z = 0.014
      overshoot_f_below_burn_z = 0.014

      which_atm_option = 'simple_photosphere'

! mass loss       
      RGB_wind_scheme = 'Reimers'
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.0003d0  

/ ! end of controls namelist