1 |
C $Header: /u/gcmpack/MITgcm/model/inc/EOS.h,v 1.4 2003/07/18 20:06:45 heimbach Exp $ |
2 |
C $Name: $ |
3 |
CBOP |
4 |
C !ROUTINE: EOS.h |
5 |
C !INTERFACE: |
6 |
C include EOS.h |
7 |
C !DESCRIPTION: \bv |
8 |
C *==========================================================* |
9 |
C | EOS.h |
10 |
C | o Header file defining coefficients for equation of state. |
11 |
C *==========================================================* |
12 |
C | The values from the model standard input file are |
13 |
C | stored into the variables held here. |
14 |
C *==========================================================* |
15 |
C \ev |
16 |
CEOP |
17 |
|
18 |
C PARAMETER SItoBar |
19 |
_RL SItoBar, SItodBar |
20 |
PARAMETER ( SItoBar = 1.D-05 ) |
21 |
PARAMETER ( SItodBar = 1.D-04 ) |
22 |
|
23 |
C Linear equation of state |
24 |
C tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
25 |
C sBeta :: Linear EOS haline contraction coefficient. |
26 |
COMMON /PARM_EOS_LIN/ tAlpha,sBeta,equationOfState |
27 |
_RL tAlpha |
28 |
_RL sBeta |
29 |
character*(6) equationOfState |
30 |
|
31 |
C Equation of State (polynomial coeffients) |
32 |
COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
33 |
_RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1) |
34 |
|
35 |
C Full Equation of State |
36 |
C eosType = 'JMD95' (Jackett and McDougall 1995, JPO) |
37 |
C eosType = 'UNESCO' (Millero et al. 1980, DSR) |
38 |
C COMMON /PARM_EOS_JMD95/ |
39 |
C eosJMDCFw :: of fresh water at pressure 0 |
40 |
C eosJMDCSw :: of sea water at pressure 0 |
41 |
C eosJMDCKFw :: of secant bulk modulus K of fresh water at pressure 0 |
42 |
C eosJMDCKSw :: of secant bulk modulus K of sea water at pressure 0 |
43 |
C eosJMDCKP :: of secant bulk modulus K at pressure p |
44 |
C eosType = 'MDJWF' (McDougall et al. 2002, JAOT, submitted) |
45 |
C COMMON /PARM_EOS_MDJWF/ |
46 |
C eosMDJWFnum :: coefficients of numerator |
47 |
C eosMDJWFden :: coefficients of denominator |
48 |
C eosType = 'TEOS10' (McDougall et al. 2011, http://www.teos-10.org) |
49 |
C Note: this eos implies that variables THETA and SALT are interpreted |
50 |
C as conservative temperature and absolute salinity |
51 |
C COMMON /PARM_TEOS10/ |
52 |
C teos :: 48 coeffiencts of numerator and denominator |
53 |
C end nonlinear equation of state |
54 |
_RL eosJMDCFw(6), eosJMDCSw(9) |
55 |
_RL eosJMDCKFw(5), eosJMDCKSw(7), eosJMDCKP(14) |
56 |
COMMON /PARM_EOS_JMD95/ |
57 |
& eosJMDCFw, eosJMDCSw, eosJMDCKFw, eosJMDCKSw, eosJMDCKP |
58 |
_RL eosMDJWFnum(0:11), eosMDJWFden(0:12) |
59 |
COMMON /PARM_EOS_MDJWF/ |
60 |
& eosMDJWFnum, eosMDJWFden |
61 |
_RL teos(48) |
62 |
COMMON /PARM_TEOS10/ |
63 |
& teos |
64 |
|
65 |
C pressure :: global absolute pressure variable needed for the |
66 |
C nonlinear equation of state |
67 |
c _RL pressure(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
68 |
c COMMON /EOS_PRESSURE/ pressure |