model/inc/EOS.h

a37a13034c Mart*0001 CBOP
                0002 C    !ROUTINE: EOS.h
                0003 C    !INTERFACE:
                0004 C    include EOS.h
                0005 C    !DESCRIPTION: \bv
                0006 C     *==========================================================*
0bfef71ac3 Jean*0007 C     | EOS.h
a37a13034c Mart*0008 C     | o Header file defining coefficients for equation of state.
                0009 C     *==========================================================*
0bfef71ac3 Jean*0010 C     | The values from the model standard input file are
                0011 C     | stored into the variables held here.
a37a13034c Mart*0012 C     *==========================================================*
                0013 C     \ev
                0014 CEOP
                0015
0bfef71ac3 Jean*0016 C     SItoBar  :: conversion factor for pressure, from Pa (SI Unit) to Bar
                0017 C     SItodBar :: conversion factor for pressure, from Pa (SI Unit) to deci Bar
31a5949e20 Mart*0018       _RL SItoBar, SItodBar
8fb3f16bf6 Patr*0019       PARAMETER ( SItoBar  = 1.D-05 )
                0020       PARAMETER ( SItodBar = 1.D-04 )
a37a13034c Mart*0021
0bfef71ac3 Jean*0022 C Shared EOS Parameter
                0023 C     eosRefP0  :: reference atmospheric pressure used in EOS formulation
                0024       COMMON /PARM_EOS_SHARED/ eosRefP0, equationOfState
                0025       _RL eosRefP0
                0026       CHARACTER*(6) equationOfState
                0027
a37a13034c Mart*0028 C Linear equation of state
                0029 C     tAlpha    :: Linear EOS thermal expansion coefficient ( 1/degree ).
                0030 C     sBeta     :: Linear EOS haline contraction coefficient.
0bfef71ac3 Jean*0031       COMMON /PARM_EOS_LIN/ tAlpha, sBeta
a37a13034c Mart*0032       _RL tAlpha
                0033       _RL sBeta
                0034
                0035 C Equation of State (polynomial coeffients)
0bfef71ac3 Jean*0036       COMMON /PARM_EOS_POLY3/ eosC,eosSig0,eosRefT,eosRefS
a37a13034c Mart*0037       _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)
                0038
                0039 C Full Equation of State
0bfef71ac3 Jean*0040 C     eosType = 'JMD95' (Jackett and McDougall 1995, JAOT)
a37a13034c Mart*0041 C     eosType = 'UNESCO' (Millero et al. 1980, DSR)
0bfef71ac3 Jean*0042 C     COMMON /PARM_EOS_JMD95/
a37a13034c Mart*0043 C     eosJMDCFw  :: of fresh water at pressure 0
                0044 C     eosJMDCSw  :: of sea water at pressure 0
                0045 C     eosJMDCKFw :: of secant bulk modulus K of fresh water at pressure 0
                0046 C     eosJMDCKSw :: of secant bulk modulus K of sea water at pressure 0
                0047 C     eosJMDCKP  :: of secant bulk modulus K at pressure p
0bfef71ac3 Jean*0048 C     eosType = 'MDJWF' (McDougall et al. 2003, JAOT)
                0049 C     COMMON /PARM_EOS_MDJWF/
31a5949e20 Mart*0050 C     eosMDJWFnum :: coefficients of numerator
                0051 C     eosMDJWFden :: coefficients of denominator
ac1c5b24e7 Mart*0052 C     eosType = 'TEOS10' (McDougall et al. 2011, http://www.teos-10.org)
                0053 C     Note: this eos implies that variables THETA and SALT are interpreted
                0054 C     as conservative temperature and absolute salinity
                0055 C     COMMON /PARM_TEOS10/
                0056 C     teos        :: 48 coeffiencts of numerator and denominator
a37a13034c Mart*0057 C     end nonlinear equation of state
                0058       _RL eosJMDCFw(6), eosJMDCSw(9)
                0059       _RL eosJMDCKFw(5), eosJMDCKSw(7), eosJMDCKP(14)
                0060       COMMON /PARM_EOS_JMD95/
                0061      &     eosJMDCFw, eosJMDCSw, eosJMDCKFw, eosJMDCKSw, eosJMDCKP
31a5949e20 Mart*0062       _RL eosMDJWFnum(0:11), eosMDJWFden(0:12)
0bfef71ac3 Jean*0063       COMMON /PARM_EOS_MDJWF/
31a5949e20 Mart*0064      &     eosMDJWFnum, eosMDJWFden
ac1c5b24e7 Mart*0065       _RL teos(48)
                0066       COMMON /PARM_TEOS10/
                0067      &     teos