/[MITgcm]/MITgcm/model/src/ini_masks_etc.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_masks_etc.F

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

revision 1.22 by jmc, Thu Jul 5 21:44:25 2001 UTC revision 1.23 by jmc, Mon Aug 27 18:42:37 2001 UTC
# Line 117  C      o Impose minimum fraction and/or Line 117  C      o Impose minimum fraction and/or
117    
118  C-  Re-calculate Reference surface position, taking into account hFacC  C-  Re-calculate Reference surface position, taking into account hFacC
119  C   initialize Total column fluid thickness and surface k index  C   initialize Total column fluid thickness and surface k index
120    C       Note: if no fluid (continent) ==> ksurf = Nr+1
121          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
122           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
123            tmpfld(I,J,bi,bj) = 0.            tmpfld(I,J,bi,bj) = 0.
124            k_surf(I,J,bi,bj) = Nr            ksurfC(I,J,bi,bj) = Nr+1
125            Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj)            Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj)
126            DO K=Nr,1,-1            DO K=Nr,1,-1
127             Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj)             Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj)
128       &                        + drF(k)*hFacC(I,J,K,bi,bj)       &                        + drF(k)*hFacC(I,J,K,bi,bj)
129             IF (hFacC(I,J,K,bi,bj).NE.0.) THEN             IF (hFacC(I,J,K,bi,bj).NE.0.) THEN
130              k_surf(I,J,bi,bj) = k              ksurfC(I,J,bi,bj) = k
131              tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1.              tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1.
132             ENDIF             ENDIF
133            ENDDO            ENDDO
# Line 240  C     Masks and reciprocals of hFac[CWS] Line 241  C     Masks and reciprocals of hFac[CWS]
241            ENDDO            ENDDO
242           ENDDO           ENDDO
243          ENDDO          ENDDO
244    C-    Calculate surface k index for interface W & S (U & V points)
245            DO J=1-Oly,sNy+Oly
246             DO I=1-Olx,sNx+Olx
247              ksurfW(I,J,bi,bj) = Nr+1
248              ksurfS(I,J,bi,bj) = Nr+1
249              DO k=Nr,1,-1
250               IF (hFacW(I,J,K,bi,bj).NE.0.) ksurfW(I,J,bi,bj) = k
251               IF (hFacS(I,J,K,bi,bj).NE.0.) ksurfS(I,J,bi,bj) = k
252              ENDDO
253             ENDDO
254            ENDDO
255    C - end bi,bj loops.
256         ENDDO         ENDDO
257        ENDDO        ENDDO
258  C     _EXCH_XYZ_R4(recip_HFacC    , myThid )  C     _EXCH_XYZ_R4(recip_HFacC    , myThid )
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22