/[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.11 by jmc, Sun Jun 25 22:31:02 2006 UTC revision 1.12 by jmc, Wed Apr 4 02:40:42 2007 UTC
# Line 36  C     === Global variables === Line 36  C     === Global variables ===
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
# Line 114  C--     Mixed layer thickness: take the Line 114  C--     Mixed layer thickness: take the
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
# Line 125  C--     Mixed layer thickness: take the Line 125  C--     Mixed layer thickness: take the
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
# Line 176  CADJ STORE empmr   = comlev1, key = iice Line 176  CADJ STORE empmr   = comlev1, key = iice
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(
# Line 239  C--  end bi,bj loop Line 245  C--  end bi,bj loop
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 )
# Line 252  C--   Exchange fields that are advected Line 257  C--   Exchange fields that are advected
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*/
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22