[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