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

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

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

revision 1.7 by cnh, Mon Jun 15 05:13:56 1998 UTC revision 1.8 by adcroft, Thu Jul 2 14:16:24 1998 UTC
# Line 85  C          Island Line 85  C          Island
85        ENDIF        ENDIF
86        _EXCH_XY_R4(    H, myThid )        _EXCH_XY_R4(    H, myThid )
87    
 C--   Cacluate 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  
            rH(i,j,bi,bj) = 0. _d 0  
           ELSE  
            rH(i,j,bi,bj) = 1. _d 0 /  h(i,j,bi,bj)  
           ENDIF  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       _EXCH_XY_R4(   rH, myThid )  
   
88  C--   Now calculate "lopping" factors hFac.  C--   Now calculate "lopping" factors hFac.
89        zG = delZ(1)*0.5 D0        zG = delZ(1)*0.5 D0
90        zFace(1) = 0        zFace(1) = 0
# Line 148  C           Set hFac tp the fraction of Line 131  C           Set hFac tp the fraction of
131        ENDDO        ENDDO
132        _EXCH_XYZ_R4(hFacC , myThid )        _EXCH_XYZ_R4(hFacC , myThid )
133    
       DO bj=myByLo(myThid), myByHi(myThid)  
        DO bi=myBxLo(myThid), myBxHi(myThid)  
         DO K=1, Nz  
          DO J=1,sNy  
           DO I=1,sNx  
            hFacW(I,J,K,bi,bj)=  
      &       MIN(hFacC(I,J,K,bi,bj),hFacC(I-1,J,K,bi,bj))  
            hFacS(I,J,K,bi,bj)=  
      &       MIN(hFacC(I,J,K,bi,bj),hFacC(I,J-1,K,bi,bj))  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       _EXCH_XYZ_R4(hFacW , myThid )  
       _EXCH_XYZ_R4(hFacS , myThid )  
   
 C--   Calculate recipricols of hFacC, hFacW and hFacS  
       DO bj=myByLo(myThid), myByHi(myThid)  
        DO bi=myBxLo(myThid), myBxHi(myThid)  
         DO K=1, Nz  
          DO J=1,sNy  
           DO I=1,sNx  
            rhFacC(I,J,K,bi,bj)=0. D0  
            if (hFacC(I,J,K,bi,bj).ne.0.)  
      &            rhFacC(I,J,K,bi,bj)=1. D0 /hFacC(I,J,K,bi,bj)  
            rhFacW(I,J,K,bi,bj)=0. D0  
            if (hFacW(I,J,K,bi,bj).ne.0.)  
      &            rhFacW(I,J,K,bi,bj)=1. D0 /hFacW(I,J,K,bi,bj)  
            rhFacS(I,J,K,bi,bj)=0. D0  
            if (hFacS(I,J,K,bi,bj).ne.0.)  
      &            rhFacS(I,J,K,bi,bj)=1. D0 /hFacS(I,J,K,bi,bj)  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
       _EXCH_XYZ_R4(rhFacC , myThid )  
       _EXCH_XYZ_R4(rhFacW , myThid )  
       _EXCH_XYZ_R4(rhFacS , myThid )  
   
134  C  C
135        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )
136  C  C
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22