/[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.13 by adcroft, Mon Jun 8 18:45:28 1998 UTC revision 1.14 by cnh, Mon Jun 8 21:43:01 1998 UTC
# Line 39  C     myThid - Thread number for this in Line 39  C     myThid - Thread number for this in
39  C     == Local variables  C     == Local variables
40  C     xA, yA                 - Per block temporaries holding face areas  C     xA, yA                 - Per block temporaries holding face areas
41  C     uTrans, vTrans, wTrans - Per block temporaries holding flow transport  C     uTrans, vTrans, wTrans - Per block temporaries holding flow transport
42  C                              o uTrans: Zonal transport  C     wVel                     o uTrans: Zonal transport
43  C                              o vTrans: Meridional transport  C                              o vTrans: Meridional transport
44  C                              o wTrans: Vertical transport  C                              o wTrans: Vertical transport
45    C                              o wVel:   Vertical velocity at upper and lower
46    C                                        cell faces.
47  C     maskC,maskUp             o maskC: land/water mask for tracer cells  C     maskC,maskUp             o maskC: land/water mask for tracer cells
48  C                              o maskUp: land/water mask for W points  C                              o maskUp: land/water mask for W points
49  C     aTerm, xTerm, cTerm    - Work arrays for holding separate terms in  C     aTerm, xTerm, cTerm    - Work arrays for holding separate terms in
# Line 68  C                          into fVerTerm Line 70  C                          into fVerTerm
70        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
71        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
72        _RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
73          _RL wVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
74        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
75        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
76        _RL aTerm (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL aTerm (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 181  C--     Set up work arrays that need val Line 184  C--     Set up work arrays that need val
184          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
185           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
186            wTrans(i,j)  = 0. _d 0            wTrans(i,j)  = 0. _d 0
187              wVel  (i,j,1) = 0. _d 0
188              wVel  (i,j,2) = 0. _d 0
189            fVerT(i,j,1) = 0. _d 0            fVerT(i,j,1) = 0. _d 0
190            fVerT(i,j,2) = 0. _d 0            fVerT(i,j,2) = 0. _d 0
191            fVerS(i,j,1) = 0. _d 0            fVerS(i,j,1) = 0. _d 0
# Line 300  C--     Initial boundary condition on ba Line 305  C--     Initial boundary condition on ba
305  C--      Get temporary terms used by tendency routines  C--      Get temporary terms used by tendency routines
306           CALL CALC_COMMON_FACTORS (           CALL CALC_COMMON_FACTORS (
307       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
308       O        xA,yA,uTrans,vTrans,wTrans,maskC,maskUp,       O        xA,yA,uTrans,vTrans,wTrans,wVel,maskC,maskUp,
309       I        myThid)       I        myThid)
310    
311  C--      Calculate the total vertical diffusivity  C--      Calculate the total vertical diffusivity
# Line 310  C--      Calculate the total vertical di Line 315  C--      Calculate the total vertical di
315       O        KappaZT,       O        KappaZT,
316       I        myThid)       I        myThid)
317    
   
318  C--      Calculate accelerations in the momentum equations  C--      Calculate accelerations in the momentum equations
319           IF ( momStepping ) THEN           IF ( momStepping ) THEN
320            CALL CALC_MOM_RHS(            CALL CALC_MOM_RHS(
321       I         bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,       I         bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
322       I         xA,yA,uTrans,vTrans,wTrans,maskC,       I         xA,yA,uTrans,vTrans,wTrans,wVel,maskC,
323       I         pH,       I         pH,
324       U         aTerm,xTerm,cTerm,mTerm,pTerm,       U         aTerm,xTerm,cTerm,mTerm,pTerm,
325       U         fZon, fMer, fVerU, fVerV,       U         fZon, fMer, fVerU, fVerV,
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22