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

Diff of /MITgcm/pkg/mom_fluxform/mom_u_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:38:31 2005 UTC
# Line 44  C  uMetricTerms         :: metric term Line 44  C  uMetricTerms         :: 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    
50        DO j=1-Oly,sNy+Oly-1        DO j=1-Oly,sNy+Oly-1
51         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
52            if (yC(i,j,bi,bj).GT. 0) then            radius = (yG(i,j,bi,bj)+yG(i,j+1,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               uMetricTerms(i,j) =               uMetricTerms(i,j) =
57       &            uFld(i,j)       &            uFld(i,j)
58       &            *0.25*(vFld(i,j) + vFld(i-1,j)       &            *0.25*(vFld(i,j) + vFld(i-1,j) +
59       &            +    vFld(i,j+1) + vFld(i-1,j+1))/yC(i,j,bi,bj)       &                   vFld(i,j+1)+vFld(i-1,j+1))/radius
60            else            ELSE
61               uMetricTerms(i,j) = 0.               uMetricTerms(i,j) = 0.
62            endif            ENDIF
63    
64         ENDDO         ENDDO
65        ENDDO        ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22