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

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

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

Revision 1.2 - (hide 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 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
  ViewVC Help
Powered by ViewVC 1.1.22