/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_fluxform.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_fluxform.F

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

revision 1.20 by afe, Thu Oct 21 21:14:02 2004 UTC revision 1.21 by jmc, Fri Oct 29 16:25:37 2004 UTC
# Line 123  C     uDudxFac, AhDudxFac, etc ... indiv Line 123  C     uDudxFac, AhDudxFac, etc ... indiv
123        _RL  rTransU(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  rTransU(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
124        _RL  rTransV(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  rTransV(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
125        _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
126        _RL viscAhD(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL viscAh_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
127        _RL viscAhZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL viscAh_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
128        _RL viscA4D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL viscA4_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
129        _RL viscA4Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL viscA4_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
130        _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
131        _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  c     _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
132        _RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
133        _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
134  C     I,J,K - Loop counters  C     I,J,K - Loop counters
# Line 265  C     Calculate velocity field "volume t Line 265  C     Calculate velocity field "volume t
265    
266        CALL MOM_CALC_KE(bi,bj,k,3,uFld,vFld,KE,myThid)        CALL MOM_CALC_KE(bi,bj,k,3,uFld,vFld,KE,myThid)
267    
268  c     IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN        IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN
269           CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,           CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,
270       O                         tension,       O                         tension,
271       I                         myThid)       I                         myThid)
272           CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,           CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,
273       O                        strain,       O                        strain,
274       I                        myThid)       I                        myThid)
275  c     ENDIF        ENDIF
276    
277  C---  First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp)  C---  First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp)
278        IF (momAdvection.AND.k.EQ.1) THEN        IF (momAdvection.AND.k.EQ.1) THEN
# Line 311  C---  Calculate vertical transports (at Line 311  C---  Calculate vertical transports (at
311  c     IF (momViscosity) THEN  c     IF (momViscosity) THEN
312  c    &  CALL MOM_CALC_VISCOSITY(bi,bj,k,  c    &  CALL MOM_CALC_VISCOSITY(bi,bj,k,
313  c    I                         uFld,vFld,  c    I                         uFld,vFld,
314  c    O                         viscAhD,viscAhZ,myThid)  c    O                         viscAh_D,viscAh_Z,myThid)
315    
316  C---- Zonal momentum equation starts here  C---- Zonal momentum equation starts here
317    
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22