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 UTCa37a13034c 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
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine. The LXR team |
|