pkg/atm_compon_interf/atm_apply_import.F

C $Header: /u/gcmpack/MITgcm/pkg/atm_compon_interf/atm_apply_import.F,v 1.3 2015/11/12 01:01:57 jmc Exp $
C $Name:  $

#include "ATM_CPL_OPTIONS.h"

CBOP
C     !ROUTINE: ATM_APPLY_IMPORT
C     !INTERFACE:
      SUBROUTINE ATM_APPLY_IMPORT(
     I               land_frc,
     U               atmSST, atmSIfrc,
     I               myTime, myIter, bi, bj, myThid )

C     !DESCRIPTION: \bv
C     *================================================================*
C     | S/R ATM_APPLY_IMPORT
C     | o Apply imported coupling data to ATM surface BC over ocean
C     *================================================================*
C     | Note: when using sea-ice, fill in Mixed layer fields instead
C     |      to be used later as SST by Atmos. + Sea-Ice
C     *================================================================*
C     \ev

C     !USES:
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "CPL_PARAMS.h"

#ifdef ALLOW_THSICE
# include "THSICE_VARS.h"
#endif

C-- Coupled to the Ocean :
#include "ATMCPL.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     land_frc :: land fraction [0-1]
C     atmSST   :: sea surface temp [K], used in ATM component
C     atmSIfrc :: sea-ice fraction [0-1], used in ATM component
C     myTime   :: Current time of simulation ( s )
C     myIter   :: Current iteration number in simulation
C     bi,bj    :: Tile index
C     myThid   :: Number of this instance of the routine
      _RS  land_frc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  atmSST  (sNx,sNy)
      _RL  atmSIfrc(sNx,sNy)
      _RL myTime
      INTEGER myIter, bi, bj, myThid
CEOP

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

#ifdef COMPONENT_MODULE
C     == Local variables ==
C     i,j    :: Loop counters
      INTEGER i,j

#ifdef ALLOW_THSICE
      IF ( useImportThSIce .AND. useThSIce ) THEN
       IF ( MOD(myIter,cplSendFrq_iter).EQ.0 ) THEN
C--   Put thSIce imported fields from the ocean component in thSIce common block
        DO j=1,sNy
         DO i=1,sNx
           iceMask   (i,j,bi,bj) = sIceFrac_cpl (i,j,bi,bj)
           iceHeight (i,j,bi,bj) = sIceThick_cpl(i,j,bi,bj)
           snowHeight(i,j,bi,bj) = sIceSnowH_cpl(i,j,bi,bj)
           Qice1     (i,j,bi,bj) = sIceQ1_cpl   (i,j,bi,bj)
           Qice2     (i,j,bi,bj) = sIceQ2_cpl   (i,j,bi,bj)
         ENDDO
        ENDDO
       ENDIF
C--   end if useImportThSIce & useThSIce
      ENDIF
#endif /* ALLOW_THSICE */

#ifdef ALLOW_THSICE
      IF ( useThSIce ) THEN
C--   Put fields from the ocean component in Mixed-layer arrays:

C-    fill in hOceMxL with Mixed-layer Depth from the ocean component
       IF ( useImportMxlD ) THEN
        DO j=1,sNy
         DO i=1,sNx
          IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
            hOceMxL(i,j,bi,bj) = ocMxlD(i,j,bi,bj)
          ENDIF
         ENDDO
        ENDDO
       ENDIF

C-    fill in tOceMxL with Sea-Surface Temp. from the ocean component
       IF ( useImportSST ) THEN
        DO j=1,sNy
         DO i=1,sNx
          IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
            tOceMxL(i,j,bi,bj) = SSTocn(i,j,bi,bj)
          ENDIF
         ENDDO
        ENDDO
       ENDIF

C-    fill in sOceMxL with Sea-Surf Salinity from the ocean component
       IF ( useImportSSS ) THEN
        DO j=1,sNy
         DO i=1,sNx
          IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
            sOceMxL(i,j,bi,bj) = SSSocn(i,j,bi,bj)
          ENDIF
         ENDDO
        ENDDO
       ENDIF

C-    fill in v2ocMxL with surf. velocity^2  from the ocean component
       IF ( useImportVsq ) THEN
        DO j=1,sNy
         DO i=1,sNx
          IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
            v2ocMxL(i,j,bi,bj) = vSqocn(i,j,bi,bj)
          ENDIF
         ENDDO
        ENDDO
       ENDIF

      ELSEIF ( useAtm_Phys ) THEN
#else /* ALLOW_THSICE */
      IF ( useAtm_Phys ) THEN
#endif /* ALLOW_THSICE */
C--   supply imported fields to Atm_Phys pkg:

       IF ( useImportSST ) THEN
        DO j=1,sNy
         DO i=1,sNx
            atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
         ENDDO
        ENDDO
       ENDIF

      ELSE
C--   supply imported fields to AIM (Phys) pkg:

       IF ( useImportSST ) THEN
        DO j=1,sNy
         DO i=1,sNx
          IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN

C--   take SST from the ocean compon where Sea-Ice fraction is zero
           IF ( atmSIfrc(i,j).EQ.0. ) THEN
            atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
           ELSEIF ( SSTocn(i,j,bi,bj).GE. -1. _d 0)  THEN
C--   take SST from the ocean compon if clearly warmer than freezing
C       then reset sea-ice fraction
            atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
            atmSIfrc(i,j) = 0.
           ENDIF

          ENDIF
         ENDDO
        ENDDO
       ENDIF

C-- if useThSIce / else / endif
      ENDIF

#endif /* COMPONENT_MODULE */

      RETURN
      END
1	jmc	1.4	C $Header: /u/gcmpack/MITgcm/pkg/atm_compon_interf/atm_apply_import.F,v 1.3 2015/11/12 01:01:57 jmc Exp $
2	jmc	1.1	C $Name: $
3
4	jmc	1.3	#include "ATM_CPL_OPTIONS.h"
5	jmc	1.1
6			CBOP
7			C !ROUTINE: ATM_APPLY_IMPORT
8			C !INTERFACE:
9	jmc	1.2	SUBROUTINE ATM_APPLY_IMPORT(
10	jmc	1.1	I land_frc,
11	jmc	1.2	U atmSST, atmSIfrc,
12	jmc	1.1	I myTime, myIter, bi, bj, myThid )
13
14			C !DESCRIPTION: \bv
15			C ================================================================
16			C \| S/R ATM_APPLY_IMPORT
17	jmc	1.2	C \| o Apply imported coupling data to ATM surface BC over ocean
18	jmc	1.1	C ================================================================
19			C \| Note: when using sea-ice, fill in Mixed layer fields instead
20	jmc	1.2	C \| to be used later as SST by Atmos. + Sea-Ice
21	jmc	1.1	C ================================================================
22			C \ev
23
24			C !USES:
25			IMPLICIT NONE
26
27			C == Global variables ===
28			#include "SIZE.h"
29			#include "EEPARAMS.h"
30			#include "PARAMS.h"
31			#include "CPL_PARAMS.h"
32
33			#ifdef ALLOW_THSICE
34			# include "THSICE_VARS.h"
35			#endif
36
37			C-- Coupled to the Ocean :
38			#include "ATMCPL.h"
39
40			C !INPUT/OUTPUT PARAMETERS:
41			C == Routine arguments ==
42			C land_frc :: land fraction [0-1]
43	jmc	1.2	C atmSST :: sea surface temp [K], used in ATM component
44			C atmSIfrc :: sea-ice fraction [0-1], used in ATM component
45	jmc	1.1	C myTime :: Current time of simulation ( s )
46			C myIter :: Current iteration number in simulation
47			C bi,bj :: Tile index
48			C myThid :: Number of this instance of the routine
49			_RS land_frc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
50	jmc	1.2	_RL atmSST (sNx,sNy)
51			_RL atmSIfrc(sNx,sNy)
52	jmc	1.1	_RL myTime
53			INTEGER myIter, bi, bj, myThid
54			CEOP
55
56			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
57
58			#ifdef COMPONENT_MODULE
59			C == Local variables ==
60			C i,j :: Loop counters
61			INTEGER i,j
62
63			#ifdef ALLOW_THSICE
64	jmc	1.4	IF ( useImportThSIce .AND. useThSIce ) THEN
65			IF ( MOD(myIter,cplSendFrq_iter).EQ.0 ) THEN
66			C-- Put thSIce imported fields from the ocean component in thSIce common block
67			DO j=1,sNy
68			DO i=1,sNx
69			iceMask (i,j,bi,bj) = sIceFrac_cpl (i,j,bi,bj)
70			iceHeight (i,j,bi,bj) = sIceThick_cpl(i,j,bi,bj)
71			snowHeight(i,j,bi,bj) = sIceSnowH_cpl(i,j,bi,bj)
72			Qice1 (i,j,bi,bj) = sIceQ1_cpl (i,j,bi,bj)
73			Qice2 (i,j,bi,bj) = sIceQ2_cpl (i,j,bi,bj)
74			ENDDO
75			ENDDO
76			ENDIF
77			C-- end if useImportThSIce & useThSIce
78			ENDIF
79			#endif /* ALLOW_THSICE */
80
81			#ifdef ALLOW_THSICE
82	jmc	1.1	IF ( useThSIce ) THEN
83			C-- Put fields from the ocean component in Mixed-layer arrays:
84
85			C- fill in hOceMxL with Mixed-layer Depth from the ocean component
86	jmc	1.2	IF ( useImportMxlD ) THEN
87	jmc	1.1	DO j=1,sNy
88			DO i=1,sNx
89			IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
90			hOceMxL(i,j,bi,bj) = ocMxlD(i,j,bi,bj)
91			ENDIF
92			ENDDO
93			ENDDO
94			ENDIF
95
96			C- fill in tOceMxL with Sea-Surface Temp. from the ocean component
97	jmc	1.2	IF ( useImportSST ) THEN
98	jmc	1.1	DO j=1,sNy
99			DO i=1,sNx
100			IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
101			tOceMxL(i,j,bi,bj) = SSTocn(i,j,bi,bj)
102			ENDIF
103			ENDDO
104			ENDDO
105			ENDIF
106
107			C- fill in sOceMxL with Sea-Surf Salinity from the ocean component
108	jmc	1.2	IF ( useImportSSS ) THEN
109	jmc	1.1	DO j=1,sNy
110			DO i=1,sNx
111			IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
112			sOceMxL(i,j,bi,bj) = SSSocn(i,j,bi,bj)
113			ENDIF
114			ENDDO
115			ENDDO
116			ENDIF
117
118			C- fill in v2ocMxL with surf. velocity^2 from the ocean component
119	jmc	1.2	IF ( useImportVsq ) THEN
120	jmc	1.1	DO j=1,sNy
121			DO i=1,sNx
122			IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
123			v2ocMxL(i,j,bi,bj) = vSqocn(i,j,bi,bj)
124			ENDIF
125			ENDDO
126			ENDDO
127			ENDIF
128
129	jmc	1.2	ELSEIF ( useAtm_Phys ) THEN
130	jmc	1.4	#else /* ALLOW_THSICE */
131	jmc	1.2	IF ( useAtm_Phys ) THEN
132	jmc	1.1	#endif /* ALLOW_THSICE */
133	jmc	1.2	C-- supply imported fields to Atm_Phys pkg:
134
135			IF ( useImportSST ) THEN
136			DO j=1,sNy
137			DO i=1,sNx
138			atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
139			ENDDO
140			ENDDO
141			ENDIF
142
143			ELSE
144			C-- supply imported fields to AIM (Phys) pkg:
145	jmc	1.1
146			IF ( useImportSST ) THEN
147			DO j=1,sNy
148			DO i=1,sNx
149			IF ( land_frc(i,j,bi,bj) .LT. 1. ) THEN
150
151			C-- take SST from the ocean compon where Sea-Ice fraction is zero
152	jmc	1.2	IF ( atmSIfrc(i,j).EQ.0. ) THEN
153			atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
154	jmc	1.1	ELSEIF ( SSTocn(i,j,bi,bj).GE. -1. _d 0) THEN
155			C-- take SST from the ocean compon if clearly warmer than freezing
156			C then reset sea-ice fraction
157	jmc	1.2	atmSST(i,j) = SSTocn(i,j,bi,bj)+celsius2K
158			atmSIfrc(i,j) = 0.
159	jmc	1.1	ENDIF
160
161			ENDIF
162			ENDDO
163			ENDDO
164			ENDIF
165
166			C-- if useThSIce / else / endif
167			ENDIF
168
169			#endif /* COMPONENT_MODULE */
170
171			RETURN
172			END