pkg/seaice/seaice_itd_sum.F

86b84a92fc Patr*0001 #include "SEAICE_OPTIONS.h"
                0002
e7af59f6fd Jean*0003 CBOP
86b84a92fc Patr*0004 C !ROUTINE: SEAICE_ITD_SUM
                0005
                0006 C !INTERFACE: ==========================================================
                0007       SUBROUTINE SEAICE_ITD_SUM(
e7af59f6fd Jean*0008      I                  bi, bj, myTime, myIter, myThid )
86b84a92fc Patr*0009
                0010 C !DESCRIPTION: \bv
                0011 C     *===========================================================*
                0012 C     | SUBROUTINE SEAICE_ITD_SUM
                0013 C     | o sum ice area and volume over all ITD categories
                0014 C     |   and write into AREA and HEFF
                0015 C     |
                0016 C     | Torge Martin, Mar. 2012, torge@mit.edu
                0017 C     *===========================================================*
                0018 C \ev
                0019
                0020 C !USES: ===============================================================
                0021       IMPLICIT NONE
                0022
                0023 C     === Global variables ===
e7af59f6fd Jean*0024 C     nITD      :: number of seaice thickness categories
86b84a92fc Patr*0025 C     AREA      :: sea ice area   total
                0026 C     AREAITD   :: sea ice area   by category
                0027 C     HEFF      :: sea ice volume total
                0028 C     HEFFITD   :: sea ice volume by category
                0029 C     HSNOW     :: snow    volume total
                0030 C     HSNOWITD  :: snow    volume by category
                0031 #include "SIZE.h"
                0032 #include "EEPARAMS.h"
e7af59f6fd Jean*0033 c#include "PARAMS.h"
86b84a92fc Patr*0034 #include "SEAICE_SIZE.h"
e7af59f6fd Jean*0035 c#include "SEAICE_PARAMS.h"
86b84a92fc Patr*0036 #include "SEAICE.h"
                0037
                0038 C !INPUT PARAMETERS: ===================================================
e7af59f6fd Jean*0039 C     bi,bj     :: Current tile indices
                0040 C     myTime    :: Current time in simulation
                0041 C     myIter    :: Current iteration number
                0042 C     myThid    :: my Thread Id number
86b84a92fc Patr*0043       _RL myTime
                0044       INTEGER bi, bj
                0045       INTEGER myIter
                0046       INTEGER myThid
e7af59f6fd Jean*0047 CEOP
86b84a92fc Patr*0048
                0049 #ifdef SEAICE_ITD
                0050 C !LOCAL VARIABLES: ====================================================
e7af59f6fd Jean*0051 C     i,j,k     :: inner loop counters
86b84a92fc Patr*0052       INTEGER i, j, k
                0053
                0054 C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0055
                0056 c      DO bj=myByLo(myThid),myByHi(myThid)
                0057 c       DO bi=myBxLo(myThid),myBxHi(myThid)
                0058 C must now be called within bi,bj loop
                0059
                0060         DO j=1-OLy,sNy+OLy
                0061          DO i=1-OLx,sNx+OLx
e7af59f6fd Jean*0062            AREA (i,j,bi,bj) = AREAITD (i,j,1,bi,bj)
                0063            HEFF (i,j,bi,bj) = HEFFITD (i,j,1,bi,bj)
                0064            HSNOW(i,j,bi,bj) = HSNOWITD(i,j,1,bi,bj)
86b84a92fc Patr*0065          ENDDO
                0066         ENDDO
                0067
                0068         DO k=2,nITD
                0069          DO j=1-OLy,sNy+OLy
                0070           DO i=1-OLx,sNx+OLx
                0071            AREA (i,j,bi,bj) = AREA (i,j,bi,bj) + AREAITD (i,j,k,bi,bj)
                0072            HEFF (i,j,bi,bj) = HEFF (i,j,bi,bj) + HEFFITD (i,j,k,bi,bj)
                0073            HSNOW(i,j,bi,bj) = HSNOW(i,j,bi,bj) + HSNOWITD(i,j,k,bi,bj)
                0074           ENDDO
                0075          ENDDO
                0076         ENDDO
                0077
                0078 c       ENDDO
                0079 c      ENDDO
                0080
                0081 C---+-|--1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0082 #endif /* SEAICE_ITD */
                0083       RETURN
b8f04b8c26 Jean*0084       END