pkg/atm2d/sum_thsice_out.F

aea46876b4 Jeff*0001 #include "ctrparam.h"
                0002 #include "ATM2D_OPTIONS.h"
                0003
                0004 C     !INTERFACE:
                0005       SUBROUTINE SUM_THSICE_OUT( myThid )
                0006 C     *==========================================================*
                0007 C     | Sum output from temp seaice step for passing to ocean.   |
                0008 C     *==========================================================*
                0009         IMPLICIT NONE
                0010
                0011 #include "ATMSIZE.h"
                0012 #include "SIZE.h"
                0013 #include "GRID.h"
                0014 #include "EEPARAMS.h"
                0015
                0016 C     === Global SeaIce Variables ===
                0017 #include "THSICE_VARS.h"
                0018
                0019 C     === Atmos/Ocean/Seaice Interface Variables ===
                0020 #include "ATM2D_VARS.h"
                0021
                0022
                0023 C     !INPUT/OUTPUT PARAMETERS:
                0024 C     === Routine arguments ===
                0025 C     myThid - Thread no. that called this routine.
                0026       INTEGER myThid
                0027
                0028 C     LOCAL VARIABLES:
                0029       INTEGER i,j
                0030
                0031       DO j=1, sNy
                0032         DO i=1,sNx
9274434acc Jean*0033
aea46876b4 Jeff*0034           IF (iceMask(i,j,1,1) .NE. 0. _d 0) THEN
                0035             sum_sHeat(i,j) = sum_sHeat(i,j) + sHeating(i,j,1,1)
                0036             sum_flxCnB(i,j)= sum_flxCnB(i,j) + flxCndBt(i,j,1,1)
                0037
                0038 C Add contribution of SW passing through ice
                0039             sum_solarnet(i,j)= sum_solarnet(i,j)
                0040      &                  + icFlxSW(i,j,1,1)*iceMask(i,j,1,1)
                0041 C zero out icFlxSw or it will be double counted in thsice_step_fwd
                0042             icFlxSW(i,j,1,1)= 0. _d 0
                0043
                0044           ENDIF
9274434acc Jean*0045
aea46876b4 Jeff*0046         ENDDO
9274434acc Jean*0047       ENDDO
aea46876b4 Jeff*0048
                0049
                0050       RETURN
                0051       END