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

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

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

revision 1.2 by jmc, Wed Jul 7 18:34:07 2004 UTC revision 1.3 by jmc, Tue Jul 12 00:44:12 2005 UTC
# Line 44  C  vMetricTerms         :: metric term Line 44  C  vMetricTerms         :: metric term
44  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
45  C  i,j                  :: loop indices  C  i,j                  :: loop indices
46        INTEGER i,j        INTEGER i,j
47          _RL radius
48  CEOP  CEOP
49    
 C jmc: using yG make more sense to me but I might be wrong ?  
50        DO j=1-Oly+1,sNy+Oly        DO j=1-Oly+1,sNy+Oly
51         DO i=1-Olx,sNx+Olx-1         DO i=1-Olx,sNx+Olx-1
52            if (yC(i,j,bi,bj).GT. 0) then            radius = (yG(i,j,bi,bj)+yG(i+1,j,bi,bj))*0.5 _d 0
53    C-   to get the same results as before:
54    c         radius = yC(i,j,bi,bj)
55              IF ( radius.GT.0. ) THEN
56               vMetricTerms(i,j) =               vMetricTerms(i,j) =
57       &            -1.*0.25*(uFld(i,j) + uFld(i+1,j) +       &           -0.25*(uFld(i,j) + uFld(i+1,j) +
58       &            uFld(i,j-1)+uFld(i+1,j-1))*       &                  uFld(i,j-1)+uFld(i+1,j-1))
59       &            0.25*(uFld(i,j) + uFld(i+1,j) +       &           *0.25*(uFld(i,j) + uFld(i+1,j) +
60       &            uFld(i,j-1)+uFld(i+1,j-1))/yC(i,j,bi,bj)       &                  uFld(i,j-1)+uFld(i+1,j-1))/radius
61            else            ELSE
62               vMetricTerms(i,j) = 0.               vMetricTerms(i,j) = 0.
63            endif            ENDIF
64         ENDDO         ENDDO
65        ENDDO        ENDDO
66    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22