Back to home page

MITgcm

 
 

    


File indexing completed on 2025-02-02 06:10:33 UTC

view on githubraw file Latest commit 701e10a9 on 2025-02-01 19:15:20 UTC
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
701e10a905 Mart*0065 
                0066 C     TEOS10 coefficients
ac1c5b24e7 Mart*0067       _RL teos(48)
701e10a905 Mart*0068 
                0069 C     Parameters in the temperature conversion code for TEOS10
                0070 C     The TEOS 10 conversion factor to go from reference salinity to
                0071 C     practical salinity (nondim)
                0072       _RL Sprac_Sref
                0073 C     The inverse of a plausible range of oceanic salinities (kg g-1)
                0074       _RL I_S0
                0075 C     The inverse of a plausible range of oceanic temperatures (degC-1)
                0076       _RL I_Ts
                0077 C     The inverse of the "specific heat" for use
                0078 C     with Conservative Temperature, as defined with TEOS10 (degC kg J-1)
                0079       _RL I_cp0
                0080 
                0081 C     The following are coefficients of contributions to conservative
                0082 C     temperature as a function of the square root of normalized
                0083 C     absolute salinity with an offset (zS) and potential temperature
                0084 C     (T) with a contribution Hab * zS**a * T**b.  The numbers here are
                0085 C     copied directly from the corresponding gsw module, but the
                0086 C     expressions here do not use the same nondimensionalization for
                0087 C     pressure or temperature as they do.
                0088 
                0089 C     Tp to Tc fit constant (degC)
                0090       _RL H00
                0091 C     Tp to Tc fit T coef. (nondim)
                0092       _RL H01
                0093 C     Tp to Tc fit T**2 coef. (degC-1)
                0094       _RL H02
                0095 C     Tp to Tc fit T**3 coef. (degC-2)
                0096       _RL H03
                0097 C     Tp to Tc fit T**4 coef. (degC-3)
                0098       _RL H04
                0099 C     Tp to Tc fit T**5 coef. (degC-4)
                0100       _RL H05
                0101 C     Tp to Tc fit T**6 coef. (degC-5)
                0102       _RL H06
                0103 C     Tp to Tc fit T**7 coef. (degC-6)
                0104       _RL H07
                0105 C     Tp to Tc fit zS**2 coef. (degC)
                0106       _RL H20
                0107 C     Tp to Tc fit zS**2 * T coef. (nondim)
                0108       _RL H21
                0109 C     Tp to Tc fit zS**2 * T**2 coef. (degC-1)
                0110       _RL H22
                0111 C     Tp to Tc fit zS**2 * T**3 coef. (degC-2)
                0112       _RL H23
                0113 C     Tp to Tc fit zS**2 * T**4 coef. (degC-3)
                0114       _RL H24
                0115 C     Tp to Tc fit zS**2 * T**5 coef. (degC-4)
                0116       _RL H25
                0117 C     Tp to Tc fit zS**2 * T**6 coef. (degC-5)
                0118       _RL H26
                0119 C     Tp to Tc fit zS**3 coef. (degC)
                0120       _RL H30
                0121 C     Tp to Tc fit zS** 3* T coef. (nondim)
                0122       _RL H31
                0123 C     Tp to Tc fit zS**3 * T**2 coef. (degC-1)
                0124       _RL H32
                0125 C     Tp to Tc fit zS**3 * T**3 coef. (degC-2)
                0126       _RL H33
                0127 C     Tp to Tc fit zS**3 * T**4 coef. (degC-3)
                0128       _RL H34
                0129 C     Tp to Tc fit zS**4 coef. (degC)
                0130       _RL H40
                0131 C     Tp to Tc fit zS**4 * T coef. (nondim)
                0132       _RL H41
                0133 C     Tp to Tc fit zS**4 * T**2 coef. (degC-1)
                0134       _RL H42
                0135 C     Tp to Tc fit zS**4 * T**3 coef. (degC-2)
                0136       _RL H43
                0137 C     Tp to Tc fit zS**4 * T**4 coef. (degC-3)
                0138       _RL H44
                0139 C     Tp to Tc fit zS**4 * T**5 coef. (degC-4)
                0140       _RL H45
                0141 C     Tp to Tc fit zS**5 coef. (degC)
                0142       _RL H50
                0143 C     Tp to Tc fit zS**6 coef. (degC)
                0144       _RL H60
                0145 C     Tp to Tc fit zS**7 coef. (degC)
                0146       _RL H70
                0147 
                0148 C     The following are coefficients in the nominator (TPNxx) or
                0149 C     denominator (TPDxx) of a simple rational expression that
                0150 C     approximately converts conservative temperature to potential
                0151 C     temperature.
                0152 C     Simple fit numerator constant (degC)
                0153       _RL TPN00
                0154 C     Simple fit numerator Sa coef. (degC ppt-1)
                0155       _RL TPN10
                0156 C     Simple fit numerator Sa**2 coef. (degC ppt-2)
                0157       _RL TPN20
                0158 C     Simple fit numerator Tc coef. (nondim)
                0159       _RL TPN01
                0160 C     Simple fit numerator Sa * Tc coef. (ppt-1)
                0161       _RL TPN11
                0162 C     Simple fit numerator Tc**2 coef. (degC-1)
                0163       _RL TPN02
                0164 C     Simple fit denominator Sa coef. (ppt-1)
                0165       _RL TPD10
                0166 C     Simple fit denominator Tc coef. (degC-1)
                0167       _RL TPD01
                0168 C     Simple fit denominator Tc**2 coef. (degC-2)
                0169       _RL TPD02
                0170 
ac1c5b24e7 Mart*0171       COMMON /PARM_TEOS10/
701e10a905 Mart*0172      &     teos,
                0173      &     Sprac_Sref, I_S0, I_Ts, I_cp0,
                0174      &     H00, H01, H02, H03, H04, H05, H06, H07,
                0175      &     H20, H21, H22, H23, H24, H25, H26,
                0176      &     H30, H31, H32, H33, H34,
                0177      &     H40, H41, H42, H43, H44, H45,
                0178      &     H50, H60, H70,
                0179      &     TPN00, TPN10, TPN20,
                0180      &     TPN01, TPN11, TPN02, TPD10, TPD01, TPD02