110 |
# endif |
# endif |
111 |
#endif |
#endif |
112 |
|
|
|
cph( |
|
|
c print *, 'ph-thsice-1 in thsice_main' |
|
|
cph) |
|
113 |
C-- Mixed layer thickness: take the 1rst layer |
C-- Mixed layer thickness: take the 1rst layer |
114 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
115 |
IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN |
IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN |
185 |
CADJ STORE qnet = comlev1, key = iicekey |
CADJ STORE qnet = comlev1, key = iicekey |
186 |
#endif |
#endif |
187 |
|
|
|
cph( |
|
|
c print *, 'ph-thsice-2 in thsice_main' |
|
|
cph) |
|
188 |
C- do sea-ice advection before getting surface fluxes |
C- do sea-ice advection before getting surface fluxes |
189 |
C Note: will inline this S/R once thSIce in Atmos. set-up is settled |
C Note: will inline this S/R once thSIce in Atmos. set-up is settled |
190 |
IF ( thSIceAdvScheme.GT.0 ) |
IF ( thSIceAdvScheme.GT.0 ) |
210 |
ENDIF |
ENDIF |
211 |
#endif |
#endif |
212 |
|
|
|
|
|
|
cph( |
|
|
c print *, 'ph-thsice-3 in thsice_main' |
|
|
cph) |
|
213 |
CALL THSICE_STEP_TEMP( |
CALL THSICE_STEP_TEMP( |
214 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
215 |
I myTime, myIter, myThid ) |
I myTime, myIter, myThid ) |
225 |
CADJ STORE snowAge = comlev1, key = iicekey |
CADJ STORE snowAge = comlev1, key = iicekey |
226 |
#endif |
#endif |
227 |
|
|
|
cph( |
|
|
c print *, 'ph-thsice-4 in thsice_main' |
|
|
cph) |
|
228 |
CALL THSICE_STEP_FWD( |
CALL THSICE_STEP_FWD( |
229 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
230 |
I prcAtm, |
I prcAtm, |
231 |
I myTime, myIter, myThid ) |
I myTime, myIter, myThid ) |
232 |
|
|
|
cph( |
|
|
c print *, 'ph-thsice-5 in thsice_main' |
|
|
cph) |
|
233 |
CALL THSICE_AVE( |
CALL THSICE_AVE( |
234 |
I bi,bj, myTime, myIter, myThid ) |
I bi,bj, myTime, myIter, myThid ) |
235 |
|
|
263 |
ENDDO |
ENDDO |
264 |
ENDDO |
ENDDO |
265 |
|
|
|
|
|
|
cph( |
|
|
c print *, 'ph-thsice-6 in thsice_main' |
|
|
cph) |
|
266 |
IF ( useSEAICE .OR. thSIceAdvScheme.GT.0 ) THEN |
IF ( useSEAICE .OR. thSIceAdvScheme.GT.0 ) THEN |
267 |
C-- Exchange fields that are advected by seaice dynamics |
C-- Exchange fields that are advected by seaice dynamics |
268 |
_EXCH_XY_R8( iceMask, myThid ) |
_EXCH_XY_R8( iceMask, myThid ) |