Matlab software for constructing 1D interior structure models based on planetary properties. Self-consistent thermodynamics are used for fluid, rock, and mineral phases. Sound speeds, attenuation, and electrical conductivities are computed as outputs.
The main code is called from an input file containing all the planetary data. Ideally, no tweaks to the main code are needed in order to change the outputs of the model.
Some calculations use Matlab's Parallel Computing package. If you don't have access to this package then parfor loops should be changed to for loops. A future version will check and do this automatically.
Calculations with seawater solutions use the Gibbs Seawater package for Matlab: http://www.teos-10.org/software.htm#1
Calculations with NH3 solutions use REFPROP and require a compiled dynamic library based on the REFPROP source code (see below for placement of files). The source can be obtained from the National Institute of Standards and Technology https://www.nist.gov/refprop Access to REFPROP functions is through python 3 using librefprop.so: https://github.com/jowr/librefprop.so The python capabilities of are employed using the included matlab code refproppy.m REFPROP version 10, expected in October 2017, will provide matlab functions and Mac modules, which may eliminate the need for the above workarounds.
Rock properties are from Perple_X: http://www.perplex.ethz.ch/ Input files were developed by Fabio Cammarano. Version 6.7.8 is currently being used.
TODOs: Modularization is not complete. Further equations of state are under development Update to work with REFPROP V10
The source files and library for REFPROP should be placed in the top-level directory under /opt, as per below $ls /opt librefprop.dylib refprop
$ls /opt/refprop fluids mixtures
$ls /opt/refprop/fluids 1BUTENE.FLD C1CC6.FLD COS.FLD D6.FLD ETHYLENE.FLD HYDROGEN.FLD MD3M.FLD MOLEATE.FLD NITROGEN.FLD PENTANE.FLD R115.FLD R124.FLD R152A.FLD R236FA.FLD RE143A.FLD WATER.FLD ACETONE.FLD C2BUTENE.FLD CYCLOHEX.FLD DECANE.FLD FLUORINE.FLD IBUTENE.FLD MD4M.FLD MPALMITA.FLD NONANE.FLD PROPANE.FLD R116.FLD R125.FLD R161.FLD R245CA.FLD RE245CB2.FLD XENON.FLD AMMONIA.FLD C3CC6.FLD CYCLOPEN.FLD DEE.FLD H2S.FLD IHEXANE.FLD MDM.FLD MSTEARAT.FLD NOVEC649.FLD PROPYLEN.FLD R12.FLD R13.FLD R21.FLD R245FA.FLD RE245FA2.FLD ARGON.FLD C4F10.FLD CYCLOPRO.FLD DMC.FLD HCL.FLD IOCTANE.FLD METHANE.FLD MXYLENE.FLD OCTANE.FLD PROPYNE.FLD R1216.FLD R134A.FLD R218.FLD R32.FLD RE347MCC.FLD BENZENE.FLD C5F12.FLD D2.FLD DME.FLD HELIUM.FLD IPENTANE.FLD METHANOL.FLD N2O.FLD ORTHOHYD.FLD PXYLENE.FLD R123.FLD R14.FLD R22.FLD R365MFC.FLD SF6.FLD BUTANE.FLD CF3I.FLD D2O.FLD EBENZENE.FLD HEPTANE.FLD ISOBUTAN.FLD MLINOLEA.FLD NEON.FLD OXYGEN.FLD R11.FLD R1233ZD.FLD R141B.FLD R227EA.FLD R40.FLD SO2.FLD C11.FLD CO.FLD D4.FLD ETHANE.FLD HEXANE.FLD KRYPTON.FLD MLINOLEN.FLD NEOPENTN.FLD OXYLENE.FLD R113.FLD R1234YF.FLD R142B.FLD R23.FLD R41.FLD T2BUTENE.FLD C12.FLD CO2.FLD D5.FLD ETHANOL.FLD HMX.BNC MD2M.FLD MM.FLD NF3.FLD PARAHYD.FLD R114.FLD R1234ZE.FLD R143A.FLD R236EA.FLD RC318.FLD TOLUENE.FLD
$ ls /opt/refprop/mixtures/ AIR.MIX HIGHN2.MIX R402A.MIX R405A.MIX R407D.MIX R409B.MIX R412A.MIX R415B.MIX R420A.MIX R422C.MIX R426A.MIX R431A.mix R436A.MIX R442A.MIX R502.MIX R508B.MIX AMARILLO.MIX NGSAMPLE.MIX R402B.MIX R406A.MIX R407E.MIX R410A.MIX R413A.MIX R416A.MIX R421A.MIX R422D.MIX R427A.MIX R432A.mix R436B.MIX R443A.MIX R503.MIX R509A.MIX EKOFISK.MIX R401A.MIX R403A.MIX R407A.MIX R407F.MIX R410B.MIX R414A.MIX R417A.MIX R421B.MIX R423A.MIX R428A.MIX R433A.mix R437A.MIX R444A.MIX R504.MIX R510A.MIX GLFCOAST.MIX R401B.MIX R403B.MIX R407B.MIX R408A.MIX R411A.MIX R414B.MIX R418A.MIX R422A.MIX R424A.MIX R429A.mix R434A.mix R438A.mix R500.MIX R507A.MIX R512A.MIX HIGHCO2.MIX R401C.MIX R404A.MIX R407C.MIX R409A.MIX R411B.MIX R415A.MIX R419A.MIX R422B.MIX R425A.MIX R430A.mix R435A.MIX R441A.MIX R501.MIX R508A.MIX