/[MITgcm]/MITgcm/pkg/ptracers/ptracers_convect.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_convect.F

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

revision 1.5 by jmc, Mon Nov 5 18:48:04 2007 UTC revision 1.6 by heimbach, Tue Jan 18 19:37:41 2011 UTC
# Line 19  C !USES: =============================== Line 19  C !USES: ===============================
19  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
20  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
21  #include "PTRACERS_FIELDS.h"  #include "PTRACERS_FIELDS.h"
22    #ifdef ALLOW_AUTODIFF_TAMC
23    # include "tamc.h"
24    # include "tamc_keys.h"
25    #endif
26    
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
28  C  bi,bj,k :: tile indices  C  bi,bj,k :: tile indices
# Line 43  CEOP Line 47  CEOP
47  C Loop over tracers  C Loop over tracers
48        DO iTracer=1,PTRACERS_numInUse        DO iTracer=1,PTRACERS_numInUse
49    
50    #ifdef ALLOW_AUTODIFF_TAMC
51              act0 = iTracer - 1
52              max0 = PTRACERS_num
53              act1 = bi - myBxLo(myThid)
54              max1 = myBxHi(myThid) - myBxLo(myThid) + 1
55              act2 = bj - myByLo(myThid)
56              max2 = myByHi(myThid) - myByLo(myThid) + 1
57              act3 = myThid - 1
58              max3 = nTx*nTy
59              act4 = ikey_dynamics - 1
60              iptrkey = (act0 + 1)
61         &                      + act1*max0
62         &                      + act2*max0*max1
63         &                      + act3*max0*max1*max2
64         &                      + act4*max0*max1*max2*max3
65              kkey = (iptrkey-1)*Nr + k
66    CADJ STORE pTracer(:,:,k,bi,bj,iTracer)
67    CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
68    #endif /* ALLOW_AUTODIFF_TAMC */
69    
70         CALL CONVECTIVELY_MIXTRACER(         CALL CONVECTIVELY_MIXTRACER(
71       I         bi,bj,k,weightA,weightB,       I         bi,bj,k,weightA,weightB,
72       U         pTracer(1-Olx,1-Oly,1,1,1,iTracer),       U         pTracer(1-Olx,1-Oly,1,1,1,iTracer),
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22