/[MITgcm]/MITgcm/pkg/thsice/thsice_main.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_main.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.13 by heimbach, Mon Apr 16 22:38:24 2007 UTC revision 1.14 by heimbach, Tue Apr 17 23:42:33 2007 UTC
# Line 105  C      to be valid at the boundaries ( d Line 105  C      to be valid at the boundaries ( d
105  CADJ STORE ocefwfx(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte  CADJ STORE ocefwfx(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte
106  CADJ STORE oceqnet(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte  CADJ STORE oceqnet(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte
107  CADJ STORE ocesflx(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte  CADJ STORE ocesflx(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte
108    # ifdef ALLOW_EXF
109    CADJ STORE qsw(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte
110    # endif
111  #endif  #endif
112    
113    cph(
114              print *, 'ph-thsice-1 in thsice_main'
115    cph)
116  C--     Mixed layer thickness: take the 1rst layer  C--     Mixed layer thickness: take the 1rst layer
117  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
118          IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN          IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
# Line 182  CADJ STORE empmr   = comlev1, key = iice Line 188  CADJ STORE empmr   = comlev1, key = iice
188  CADJ STORE qnet    = comlev1, key = iicekey  CADJ STORE qnet    = comlev1, key = iicekey
189  #endif  #endif
190    
191    cph(
192              print *, 'ph-thsice-2 in thsice_main'
193    cph)
194  C-      do sea-ice advection before getting surface fluxes  C-      do sea-ice advection before getting surface fluxes
195  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
196          IF ( thSIceAdvScheme.GT.0 )          IF ( thSIceAdvScheme.GT.0 )
# Line 208  C Note: will inline this S/R once thSIce Line 217  C Note: will inline this S/R once thSIce
217  #endif  #endif
218    
219    
220    cph(
221              print *, 'ph-thsice-3 in thsice_main'
222    cph)
223          CALL THSICE_STEP_TEMP(          CALL THSICE_STEP_TEMP(
224       I                     bi, bj, iMin, iMax, jMin, jMax,       I                     bi, bj, iMin, iMax, jMin, jMax,
225       I                     myTime, myIter, myThid )       I                     myTime, myIter, myThid )
# Line 223  CADJ STORE Qice2   = comlev1, key = iice Line 235  CADJ STORE Qice2   = comlev1, key = iice
235  CADJ STORE snowAge = comlev1, key = iicekey  CADJ STORE snowAge = comlev1, key = iicekey
236  #endif  #endif
237    
238    cph(
239              print *, 'ph-thsice-4 in thsice_main'
240    cph)
241          CALL THSICE_STEP_FWD(          CALL THSICE_STEP_FWD(
242       I                     bi, bj, iMin, iMax, jMin, jMax,       I                     bi, bj, iMin, iMax, jMin, jMax,
243       I                     prcAtm,       I                     prcAtm,
244       I                     myTime, myIter, myThid )       I                     myTime, myIter, myThid )
245    
246    cph(
247              print *, 'ph-thsice-5 in thsice_main'
248    cph)
249          CALL THSICE_AVE(          CALL THSICE_AVE(
250       I                     bi,bj, myTime, myIter, myThid )       I                     bi,bj, myTime, myIter, myThid )
251    
# Line 262  C--  end bi,bj loop Line 280  C--  end bi,bj loop
280        ENDDO        ENDDO
281    
282    
283    cph(
284              print *, 'ph-thsice-6 in thsice_main'
285    cph)
286        IF ( useSEAICE .OR. thSIceAdvScheme.GT.0 ) THEN        IF ( useSEAICE .OR. thSIceAdvScheme.GT.0 ) THEN
287  C--   Exchange fields that are advected by seaice dynamics  C--   Exchange fields that are advected by seaice dynamics
288          _EXCH_XY_R8( iceMask, myThid )          _EXCH_XY_R8( iceMask, myThid )
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22