/[MITgcm]/MITgcm/model/src/seawater.F
ViewVC logotype

Diff of /MITgcm/model/src/seawater.F

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

revision 1.2 by jmc, Sat Aug 9 00:51:34 2008 UTC revision 1.3 by jmc, Sat Sep 6 17:42:27 2008 UTC
# Line 19  C--              used by sw_ptmp Line 19  C--              used by sw_ptmp
19  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
20  C     *==========================================================*  C     *==========================================================*
21  C     | o SUBROUTINE FIND_RHO_SCALAR  C     | o SUBROUTINE FIND_RHO_SCALAR
22  C     |   Calculates [rho(S,T,p)-rhoConst]  C     |   Calculates rho(S,T,p)
23  C     *==========================================================*  C     *==========================================================*
24  C     \ev  C     \ev
25    
# Line 81  C     issue a warning Line 81  C     issue a warning
81           rholoc = rhoNil*(           rholoc = rhoNil*(
82       &                      sBeta *(sLoc-sRef(1))       &                      sBeta *(sLoc-sRef(1))
83       &                     -tAlpha*(tLoc-tRef(1))       &                     -tAlpha*(tLoc-tRef(1))
84       &                   ) + (rhoNil-rhoConst)       &                   ) + rhoNil
85  c        rhoLoc = 0. _d  0  c        rhoLoc = 0. _d  0
86    
87        ELSEIF (equationOfState.EQ.'POLY3') THEN        ELSEIF (equationOfState.EQ.'POLY3') THEN
# Line 175  C     secant bulk modulus of sea water a Line 175  C     secant bulk modulus of sea water a
175           bulkMod = bMfresh + bMsalt + bMpres           bulkMod = bMfresh + bMsalt + bMpres
176    
177  C     density of sea water at pressure p  C     density of sea water at pressure p
178           rhoLoc = rhoP0/(1. _d 0 - p1/bulkMod) - rhoConst           rhoLoc = rhoP0/(1. _d 0 - p1/bulkMod)
179    
180        ELSEIF ( equationOfState.EQ.'MDJWF' ) THEN        ELSEIF ( equationOfState.EQ.'MDJWF' ) THEN
181    
# Line 207  C     density of sea water at pressure p Line 207  C     density of sea water at pressure p
207    
208           rhoDen = 1.0/(epsln+den)           rhoDen = 1.0/(epsln+den)
209    
210           rhoLoc = rhoNum*rhoDen - rhoConst           rhoLoc = rhoNum*rhoDen
211    
212        ELSEIF( equationOfState .EQ. 'IDEALG' ) THEN        ELSEIF( equationOfState .EQ. 'IDEALG' ) THEN
213  C  C
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22