/[MITgcm]/MITgcm/model/src/timestep_tracer.F
ViewVC logotype

Diff of /MITgcm/model/src/timestep_tracer.F

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

revision 1.14 by jahn, Fri Jun 26 23:10:09 2009 UTC revision 1.15 by jmc, Fri May 20 16:34:19 2011 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CBOP  CBOP
8  C     !ROUTINE: TIMESTEP_TRACER  C     !ROUTINE: TIMESTEP_TRACER
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE TIMESTEP_TRACER(        SUBROUTINE TIMESTEP_TRACER(
11       I                     bi, bj, iMin, iMax, jMin, jMax, K,       I                     bi, bj, iMin, iMax, jMin, jMax, k,
12       I                     advection_scheme, deltaTloc,       I                     advection_scheme, deltaTloc,
13       I                     tracer, gTracer,       I                     tracer, gTracer,
14       I                     myIter, myThid )       I                     myIter, myThid )
15  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
16  C     *==========================================================*  C     *==========================================================*
17  C     | S/R TIMESTEP_TRACER                                        C     | S/R TIMESTEP_TRACER
18  C     | o Step model tracer field forward in time                  C     | o Step model tracer field forward in time
19  C     *==========================================================*  C     *==========================================================*
20  C     \ev  C     \ev
21    
# Line 24  C     == Global variables === Line 25  C     == Global variables ===
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
28    #include "GRID.h"
29    
30  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
31  C     == Routine Arguments ==  C     == Routine Arguments ==
32        INTEGER bi,bj,iMin,iMax,jMin,jMax,K        INTEGER bi,bj,iMin,iMax,jMin,jMax,k
33        INTEGER advection_scheme        INTEGER advection_scheme
34        _RL     deltaTloc        _RL     deltaTloc
35        _RL  tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
# Line 41  CEOP Line 43  CEOP
43    
44  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
45    
46  C-    Extrapolate tendencies to mid-time level (n+1/2)  C-    Step tracer forward in time and store provisional value in gTracer array
47        DO j=jMin,jMax  c     DO j=jMin,jMax
48         DO i=iMin,iMax  c      DO i=iMin,iMax
49          DO j=1-Oly,sNy+Oly
50           DO i=1-Olx,sNx+Olx
51           gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj)           gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj)
52       &            +deltaTloc*gTracer(i,j,k,bi,bj)       &            +deltaTloc*gTracer(i,j,k,bi,bj)
53    #ifdef ALLOW_OBCS
54         &                      *maskInC(i,j,bi,bj)
55    #endif
56         ENDDO         ENDDO
57        ENDDO        ENDDO
58    
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22