[mesa-users] Interacting with EOS module

Ehsan Moravveji e.moravveji at gmail.com
Thu Jun 12 10:24:04 EDT 2014


Thanks Ed for your reply. Now it's solved.
However, I came up with another question:

In /eos/private/eosPT_eval.f, 
            d_dabar_c_TRho = 0d0
            d_dzbar_c_TRho = 0d0
for the case where a control parameter, alfa, is zero. I remember from some old emails of Bill that alfa and beta control merging different EOS tables.
Now, I would like to know why these two differentials are hardcoded to zero, and how to weave the case to where they are explicitly evaluated?
Any useful &star_job variable to change?

Cheers
Ehsan.



On Jun 12, 2014, at 3:13 PM, Ed Brown wrote:

> Hi Ehsan,
> The indices for particular entries in d_dabar_const_TRho(:) are defined in lines 41-72 of eos_def.f, I believe.  As for errors with the code, I notice you define res(:), etc. as a pointer, but in eos_lib.f they are defined as assumed shape arrays—perhaps that is a problem?
> 
> Best,
> Ed
> 
> 
> 
> On Jun 12, 2014, at 14:56 , Ehsan Moravveji <e.moravveji at gmail.com> wrote:
> 
>> Hi all,
>> I have a difficulty and a question regarding the public EOS.
>> 
>> - First the question:
>> I would like to calculate \partial P / \partial \bar{A} at constant (T,Rho). I see from mesa/eos/public/eos_lib.f that eosPT_get returns e.g. d_dabar_const_TRho(:). How can I use it to get my desired derivative (assuming the difficulty below is resolved)?
>> 
>> - And the difficulty:
>> In my inlist, I set 
>>    change_lnPgas_flag = .true.
>>    new_lnPgas_flag = .true.
>> and also put a short subroutine together to interact with the EOS and query for eosPT_get(). However, I receive a segmentation fault during the runtime. 
>> I appreciate if someone helps me finding the problem. I have tried hard to preserve argument shapes, but something is missing for sure.
>> I have attached my short subroutine.
>> 
>> Best regards
>> Ehsan.
>> <code.f>------------------------------------------------------------------------------
>> HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions
>> Find What Matters Most in Your Big Data with HPCC Systems
>> Open Source. Fast. Scalable. Simple. Ideal for Dirty Data.
>> Leverages Graph Analysis for Fast Processing & Easy Data Exploration
>> http://p.sf.net/sfu/hpccsystems_______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
> 





More information about the Mesa-users mailing list