/[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.30 by jmc, Fri Aug 1 04:03:54 2003 UTC revision 1.33 by mlosch, Tue Feb 7 11:47:49 2006 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CBOP  CBOP
# Line 148  CADJ STORE sFld (:,:,k,bi,bj) = comlev1_ Line 149  CADJ STORE sFld (:,:,k,bi,bj) = comlev1_
149          CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k,          CALL FIND_RHO( bi, bj, iMin, iMax, jMin, jMax, k, k,
150       &                 tFld, sFld,       &                 tFld, sFld,
151       &                 alphaRho, myThid)       &                 alphaRho, myThid)
152    #ifdef ALLOW_SHELFICE
153    C     mask rho, so that there is no contribution of phiHyd from
154    C     overlying shelfice (whose density we do not know)
155            IF ( useShelfIce ) THEN
156             DO j=jMin,jMax
157              DO i=iMin,iMax
158               alphaRho(i,j) = alphaRho(i,j)*maskC(i,j,k,bi,bj)
159              ENDDO
160             ENDDO
161            ENDIF
162    #endif /* ALLOW_SHELFICE */
163    
164    #ifdef ALLOW_DIAGNOSTICS
165            IF ( useDiagnostics )
166         &   CALL DIAGNOSTICS_FILL(alphaRho,'RHOAnoma',k,1,2,bi,bj,myThid)
167    #endif
168    
169  C Quasi-hydrostatic terms are added in as if they modify the buoyancy  C Quasi-hydrostatic terms are added in as if they modify the buoyancy
170          IF (quasiHydrostatic) THEN          IF (quasiHydrostatic) THEN
# Line 233  CADJ STORE sFld (:,:,k,bi,bj) = comlev1_ Line 250  CADJ STORE sFld (:,:,k,bi,bj) = comlev1_
250  CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE alphaRho (:,:) = comlev1_bibj_k, key=kkey, byte=isbyte
251  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
252    
253    #ifdef ALLOW_DIAGNOSTICS
254            IF ( useDiagnostics )
255         &   CALL DIAGNOSTICS_FILL(alphaRho,'RHOAnoma',k,1,2,bi,bj,myThid)
256    #endif
257    
258  C--     Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst  C--     Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst
259          DO j=jMin,jMax          DO j=jMin,jMax
260            DO i=iMin,iMax            DO i=iMin,iMax
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22