/[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.8 by cnh, Tue Sep 29 18:50:57 1998 UTC revision 1.9 by cnh, Wed Oct 28 03:11:37 1998 UTC
# Line 59  C     Calculate lopping factor hFacC Line 59  C     Calculate lopping factor hFacC
59            DO I=1,sNx            DO I=1,sNx
60  C          Round depths within a small fraction of layer depth to that  C          Round depths within a small fraction of layer depth to that
61  C          layer depth.  C          layer depth.
62             IF ( ABS(H(I,J,bi,bj)-rF(K)) .LT. 1. _d -6*ABS(rF(K)) .AND.             IF ( ABS(H(I,J,bi,bj)-rF(K)) .LT.
63       &          ABS(H(I,J,bi,bj)-rF(K)) .LT. 1. _d -6*ABS(H(I,J,bi,bj)) )THEN       &          1. _d -6*ABS(rF(K)) .AND.
64         &          ABS(H(I,J,bi,bj)-rF(K)) .LT.
65         &          1. _d -6*ABS(H(I,J,bi,bj)) )THEN
66              H(I,J,bi,bj) = rF(K)              H(I,J,bi,bj) = rF(K)
67             ENDIF             ENDIF
68             IF     ( H(I,J,bi,bj)*rkFac .GE. rF(K)*rkFac ) THEN             IF     ( H(I,J,bi,bj)*rkFac .GE. rF(K)*rkFac ) THEN
# Line 94  C Impose minimum size (dimensional) Line 96  C Impose minimum size (dimensional)
96               hFacC(I,J,K,bi,bj)=hFacMinDr*recip_drF(k)               hFacC(I,J,K,bi,bj)=hFacMinDr*recip_drF(k)
97              ENDIF              ENDIF
98             ENDIF             ENDIF
99             depthInK(i,j,bi,bj) = depthInK(i,j,bi,bj) + hFacC(i,j,k,bi,bj)             depthInK(i,j,bi,bj) = depthInK(i,j,bi,bj)
100         &                          +hFacC(i,j,k,bi,bj)
101            ENDDO            ENDDO
102           ENDDO           ENDDO
103          ENDDO          ENDDO
# Line 106  C Impose minimum size (dimensional) Line 109  C Impose minimum size (dimensional)
109        _EXCH_XY_R4( depthInK, myThid )        _EXCH_XY_R4( depthInK, myThid )
110    
111        IF ( myThid .EQ. 1 ) WRITE(0,*) 'DDDD'        IF ( myThid .EQ. 1 ) WRITE(0,*) 'DDDD'
112        CALL PLOT_FIELD_XYRS( depthInK, 'Model Depths K Index' , 1, myThid )        CALL PLOT_FIELD_XYRS( depthInK,
113         & 'Model Depths K Index' , 1, myThid )
114    
115  C     hFacW and hFacS (at U and V points)  C     hFacW and hFacS (at U and V points)
116        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22