/[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.17 by adcroft, Wed Dec 9 16:11:52 1998 UTC revision 1.20 by adcroft, Mon Mar 27 22:25:44 2000 UTC
# Line 35  CEndOfInterface Line 35  CEndOfInterface
35  C     == Local variables ==  C     == Local variables ==
36  C     iG, jG - Global coordinate index  C     iG, jG - Global coordinate index
37  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
38  C     I,J,K  C     I,J
39  C     phi    - total depth of model  C     phi    - total depth of model
40  C     oldPrec - Temporary used in controlling binary input dataset precision  C     oldPrec - Temporary used in controlling binary input dataset precision
41        INTEGER iG, jG        INTEGER iG, jG
42        INTEGER bi, bj        INTEGER bi, bj
43        INTEGER  I,  J, K        INTEGER  I,  J
       INTEGER oldPrec  
44        _RL phi        _RL phi
45    
46        _BARRIER        _BARRIER
# Line 60  C      realistic problem. Line 59  C      realistic problem.
59  C          Default depth of full domain  C          Default depth of full domain
60             H(i,j,bi,bj) = phi             H(i,j,bi,bj) = phi
61  C          Test for eastern edge  C          Test for eastern edge
62             IF ( iG .EQ. nX ) H(i,j,bi,bj) = 0.  C          IF ( iG .EQ. nX ) H(i,j,bi,bj) = 0.
63  C          Test for northern edge  C          Test for northern edge
64             IF ( jG .EQ. nY ) H(i,j,bi,bj) = 0.  C          IF ( jG .EQ. nY ) H(i,j,bi,bj) = 0.
 C          Island  
            IF ( iG .EQ. 1 .AND.  
      &          jG .EQ. 24 ) H(i,j,bi,bj) = 0.75*phi  
65            ENDDO            ENDDO
66           ENDDO           ENDDO
67          ENDDO          ENDDO
68         ENDDO         ENDDO
69        ELSE        ELSE
70         _BEGIN_MASTER( myThid )         _BEGIN_MASTER( myThid )
71    C Read the bathymetry using the mid-level I/O pacakage read_write_rec
72  CcnhDebugStarts  C The 0 is the "iteration" argument. The 1 is the record number.
73  C       Force 64-bit IO         CALL READ_REC_XY_RS( bathyFile, H, 1, 0, myThid )
74          oldPrec        = readBinaryPrec  C Read the bathymetry using the mid-level I/O pacakage read_write_fld
75          readBinaryPrec = precFloat64  C The 0 is the "iteration" argument. The ' ' is an empty suffix
76  CcnhDEbugEnds  C      CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )
77         CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )  C Read the bathymetry using the low-level I/O package
78  CcnhDebugStarts  C      CALL MDSREADFIELD( bathyFile, readBinaryPrec,
79         readBinaryPrec = oldPrec  C    &                    'RS', 1, H, 1, myThid )
 CcnhdDebugEnds  
   
80         _END_MASTER(myThid)         _END_MASTER(myThid)
81        ENDIF        ENDIF
82                
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22