/[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.28 by heimbach, Sat Oct 16 19:22:34 2010 UTC revision 1.29 by jmc, Thu Oct 21 19:03:31 2010 UTC
# Line 240  CADJ STORE iceMask(:,:,bi,bj) = comlev1_ Line 240  CADJ STORE iceMask(:,:,bi,bj) = comlev1_
240       O          frw2oc, fsalt,       O          frw2oc, fsalt,
241       I          myTime, myIter, myThid )       I          myTime, myIter, myThid )
242    
243    #ifdef ALLOW_AUTODIFF_TAMC
244    CADJ STORE iceMask(:,:,bi,bj) = comlev1_bibj,key=iicekey,byte=isbyte
245    CADJ STORE fsalt(:,:)  = comlev1_bibj,key=iicekey,byte=isbyte
246    CADJ STORE flx2oc(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
247    CADJ STORE frw2oc(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
248    #endif
249  C--    Net fluxes :  C--    Net fluxes :
250        DO j = jMin, jMax        DO j = jMin, jMax
251         DO i = iMin, iMax         DO i = iMin, iMax
252  #ifdef ALLOW_AUTODIFF_TAMC  c#ifdef ALLOW_AUTODIFF_TAMC
253            ikey_1 = i  c         ikey_1 = i
254       &         + sNx*(j-1)  c    &         + sNx*(j-1)
255       &         + sNx*sNy*act1  c    &         + sNx*sNy*act1
256       &         + sNx*sNy*max1*act2  c    &         + sNx*sNy*max1*act2
257       &         + sNx*sNy*max1*max2*act3  c    &         + sNx*sNy*max1*max2*act3
258       &         + sNx*sNy*max1*max2*max3*act4  c    &         + sNx*sNy*max1*max2*max3*act4
259  #endif /* ALLOW_AUTODIFF_TAMC */  c#endif /* ALLOW_AUTODIFF_TAMC */
260  C--  c#ifdef ALLOW_AUTODIFF_TAMC
261  #ifdef ALLOW_AUTODIFF_TAMC  cCADJ STORE  icemask(i,j,bi,bj) = comlev1_thsice_1, key=ikey_1
262  CADJ STORE  icemask(i,j,bi,bj) = comlev1_thsice_1, key=ikey_1  c#endif
 #endif  
263          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN          IF (iceMask(i,j,bi,bj).GT.0. _d 0) THEN
264  C-     weighted average net fluxes:  C-     weighted average net fluxes:
265  #ifdef ALLOW_AUTODIFF_TAMC  c#ifdef ALLOW_AUTODIFF_TAMC
266  CADJ STORE  fsalt(i,j) = comlev1_thsice_1, key=ikey_1  cCADJ STORE  fsalt(i,j) = comlev1_thsice_1, key=ikey_1
267  CADJ STORE  flx2oc(i,j) = comlev1_thsice_1, key=ikey_1  cCADJ STORE  flx2oc(i,j) = comlev1_thsice_1, key=ikey_1
268  CADJ STORE  frw2oc(i,j) = comlev1_thsice_1, key=ikey_1  cCADJ STORE  frw2oc(i,j) = comlev1_thsice_1, key=ikey_1
269  CADJ STORE  icemask(i,j,bi,bj) = comlev1_thsice_1, key=ikey_1  cCADJ STORE  icemask(i,j,bi,bj) = comlev1_thsice_1, key=ikey_1
270  #endif  c#endif
271            icFrac = iceMask(i,j,bi,bj)            icFrac = iceMask(i,j,bi,bj)
272            opFrac= 1. _d 0-icFrac            opFrac= 1. _d 0-icFrac
273  #ifdef ALLOW_ATM2D  #ifdef ALLOW_ATM2D
# Line 314  C------- Line 319  C-------
319       I          myTime, myIter, myThid )       I          myTime, myIter, myThid )
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=iicekey,byte=isbyte
324  #endif  #endif
325        DO j = jMin, jMax        DO j = jMin, jMax
# Line 361  C--    Net fluxes : Line 366  C--    Net fluxes :
366        ENDDO        ENDDO
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=iicekey,byte=isbyte
371  # endif  # endif
372        DO j = jMin, jMax        DO j = jMin, jMax
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22