/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_advect_thickness.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_advect_thickness.F

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

revision 1.3 by dgoldberg, Thu Jul 26 16:13:18 2012 UTC revision 1.4 by dgoldberg, Tue Sep 4 21:11:44 2012 UTC
# Line 36  C     === Global variables === Line 36  C     === Global variables ===
36        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
37        _RL thick_bd        _RL thick_bd
38        _RL SLOPE_LIMITER        _RL SLOPE_LIMITER
39        _RL sec_per_year, time_step_loc        _RL sec_per_year, time_step_loc, MR
40        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
41        external SLOPE_LIMITER        external SLOPE_LIMITER
42    
# Line 125  CADJ STORE streamice_hmask  = comlev1, k Line 125  CADJ STORE streamice_hmask  = comlev1, k
125        CALL STREAMICE_ADV_FRONT ( myThid, time_step_loc )        CALL STREAMICE_ADV_FRONT ( myThid, time_step_loc )
126    
127    
128    ! NOW WE APPLY MELT RATES !!
129    ! THIS MAY BE MOVED TO A SEPARATE SUBROUTINE
130    
131        _EXCH_XY_RL( H_streamice, myThid )        DO bj=myByLo(myThid),myByHi(myThid)
132        _EXCH_XY_RL( area_shelf_streamice, myThid )         DO bi=myBxLo(myThid),myBxHi(myThid)
133        _EXCH_XY_RL( STREAMICE_hmask, myThid )          DO j=1-OLy,sNy+OLy
134             DO i=1-OLx,sNx+OLx
135               IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0 .or.
136         &       STREAMICE_hmask(i,j,bi,bj).eq.2.0) THEN
137                 MR = (1.-float_frac_streamice(i,j,bi,bj)) *
138         &             BDOT_STREAMICE(i,j,bi,bj)
139                 IF (MR * time_step_loc .gt.
140         &         H_streamice (i,j,bi,bj)) THEN
141                    H_streamice (i,j,bi,bj) = 0. _d 0
142                    STREAMICE_hmask(i,j,bi,bj) = 0.
143                 ELSE
144                   H_streamice (i,j,bi,bj) =
145         &          H_streamice (i,j,bi,bj) - MR
146                 ENDIF
147               ENDIF
148             ENDDO
149            ENDDO
150           ENDDO
151          ENDDO    
152    
153                
154        WRITE(msgBuf,'(A)') 'END STREAMICE_ADVECT_THICKNESS'        WRITE(msgBuf,'(A)') 'END STREAMICE_ADVECT_THICKNESS'
155         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22