/[MITgcm]/MITgcm/model/inc/EOS.h
ViewVC logotype

Contents of /MITgcm/model/inc/EOS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.2 - (show annotations) (download)
Thu Sep 5 20:49:33 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, checkpoint46f_post, checkpoint48e_post, checkpoint46l_pre, checkpoint48b_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint46j_pre, checkpoint47g_post, checkpoint46j_post, checkpoint46k_post, checkpoint48a_post, checkpoint47j_post, branch-exfmods-tag, checkpoint48c_post, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46g_post, checkpoint47f_post, checkpoint46i_post, checkpoint46e_post, checkpoint47, checkpoint48, checkpoint46h_post, checkpoint47h_post
Branch point for: branch-exfmods-curt
Changes since 1.1: +12 -4 lines
File MIME type: text/plain 
o Added new equation of state -> MDJWF
  - EOS of McDougall et al., 2002, JAOT, submitted
  - caveat: the equation of state is only valid for a smaller (more
    realistic?) range of values than JMD95P/Z and UNESCO
  - added masks to the calculation of pressure in store_pressure
  - added more check values for density in check_eos (ini_eos.F), some of
    the old check values are out of the range of the MDJWF-eos, so don't
    expect perfect matches for those
1 C $Header: /u/gcmpack/MITgcm/model/inc/EOS.h,v 1.1 2002/08/07 16:55:52 mlosch 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 = '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 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 _RL eosMDJWFnum(0:11), eosMDJWFden(0:12)
54 COMMON /PARM_EOS_MDJWF/
55 & eosMDJWFnum, eosMDJWFden
56
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
  ViewVC Help
Powered by ViewVC 1.1.22