/[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.38 by jahn, Fri Jun 26 23:10:10 2009 UTC revision 1.42 by heimbach, Tue Jan 18 22:58:20 2011 UTC
# Line 128  CADJ STORE pTracer(:,:,k,bi,bj,iTracer) Line 128  CADJ STORE pTracer(:,:,k,bi,bj,iTracer)
128  CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte  CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
129  CADJ STORE gpTrNm1(:,:,k,bi,bj,iTracer)  CADJ STORE gpTrNm1(:,:,k,bi,bj,iTracer)
130  CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte  CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
131    # ifdef NONLIN_FRSURF
132    CADJ STORE gpTr(:,:,k,bi,bj,iTracer)
133    CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
134    CADJ STORE rFlx(:,:,:,iTracer)
135    CADJ &      = comlev1_bibj_k_ptracers, key=kkey, byte=isbyte
136    # endif
137  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
138    
139  C Calculate active tracer tendencies (gPtr) due to internal processes  C Calculate active tracer tendencies (gPtr) due to internal processes
# Line 173  C  prevent gPtr from being replaced by t Line 179  C  prevent gPtr from being replaced by t
179  C       compute iter at beginning of ptracer time step  C       compute iter at beginning of ptracer time step
180  #ifdef ALLOW_LONGSTEP  #ifdef ALLOW_LONGSTEP
181          iterNb = myIter - LS_nIter + 1          iterNb = myIter - LS_nIter + 1
182          IF (LS_staggerTimeStep) iterNb = myIter - LS_nIter          IF (LS_whenToSample.GE.2) iterNb = myIter - LS_nIter
183  #else  #else
184          iterNb = myIter          iterNb = myIter
185          IF (staggerTimeStep) iterNb = myIter - 1          IF (staggerTimeStep) iterNb = myIter - 1
# Line 220  C Integrate forward in time, storing in Line 226  C Integrate forward in time, storing in
226       I                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),       I                        gPtr(1-Olx,1-Oly,1,1,1,iTracer),
227       I                        myIter,myThid )       I                        myIter,myThid )
228    
 #ifdef   ALLOW_OBCS  
 C Apply open boundary conditions  
          IF (useOBCS) THEN  
           CALL OBCS_APPLY_PTRACER(  
      I         bi, bj, k, iTracer,  
      U         gPtr(1-Olx,1-Oly,k,bi,bj,iTracer),  
      I         myThid )  
          ENDIF  
 #endif   /* ALLOW_OBCS */  
229  C end of tracer loop  C end of tracer loop
230        ENDDO        ENDDO
231    
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.42
  ViewVC Help
Powered by ViewVC 1.1.22