--- MITgcm/pkg/ptracers/ptracers_advection.F	2005/10/10 05:53:48	1.7
+++ MITgcm/pkg/ptracers/ptracers_advection.F	2006/04/04 14:52:43	1.8
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/ptracers/Attic/ptracers_advection.F,v 1.7 2005/10/10 05:53:48 mlosch Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/ptracers/Attic/ptracers_advection.F,v 1.8 2006/04/04 14:52:43 heimbach Exp $
 C $Name:  $
 
 #include "PTRACERS_OPTIONS.h"
@@ -22,6 +22,10 @@
 #include "PTRACERS_SIZE.h"
 #include "PTRACERS.h"
 #include "GAD.h"
+#ifdef ALLOW_AUTODIFF_TAMC
+# include "tamc.h"
+# include "tamc_keys.h"
+#endif
 
 C !INPUT PARAMETERS: ===================================================
 C  bi,bj                :: tile indices
@@ -50,6 +54,28 @@
 C Loop over tracers
       DO iTracer=1,PTRACERS_numInUse
 
+#ifdef ALLOW_AUTODIFF_TAMC
+          act0 = iTracer - 1
+          max0 = PTRACERS_num
+          act1 = bi - myBxLo(myThid)
+          max1 = myBxHi(myThid) - myBxLo(myThid) + 1
+          act2 = bj - myByLo(myThid)
+          max2 = myByHi(myThid) - myByLo(myThid) + 1
+          act3 = myThid - 1
+          max3 = nTx*nTy
+          act4 = ikey_dynamics - 1
+          iptrkey = (act0 + 1)
+     &                      + act1*max0
+     &                      + act2*max0*max1
+     &                      + act3*max0*max1*max2
+     &                      + act4*max0*max1*max2*max3
+#endif /* ALLOW_AUTODIFF_TAMC */
+
+#ifdef ALLOW_AUTODIFF_TAMC
+CADJ STORE ptracer(:,:,:,bi,bj,iTracer)
+CADJ &      = comlev1_bibj_ptracers, key=iptrkey, byte=isbyte
+#endif /* ALLOW_AUTODIFF_TAMC */
+ 
        IF (
      &     multiDimAdvection .AND.
      &     PTRACERS_advScheme(iTracer).NE.ENUM_CENTERED_2ND .AND.