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

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

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

revision 1.7 by jmc, Sat Jan 2 23:42:51 2010 UTC revision 1.8 by jmc, Sun Nov 14 23:32:15 2010 UTC
# Line 14  C !INTERFACE: ========================== Line 14  C !INTERFACE: ==========================
14  C !DESCRIPTION:  C !DESCRIPTION:
15  C     Calls the implicit vertical advection/diffusion routine  C     Calls the implicit vertical advection/diffusion routine
16  C     for each passive tracer.  C     for each passive tracer.
17    C     And apply open boundary conditions for each passive tracer
18    C Note: would be better to apply OBC in a dedicated S/R ;
19    C     can be done here assuming no other contribution modify passive
20    C     tracer after implicit vertical diffus/advect is applied.
21    
22  C !USES: ===============================================================  C !USES: ===============================================================
23        IMPLICIT NONE        IMPLICIT NONE
# Line 56  C  iTracer :: tracer index Line 60  C  iTracer :: tracer index
60        INTEGER iTracer        INTEGER iTracer
61        INTEGER iMin,iMax,jMin,jMax        INTEGER iMin,iMax,jMin,jMax
62        INTEGER GAD_TR        INTEGER GAD_TR
 #ifdef ALLOW_OBCS  
       INTEGER k  
 #endif  
63  CEOP  CEOP
64    
65        iMin=0        iMin=0
# Line 129  CADJ &     key=iptrkey, byte=isbyte Line 130  CADJ &     key=iptrkey, byte=isbyte
130       I         GAD_TR, kappaRk, recip_HFacC,       I         GAD_TR, kappaRk, recip_HFacC,
131       U         gPtr(1-Olx,1-Oly,1,1,1,iTracer),       U         gPtr(1-Olx,1-Oly,1,1,1,iTracer),
132       I         myThid )       I         myThid )
133            ENDIF
134    
135  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
136  C--      Apply open boundary conditions  C--     Apply open boundary conditions
137          IF ( useOBCS ) THEN          IF ( useOBCS ) THEN
          DO k=1,Nr  
138            CALL OBCS_APPLY_PTRACER(            CALL OBCS_APPLY_PTRACER(
139       I         bi, bj, k, iTracer,       I         bi, bj, 0, iTracer,
140       U         gPtr(1-Olx,1-Oly,k,bi,bj,iTracer),       U         gPtr(1-Olx,1-Oly,1,bi,bj,iTracer),
141       I         myThid )       I         myThid )
          ENDDO  
142          ENDIF          ENDIF
143  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
144    
         ENDIF  
145  C End of tracer loop  C End of tracer loop
146        ENDDO        ENDDO
147    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22