/[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.29 by jmc, Wed Oct 12 01:10:10 2005 UTC revision 1.32 by mlosch, Tue Feb 7 11:46:18 2006 UTC
# Line 166  C     Initialise intermediate terms Line 166  C     Initialise intermediate terms
166          tension(i,j)= 0.          tension(i,j)= 0.
167          guDiss(i,j) = 0.          guDiss(i,j) = 0.
168          gvDiss(i,j) = 0.          gvDiss(i,j) = 0.
169    #ifdef ALLOW_AUTODIFF_TAMC
170            vort3(i,j)   = 0. _d 0
171            strain(i,j)  = 0. _d 0
172            tension(i,j) = 0. _d 0
173    #endif
174         ENDDO         ENDDO
175        ENDDO        ENDDO
176    
# Line 349  C--   Tendency is minus divergence of th Line 354  C--   Tendency is minus divergence of th
354    
355  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
356  C-- account for 3.D divergence of the flow in rStar coordinate:  C-- account for 3.D divergence of the flow in rStar coordinate:
357    # ifndef DISABLE_RSTAR_CODE
358          IF ( select_rStar.GT.0 ) THEN          IF ( select_rStar.GT.0 ) THEN
359           DO j=jMin,jMax           DO j=jMin,jMax
360            DO i=iMin,iMax            DO i=iMin,iMax
# Line 366  C-- account for 3.D divergence of the fl Line 372  C-- account for 3.D divergence of the fl
372            ENDDO            ENDDO
373           ENDDO           ENDDO
374          ENDIF          ENDIF
375    # endif /* DISABLE_RSTAR_CODE */
376  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
377    
378        ELSE        ELSE
# Line 453  C-    No-slip BCs impose a drag at botto Line 460  C-    No-slip BCs impose a drag at botto
460           ENDDO           ENDDO
461          ENDIF          ENDIF
462    
463    #ifdef ALLOW_SHELFICE
464            IF (useShelfIce) THEN
465             CALL SHELFICE_U_DRAG(bi,bj,k,uFld,KE,KappaRU,vF,myThid)
466             DO j=jMin,jMax
467              DO i=iMin,iMax
468               gUdiss(i,j) = gUdiss(i,j) + vF(i,j)
469              ENDDO
470             ENDDO
471            ENDIF
472    #endif /* ALLOW_SHELFICE */
473    
474  C-    endif momViscosity  C-    endif momViscosity
475        ENDIF        ENDIF
476    
# Line 537  C--   Tendency is minus divergence of th Line 555  C--   Tendency is minus divergence of th
555    
556  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
557  C-- account for 3.D divergence of the flow in rStar coordinate:  C-- account for 3.D divergence of the flow in rStar coordinate:
558    # ifndef DISABLE_RSTAR_CODE
559          IF ( select_rStar.GT.0 ) THEN          IF ( select_rStar.GT.0 ) THEN
560           DO j=jMin,jMax           DO j=jMin,jMax
561            DO i=iMin,iMax            DO i=iMin,iMax
# Line 554  C-- account for 3.D divergence of the fl Line 573  C-- account for 3.D divergence of the fl
573            ENDDO            ENDDO
574           ENDDO           ENDDO
575          ENDIF          ENDIF
576    # endif /* DISABLE_RSTAR_CODE */
577  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
578    
579        ELSE        ELSE
# Line 615  C--   Tendency is minus divergence of th Line 635  C--   Tendency is minus divergence of th
635  #endif  #endif
636    
637  C-- No-slip and drag BCs appear as body forces in cell abutting topography  C-- No-slip and drag BCs appear as body forces in cell abutting topography
638        IF (no_slip_sides) THEN          IF (no_slip_sides) THEN
639  C-     No-slip BCs impose a drag at walls...  C-     No-slip BCs impose a drag at walls...
640           CALL MOM_V_SIDEDRAG(           CALL MOM_V_SIDEDRAG(
641       I        bi,bj,k,       I        bi,bj,k,
# Line 640  C-    No-slip BCs impose a drag at botto Line 660  C-    No-slip BCs impose a drag at botto
660           ENDDO           ENDDO
661          ENDIF          ENDIF
662    
663    #ifdef ALLOW_SHELFICE
664            IF (useShelfIce) THEN
665             CALL SHELFICE_V_DRAG(bi,bj,k,vFld,KE,KappaRU,vF,myThid)
666             DO j=jMin,jMax
667              DO i=iMin,iMax
668               gvDiss(i,j) = gvDiss(i,j) + vF(i,j)
669              ENDDO
670             ENDDO
671            ENDIF
672    #endif /* ALLOW_SHELFICE */
673    
674  C-    endif momViscosity  C-    endif momViscosity
675        ENDIF        ENDIF
676    
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22