[Mesa-users] Change_E_flag

[Francis Timmes]

> Could you please suggest what would be the solution for the problem: Fortran runtime error: Cannot match namelist object name use_hydrodynamic_surface_bcs ?


http://mesa.sourceforge.net/controls_defaults.html#use_compression_outer_BC

fxt





> On Jun 25, 2018, at 12:05 PM, Subedi, Shiv Kumar <ss383615 at ohio.edu> wrote:
> 
> 
> 
> Hello Bill,
> 
> 
> I'd some homework for the "Change_E_flag" and "use_hydrodynamic_surface_BCs" in two mesa versions (mesa-r7624 and mesa-r10108). Here's is what I found.
> 
>  For mesa-r7624:
> grep -r "use_hydrodynamic_surface_BCs" test_suite/
> test_suite/ccsn/inlist_edep:      !use_hydrodynamic_surface_BCs = .true.
> i.e. "use_hydrodynamic_surface_BCs" is present in the inlist_edep of "ccsn" in test_suite.
> 
> Also,
> 
> grep -r "change_E_flag" test_suite/
> test_suite/envelope_shock/inlist_start_shock:      change_E_flag = .true.
> test_suite/ccsn/inlist_edep:      change_E_flag = .true.
> i.e. "change_E_flag" is present in the inlist_edep of "ccsn" and the initial_start_shock of "envelope_shock" of test_suite.
> 
> For mesa-r10108:
> Both the directories "envelope_shock" and "ccsn" are absent in the test_suite and the above variables are not present in any of other current directories in the test_suite.
> 
> 
> MESA INSTRUMENT PAPER 3
> I've been using the inlists of Mesa instrument paper 3 for my work on ccsn with mesa-r7624. In the "inlist_edep" of the inlists provided, above variables come into play as follows.
> 
> &star_job
> change_E_flag = .true.
>  new_E_flag = .true.
> 
> &controls
> ! equation controls
>       use_ODE_var_eqn_pairing = .true.
>       use_hydrodynamic_surface_BCs = .true.
>       use_dvdt_form_of_momentum_eqn = .true.
>       use_dPrad_dm_form_of_T_gradient_eqn = .true.
> 
>       use_dedt_form_of_energy_eqn = .true.
> 
> 
> MESA INSTRUMENT PAPER 4
> I then looked into the "inlist_edep" of the inlists provided for "ccsn" in Mesa instrument paper 4 (mesa-r10108) for all the cases of Core-collapse SN models.
> There is no "E_flag" in the "star_job" section and the "controls" section has following variables.
> 
> &controls 
>   ! equation controls
>       use_ODE_var_eqn_pairing = .true.      
>       use_dPrad_dm_form_of_T_gradient_eqn = .true.
>       use_dedt_form_of_energy_eqn = .true.
>       
>       !use_momentum_outer_BC = .true.
>       use_fixed_vsurf_outer_BC = .true.
>       fixed_vsurf = 0d0
> 
> i.e. there is no "use_hydrodynamic_surface_BCs" in the equation controls here, and comparing above I cannot say for 100% what variable replaces the " use_hydrodynamic_surface_BCs".
> 
> 
> Also, there is no mention of the above two variables in the controls defaults section in MESA web page 
> http://mesa.sourceforge.net/controls_defaults.html#use_compression_outer_BC.
> 
> 
> 
> Could you please suggest what would be the solution for the problem: Fortran runtime error: Cannot match namelist object name use_hydrodynamic_surface_bcs ?
> Also what is best solution if any of such problem arises for a variable which has no definition in the defaults files while upgrading a new mesa version?
> 
> Best,
> Shiv Subedi
> Graduate Student- Ohio University
> 
> 
> 
> 
> 
> 
> 
> 
> 
> From: Bill Paxton <paxton at kitp.ucsb.edu>
> Sent: Sunday, June 24, 2018 11:18 AM
> To: Subedi, Shiv Kumar
> Cc: Frank Timmes; mesa-users at lists.mesastar.org
> Subject: Re: [Mesa-users] Change_E_flag
>  
> time for some detective work.   look in star/test_suite from the old version to find cases that use the old control name (use_hydrodynamic_surface_BCs) in their inlists.
> 
> then look at the corresponding inlists in the newer version to see what has taken its place.
> 
> let us know what you discover.
> 
> bill
> 
> 
> 
> 
>> On Jun 24, 2018, at 4:13 PM, Subedi, Shiv Kumar via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>> 
>> !### use_hydrodynamic_surface_BCs
>>          
>>          ! 1) gradient of compression vanishes at surface
>>             ! see Grott, Chernigovski, Glatzel, 2005.
>>             ! d_dm(d_dm(r^2*v)) = 0 at surface
>>             ! by continuity, this is d_dm(d_dt(1/rho)) = 0 at surface
>>             ! finite volume form is
>>             ! (1/rho(1) - 1/rho_start(1)) = (1/rho(2) - 1/rho_start(2))
>>             ! this BC determines the density for surface cell.
>>          
>>          ! 2) black body radiation at surface
>>             ! L = 4*pi*R^2*boltz_sig*T^4
>>             ! this BC determines the temperature for the surface cell given R and L.
>>                   
>>       use_hydrodynamic_surface_BCs = .false.
>> 
>> However controls.defaults in mesa-r10108 doesnot have any description of the above variable.