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

Diff of /MITgcm/pkg/mom_vecinv/mom_vecinv.F

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

revision 1.8 by heimbach, Fri Oct 10 23:00:01 2003 UTC revision 1.9 by heimbach, Sun Oct 26 01:01:23 2003 UTC
# Line 115  C     xxxFac - On-off tracer parameters Line 115  C     xxxFac - On-off tracer parameters
115        _RL  phyFac        _RL  phyFac
116        _RL  vForcFac        _RL  vForcFac
117        _RL  mtFacV        _RL  mtFacV
       INTEGER km1,kp1  
118        _RL wVelBottomOverride        _RL wVelBottomOverride
119        LOGICAL bottomDragTerms        LOGICAL bottomDragTerms
120        _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 123  C     xxxFac - On-off tracer parameters Line 122  C     xxxFac - On-off tracer parameters
122        _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
123        _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
124    
125        km1=MAX(1,k-1)  #ifdef ALLOW_AUTODIFF_TAMC
126        kp1=MIN(Nr,k+1)  C--   only the kDown part of fverU/V is set in this subroutine
127    C--   the kUp is still required
128    C--   In the case of mom_fluxform Kup is set as well
129    C--   (at least in part)
130          fVerU(1,1,kUp) = fVerU(1,1,kUp)
131          fVerV(1,1,kUp) = fVerV(1,1,kUp)
132    #endif
133    
134        rVelMaskOverride=1.        rVelMaskOverride=1.
135        IF ( k .EQ. 1 ) rVelMaskOverride=freeSurfFac        IF ( k .EQ. 1 ) rVelMaskOverride=freeSurfFac
136        wVelBottomOverride=1.        wVelBottomOverride=1.
# Line 152  C     Initialise intermediate terms Line 158  C     Initialise intermediate terms
158  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
159          strain(i,j)  = 0. _d 0          strain(i,j)  = 0. _d 0
160          tension(i,j) = 0. _d 0          tension(i,j) = 0. _d 0
         fVerU(i,j,1) = 0. _d 0  
         fVerU(i,j,2) = 0. _d 0  
         fVerV(i,j,1) = 0. _d 0  
         fVerV(i,j,2) = 0. _d 0  
161  #endif  #endif
162         ENDDO         ENDDO
163        ENDDO        ENDDO
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22