[Mesa-users] Change_E_flag

[Bill Paxton]

It is great that you are going to jump forward to the current public release of mesa (10398) to take advantage of the improvements that have been made in the treatment of ccsn’s.  The test_suite cases “example_make_pre_ccsn” and “example_ccsn_IIp” are up-to-date and should give you a good start for migrating your work from 7624.   You might do better if you think in terms of modifying the inlists in the current example_ccsn_IIp test case to do what you want rather than trying to modify your old inlists to work in the new version.  That might seem like more effort, but I think it is the best way to catch up with what has happened in mesa in the last 3 years (7624 was released in June 2015). 

good luck
bill







> 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 <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 <mailto: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 <mailto: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 <mailto: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.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180625/94a9d372/attachment.html>