[mesa-users] Interacting with EOS module

[Ed Brown]

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