pkg/shelfice/shelfice_mask_seaice.F

148dd84005 jm-c 0001 #include "SHELFICE_OPTIONS.h"
                0002
                0003 CBOP
                0004 C !ROUTINE: SHELFICE_MASK_SEAICE
                0005
                0006 C !INTERFACE: ==========================================================
                0007       SUBROUTINE SHELFICE_MASK_SEAICE(
                0008      U                    SIfield,
                0009      I                    bi, bj, myIter, myThid )
                0010
                0011 C !DESCRIPTION:
                0012 C  Mask seaice field (mask) to prevent seaice to spread
                0013 C         in grid-cells that are occupied by ice-shelf
                0014
                0015 C !USES: ===============================================================
                0016       IMPLICIT NONE
                0017 #include "SIZE.h"
                0018 #include "EEPARAMS.h"
                0019 #include "PARAMS.h"
                0020 #include "SHELFICE.h"
                0021
                0022 C !INPUT/OUTPUT PARAMETERS: ============================================
                0023 C     SIfield     :: seaice field @ grid-cell center
                0024 C     bi, bj      :: Current tile indices
                0025 C     myIter      :: Current iteration number in simulation
                0026 C     myThid      :: my Thread Id number
                0027       _RL SIfield(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0028       INTEGER bi, bj
                0029       INTEGER myIter
                0030       INTEGER myThid
                0031
                0032 C !LOCAL VARIABLES: ====================================================
                0033 C     i, j        :: Loop counters
                0034       INTEGER i, j
                0035 CEOP
                0036 C     msgBuf      :: Informational/error message buffer
                0037 c     CHARACTER*(MAX_LEN_MBUF) msgBuf
                0038
                0039 c     DO bj=myByLo(myThid),myByHi(myThid)
                0040 c      DO bi=myBxLo(myThid),myBxHi(myThid)
                0041
                0042         DO j=1-OLy,sNy+OLy
                0043          DO i=1-OLx,sNx+OLx
                0044           IF ( kTopC(i,j,bi,bj).NE.0 ) SIfield(i,j,bi,bj) = zeroRL
                0045          ENDDO
                0046         ENDDO
                0047
                0048 c      ENDDO
                0049 c     ENDDO
                0050
                0051       RETURN
                0052       END