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

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

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

revision 1.44 by jmc, Fri Apr 4 20:54:11 2014 UTC revision 1.45 by jmc, Mon May 12 01:27:47 2014 UTC
# Line 91  C     pKappaC    :: (p/Po)^kappa at cell Line 91  C     pKappaC    :: (p/Po)^kappa at cell
91        _RL pKappaU (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL pKappaU (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
92        _RL pKappaC, locDepth, fullDepth        _RL pKappaC, locDepth, fullDepth
93  #endif /* DISABLE_SIGMA_CODE */  #endif /* DISABLE_SIGMA_CODE */
94        _RL dRlocM,dRlocP, ddRloc, locAlpha        _RL thetaRef, locAlpha
95          _RL dRlocM,dRlocP, ddRloc
96        _RL ddPIm, ddPIp, rec_dRm, rec_dRp        _RL ddPIm, ddPIp, rec_dRm, rec_dRp
97        _RL surfPhiFac        _RL surfPhiFac
98        LOGICAL useDiagPhiRlow, addSurfPhiAnom        LOGICAL useDiagPhiRlow, addSurfPhiAnom
# Line 428  C       The ideal gas law is used implic Line 429  C       The ideal gas law is used implic
429  C       the specific volume, analogous to the oceanic case.  C       the specific volume, analogous to the oceanic case.
430    
431  C--     virtual potential temperature anomaly (including water vapour effect)  C--     virtual potential temperature anomaly (including water vapour effect)
432            IF ( select_rStar.GE.1 .OR. selectSigmaCoord.GE.1 ) THEN
433    C-      isothermal (theta=const) reference state
434              thetaRef = thetaConst
435            ELSE
436    C-      horizontally uniform (tRef) reference state
437              thetaRef = tRef(k)
438            ENDIF
439          DO j=jMin,jMax          DO j=jMin,jMax
440           DO i=iMin,iMax           DO i=iMin,iMax
441            alphaRho(i,j) = ( tFld(i,j,k,bi,bj)            alphaRho(i,j) = ( tFld(i,j,k,bi,bj)
442       &                      *( sFld(i,j,k,bi,bj)*atm_Rq + oneRL )       &                      *( sFld(i,j,k,bi,bj)*atm_Rq + oneRL )
443       &                    - tRef(k) )*maskC(i,j,k,bi,bj)       &                    - thetaRef )*maskC(i,j,k,bi,bj)
444           ENDDO           ENDDO
445          ENDDO          ENDDO
446    
Legend:
Removed from v.1.44  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22