/[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.12 by cnh, Sat Aug 22 17:51:08 1998 UTC revision 1.17 by adcroft, Wed Dec 9 16:11:52 1998 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 36  C     iG, jG - Global coordinate index Line 37  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,K
39  C     phi    - total depth of model  C     phi    - total depth of model
40    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, K
44          INTEGER oldPrec
45        _RL phi        _RL phi
46    
47        _BARRIER        _BARRIER
# Line 69  C          Island Line 72  C          Island
72         ENDDO         ENDDO
73        ELSE        ELSE
74         _BEGIN_MASTER( myThid )         _BEGIN_MASTER( myThid )
75    
76    CcnhDebugStarts
77    C       Force 64-bit IO
78            oldPrec        = readBinaryPrec
79            readBinaryPrec = precFloat64
80    CcnhDEbugEnds
81         CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )         CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )
82    CcnhDebugStarts
83           readBinaryPrec = oldPrec
84    CcnhdDebugEnds
85    
86         _END_MASTER(myThid)         _END_MASTER(myThid)
87        ENDIF        ENDIF
88          
89        _EXCH_XY_R4(    H, myThid )        _EXCH_XY_R4(    H, myThid )
90  C  C
91        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22