jscott/pkg_atm2d/sum_seaice_out.F

#include "ctrparam.h"
#include "ATM2D_OPTIONS.h"

C     !INTERFACE:
      SUBROUTINE SUM_SEAICE_OUT( myThid )
C     *==========================================================*
C     | |
c     | |
C     *==========================================================*
        IMPLICIT NONE

#include "ATMSIZE.h"
#include "SIZE.h"
#include "GRID.h"
#include "EEPARAMS.h"

C     === Global SeaIce Variables ===
#include "THSICE_VARS.h"

C     === Atmos/Ocean/Seaice Interface Variables ===
#include "ATM2D_VARS.h"


C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid - Thread no. that called this routine.
      INTEGER myThid

C     LOCAL VARIABLES:
      INTEGER i,j

      DO j=1, sNy
        DO i=1,sNx
        
          IF (iceMask(i,j,1,1) .NE. 0. _d 0) THEN
            sum_sHeat(i,j) = sum_sHeat(i,j) + sHeating(i,j,1,1)
            sum_flxCnB(i,j)= sum_flxCnB(i,j) + flxCndBt(i,j,1,1)

C Add contribution of SW passing through ice
C note sign convention of sum_solarnet is +=up, opposite of ice model
            sum_solarnet(i,j)= sum_solarnet(i,j)
     &                  - icFlxSW(i,j,1,1)*iceMask(i,j,1,1)
          ENDIF
          
        ENDDO
      ENDDO  


      RETURN
      END
1	#include "ctrparam.h"
2	#include "ATM2D_OPTIONS.h"
3
4	C !INTERFACE:
5	SUBROUTINE SUM_SEAICE_OUT( myThid )
6	C ==========================================================
7	C \| \|
8	c \| \|
9	C ==========================================================
10	IMPLICIT NONE
11
12	#include "ATMSIZE.h"
13	#include "SIZE.h"
14	#include "GRID.h"
15	#include "EEPARAMS.h"
16
17	C === Global SeaIce Variables ===
18	#include "THSICE_VARS.h"
19
20	C === Atmos/Ocean/Seaice Interface Variables ===
21	#include "ATM2D_VARS.h"
22
23
24	C !INPUT/OUTPUT PARAMETERS:
25	C === Routine arguments ===
26	C myThid - Thread no. that called this routine.
27	INTEGER myThid
28
29	C LOCAL VARIABLES:
30	INTEGER i,j
31
32	DO j=1, sNy
33	DO i=1,sNx
34
35	IF (iceMask(i,j,1,1) .NE. 0. _d 0) THEN
36	sum_sHeat(i,j) = sum_sHeat(i,j) + sHeating(i,j,1,1)
37	sum_flxCnB(i,j)= sum_flxCnB(i,j) + flxCndBt(i,j,1,1)
38
39	C Add contribution of SW passing through ice
40	C note sign convention of sum_solarnet is +=up, opposite of ice model
41	sum_solarnet(i,j)= sum_solarnet(i,j)
42	& - icFlxSW(i,j,1,1)*iceMask(i,j,1,1)
43	ENDIF
44
45	ENDDO
46	ENDDO
47
48
49	RETURN
50	END