/[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.2 by adcroft, Thu Jul 2 15:46:21 1998 UTC revision 1.3 by adcroft, Wed Jul 29 18:33:47 1998 UTC
# Line 41  C         Inverse of depth Line 41  C         Inverse of depth
41            IF ( h(i,j,bi,bj) .EQ. 0. _d 0 ) THEN            IF ( h(i,j,bi,bj) .EQ. 0. _d 0 ) THEN
42             rH(i,j,bi,bj) = 0. _d 0             rH(i,j,bi,bj) = 0. _d 0
43            ELSE            ELSE
44             rH(i,j,bi,bj) = 1. _d 0 /  H(i,j,bi,bj)             rH(i,j,bi,bj) = 1. _d 0 /  abs( H(i,j,bi,bj) )
45            ENDIF            ENDIF
46           ENDDO           ENDDO
47          ENDDO          ENDDO
# Line 55  C     Calculate lopping factor hFacC Line 55  C     Calculate lopping factor hFacC
55          DO K=1, Nz          DO K=1, Nz
56           DO J=1,sNy           DO J=1,sNy
57            DO I=1,sNx            DO I=1,sNx
58             IF     ( H(I,J,bi,bj) .LE. zFace(K) ) THEN             IF     ( H(I,J,bi,bj) .GE. zFace(K) ) THEN
59  C           Below base of domain  C           Top of cell is below base of domain
60              hFacC(I,J,K,bi,bj) = 0.              hFacC(I,J,K,bi,bj) = 0.
61             ELSEIF ( H(I,J,bi,bj) .GT. zFace(K+1) ) THEN             ELSEIF ( H(I,J,bi,bj) .LE. zFace(K+1) ) THEN
62  C           Base of domain is below this cell  C           Base of domain is below bottom of this cell
63              hFacC(I,J,K,bi,bj) = 1.              hFacC(I,J,K,bi,bj) = 1.
64             ELSE             ELSE
65  C           Base of domain is in this cell  C           Base of domain is in this cell
66  C           Set hFac tp the fraction of the cell that is open.  C           Set hFac to the fraction of the cell that is open.
67              hFacC(I,J,K,bi,bj) =              hFacC(I,J,K,bi,bj) = (zFace(K)-H(I,J,bi,bj))*rdzF(K)
68       &        (zFace(K)-H(I,J,bi,bj))/(zFace(K)-zFace(K+1))             ENDIF
69    C Impose minimum fraction
70               IF (hFacC(I,J,K,bi,bj).LT.hFacMin) THEN
71                IF (hFacC(I,J,K,bi,bj).LT.hFacMin*0.5) THEN
72                 hFacC(I,J,K,bi,bj)=0.
73                ELSE
74                 hFacC(I,J,K,bi,bj)=hFacMin
75                ENDIF
76               ENDIF
77    C Impose minimum size (dimensional)
78               IF (dzF(k)*hFacC(I,J,K,bi,bj).LT.hFacMinDz) THEN
79                IF (dzF(k)*hFacC(I,J,K,bi,bj).LT.hFacMinDz*0.5) THEN
80                 hFacC(I,J,K,bi,bj)=0.
81                ELSE
82                 hFacC(I,J,K,bi,bj)=hFacMinDz*rDzF(k)
83                ENDIF
84             ENDIF             ENDIF
85            ENDDO            ENDDO
86           ENDDO           ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22