/[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.3 by jmc, Sat Jun 15 03:37:44 2002 UTC revision 1.7 by heimbach, Fri Jun 27 01:56:17 2003 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "PTRACERS_OPTIONS.h"  #include "PTRACERS_OPTIONS.h"
5    cswdptr -- add ---
6    #ifdef ALLOW_GCHEM
7    # include "GCHEM_OPTIONS.h"
8    #endif
9    cswdptr -- end add ---
10    
11  CBOP  CBOP
12  C !ROUTINE: PTRACERS_INTEGERATE  C !ROUTINE: PTRACERS_INTEGRATE
13    
14  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
15        SUBROUTINE PTRACERS_INTEGERATE(        SUBROUTINE PTRACERS_INTEGRATE(
16       I                               bi,bj,k,       I                               bi,bj,k,
17       I                               xA,yA,uTrans,vTrans,rTrans,maskUp,       I                               xA,yA,uTrans,vTrans,rTrans,maskUp,
18       X                               KappaRtr,       X                               KappaRtr,
# Line 21  C !USES: =============================== Line 26  C !USES: ===============================
26        IMPLICIT NONE        IMPLICIT NONE
27  #include "SIZE.h"  #include "SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29    #include "PARAMS.h"
30  #include "PTRACERS.h"  #include "PTRACERS.h"
31  #include "GAD.h"  #include "GAD.h"
32    
# Line 114  C (advection, [explicit] diffusion, para Line 120  C (advection, [explicit] diffusion, para
120       I                     myThid )       I                     myThid )
121    
122  C External forcing term(s)  C External forcing term(s)
123    cswdptr - add--
124    #ifndef PTRACERS_SEPERATE_FORCING
125    cswdptr - end add ---
126         IF ( forcing_In_AB )         IF ( forcing_In_AB )
127       &   CALL PTRACERS_FORCING(       &   CALL PTRACERS_FORCING(
128       I                        bi,bj,k,iTracer,       I                        bi,bj,k,iTracer,
129       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),
130       I                        myIter,myTime,myThid)       I                        myIter,myTime,myThid)
131    cswdptr --add---
132    #endif
133    cswdptr -- end add ---
134    
135  C If using Adams-Bashforth II, then extrapolate tendancies  C If using Adams-Bashforth II, then extrapolate tendancies
136        IF ( PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND        IF ( PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND
# Line 132  C If using Adams-Bashforth II, then extr Line 144  C If using Adams-Bashforth II, then extr
144        ENDIF        ENDIF
145    
146  C External forcing term(s)  C External forcing term(s)
147    cswdptr - add--
148    #ifndef PTRACERS_SEPERATE_FORCING
149    cswdptr - end add ---
150         IF ( .NOT.forcing_In_AB )         IF ( .NOT.forcing_In_AB )
151       &   CALL PTRACERS_FORCING(       &   CALL PTRACERS_FORCING(
152       I                        bi,bj,k,iTracer,       I                        bi,bj,k,iTracer,
153       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),
154       I                        myIter,myTime,myThid)       I                        myIter,myTime,myThid)
155    cswdptr - add--
156    #endif
157    cswdptr -- end add ---
158    
159  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
160  C Account for change in level thickness  C Account for change in level thickness
# Line 145  C Account for change in level thickness Line 163  C Account for change in level thickness
163       I                          bi,bj,K,       I                          bi,bj,K,
164       U                          gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U                          gPtr(1-Olx,1-Oly,1,1,1,iTracer),
165       I                          myThid )       I                          myThid )
166            IF ( PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND
167         &   .OR.PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD
168         &   .OR.PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_4TH )
169         &  CALL FREESURF_RESCALE_G(
170         I                          bi,bj,K,
171         U                          gPtrNm1(1-Olx,1-Oly,1,1,1,iTracer),
172         I                          myThid )
173        ENDIF        ENDIF
174  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
175    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22