/[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.10 by adcroft, Wed Jul 15 22:20:02 1998 UTC revision 1.16 by cnh, Fri Nov 6 22:44:47 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 36  C     iG, jG - Global coordinate index Line 36  C     iG, jG - Global coordinate index
36  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
37  C     I,J,K  C     I,J,K
38  C     phi    - total depth of model  C     phi    - total depth of model
39    C     oldPrec - Temporary used in controlling binary input dataset precision
40        INTEGER iG, jG        INTEGER iG, jG
41        INTEGER bi, bj        INTEGER bi, bj
42        INTEGER  I,  J, K        INTEGER  I,  J, K
43          INTEGER oldPrec
44        _RL phi        _RL phi
45    
46        _BARRIER        _BARRIER
# Line 47  C      Set up a flat bottom box with dou Line 49  C      Set up a flat bottom box with dou
49  C      H is the basic variable from which other terms are derived. It  C      H is the basic variable from which other terms are derived. It
50  C      is the term that would be set from an external file for a  C      is the term that would be set from an external file for a
51  C      realistic problem.  C      realistic problem.
52         phi = zFace(nZ+1)         phi = rF(Nr+1)
53         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
54          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
55           DO j=1,sNy           DO j=1,sNy
# Line 69  C          Island Line 71  C          Island
71         ENDDO         ENDDO
72        ELSE        ELSE
73         _BEGIN_MASTER( myThid )         _BEGIN_MASTER( myThid )
74    
75    CcnhDebugStarts
76    C       Force 64-bit IO
77            oldPrec        = readBinaryPrec
78            readBinaryPrec = precFloat64
79    CcnhDEbugEnds
80         CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )         CALL READ_FLD_XY_RS( bathyFile, ' ', H, 0, myThid )
81    CcnhDebugStarts
82           readBinaryPrec = oldPrec
83    CcnhdDebugEnds
84    
85         _END_MASTER(myThid)         _END_MASTER(myThid)
 C--    Reverse sign if necessary (We assume H>0)  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            H(i,j,bi,bj)=-H(i,j,bi,bj)  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
86        ENDIF        ENDIF
87          
88        _EXCH_XY_R4(    H, myThid )        _EXCH_XY_R4(    H, myThid )
89  C  C
90        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )        CALL PLOT_FIELD_XYRS( H, 'Model depths' , 1, myThid )
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22