/[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.1 by adcroft, Mon Mar 4 19:01:29 2002 UTC revision 1.7 by jmc, Thu Mar 8 17:13:31 2012 UTC
# Line 10  C !INTERFACE: ========================== Line 10  C !INTERFACE: ==========================
10        SUBROUTINE PTRACERS_CONVECT( bi,bj,k,weightA,weightB,myThid )        SUBROUTINE PTRACERS_CONVECT( bi,bj,k,weightA,weightB,myThid )
11    
12  C !DESCRIPTION:  C !DESCRIPTION:
13  C     Calculates tendancy for passive tracers and integrates forward  C     do passive tracers convection
 C     in time.  
14    
15  C !USES: ===============================================================  C !USES: ===============================================================
16        IMPLICIT NONE        IMPLICIT NONE
17  #include "SIZE.h"  #include "SIZE.h"
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PTRACERS.h"  #include "PTRACERS_SIZE.h"
20    #include "PTRACERS_PARAMS.h"
21    #include "PTRACERS_START.h"
22    #include "PTRACERS_FIELDS.h"
23    #ifdef ALLOW_AUTODIFF_TAMC
24    # include "tamc.h"
25    # include "tamc_keys.h"
26    #endif
27    
28  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
29  C  bi,bj,k :: tile indices  C  bi,bj,k :: tile indices
# Line 41  CEOP Line 47  CEOP
47    
48  C Loop over tracers  C Loop over tracers
49        DO iTracer=1,PTRACERS_numInUse        DO iTracer=1,PTRACERS_numInUse
50           IF ( PTRACERS_StepFwd(iTracer) ) THEN
51    
52         CALL CONVECTIVELY_MIXTRACER(  #ifdef ALLOW_AUTODIFF_TAMC
53       I         bi,bj,k,weightA,weightB,            act0 = iTracer - 1
54       U         pTracer(1-Olx,1-Oly,1,bi,bj,iTracer),            max0 = PTRACERS_num
55       I         myThid)            act1 = bi - myBxLo(myThid)
56              max1 = myBxHi(myThid) - myBxLo(myThid) + 1
57              act2 = bj - myByLo(myThid)
58              max2 = myByHi(myThid) - myByLo(myThid) + 1
59              act3 = myThid - 1
60              max3 = nTx*nTy
61              act4 = ikey_dynamics - 1
62              iptrkey = (act0 + 1)
63         &                      + act1*max0
64         &                      + act2*max0*max1
65         &                      + act3*max0*max1*max2
66         &                      + act4*max0*max1*max2*max3
67              kkey = (iptrkey-1)*Nr + k
68    CADJ STORE pTracer(:,:,k,bi,bj,iTracer)
69    CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
70    #endif /* ALLOW_AUTODIFF_TAMC */
71    
72             CALL CONVECTIVELY_MIXTRACER(
73         I           bi,bj,k,weightA,weightB,
74         U           pTracer(1-OLx,1-OLy,1,1,1,iTracer),
75         I           myThid )
76    
77  C End of tracer loop  C End of tracer loop
78           ENDIF
79        ENDDO        ENDDO
80    
81  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22