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

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

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

revision 1.30 by baylor, Mon Oct 10 15:56:34 2005 UTC revision 1.31 by jmc, Tue Oct 11 00:09:02 2005 UTC
# Line 34  C     == Global variables == Line 34  C     == Global variables ==
34  #include "PARAMS.h"  #include "PARAMS.h"
35  #include "EOS.h"  #include "EOS.h"
36  #include "GRID.h"  #include "GRID.h"
 #include "SURFACE.h"  
37    
38  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
39  C     == Routine arguments ==  C     == Routine arguments ==
# Line 50  C     kRef :: Pressure reference level Line 49  C     kRef :: Pressure reference level
49    
50  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
51  C     == Local variables ==  C     == Local variables ==
52        INTEGER i,j,kr        INTEGER i,j
53        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig,dRho        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig,dRho
54        _RL locPres(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL locPres(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
55        _RL rhoP0  (1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL rhoP0  (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
# Line 77  CEOP Line 76  CEOP
76    
77        IF (equationOfState.EQ.'LINEAR') THEN        IF (equationOfState.EQ.'LINEAR') THEN
78    
79         IF ( usingPCoords ) THEN  C ***NOTE***
80          kr        = Nr  C In the linear EOS, to make the static stability calculation meaningful
81         ELSE  C we alway calculate the perturbation with respect to the surface level.
82          kr        = 1  C **********
83         ENDIF         refTemp=tRef(kRef)
84         refTemp=tRef(kr)         refSalt=sRef(kRef)
85         refSalt=sRef(kr)  
86         dRho = rhoNil-rhoConst         dRho = rhoNil-rhoConst
87    
88         DO j=jMin,jMax         DO j=jMin,jMax
# Line 94  CEOP Line 93  CEOP
93       &        + dRho       &        + dRho
94          ENDDO          ENDDO
95         ENDDO         ENDDO
96          
97        ELSEIF (equationOfState.EQ.'POLY3') THEN        ELSEIF (equationOfState.EQ.'POLY3') THEN
98    
99         refTemp=eosRefT(kRef)         refTemp=eosRefT(kRef)
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22