/[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.26 by cnh, Wed Sep 26 18:09:15 2001 UTC revision 1.27 by jmc, Fri Nov 22 03:01:18 2002 UTC
# Line 35  C     === Global variables === Line 35  C     === Global variables ===
35  #include "EEPARAMS.h"  #include "EEPARAMS.h"
36  #include "PARAMS.h"  #include "PARAMS.h"
37  #include "GRID.h"  #include "GRID.h"
38    #include "SURFACE.h"
39    
40  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
41  C     == Routine arguments ==  C     == Routine arguments ==
# Line 43  C     myThid -  Number of this instance Line 44  C     myThid -  Number of this instance
44  CEndOfInterface  CEndOfInterface
45    
46  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
 C     == Local variables in common ==  
 C     Hloc  - Temporary array used to read surface topography  
 C             has to be in common for multi threading      
       COMMON / LOCAL_INI_DEPTHS / Hloc  
       _RS Hloc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
47  C     == Local variables ==  C     == Local variables ==
48  C     iG, jG - Global coordinate index  C     iG, jG - Global coordinate index
49  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
# Line 78  C------ Line 74  C------
74           DO i=1-Olx,sNx+Olx           DO i=1-Olx,sNx+Olx
75            R_low(i,j,bi,bj) = 0.            R_low(i,j,bi,bj) = 0.
76            Ro_surf(i,j,bi,bj) = 0.            Ro_surf(i,j,bi,bj) = 0.
77              topoZ(i,j,bi,bj) = 0.
78           ENDDO           ENDDO
79          ENDDO          ENDDO
80         ENDDO         ENDDO
# Line 154  C------ read from file: Line 151  C------ read from file:
151    
152  C- read surface topography (in m) from topoFile (case topoFile.NE.' '):  C- read surface topography (in m) from topoFile (case topoFile.NE.' '):
153          _BEGIN_MASTER( myThid )          _BEGIN_MASTER( myThid )
154          CALL READ_REC_XY_RS( topoFile, Hloc, 1, 0, myThid )          CALL READ_REC_XY_RS( topoFile, topoZ, 1, 0, myThid )
155          _END_MASTER(myThid)          _END_MASTER(myThid)
156          _BARRIER          _BARRIER
157    
# Line 164  C   Convert Surface Geopotential to (ref Line 161  C   Convert Surface Geopotential to (ref
161  C   according to Tref profile, using same discretisation as in calc_phi_hyd  C   according to Tref profile, using same discretisation as in calc_phi_hyd
162  C----  C----
163  c         _BEGIN_MASTER( myThid )  c         _BEGIN_MASTER( myThid )
164  c         CALL WRITE_FLD_XY_RS( 'topo_Z',' ',Hloc,0,myThid)  c         CALL WRITE_FLD_XY_RS( 'topo_Z',' ',topoZ,0,myThid)
165  c         _END_MASTER(myThid)  c         _END_MASTER(myThid)
166    
167            CALL INI_P_GROUND( Hloc, Ro_surf, myThid )            CALL INI_P_GROUND( topoZ, Ro_surf, myThid )
168    
169            _BARRIER            _BARRIER
170            _BEGIN_MASTER( myThid )            _BEGIN_MASTER( myThid )
# Line 181  C   Direct Transfer to Ro_surf : Line 178  C   Direct Transfer to Ro_surf :
178             DO bi = myBxLo(myThid), myBxHi(myThid)             DO bi = myBxLo(myThid), myBxHi(myThid)
179              DO j=1,sNy              DO j=1,sNy
180               DO i=1,sNx               DO i=1,sNx
181                Ro_surf(i,j,bi,bj) = Hloc(i,j,bi,bj)                Ro_surf(i,j,bi,bj) = topoZ(i,j,bi,bj)
182               ENDDO               ENDDO
183              ENDDO              ENDDO
184             ENDDO             ENDDO
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22