/[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.29 by jmc, Thu Oct 21 19:03:31 2010 UTC revision 1.30 by gforget, Fri Dec 17 04:00:14 2010 UTC
# Line 116  C---+----1----+----2----+----3----+----4 Line 116  C---+----1----+----2----+----3----+----4
116        act3 = myThid - 1        act3 = myThid - 1
117        max3 = nTx*nTy        max3 = nTx*nTy
118        act4 = ikey_dynamics - 1        act4 = ikey_dynamics - 1
119        iicekey = (act1 + 1) + act2*max1        ticekey = (act1 + 1) + act2*max1
120       &                     + act3*max1*max2       &                     + act3*max1*max2
121       &                     + act4*max1*max2*max3       &                     + act4*max1*max2*max3
122  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
# Line 141  C-    Initialise Line 141  C-    Initialise
141        snowFac = thSIce_deltaT/(rhos*hNewSnowAge)        snowFac = thSIce_deltaT/(rhos*hNewSnowAge)
142    
143  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
144  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
145  CADJ STORE iceheight(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE iceheight(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
146  CADJ STORE icfrwatm(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE icfrwatm(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
147  CADJ STORE qice1(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE qice1(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
148  CADJ STORE qice2(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE qice2(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
149  CADJ STORE snowheight(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE snowheight(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
150  #endif  #endif
151        DO j = jMin, jMax        DO j = jMin, jMax
152         DO i = iMin, iMax         DO i = iMin, iMax
# Line 223  C------- Line 223  C-------
223        ENDDO        ENDDO
224    
225  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
226  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
227  #endif  #endif
228    
229        CALL THSICE_CALC_THICKN(        CALL THSICE_CALC_THICKN(
# Line 241  CADJ STORE iceMask(:,:,bi,bj) = comlev1_ Line 241  CADJ STORE iceMask(:,:,bi,bj) = comlev1_
241       I          myTime, myIter, myThid )       I          myTime, myIter, myThid )
242    
243  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
244  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=ticekey,byte=isbyte
245  CADJ STORE fsalt(:,:)  = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE fsalt(:,:)  = comlev1_bibj,key=ticekey,byte=isbyte
246  CADJ STORE flx2oc(:,:) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE flx2oc(:,:) = comlev1_bibj,key=ticekey,byte=isbyte
247  CADJ STORE frw2oc(:,:) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE frw2oc(:,:) = comlev1_bibj,key=ticekey,byte=isbyte
248  #endif  #endif
249  C--    Net fluxes :  C--    Net fluxes :
250        DO j = jMin, jMax        DO j = jMin, jMax
# Line 320  C------- Line 320  C-------
320    
321  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
322  CADJ STORE snowHeight(:,:,bi,bj) =  CADJ STORE snowHeight(:,:,bi,bj) =
323  CADJ &     comlev1_bibj,key=iicekey,byte=isbyte  CADJ &     comlev1_bibj,key=ticekey,byte=isbyte
324  #endif  #endif
325        DO j = jMin, jMax        DO j = jMin, jMax
326         DO i = iMin, iMax         DO i = iMin, iMax
# Line 367  C--    Net fluxes : Line 367  C--    Net fluxes :
367    
368  # ifdef ALLOW_AUTODIFF_TAMC  # ifdef ALLOW_AUTODIFF_TAMC
369  CADJ STORE snowHeight(:,:,bi,bj) =  CADJ STORE snowHeight(:,:,bi,bj) =
370  CADJ &     comlev1_bibj,key=iicekey,byte=isbyte  CADJ &     comlev1_bibj,key=ticekey,byte=isbyte
371  # endif  # endif
372        DO j = jMin, jMax        DO j = jMin, jMax
373         DO i = iMin, iMax         DO i = iMin, iMax
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22