[mesa-users] Interacting with EOS module

[Ehsan Moravveji]

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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: code.f
Type: application/octet-stream
Size: 1809 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140612/bbbd37e8/attachment.obj>