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

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

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

revision 1.28 by cnh, Thu Aug 20 19:25:05 1998 UTC revision 1.29 by cnh, Thu Aug 20 19:26:40 1998 UTC
# Line 65  C     buoyK, buoyKM1 - Buoyancy at curre Line 65  C     buoyK, buoyKM1 - Buoyancy at curre
65  C     phiHyd         - Hydrostatic part of the potential phi.  C     phiHyd         - Hydrostatic part of the potential phi.
66  C                      In z coords phiHyd is the hydrostatic pressure anomaly  C                      In z coords phiHyd is the hydrostatic pressure anomaly
67  C                      In p coords phiHyd is the geopotential surface height anomaly.  C                      In p coords phiHyd is the geopotential surface height anomaly.
68    C     etaSurfX, etaSurfY - Holds surface elevation gradient in X and Y.
69  C     iMin, iMax - Ranges and sub-block indices on which calculations  C     iMin, iMax - Ranges and sub-block indices on which calculations
70  C     jMin, jMax   are applied.  C     jMin, jMax   are applied.
71  C     bi, bj  C     bi, bj
# Line 97  C                          into fVerTerm Line 98  C                          into fVerTerm
98        _RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
99        _RL buoyK (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL buoyK (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
100        _RL rhotmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rhotmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
101        _RL pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL etaSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
102        _RL pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL etaSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
103        _RL K13   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)        _RL K13   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)
104        _RL K23   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)        _RL K23   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)
105        _RL K33   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)        _RL K33   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz)
# Line 229  C--     Set up work arrays that need val Line 230  C--     Set up work arrays that need val
230  C--     Calculate gradient of surface pressure  C--     Calculate gradient of surface pressure
231          CALL CALC_GRAD_ETA_SURF(          CALL CALC_GRAD_ETA_SURF(
232       I       bi,bj,iMin,iMax,jMin,jMax,       I       bi,bj,iMin,iMax,jMin,jMax,
233       O       pSurfX,pSurfY,       O       etaSurfX,etaSurfY,
234       I       myThid)       I       myThid)
235    
236  C--     Update fields in top level according to tendency terms  C--     Update fields in top level according to tendency terms
237          CALL CORRECTION_STEP(          CALL CORRECTION_STEP(
238       I       bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myTime,myThid)       I       bi,bj,iMin,iMax,jMin,jMax,K,etaSurfX,etaSurfY,myTime,myThid)
239    
240          IF ( .NOT. BOTTOM_LAYER ) THEN          IF ( .NOT. BOTTOM_LAYER ) THEN
241  C--      Update fields in layer below according to tendency terms  C--      Update fields in layer below according to tendency terms
242           CALL CORRECTION_STEP(           CALL CORRECTION_STEP(
243       I        bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)       I        bi,bj,iMin,iMax,jMin,jMax,K+1,etaSurfX,etaSurfY,myTime,myThid)
244          ENDIF          ENDIF
245    
246  C--     Density of 1st level (below W(1)) reference to level 1  C--     Density of 1st level (below W(1)) reference to level 1
# Line 285  C--     Integrate hydrostatic balance fo Line 286  C--     Integrate hydrostatic balance fo
286           IF ( .NOT. BOTTOM_LAYER ) THEN           IF ( .NOT. BOTTOM_LAYER ) THEN
287  C--       Update fields in layer below according to tendency terms  C--       Update fields in layer below according to tendency terms
288            CALL CORRECTION_STEP(            CALL CORRECTION_STEP(
289       I         bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)       I         bi,bj,iMin,iMax,jMin,jMax,K+1,etaSurfX,etaSurfY,myTime,myThid)
290           ENDIF           ENDIF
291    
292  C--      Density of K level (below W(K)) reference to K level  C--      Density of K level (below W(K)) reference to K level
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22