[mesa-users] Metal mass fractions from chem_def

[Bill Paxton]

Hi Josiah,

Perhaps we need to publicize the existence of mesa/kap/test/src/sample_kap.f and similar files.

It is intended to provide a simple to use temple with all of the issues of what to initialize and in what order illustrated explicitly.

For example, in includes the following lines that do the job including the call on crlibm:

      ! initialization and setup

      my_mesa_dir = '../..'         
      call const_init(my_mesa_dir,ierr)     
   	if (ierr /= 0) then
   	   write(*,*) 'const_init failed'
   	   stop 1
   	end if        
         
      call crlibm_init

      call chem_init('isotopes.data', ierr)
      call kap_init( &
         kappa_file_prefix, kappa_CO_prefix, kappa_lowT_prefix, &
         kappa_blend_logT_upper_bdy, kappa_blend_logT_lower_bdy, &
         type2_logT_lower_bdy, use_cache, '', ierr) 
      if(ierr/=0) stop 'problem in kap_init'


-Bill





On May 20, 2015, at 5:32 PM, Josiah Schwab wrote:

> Hi Ehsan,
> 
>> Perhaps, crlibm_init() can go into chem_init() to ensure correct use
>> of exponentials, otherwise, a user — who is not very updated with all
>> design changes and updates — may miss this tiny little detail, and do
>> computations wrong. Consider this, please, as a suggestion only.
> 
> You should call crlibm_init before using other MESA modules too, so I
> don't think it makes sense to directly include it in chem_init.
> 
> I do agree it would be nice if the code didn't let you do what you did
> without complaining, but I don't know if that is feasible.  At least it
> gives obviously bogus results and not subtly wrong ones :)
> 
> Josiah
> 
> 
> ------------------------------------------------------------------------------
> One dashboard for servers and applications across Physical-Virtual-Cloud 
> Widest out-of-the-box monitoring support with 50+ applications
> Performance metrics, stats and reports that give you Actionable Insights
> Deep dive visibility with transaction tracing using APM Insight.
> http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users