[mesa-users] About eos and kap module

[无心]

Hi Josiah,
  
 Thank you for your tip and now i basically can modify those sample programs to get my result. But i still have one question: 
  
 in sample_kap.f, calling kap_get_Type1 function, there is a parameter zbar, you give the init value zbar=1.5 and comment says : needed for electron conduction at high rho. I do not know whether this parameter should be modified in my program or it fixed to 1.5 is ok? If not ok, what range for this parameter and how should i estimate the init value?
  
 and similar as above, in sample_eos.f eosDT_get function there are zbar and abar parameters. In this program they are get from Init_Compositon function, i think in most cases i need not modify this function to get other value in addition to T,rho,X,Z, is that right?
  
 Thank you for your help in advance
  
 Chen DanDan
  

 

 ------------------ Original ------------------

Hi DanDan,

The answers to both these questions are contained in the sample_eos.f
code.  You need to read the code and try and understand it.  There's no
way around that.  It is only 233 lines, so that's doable a few hours.

Here are some questions to think about while you read it.

> Could you tell me which function can i call in order to get gamma1 and
> gamm3 parameter from T and Rho input?

What does the comment on line 84 mean?

  ! the indices for the results are defined in eos_def.f

How is the code retrieving the results to display?  What sorts of
things are in eos/public/eos_def.f?

> And i still find there are eosDE,eosDT,eosPT,eosQT, could you tell me
> the difference of them?

There is a call to eos_DT on line 85.  What does the preceding comment
say?  There is a call to eos_PT on line 109.  How does this look
different than the call to eos_DT?

Best,
Josiah

P.S. Line numbers are from r7624.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150903/27e5ceb3/attachment.html>