1 |
mlosch |
1.1 |
C $Header: $ |
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 |
20 |
|
|
PARAMETER ( SItoBar = 1. _d -05 ) |
21 |
|
|
|
22 |
|
|
C Linear equation of state |
23 |
|
|
C tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
24 |
|
|
C sBeta :: Linear EOS haline contraction coefficient. |
25 |
|
|
COMMON /PARM_EOS_LIN/ tAlpha,sBeta,equationOfState |
26 |
|
|
_RL tAlpha |
27 |
|
|
_RL sBeta |
28 |
|
|
character*(6) equationOfState |
29 |
|
|
|
30 |
|
|
C Equation of State (polynomial coeffients) |
31 |
|
|
COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
32 |
|
|
_RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1) |
33 |
|
|
|
34 |
|
|
C Full Equation of State |
35 |
|
|
C eosType = 'JMD95' (Jackett and McDougall 1995, JPO) |
36 |
|
|
C eosType = 'UNESCO' (Millero et al. 1980, DSR) |
37 |
|
|
C COMMON /PARM_EOS_JMD95/ |
38 |
|
|
C eosJMDCFw :: of fresh water at pressure 0 |
39 |
|
|
C eosJMDCSw :: of sea water at pressure 0 |
40 |
|
|
C eosJMDCKFw :: of secant bulk modulus K of fresh water at pressure 0 |
41 |
|
|
C eosJMDCKSw :: of secant bulk modulus K of sea water at pressure 0 |
42 |
|
|
C eosJMDCKP :: of secant bulk modulus K at pressure p |
43 |
|
|
C end nonlinear equation of state |
44 |
|
|
_RL eosJMDCFw(6), eosJMDCSw(9) |
45 |
|
|
_RL eosJMDCKFw(5), eosJMDCKSw(7), eosJMDCKP(14) |
46 |
|
|
COMMON /PARM_EOS_JMD95/ |
47 |
|
|
& eosJMDCFw, eosJMDCSw, eosJMDCKFw, eosJMDCKSw, eosJMDCKP |
48 |
|
|
|
49 |
|
|
C pressure :: global absolute pressure variable needed for the |
50 |
|
|
C nonlinear equation of state |
51 |
|
|
_RL pressure(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
52 |
|
|
COMMON /EOS_PRESSURE/ pressure |