/[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.44 by jmc, Thu Dec 17 23:47:06 2009 UTC revision 1.45 by jmc, Sat Jan 16 23:00:16 2010 UTC
# Line 137  C       Note: if no fluid (continent) == Line 137  C       Note: if no fluid (continent) ==
137           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
138            tmpfld(I,J,bi,bj) = 0.            tmpfld(I,J,bi,bj) = 0.
139            ksurfC(I,J,bi,bj) = Nr+1            ksurfC(I,J,bi,bj) = Nr+1
140            maskH(I,J,bi,bj)  = 0.  c         maskH(I,J,bi,bj)  = 0.
141            Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj)            Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj)
142            DO K=Nr,1,-1            DO K=Nr,1,-1
143             Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj)             Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj)
144       &                        + drF(k)*hFacC(I,J,K,bi,bj)       &                        + drF(k)*hFacC(I,J,K,bi,bj)
145             IF (hFacC(I,J,K,bi,bj).NE.0.) THEN             IF (hFacC(I,J,K,bi,bj).NE.0.) THEN
146              ksurfC(I,J,bi,bj) = k              ksurfC(I,J,bi,bj) = k
147              maskH(I,J,bi,bj)  = 1.  c           maskH(I,J,bi,bj)  = 1.
148              tmpfld(I,J,bi,bj) = tmpfld(I,J,bi,bj) + 1.              tmpfld(I,J,bi,bj) = tmpfld(I,J,bi,bj) + 1.
149             ENDIF             ENDIF
150            ENDDO            ENDDO
Legend:
Removed from v.1.44  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22