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

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

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

revision 1.20 by heimbach, Mon Apr 16 22:38:24 2007 UTC revision 1.21 by heimbach, Tue Apr 17 23:42:33 2007 UTC
# Line 129  C-    Initialise Line 129  C-    Initialise
129    
130        ageFac = 1. _d 0 - thSIce_deltaT/snowAgTime        ageFac = 1. _d 0 - thSIce_deltaT/snowAgTime
131        snowFac = thSIce_deltaT/(rhos*hNewSnowAge)        snowFac = thSIce_deltaT/(rhos*hNewSnowAge)
132    
133    #ifdef ALLOW_AUTODIFF_TAMC
134    CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj, key=iicekey, byte=isbyte
135    #endif
136        DO j = jMin, jMax        DO j = jMin, jMax
137         DO i = iMin, iMax         DO i = iMin, iMax
138          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN
# Line 287  C------- Line 291  C-------
291       O          flx2oc, frw2oc, fsalt,       O          flx2oc, frw2oc, fsalt,
292       I          myTime, myIter, myThid )       I          myTime, myIter, myThid )
293    
294    #ifdef ALLOW_AUTODIFF_TAMC
295    CADJ STORE snowHeight(:,:,bi,bj) =
296    CADJ &     comlev1_bibj, key=iicekey, byte=isbyte
297    #endif
298        DO j = jMin, jMax        DO j = jMin, jMax
299         DO i = iMin, iMax         DO i = iMin, iMax
300          IF (frzmltMxL(i,j).GT.0. _d 0) THEN          IF (frzmltMxL(i,j).GT.0. _d 0) THEN
# Line 321  C--    Net fluxes : Line 329  C--    Net fluxes :
329            Qice1(i,j,bi,bj)    = Lfresh            Qice1(i,j,bi,bj)    = Lfresh
330            Qice2(i,j,bi,bj)    = Lfresh            Qice2(i,j,bi,bj)    = Lfresh
331          ENDIF          ENDIF
332           ENDDO
333          ENDDO
334    
335  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
336    # ifdef ALLOW_AUTODIFF_TAMC
337    CADJ STORE snowHeight(:,:,bi,bj) =
338    CADJ &     comlev1_bibj, key=iicekey, byte=isbyte
339    # endif
340          DO j = jMin, jMax
341           DO i = iMin, iMax
342  C--     Compute Sea-Ice Loading (= mass of sea-ice + snow / area unit)  C--     Compute Sea-Ice Loading (= mass of sea-ice + snow / area unit)
343          sIceLoad(i,j,bi,bj) = ( snowHeight(i,j,bi,bj)*rhos          sIceLoad(i,j,bi,bj) = ( snowHeight(i,j,bi,bj)*rhos
344       &                         + iceHeight(i,j,bi,bj)*rhoi       &                         + iceHeight(i,j,bi,bj)*rhoi
345       &                        )*iceMask(i,j,bi,bj)       &                        )*iceMask(i,j,bi,bj)
 #endif  
   
346         ENDDO         ENDDO
347        ENDDO        ENDDO
348    #endif
349    
350        IF ( thSIceAdvScheme.GT.0 ) THEN        IF ( thSIceAdvScheme.GT.0 ) THEN
351  C--   note: those fluxes should to be added directly to Qnet, EmPmR & saltFlux  C--   note: those fluxes should to be added directly to Qnet, EmPmR & saltFlux
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22