1 |
mlosch |
1.2 |
C $Header: /u/gcmpack/MITgcm/model/inc/EOS.h,v 1.1 2002/08/07 16:55:52 mlosch Exp $ |
2 |
|
|
C $Name: $ |
3 |
mlosch |
1.1 |
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 |
mlosch |
1.2 |
_RL SItoBar, SItodBar |
20 |
|
|
PARAMETER ( SItoBar = 1. _d -05 ) |
21 |
|
|
PARAMETER ( SItodBar = 1. _d -04 ) |
22 |
mlosch |
1.1 |
|
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 |
mlosch |
1.2 |
C eosType = 'DJWF02' (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 |
mlosch |
1.1 |
C end nonlinear equation of state |
49 |
|
|
_RL eosJMDCFw(6), eosJMDCSw(9) |
50 |
|
|
_RL eosJMDCKFw(5), eosJMDCKSw(7), eosJMDCKP(14) |
51 |
|
|
COMMON /PARM_EOS_JMD95/ |
52 |
|
|
& eosJMDCFw, eosJMDCSw, eosJMDCKFw, eosJMDCKSw, eosJMDCKP |
53 |
mlosch |
1.2 |
_RL eosMDJWFnum(0:11), eosMDJWFden(0:12) |
54 |
|
|
COMMON /PARM_EOS_MDJWF/ |
55 |
|
|
& eosMDJWFnum, eosMDJWFden |
56 |
mlosch |
1.1 |
|
57 |
|
|
C pressure :: global absolute pressure variable needed for the |
58 |
|
|
C nonlinear equation of state |
59 |
|
|
_RL pressure(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
60 |
|
|
COMMON /EOS_PRESSURE/ pressure |