/[MITgcm]/MITgcm/pkg/my82/my82_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/my82/my82_calc.F

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

revision 1.3 by jmc, Wed Apr 6 18:43:11 2005 UTC revision 1.4 by mlosch, Mon May 30 07:41:45 2005 UTC
# Line 80  CEOP Line 80  CEOP
80  C     Initialize local fields  C     Initialize local fields
81        DO J=1-Oly,sNy+Oly        DO J=1-Oly,sNy+Oly
82         DO I=1-Olx,sNx+Olx         DO I=1-Olx,sNx+Olx
83          GH(I,J) = 0.          GH(I,J) = 0. _d 0
84          GM(I,J) = 0.          GM(I,J) = 0. _d 0
85         ENDDO         ENDDO
86        ENDDO        ENDDO
87        DO K = 1, Nr        DO K = 1, Nr
88         DO J=1-Oly,sNy+Oly         DO J=1-Oly,sNy+Oly
89          DO I=1-Olx,sNx+Olx          DO I=1-Olx,sNx+Olx
90           SH(I,J,K) = 0.           SH(I,J,K)  = 0. _d 0
91           SM(I,J,K) = 0.           SM(I,J,K)  = 0. _d 0
92           tke(I,J,K) = 0.           tke(I,J,K) = 0. _d 0
93          ENDDO          ENDDO
94         ENDDO             ENDDO    
95        ENDDO        ENDDO
# Line 129  C     re-initilialize GM and GH for abus Line 129  C     re-initilialize GM and GH for abus
129  C     the meaning of shear and negative  buoyancy frequency  C     the meaning of shear and negative  buoyancy frequency
130        DO J=jMin,jMax        DO J=jMin,jMax
131         DO I=iMin,iMax         DO I=iMin,iMax
132          GH(I,J) = 0.          GH(I,J) = 0. _d 0
133          GM(I,J) = 0.          GM(I,J) = 0. _d 0
134         ENDDO         ENDDO
135        ENDDO        ENDDO
136  C     Find boundary length scale from energy weighted mean.  C     Find boundary length scale from energy weighted mean.
# Line 149  C     end of second k-loop Line 149  C     end of second k-loop
149  C     compute boundary length scale MYhbl  C     compute boundary length scale MYhbl
150        DO J=jMin,jMax        DO J=jMin,jMax
151         DO I=iMin,iMax         DO I=iMin,iMax
152          IF ( GH(I,J) .EQ. 0. ) THEN          IF ( GH(I,J) .EQ. 0. _d 0 ) THEN
153           MYhbl(I,J,bi,bj) = 0.           MYhbl(I,J,bi,bj) = 0. _d 0
154          ELSE          ELSE
155           MYhbl(I,J,bi,bj) = -GM(I,J)/GH(I,J)*MYhblScale           MYhbl(I,J,bi,bj) = -GM(I,J)/GH(I,J)*MYhblScale
156          ENDIF          ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22