/[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.26 by cnh, Wed Aug 19 16:20:49 1998 UTC revision 1.27 by cnh, Thu Aug 20 16:03:15 1998 UTC
# Line 75  C                          into fVerTerm Line 75  C                          into fVerTerm
75        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
76        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
77        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
78        _RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
79        _RL wVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL rVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
80        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
81        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
82        _RL aTerm (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL aTerm (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 127  C       "Calculation of Gs" Line 127  C       "Calculation of Gs"
127  C       ===================  C       ===================
128  C       This is where all the accelerations and tendencies (ie.  C       This is where all the accelerations and tendencies (ie.
129  C       physics, parameterizations etc...) are calculated  C       physics, parameterizations etc...) are calculated
130  C         w = sum_z ( div. u[n] )  C         rVel = sum_r ( div. u[n] )
131  C         rho = rho ( theta[n], salt[n] )  C         rho = rho ( theta[n], salt[n] )
132    C         b   = b(rho, theta)
133  C         K31 = K31 ( rho )  C         K31 = K31 ( rho )
134  C         Gu[n] = Gu( u[n], v[n], w, rho, Ph, ... )  C         Gu[n] = Gu( u[n], v[n], rVel, b, ... )
135  C         Gv[n] = Gv( u[n], v[n], w, rho, Ph, ... )  C         Gv[n] = Gv( u[n], v[n], rVel, b, ... )
136  C         Gt[n] = Gt( theta[n], u[n], v[n], w, K31, ... )  C         Gt[n] = Gt( theta[n], u[n], v[n], rVel, K31, ... )
137  C         Gs[n] = Gs( salt[n], u[n], v[n], w, K31, ... )  C         Gs[n] = Gs( salt[n], u[n], v[n], rVel, K31, ... )
138  C  C
139  C       "Time-stepping" or "Prediction"  C       "Time-stepping" or "Prediction"
140  C       ================================  C       ================================
# Line 198  C     uninitialised but inert locations. Line 199  C     uninitialised but inert locations.
199  C--     Set up work arrays that need valid initial values  C--     Set up work arrays that need valid initial values
200          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
201           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
202            wTrans(i,j)  = 0. _d 0            rTrans(i,j)   = 0. _d 0
203            wVel  (i,j,1) = 0. _d 0            rVel  (i,j,1) = 0. _d 0
204            wVel  (i,j,2) = 0. _d 0            rVel  (i,j,2) = 0. _d 0
205            fVerT(i,j,1) = 0. _d 0            fVerT(i,j,1)  = 0. _d 0
206            fVerT(i,j,2) = 0. _d 0            fVerT(i,j,2)  = 0. _d 0
207            fVerS(i,j,1) = 0. _d 0            fVerS(i,j,1)  = 0. _d 0
208            fVerS(i,j,2) = 0. _d 0            fVerS(i,j,2)  = 0. _d 0
209            fVerU(i,j,1) = 0. _d 0            fVerU(i,j,1)  = 0. _d 0
210            fVerU(i,j,2) = 0. _d 0            fVerU(i,j,2)  = 0. _d 0
211            fVerV(i,j,1) = 0. _d 0            fVerV(i,j,1)  = 0. _d 0
212            fVerV(i,j,2) = 0. _d 0            fVerV(i,j,2)  = 0. _d 0
213            pH(i,j,1) = 0. _d 0            phiHyd(i,j,1) = 0. _d 0
214            K13(i,j,1) = 0. _d 0            K13(i,j,1)    = 0. _d 0
215            K23(i,j,1) = 0. _d 0            K23(i,j,1)    = 0. _d 0
216            K33(i,j,1) = 0. _d 0            K33(i,j,1)    = 0. _d 0
217            KapGM(i,j) = GMkbackground            KapGM(i,j)    = GMkbackground
218           ENDDO           ENDDO
219          ENDDO          ENDDO
220    
# Line 240  C--      Update fields in layer below ac Line 241  C--      Update fields in layer below ac
241           CALL CORRECTION_STEP(           CALL CORRECTION_STEP(
242       I        bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)       I        bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)
243          ENDIF          ENDIF
244    
245  C--     Density of 1st level (below W(1)) reference to level 1  C--     Density of 1st level (below W(1)) reference to level 1
246          CALL FIND_RHO(          CALL FIND_RHO(
247       I     bi, bj, iMin, iMax, jMin, jMax, K, K, eosType,       I     bi, bj, iMin, iMax, jMin, jMax, K, K, eosType,
# Line 257  C        and mix as needed. Line 259  C        and mix as needed.
259           CALL CONVECT(           CALL CONVECT(
260       I       bi,bj,iMin,iMax,jMin,jMax,K+1,rhoKm1,rhoKp1,       I       bi,bj,iMin,iMax,jMin,jMax,K+1,rhoKm1,rhoKp1,
261       I       myTime,myIter,myThid)       I       myTime,myIter,myThid)
262    
263  C--      Recompute density after mixing  C--      Recompute density after mixing
264           CALL FIND_RHO(           CALL FIND_RHO(
265       I      bi, bj, iMin, iMax, jMin, jMax, K, K, eosType,       I      bi, bj, iMin, iMax, jMin, jMax, K, K, eosType,
# Line 284  C--       Update fields in layer below a Line 287  C--       Update fields in layer below a
287            CALL CORRECTION_STEP(            CALL CORRECTION_STEP(
288       I         bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)       I         bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)
289           ENDIF           ENDIF
290    
291  C--      Density of K level (below W(K)) reference to K level  C--      Density of K level (below W(K)) reference to K level
292           CALL FIND_RHO(           CALL FIND_RHO(
293       I      bi, bj, iMin, iMax, jMin, jMax,  K, K, eosType,       I      bi, bj, iMin, iMax, jMin, jMax,  K, K, eosType,
294       O      rhoK,       O      rhoK,
295       I      myThid )       I      myThid )
296    
297           IF ( .NOT. BOTTOM_LAYER ) THEN           IF ( .NOT. BOTTOM_LAYER ) THEN
298  C--       Check static stability with layer below  C--       Check static stability with layer below and mix as needed.
299  C         and mix as needed.  C--       Density of K+1 level (below W(K+1)) reference to K level.
 C--       Density of K+1 level (below W(K+1)) reference to K level  
300            CALL FIND_RHO(            CALL FIND_RHO(
301       I       bi, bj, iMin, iMax, jMin, jMax,  K+1, K, eosType,       I       bi, bj, iMin, iMax, jMin, jMax,  K+1, K, eosType,
302       O       rhoKp1,       O       rhoKp1,
# Line 330  C--      Calculate iso-neutral slopes fo Line 334  C--      Calculate iso-neutral slopes fo
334       I             myThid )       I             myThid )
335           DO J=jMin,jMax           DO J=jMin,jMax
336            DO I=iMin,iMax            DO I=iMin,iMax
337             rhoKm1(I,J) =rhoK(I,J)             rhoKm1 (I,J) = rhoK(I,J)
338             buoyKm1(I,J)=buoyK(I,J)             buoyKm1(I,J) = buoyK(I,J)
339            ENDDO            ENDDO
340           ENDDO           ENDDO
341    
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22