/[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.11 by jmc, Tue Feb 18 15:36:45 2003 UTC revision 1.18 by adcroft, Mon May 24 20:03:49 2004 UTC
# Line 25  C where ${\bf v}=(u,v,w)$ and $\tau$, th Line 25  C where ${\bf v}=(u,v,w)$ and $\tau$, th
25  C stresses as well as internal viscous stresses.  C stresses as well as internal viscous stresses.
26  CEOI  CEOI
27    
28  #include "CPP_OPTIONS.h"  #include "MOM_FLUXFORM_OPTIONS.h"
29    
30  CBOP  CBOP
31  C !ROUTINE: MOM_FLUXFORM  C !ROUTINE: MOM_FLUXFORM
# Line 122  C     uDudxFac, AhDudxFac, etc ... indiv Line 122  C     uDudxFac, AhDudxFac, etc ... indiv
122        _RL  vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
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)
126          _RL viscAhD(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
127          _RL viscAhZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
128          _RL viscA4D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
129          _RL viscA4Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
130          _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
131          _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
132          _RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
133          _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
134  C     I,J,K - Loop counters  C     I,J,K - Loop counters
135  C     rVelMaskOverride - Factor for imposing special surface boundary conditions  C     rVelMaskOverride - Factor for imposing special surface boundary conditions
136  C                        ( set according to free-surface condition ).  C                        ( set according to free-surface condition ).
# Line 155  C     xxxFac - On-off tracer parameters Line 164  C     xxxFac - On-off tracer parameters
164        INTEGER km1,kp1        INTEGER km1,kp1
165        _RL wVelBottomOverride        _RL wVelBottomOverride
166        LOGICAL bottomDragTerms        LOGICAL bottomDragTerms
       _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
167  CEOP  CEOP
168    
169        km1=MAX(1,k-1)        km1=MAX(1,k-1)
# Line 179  C     Initialise intermediate terms Line 187  C     Initialise intermediate terms
187          fMer(i,j) = 0.          fMer(i,j) = 0.
188          rTransU(i,j) = 0.          rTransU(i,j) = 0.
189          rTransV(i,j) = 0.          rTransV(i,j) = 0.
190            strain(i,j) = 0.
191            tension(i,j) = 0.
192         ENDDO         ENDDO
193        ENDDO        ENDDO
194    
# Line 253  C     Calculate velocity field "volume t Line 263  C     Calculate velocity field "volume t
263         ENDDO         ENDDO
264        ENDDO        ENDDO
265    
266        CALL MOM_CALC_KE(bi,bj,k,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
269             CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,
270         O                         tension,
271         I                         myThid)
272             CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,
273         O                        strain,
274         I                        myThid)
275    c     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 289  C---  Calculate vertical transports (at Line 308  C---  Calculate vertical transports (at
308       I                       myTime, myIter, myThid)       I                       myTime, myIter, myThid)
309        ENDIF        ENDIF
310    
311    c     IF (momViscosity) THEN
312    c    &  CALL MOM_CALC_VISCOSITY(bi,bj,k,
313    c    I                         uFld,vFld,
314    c    O                         viscAhD,viscAhZ,myThid)
315    
316  C---- Zonal momentum equation starts here  C---- Zonal momentum equation starts here
317    
# Line 409  C-    No-slip BCs impose a drag at botto Line 432  C-    No-slip BCs impose a drag at botto
432         ENDDO         ENDDO
433        ENDIF        ENDIF
434    
435  C--   Forcing term  C--   Forcing term (moved to timestep.F)
436        IF (momForcing)  c     IF (momForcing)
437       &  CALL EXTERNAL_FORCING_U(  c    &  CALL EXTERNAL_FORCING_U(
438       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
439       I     myTime,myThid)  c    I     myTime,myThid)
440    
441  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
442        IF (useNHMTerms) THEN        IF (useNHMTerms) THEN
# Line 561  C-    No-slip BCs impose a drag at botto Line 584  C-    No-slip BCs impose a drag at botto
584         ENDDO         ENDDO
585        ENDIF        ENDIF
586    
587  C--   Forcing term  C--   Forcing term (moved to timestep.F)
588        IF (momForcing)  c     IF (momForcing)
589       & CALL EXTERNAL_FORCING_V(  c    & CALL EXTERNAL_FORCING_V(
590       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
591       I     myTime,myThid)  c    I     myTime,myThid)
592    
593  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
594        IF (useNHMTerms) THEN        IF (useNHMTerms) THEN
# Line 595  C--   Set dv/dt on boundaries to zero Line 618  C--   Set dv/dt on boundaries to zero
618    
619  C--   Coriolis term  C--   Coriolis term
620  C     Note. As coded here, coriolis will not work with "thin walls"  C     Note. As coded here, coriolis will not work with "thin walls"
621  #ifdef INCLUDE_CD_CODE  c     IF (useCDscheme) THEN
622        CALL MOM_CDSCHEME(bi,bj,k,dPhiHydX,dPhiHydY,myThid)  c       CALL MOM_CDSCHEME(bi,bj,k,dPhiHydX,dPhiHydY,myThid)
623  #else  c     ELSE
624        CALL MOM_U_CORIOLIS(bi,bj,k,vFld,cf,myThid)        IF (.NOT.useCDscheme) THEN
625        DO j=jMin,jMax          CALL MOM_U_CORIOLIS(bi,bj,k,vFld,cf,myThid)
626         DO i=iMin,iMax          DO j=jMin,jMax
627          gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j)           DO i=iMin,iMax
628         ENDDO            gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j)
629        ENDDO           ENDDO
630        CALL MOM_V_CORIOLIS(bi,bj,k,uFld,cf,myThid)          ENDDO
631        DO j=jMin,jMax          CALL MOM_V_CORIOLIS(bi,bj,k,uFld,cf,myThid)
632         DO i=iMin,iMax          DO j=jMin,jMax
633          gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+fvFac*cf(i,j)           DO i=iMin,iMax
634         ENDDO            gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+fvFac*cf(i,j)
635        ENDDO           ENDDO
636  #endif /* INCLUDE_CD_CODE */          ENDDO
637          ENDIF
638    
639        IF (nonHydrostatic.OR.quasiHydrostatic) THEN        IF (nonHydrostatic.OR.quasiHydrostatic) THEN
640         CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid)         CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid)
641         DO j=jMin,jMax         DO j=jMin,jMax
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22