/[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.3 by jmc, Tue Jul 12 00:38:31 2005 UTC revision 1.4 by jmc, Sat Aug 11 19:01:35 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "MOM_FLUXFORM_OPTIONS.h"
5    
6  CBOP  CBOP
7  C !ROUTINE: MOM_U_METRIC_SPHERE  C !ROUTINE: MOM_U_METRIC_CYLINDER
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE MOM_U_METRIC_CYLINDER(        SUBROUTINE MOM_U_METRIC_CYLINDER(
11       I        bi,bj,k,       I        bi,bj,k,
12       I        uFld, vFld,       I        uFld, vFld,
13       O        uMetricTerms,       O        uMetricTerms,
14       I        myThid)       I        myThid )
15    
16  C !DESCRIPTION:  C !DESCRIPTION:
17  C Calculates the zonal metric term due to sphericity:  C Calculates the zonal metric term due to cylinder curvature:
18  C \begin{equation}  C \begin{equation}
19  C -\frac{v^2}{r}      C +\frac{u v}{r}
20  C \end{equation}  C \end{equation}
21    
22  C !USES: ===============================================================  C !USES: ===============================================================
# Line 47  C  i,j                  :: loop indices Line 47  C  i,j                  :: loop indices
47        _RL radius        _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            radius = (yG(i,j,bi,bj)+yG(i,j+1,bi,bj))*0.5 _d 0            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:  C-   to get the same results as before:
54  c         radius = yC(i,j,bi,bj)  c         radius = yC(i,j,bi,bj)
55            IF ( radius.GT.0. ) THEN            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))/radius       &                   vFld(i,j+1)+vFld(i-1,j+1))/radius
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22