jscott/pkg_atm2d/pass_seaice_fluxes.F

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

C     !INTERFACE:
      SUBROUTINE PASS_SEAICE_FLUXES( myThid )
C     *==========================================================*
C     | Called after each atm time step - put Qneti, dF/dT, solar|
c     | for seaice to use, accumlate precip fluxes.              |
C     *==========================================================*
        IMPLICIT NONE

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

C     === Global SeaIce Variables ===
#include "THSICE_PARAMS.h"
#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
        
          icFlxAtm(i,j,1,1) = -qneti_2D(i,j) 
          icFlxSw(i,j,1,1) = - solarinc_2D(i,j) !at this point, incident SW
          dFdT(i,j,1,1) = - dFdT_ice_2D(i,j)

C         Note converion to ice model sign convention below
          sum_prcIce(i,j) = sum_prcIce(i,j) -  precipi_2D(i,j)*rhofw   
          sum_evapIce(i,j) = sum_evapIce(i,j) -  evapi_2D(i,j)*rhofw   
          IF (Tair_2D(i,j) .LE. 0. _d 0)    
     &      sum_snowPrc(i,j) = sum_snowPrc(i,j) -  precipi_2D(i,j)*rhofw
 
C          IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN
C            PRINT *,'Pass seaice fluxes at, T :',i,j,Tair_2D(i,j)
C            PRINT *,precipi_2D(i,j),sum_snowPrc(i,j)
C          ENDIF

        ENDDO
      ENDDO  


      RETURN
      END
1	jscott	1.1	#include "ctrparam.h"
2			#include "ATM2D_OPTIONS.h"
3
4			C !INTERFACE:
5			SUBROUTINE PASS_SEAICE_FLUXES( myThid )
6			C ==========================================================
7			C \| Called after each atm time step - put Qneti, dF/dT, solar\|
8			c \| for seaice to use, accumlate precip fluxes. \|
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_PARAMS.h"
19			#include "THSICE_VARS.h"
20
21			C === Atmos/Ocean/Seaice Interface Variables ===
22			#include "ATM2D_VARS.h"
23
24
25			C !INPUT/OUTPUT PARAMETERS:
26			C === Routine arguments ===
27			C myThid - Thread no. that called this routine.
28			INTEGER myThid
29
30			C LOCAL VARIABLES:
31			INTEGER i,j
32
33			DO j=1, sNy
34			DO i=1, sNx
35
36			icFlxAtm(i,j,1,1) = -qneti_2D(i,j)
37			icFlxSw(i,j,1,1) = - solarinc_2D(i,j) !at this point, incident SW
38			dFdT(i,j,1,1) = - dFdT_ice_2D(i,j)
39
40	jscott	1.3	C Note converion to ice model sign convention below
41	jscott	1.1	sum_prcIce(i,j) = sum_prcIce(i,j) - precipi_2D(i,j)*rhofw
42			sum_evapIce(i,j) = sum_evapIce(i,j) - evapi_2D(i,j)*rhofw
43	jscott	1.2	IF (Tair_2D(i,j) .LE. 0. _d 0)
44	jscott	1.3	& sum_snowPrc(i,j) = sum_snowPrc(i,j) - precipi_2D(i,j)*rhofw
45
46	jscott	1.2	C IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN
47	jscott	1.1	C PRINT *,'Pass seaice fluxes at, T :',i,j,Tair_2D(i,j)
48			C PRINT *,precipi_2D(i,j),sum_snowPrc(i,j)
49			C ENDIF
50
51			ENDDO
52			ENDDO
53
54
55			RETURN
56			END