/[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.40 by mlosch, Thu Apr 24 08:23:00 2008 UTC revision 1.41 by mlosch, Wed Sep 10 09:19:22 2008 UTC
# Line 37  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"
 #ifdef ALLOW_SHELFICE  
 # include "SHELFICE.h"  
 #endif /* ALLOW_SHELFICE */  
40  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
41  # include "MNC_PARAMS.h"  # include "MNC_PARAMS.h"
42  #endif  #endif
# Line 82  C------ Line 79  C------
79            R_low(i,j,bi,bj)   = 0. _d 0            R_low(i,j,bi,bj)   = 0. _d 0
80            Ro_surf(i,j,bi,bj) = 0. _d 0            Ro_surf(i,j,bi,bj) = 0. _d 0
81            topoZ(i,j,bi,bj)   = 0. _d 0            topoZ(i,j,bi,bj)   = 0. _d 0
 #ifdef ALLOW_SHELFICE  
           R_shelfIce(i,j,bi,bj) = 0. _d 0  
 #endif /* ALLOW_SHELFICE */  
82           ENDDO           ENDDO
83          ENDDO          ENDDO
84         ENDDO         ENDDO
# Line 270  c     CALL WRITE_FLD_XY_RS('Ro_surf',' ' Line 264  c     CALL WRITE_FLD_XY_RS('Ro_surf',' '
264  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
265  c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
266        IF ( useShelfIce ) THEN        IF ( useShelfIce ) THEN
267  C------         CALL SHELFICE_INI_DEPTH( myThid )
 C   4) Set R_shelfIce = the Lower (in r sense) boundary of floating shelfice :  
 C------  
       IF (usingPCoords .OR. shelfIceFile .EQ. ' ') THEN  
 C- e.g., atmosphere : R_low = Top of atmosphere  
 C-            ocean : R_low = Bottom  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            R_shelfIce(i,j,bi,bj) = 0. _d 0  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ELSE  
 C Read the shelfIce draught using the mid-level I/O pacakage read_write_rec  
 C The 0 is the "iteration" argument. The 1 is the record number.  
         CALL READ_REC_XY_RS( shelfIceFile, R_shelfIce, 1, 0, myThid )  
 C Read the shelfIce draught using the mid-level I/O pacakage read_write_fld  
 C The 0 is the "iteration" argument. The ' ' is an empty suffix  
 C        CALL READ_FLD_XY_RS( shelfIceFile, ' ', R_shelfIce, 0, myThid )  
 c       CALL READ_FLD_XY_RS( bathyFile, ' ', R_low, 0, myThid )  
 C Read the selfIce draught using the low-level I/O package  
 c       CALL MDSREADFIELD( shelfIceFile, readBinaryPrec,  
 c    &                     'RS', 1, R_selfIce, 1, myThid )  
   
       ENDIF  
 C- end setup R_shelfIce in the interior  
   
 C- fill in the overlap (+ BARRIER):  
       _EXCH_XY_R4(R_shelfIce, myThid )  
   
 c     CALL PLOT_FIELD_XYRS(R_selfIce,'Shelf ice draught (ini_depths)',  
 c    &     1,myThid)  
 CML      CALL WRITE_FLD_XY_RS( 'R_shelfIce' ,' ', R_shelfIce, 0,myThid)  
268        ENDIF        ENDIF
269  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
270    
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22