/[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.38 by cnh, Fri Nov 6 22:44:45 1998 UTC revision 1.41 by adcroft, Mon May 3 21:45:57 1999 UTC
# Line 20  C     | ===== Line 20  C     | =====
20  C     | C*P* comments indicating place holders for which code is |  C     | C*P* comments indicating place holders for which code is |
21  C     |      presently being developed.                          |  C     |      presently being developed.                          |
22  C     \==========================================================/  C     \==========================================================/
23          IMPLICIT NONE
24    
25  C     == Global variables ===  C     == Global variables ===
26  #include "SIZE.h"  #include "SIZE.h"
# Line 251  C--     Update fields in top level accor Line 252  C--     Update fields in top level accor
252          CALL CORRECTION_STEP(          CALL CORRECTION_STEP(
253       I       bi,bj,iMin,iMax,jMin,jMax,K,       I       bi,bj,iMin,iMax,jMin,jMax,K,
254       I       etaSurfX,etaSurfY,myTime,myThid)       I       etaSurfX,etaSurfY,myTime,myThid)
255            IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K, myThid )
256          IF ( .NOT. BOTTOM_LAYER ) THEN          IF ( .NOT. BOTTOM_LAYER ) THEN
257  C--      Update fields in layer below according to tendency terms  C--      Update fields in layer below according to tendency terms
258           CALL CORRECTION_STEP(           CALL CORRECTION_STEP(
259       I        bi,bj,iMin,iMax,jMin,jMax,K+1,       I        bi,bj,iMin,iMax,jMin,jMax,K+1,
260       I        etaSurfX,etaSurfY,myTime,myThid)       I        etaSurfX,etaSurfY,myTime,myThid)
261             IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid )
262          ENDIF          ENDIF
263  #endif  #endif
264  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 308  C--       Update fields in layer below a Line 311  C--       Update fields in layer below a
311            CALL CORRECTION_STEP(            CALL CORRECTION_STEP(
312       I         bi,bj,iMin,iMax,jMin,jMax,K+1,       I         bi,bj,iMin,iMax,jMin,jMax,K+1,
313       I         etaSurfX,etaSurfY,myTime,myThid)       I         etaSurfX,etaSurfY,myTime,myThid)
314              IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid )
315           ENDIF           ENDIF
316  #endif  #endif
317  C--      Density of K level (below W(K)) reference to K level  C--      Density of K level (below W(K)) reference to K level
# Line 428  C--      Prediction step (step forward a Line 432  C--      Prediction step (step forward a
432           CALL TIMESTEP(           CALL TIMESTEP(
433       I       bi,bj,iMin,iMax,jMin,jMax,K,       I       bi,bj,iMin,iMax,jMin,jMax,K,
434       I       myThid)       I       myThid)
435    C--      Apply open boundary conditions
436             IF (openBoundaries) CALL APPLY_OBCS2( bi, bj, K, myThid )
437    C--      Freeze water
438             IF (allowFreezing)
439         &   CALL FREEZE( bi, bj, iMin, iMax, jMin, jMax, K, myThid )
440  C--      Diagnose barotropic divergence of predicted fields  C--      Diagnose barotropic divergence of predicted fields
441           CALL CALC_DIV_GHAT(           CALL CALC_DIV_GHAT(
442       I       bi,bj,iMin,iMax,jMin,jMax,K,       I       bi,bj,iMin,iMax,jMin,jMax,K,
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22