pkg/steep_icecavity/stic_init_depths.F

00f81e6785 Ou W*0001 #include "STIC_OPTIONS.h"
                0002
                0003 CBOP
                0004 C     !ROUTINE: STIC_INIT_DEPTHS
                0005 C     !INTERFACE:
                0006       SUBROUTINE STIC_INIT_DEPTHS(
                0007      I                    myThid )
                0008 C     !DESCRIPTION: \bv
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE STIC_INIT_DEPTHS
                0011 C     | o Modify ocean upper boundary position according to
                0012 C     |   ice-shelf topography
                0013 C     *==========================================================*
                0014 C     \ev
                0015
                0016 C     !USES:
                0017       IMPLICIT NONE
                0018 C     === Global variables ===
                0019 #include "SIZE.h"
                0020 #include "EEPARAMS.h"
                0021 #include "PARAMS.h"
                0022 #include "STIC.h"
                0023
                0024 C     !INPUT/OUTPUT PARAMETERS:
                0025 C     == Routine arguments ==
                0026 C     myThid    :: my Thread Id number
                0027       INTEGER myThid
                0028
                0029 C     !LOCAL VARIABLES:
                0030 C     == Local variables ==
                0031 C     bi, bj    :: tile indices
                0032 C     i, j      :: Loop counters
                0033       INTEGER bi, bj
                0034       INTEGER i, j
                0035 CEOP
                0036
                0037 C--   Initialize R_stic
                0038       DO bj = myByLo(myThid), myByHi(myThid)
                0039        DO bi = myBxLo(myThid), myBxHi(myThid)
                0040         DO j=1-OLy,sNy+OLy
                0041          DO i=1-OLx,sNx+OLx
                0042           R_stic     (i,j,bi,bj) = 0. _d 0
                0043           sticfLength(i,j,bi,bj) = 0. _d 0
                0044          ENDDO
                0045         ENDDO
                0046        ENDDO
                0047       ENDDO
                0048
                0049 C---- ICEFRONT BEGIN
                0050       IF ( STIClengthFile .NE. ' ' ) THEN
                0051        CALL READ_FLD_XY_RS( STIClengthFile, ' ',
                0052      &      sticfLength, 0, myThid )
                0053        _EXCH_XY_RS( sticfLength, myThid )
                0054       ENDIF
                0055
                0056       IF ( STICdepthFile .NE. ' ' ) THEN
                0057        CALL READ_FLD_XY_RS( STICdepthFile, ' ',
                0058      &      R_stic, 0, myThid )
                0059        _EXCH_XY_RS( R_stic, myThid )
                0060       ENDIF
                0061
                0062 C     Make sure that R_stic is positive
                0063       DO bj = myByLo(myThid), myByHi(myThid)
                0064        DO bi = myBxLo(myThid), myBxHi(myThid)
                0065         DO j = 1-OLy, sNy+OLy
                0066          DO i = 1-OLx, sNx+OLx
                0067           R_stic(i,j,bi,bj) = ABS(R_stic(i,j,bi,bj))
                0068          ENDDO
                0069         ENDDO
                0070        ENDDO
                0071       ENDDO
                0072
                0073       RETURN
                0074       END