36 |
|
|
37 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
38 |
C === Routine arguments === |
C === Routine arguments === |
39 |
C myIter :: iteration counter for this thread |
C myTime :: Current time in simulation (s) |
40 |
C myTime :: time counter for this thread |
C myIter :: Current iteration number |
41 |
C myThid :: thread number for this instance of the routine. |
C myThid :: My Thread Id. number |
42 |
_RL myTime |
_RL myTime |
43 |
INTEGER myIter |
INTEGER myIter |
44 |
INTEGER myThid |
INTEGER myThid |
45 |
CEOP |
CEOP |
114 |
IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN |
IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN |
115 |
hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj) |
hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj) |
116 |
ELSE |
ELSE |
117 |
hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj) |
hOceMxL(i,j,bi,bj) = drF(1)*hFacC(i,j,1,bi,bj) |
118 |
ENDIF |
ENDIF |
119 |
ENDDO |
ENDDO |
120 |
ENDDO |
ENDDO |
125 |
#endif /* NONLIN_FRSURF */ |
#endif /* NONLIN_FRSURF */ |
126 |
DO j = jMin, jMax |
DO j = jMin, jMax |
127 |
DO i = iMin, iMax |
DO i = iMin, iMax |
128 |
hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj) |
hOceMxL(i,j,bi,bj) = drF(1)*hFacC(i,j,1,bi,bj) |
129 |
ENDDO |
ENDDO |
130 |
ENDDO |
ENDDO |
131 |
ENDIF |
ENDIF |
176 |
CADJ STORE qnet = comlev1, key = iicekey |
CADJ STORE qnet = comlev1, key = iicekey |
177 |
#endif |
#endif |
178 |
|
|
179 |
|
C- do sea-ice advection before getting surface fluxes |
180 |
|
C Note: will inline this S/R once thSIce in Atmos. set-up is settled |
181 |
|
IF ( thSIceAdvScheme.GT.0 ) |
182 |
|
& CALL THSICE_DO_ADVECT( |
183 |
|
I bi,bj, myTime, myIter, myThid ) |
184 |
|
|
185 |
#ifdef ALLOW_BULK_FORCE |
#ifdef ALLOW_BULK_FORCE |
186 |
IF ( useBulkforce ) THEN |
IF ( useBulkforce ) THEN |
187 |
CALL THSICE_GET_PRECIP( |
CALL THSICE_GET_PRECIP( |
245 |
ENDDO |
ENDDO |
246 |
|
|
247 |
|
|
248 |
#ifdef ALLOW_SEAICE |
IF ( useSEAICE .OR. thSIceAdvScheme.GT.0 ) THEN |
|
IF ( useSEAICE ) THEN |
|
249 |
C-- Exchange fields that are advected by seaice dynamics |
C-- Exchange fields that are advected by seaice dynamics |
250 |
_EXCH_XY_R8( iceMask, myThid ) |
_EXCH_XY_R8( iceMask, myThid ) |
251 |
_EXCH_XY_R8( iceHeight, myThid ) |
_EXCH_XY_R8( iceHeight, myThid ) |
257 |
& _EXCH_XY_RS( sIceLoad, myThid ) |
& _EXCH_XY_RS( sIceLoad, myThid ) |
258 |
#endif |
#endif |
259 |
ENDIF |
ENDIF |
|
#endif /* ALLOW_SEAICE */ |
|
260 |
|
|
261 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
262 |
#endif /*ALLOW_THSICE*/ |
#endif /*ALLOW_THSICE*/ |