/[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.34 by heimbach, Wed Jan 4 12:25:40 2006 UTC revision 1.35 by mlosch, Tue Feb 7 11:47:48 2006 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CBOP  CBOP
# Line 36  C     === Global variables === Line 37  C     === Global variables ===
37  #include "PARAMS.h"  #include "PARAMS.h"
38  #include "GRID.h"  #include "GRID.h"
39  #include "SURFACE.h"  #include "SURFACE.h"
40    #ifdef ALLOW_SHELFICE
41    # include "SHELFICE.h"
42    #endif /* ALLOW_SHELFICE */
43    
44  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
45  C     == Routine arguments ==  C     == Routine arguments ==
# Line 243  c     _BEGIN_MASTER( myThid ) Line 247  c     _BEGIN_MASTER( myThid )
247  c     CALL WRITE_FLD_XY_RS('Ro_surf',' ',Ro_surf,0,myThid)  c     CALL WRITE_FLD_XY_RS('Ro_surf',' ',Ro_surf,0,myThid)
248  c     _END_MASTER(myThid)  c     _END_MASTER(myThid)
249    
250    #ifdef ALLOW_SHELFICE
251    c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
252          IF ( useShelfIce ) THEN
253    C------
254    C   4) Set R_shelfIce = the Lower (in r sense) boundary of floating shelfice :
255    C------
256          IF (usingPCoords .OR. shelfIceFile .EQ. ' ') THEN
257    C- e.g., atmosphere : R_low = Top of atmosphere
258    C-            ocean : R_low = Bottom
259           DO bj = myByLo(myThid), myByHi(myThid)
260            DO bi = myBxLo(myThid), myBxHi(myThid)
261             DO j=1,sNy
262              DO i=1,sNx
263               R_shelfIce(i,j,bi,bj) = 0. _d 0
264              ENDDO
265             ENDDO
266            ENDDO
267           ENDDO
268          ELSE
269            _BEGIN_MASTER( myThid )
270    C Read the shelfIce draught using the mid-level I/O pacakage read_write_rec
271    C The 0 is the "iteration" argument. The 1 is the record number.
272            CALL READ_REC_XY_RS( shelfIceFile, R_shelfIce, 1, 0, myThid )
273    C Read the shelfIce draught using the mid-level I/O pacakage read_write_fld
274    C The 0 is the "iteration" argument. The ' ' is an empty suffix
275    C        CALL READ_FLD_XY_RS( shelfIceFile, ' ', R_shelfIce, 0, myThid )
276    c       CALL READ_FLD_XY_RS( bathyFile, ' ', R_low, 0, myThid )
277    C Read the selfIce draught using the low-level I/O package
278    c       CALL MDSREADFIELD( shelfIceFile, readBinaryPrec,
279    c    &                     'RS', 1, R_selfIce, 1, myThid )
280            _END_MASTER(myThid)
281    
282          ENDIF
283    C- end setup R_shelfIce in the interior
284    
285    C- fill in the overlap :
286          _EXCH_XY_R4(R_shelfIce, myThid )
287    
288    c     CALL PLOT_FIELD_XYRS(R_selfIce,'Shelf ice draught (ini_depths)',
289    c    &     1,myThid)
290    CML      _BEGIN_MASTER( myThid )
291    CML      CALL WRITE_FLD_XY_RS( 'R_shelfIce' ,' ', R_shelfIce, 0,myThid)
292    CML      _END_MASTER(myThid)
293          ENDIF
294    #endif /* ALLOW_SHELFICE */
295    
296        RETURN        RETURN
297        END        END
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22