/[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.9 by cnh, Wed Oct 28 03:11:37 1998 UTC revision 1.10 by cnh, Mon Nov 2 03:34:12 1998 UTC
# Line 136  C     Masks and reciprocals of hFac[CWS] Line 136  C     Masks and reciprocals of hFac[CWS]
136          DO K=1,Nr          DO K=1,Nr
137           DO J=1,sNy           DO J=1,sNy
138            DO I=1,sNx            DO I=1,sNx
139             IF (HFacC(I,J,K,bi,bj) .NE. 0. D0 ) THEN             IF (HFacC(I,J,K,bi,bj) .NE. 0. _d 0 ) THEN
140              recip_HFacC(I,J,K,bi,bj) = 1. D0 / HFacC(I,J,K,bi,bj)              recip_HFacC(I,J,K,bi,bj) = 1. _d 0 / HFacC(I,J,K,bi,bj)
141             ELSE             ELSE
142              recip_HFacC(I,J,K,bi,bj) = 0. D0              recip_HFacC(I,J,K,bi,bj) = 0. _d 0
143             ENDIF             ENDIF
144             IF (HFacW(I,J,K,bi,bj) .NE. 0. D0 ) THEN             IF (HFacW(I,J,K,bi,bj) .NE. 0. _d 0 ) THEN
145              recip_HFacW(I,J,K,bi,bj) = 1. D0 / HFacW(I,J,K,bi,bj)              recip_HFacW(I,J,K,bi,bj) = 1. _d 0 / HFacW(I,J,K,bi,bj)
146              maskW(I,J,K,bi,bj) = 1. D0              maskW(I,J,K,bi,bj) = 1. _d 0
147             ELSE             ELSE
148              recip_HFacW(I,J,K,bi,bj) = 0. D0              recip_HFacW(I,J,K,bi,bj) = 0. _d 0
149              maskW(I,J,K,bi,bj) = 0.0 D0              maskW(I,J,K,bi,bj) = 0.0 _d 0
150             ENDIF             ENDIF
151             IF (HFacS(I,J,K,bi,bj) .NE. 0. D0 ) THEN             IF (HFacS(I,J,K,bi,bj) .NE. 0. _d 0 ) THEN
152              recip_HFacS(I,J,K,bi,bj) = 1. D0 / HFacS(I,J,K,bi,bj)              recip_HFacS(I,J,K,bi,bj) = 1. _d 0 / HFacS(I,J,K,bi,bj)
153              maskS(I,J,K,bi,bj) = 1. D0              maskS(I,J,K,bi,bj) = 1. _d 0
154             ELSE             ELSE
155              recip_HFacS(I,J,K,bi,bj) = 0. D0              recip_HFacS(I,J,K,bi,bj) = 0. _d 0
156              maskS(I,J,K,bi,bj) = 0. D0              maskS(I,J,K,bi,bj) = 0. _d 0
157             ENDIF             ENDIF
158            ENDDO            ENDDO
159           ENDDO           ENDDO
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22