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

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

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

revision 1.6 by jmc, Thu Sep 17 18:21:55 2009 UTC revision 1.7 by jmc, Thu Sep 17 21:16:21 2009 UTC
# Line 53  C-    used for side-drag: Line 53  C-    used for side-drag:
53  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
54    
55         IF (deltaTmom.NE.0.) THEN         IF (deltaTmom.NE.0.) THEN
56          recip_dt=1. _d 0/deltaTmom          recip_dt = 1. _d 0/deltaTmom
57         ELSE         ELSE
58          recip_dt=0. _d 0          recip_dt = 0. _d 0
59         ENDIF         ENDIF
60    
61         DO bj=myByLo(myThid), myByHi(myThid)         DO bj=myByLo(myThid), myByHi(myThid)
62          DO bi=myBxLo(myThid), myBxHi(myThid)          DO bi=myBxLo(myThid), myBxHi(myThid)
63           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
64            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
65             IF (useAreaViscLength) THEN              L2_D(i,j,bi,bj) = rA(i,j,bi,bj)
66              L2_D(i,j,bi,bj)=rA(i,j,bi,bj)            ENDDO
67             ELSE           ENDDO
68              L2_D(i,j,bi,bj)=2. _d 0           IF ( .NOT.useAreaViscLength ) THEN
69       &           /((recip_DXF(i,j,bi,bj)**2+recip_DYF(i,j,bi,bj)**2))            DO j=1-Oly,sNy+Oly
70             ENDIF             DO i=1-Olx,sNx+Olx
71             L3_D(i,j,bi,bj)=L2_D(i,j,bi,bj)**1.5              IF ( recip_dxF(i,j,bi,bj).NE.0. .OR.
72             L4rdt_D(i,j,bi,bj)=0.03125 _d 0*recip_dt*L2_D(i,j,bi,bj)**2       &           recip_dyF(i,j,bi,bj).NE.0. ) THEN
73                  L2_D(i,j,bi,bj) = 2. _d 0
74         &           /((recip_dxF(i,j,bi,bj)**2+recip_dyF(i,j,bi,bj)**2))
75                ENDIF
76               ENDDO
77              ENDDO
78             ENDIF
79             DO j=1-Oly,sNy+Oly
80              DO i=1-Olx,sNx+Olx
81                L3_D(i,j,bi,bj) = L2_D(i,j,bi,bj)**1.5
82                L4rdt_D(i,j,bi,bj) = 0.03125 _d 0*recip_dt
83         &                       *L2_D(i,j,bi,bj)**2
84            ENDDO            ENDDO
85           ENDDO           ENDDO
86          ENDDO          ENDDO
# Line 79  C-    used for side-drag: Line 90  C-    used for side-drag:
90          DO bi=myBxLo(myThid), myBxHi(myThid)          DO bi=myBxLo(myThid), myBxHi(myThid)
91           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
92            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
93             IF (useAreaViscLength) THEN              L2_Z(i,j,bi,bj) = rAz(i,j,bi,bj)
94              L2_Z(i,j,bi,bj)=rAz(i,j,bi,bj)            ENDDO
95             ELSE           ENDDO
96              L2_Z(i,j,bi,bj)=2. _d 0           IF ( .NOT.useAreaViscLength ) THEN
97       &           /((recip_DXV(I,J,bi,bj)**2+recip_DYU(I,J,bi,bj)**2))            DO j=1-Oly,sNy+Oly
98             ENDIF             DO i=1-Olx,sNx+Olx
99             L3_Z(i,j,bi,bj)=L2_Z(i,j,bi,bj)**1.5              IF ( recip_dxV(i,j,bi,bj).NE.0. .OR.
100             L4rdt_Z(i,j,bi,bj)=0.03125 _d 0*recip_dt*L2_Z(i,j,bi,bj)**2       &           recip_dyU(i,j,bi,bj).NE.0. ) THEN
101                  L2_Z(i,j,bi,bj) = 2. _d 0
102         &           /((recip_dxV(i,j,bi,bj)**2+recip_dyU(i,j,bi,bj)**2))
103                ENDIF
104               ENDDO
105              ENDDO
106             ENDIF
107             DO j=1-Oly,sNy+Oly
108              DO i=1-Olx,sNx+Olx
109                L3_Z(i,j,bi,bj) = L2_Z(i,j,bi,bj)**1.5
110                L4rdt_Z(i,j,bi,bj) = 0.03125 _d 0*recip_dt
111         &                       *L2_Z(i,j,bi,bj)**2
112            ENDDO            ENDDO
113           ENDDO           ENDDO
114          ENDDO          ENDDO
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22