/[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.15 by adcroft, Mon Mar 22 15:54:04 1999 UTC revision 1.16 by adcroft, Wed May 5 18:32:34 1999 UTC
# Line 37  C     I,J,K Line 37  C     I,J,K
37        _RL hFacUpper,hFacLower        _RL hFacUpper,hFacLower
38  #endif  #endif
39    
 C     Calculate quantities derived from XY depth map  
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO J=1,sNy  
          DO I=1,sNx  
 C         Inverse of depth  
           IF ( h(i,j,bi,bj) .EQ. 0. _d 0 ) THEN  
            recip_H(i,j,bi,bj) = 0. _d 0  
           ELSE  
            recip_H(i,j,bi,bj) = 1. _d 0 /  abs( H(i,j,bi,bj) )  
           ENDIF  
           depthInK(i,j,bi,bj) = 0.  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       _EXCH_XY_R4(   recip_H, myThid )  
   
40  C     Calculate lopping factor hFacC  C     Calculate lopping factor hFacC
41        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
42         DO bi=myBxLo(myThid), myBxHi(myThid)         DO bi=myBxLo(myThid), myBxHi(myThid)
# Line 110  C Impose minimum size (dimensional) Line 92  C Impose minimum size (dimensional)
92        CALL PLOT_FIELD_XYRS( depthInK,        CALL PLOT_FIELD_XYRS( depthInK,
93       & 'Model Depths K Index' , 1, myThid )       & 'Model Depths K Index' , 1, myThid )
94    
95    C Re-calculate depth of ocean, taking into account hFacC
96          DO bj=myByLo(myThid), myByHi(myThid)
97           DO bi=myBxLo(myThid), myBxHi(myThid)
98            DO J=1,sNy
99             DO I=1,sNx
100              H(I,J,bi,bj)=0.
101              DO K=1,Nr
102               H(I,J,bi,bj)=H(I,J,bi,bj)-
103         &       rkFac*drF(k)*hFacC(I,J,K,bi,bj)
104              ENDDO
105             ENDDO
106            ENDDO
107           ENDDO
108          ENDDO
109          _EXCH_XY_R4(H , myThid )
110          CALL WRITE_FLD_XY_RS( 'Depth',' ',H,0,myThid)
111    C     CALL MDSWRITEFIELD( 'Depth', writeBinaryPrec, .TRUE.,
112    C    &                    'RS', 1, H, 1, -1, myThid )
113    
114    C     Calculate quantities derived from XY depth map
115          DO bj = myByLo(myThid), myByHi(myThid)
116           DO bi = myBxLo(myThid), myBxHi(myThid)
117            DO J=1,sNy
118             DO I=1,sNx
119    C         Inverse of depth
120              IF ( h(i,j,bi,bj) .EQ. 0. _d 0 ) THEN
121               recip_H(i,j,bi,bj) = 0. _d 0
122              ELSE
123               recip_H(i,j,bi,bj) = 1. _d 0 /  abs( H(i,j,bi,bj) )
124              ENDIF
125              depthInK(i,j,bi,bj) = 0.
126             ENDDO
127            ENDDO
128           ENDDO
129          ENDDO
130          _EXCH_XY_R4(   recip_H, myThid )
131    
132  C     hFacW and hFacS (at U and V points)  C     hFacW and hFacS (at U and V points)
133        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
134         DO bi=myBxLo(myThid), myBxHi(myThid)         DO bi=myBxLo(myThid), myBxHi(myThid)
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22