/[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.15 by jmc, Thu Apr 19 02:59:29 2007 UTC revision 1.16 by heimbach, Thu Apr 19 13:22:36 2007 UTC
# Line 110  CADJ STORE qsw(:,:,bi,bj) = comlev1_bibj Line 110  CADJ STORE qsw(:,:,bi,bj) = comlev1_bibj
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
# Line 188  CADJ STORE empmr   = comlev1, key = iice Line 185  CADJ STORE empmr   = comlev1, key = iice
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 )
# Line 216  C Note: will inline this S/R once thSIce Line 210  C Note: will inline this S/R once thSIce
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 )
# Line 235  CADJ STORE Qice2   = comlev1, key = iice Line 225  CADJ STORE Qice2   = comlev1, key = iice
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    
# Line 279  C--  end bi,bj loop Line 263  C--  end bi,bj loop
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 )

Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16

  ViewVC Help
Powered by ViewVC 1.1.22