/[MITgcm]/MITgcm/pkg/mom_common/mom_v_metric_nh.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_v_metric_nh.F

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

revision 1.2 by jmc, Tue Dec 5 05:27:38 2006 UTC revision 1.3 by jmc, Mon Mar 12 23:54:26 2007 UTC
# Line 43  C  vMetricTerms         :: metric term Line 43  C  vMetricTerms         :: metric term
43    
44  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
45  C  i,j                  :: loop indices  C  i,j                  :: loop indices
46        INTEGER i,j,Kp1        INTEGER i,j,kp1
47        _RL wVelBottomOverride        _RL wVelBottomOverride
48  CEOP  CEOP
49    
50        Kp1=min(k+1,Nr)        kp1=min(k+1,Nr)
51        wVelBottomOverride=1.        wVelBottomOverride=1.
52        IF (k.EQ.Nr) wVelBottomOverride=0.        IF (k.EQ.Nr) wVelBottomOverride=0.
53    
54        DO j=1-Olx+1,sNy+Oly        DO j=1-Olx+1,sNy+Oly
55         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
56           vMetricTerms(i,j) = vFld(i,j)*recip_rSphere*recip_deepFacC(k)           vMetricTerms(i,j) = vFld(i,j)*recip_rSphere*recip_deepFacC(k)
57       &    *0.25*(wVelBottomOverride*       &    *0.25*( (wFld(i,j-1,kp1,bi,bj)+wFld(i,j,kp1,bi,bj))
58       &           (wFld(i,j,kp1,bi,bj)+wFld(i,j-1,kp1,bi,bj))       &            *rVel2wUnit(kp1)*wVelBottomOverride
59       &           +wFld(i,j, k ,bi,bj)+wFld(i,j-1, k ,bi,bj)       &          + (wFld(i,j-1, k ,bi,bj)+wFld(i,j, k ,bi,bj))
60       &          )*recip_horiVertRatio*gravitySign       &            *rVel2wUnit( k )
61         &          )*gravitySign
62         ENDDO         ENDDO
63        ENDDO        ENDDO
64    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22