/[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.14 by cnh, Tue Sep 29 18:50:57 1998 UTC revision 1.19 by adcroft, Wed May 5 18:32:34 1999 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_DEPTHS( myThid )        SUBROUTINE INI_DEPTHS( myThid )
# Line 11  C     |================================= Line 11  C     |=================================
11  C     | The depths of the bottom of the model is specified in    |  C     | The depths of the bottom of the model is specified in    |
12  C     | terms of an XY map with one depth for each column of     |  C     | terms of an XY map with one depth for each column of     |
13  C     | grid cells. Depths do not have to coincide with the      |  C     | grid cells. Depths do not have to coincide with the      |
14  C     | model levels. The model's lopping algorithm makes it     |  C     | model levels. The model lopping algorithm makes it       |
15  C     | possible to represent arbitrary depths.                  |  C     | possible to represent arbitrary depths.                  |
16  C     | The mode depths map also influences the models topology  |  C     | The mode depths map also influences the models topology  |
17  C     | By default the model domain wraps around in X and Y.     |  C     | By default the model domain wraps around in X and Y.     |
# Line 19  C     | This default doubly periodic top Line 19  C     | This default doubly periodic top
19  C     | if a depth map is defined which closes off all wrap      |  C     | if a depth map is defined which closes off all wrap      |
20  C     | around flow.                                             |  C     | around flow.                                             |
21  C     \==========================================================/  C     \==========================================================/
22          IMPLICIT NONE
23    
24  C     === Global variables ===  C     === Global variables ===
25  #include "SIZE.h"  #include "SIZE.h"
# Line 40  C     oldPrec - Temporary used in contro Line 41  C     oldPrec - Temporary used in contro
41        INTEGER iG, jG        INTEGER iG, jG
42        INTEGER bi, bj        INTEGER bi, bj
43        INTEGER  I,  J, K        INTEGER  I,  J, K
       INTEGER oldPrec  
44        _RL phi        _RL phi
45    
46        _BARRIER        _BARRIER
# Line 59  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.14  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22