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

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

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

revision 1.7 by heimbach, Fri Jun 27 01:56:17 2003 UTC revision 1.10 by stephd, Mon Oct 6 20:25:54 2003 UTC
# Line 70  C  iMin,iMax,jMin,jMax  :: loop ranges Line 70  C  iMin,iMax,jMin,jMax  :: loop ranges
70  C  kUp,kDown            :: toggle indices for even/odd level fluxes  C  kUp,kDown            :: toggle indices for even/odd level fluxes
71  C  km1                  :: =min(1,k-1)  C  km1                  :: =min(1,k-1)
72  C  rFlx                 :: vertical flux  C  rFlx                 :: vertical flux
73    C  GAD_TR               :: passive tracer id (GAD_TR1+iTracer-1)
74        INTEGER i,j,iTracer        INTEGER i,j,iTracer
75        INTEGER iMin,iMax,jMin,jMax        INTEGER iMin,iMax,jMin,jMax
76        INTEGER kUp,kDown,km1        INTEGER kUp,kDown,km1
77        _RL rFlx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,2,PTRACERS_num)        _RL rFlx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,2,PTRACERS_num)
78          INTEGER GAD_TR
79        LOGICAL calcAdvection        LOGICAL calcAdvection
80  CEOP  CEOP
81    
# Line 106  C (advection, [explicit] diffusion, para Line 108  C (advection, [explicit] diffusion, para
108       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND
109       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD
110       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_4TH       &      .OR. PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_4TH
111           GAD_TR = GAD_TR1 + iTracer - 1
112           CALL GAD_CALC_RHS(           CALL GAD_CALC_RHS(
113       I                     bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,       I                     bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
114       I                     xA,yA,uTrans,vTrans,rTrans,maskUp,       I                     xA,yA,uTrans,vTrans,rTrans,maskUp,
# Line 113  C (advection, [explicit] diffusion, para Line 116  C (advection, [explicit] diffusion, para
116       I                     PTRACERS_diffK4(iTracer),       I                     PTRACERS_diffK4(iTracer),
117       I                     KappaRtr,       I                     KappaRtr,
118       I                     pTracer(1-Olx,1-Oly,1,1,1,iTracer),       I                     pTracer(1-Olx,1-Oly,1,1,1,iTracer),
119       I                     GAD_TR1,       I                     GAD_TR,
120       I                     PTRACERS_advScheme(iTracer),calcAdvection,       I                     PTRACERS_advScheme(iTracer),calcAdvection,
121       U                     rFlx(1-Olx,1-Oly,1,iTracer),       U                     rFlx(1-Olx,1-Oly,1,iTracer),
122       U                     gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                     gPtr(1-Olx,1-Oly,1,1,1,iTracer),
123       I                     myThid )       I                     myThid )
124    
125  C External forcing term(s)  C External forcing term(s)
126  cswdptr - add--  cswdptr --add --
127  #ifndef PTRACERS_SEPERATE_FORCING  #ifdef ALLOW_GCHEM
128    #ifndef PTRACERS_SEPARATE_FORCING
129           IF ( forcing_In_AB )
130         &    CALL GCHEM_FORCING_INT(
131         I                        bi,bj,iMin,iMax,jMin,jMax,k,
132         I                        iTracer,
133         I                        myTime,myIter, myThid)
134    #endif
135    #else
136  cswdptr - end add ---  cswdptr - end add ---
137         IF ( forcing_In_AB )         IF ( forcing_In_AB )
138       &   CALL PTRACERS_FORCING(       &   CALL PTRACERS_FORCING(
139       I                        bi,bj,k,iTracer,       I                        bi,bj,iMin,iMax,jMin,jMax,k,iTracer,
140       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),
141       I                        myIter,myTime,myThid)       I                        myIter,myTime,myThid)
142  cswdptr --add---  cswdptr --add---
# Line 145  C If using Adams-Bashforth II, then extr Line 156  C If using Adams-Bashforth II, then extr
156    
157  C External forcing term(s)  C External forcing term(s)
158  cswdptr - add--  cswdptr - add--
159  #ifndef PTRACERS_SEPERATE_FORCING  #ifdef ALLOW_GCHEM
160    #ifndef PTRACERS_SEPARATE_FORCING
161           IF ( .NOT.forcing_In_AB )
162         &    CALL GCHEM_FORCING_INT(
163         I                        bi,bj,iMin,iMax,jMin,jMax,k,
164         I                        iTracer,
165         I                        myTime,myIter, myThid)
166    #endif
167    #else
168  cswdptr - end add ---  cswdptr - end add ---
169         IF ( .NOT.forcing_In_AB )         IF ( .NOT.forcing_In_AB )
170       &   CALL PTRACERS_FORCING(       &   CALL PTRACERS_FORCING(
171       I                        bi,bj,k,iTracer,       I                        bi,bj,iMin,iMax,jMin,jMax,k,iTracer,
172       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),
173       I                        myIter,myTime,myThid)       I                        myIter,myTime,myThid)
174  cswdptr - add--  cswdptr - add--
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22