model/src/timestep_tracer.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/timestep_tracer.F,v 1.6 2001/09/04 14:53:10 adcroft Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

      SUBROUTINE TIMESTEP_TRACER(
     I                     bi, bj, iMin, iMax, jMin, jMax, K,
     I                     advection_scheme,
     I                     tracer, gTracer,
     O                     gTrNm1,
     I                     myIter, myThid )
C     /==========================================================\
C     | S/R TIMESTEP_TRACER                                      |
C     | o Step model tracer field forward in time                |
C     \==========================================================/
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GAD.h"
#include "GRID.h"
#include "SURFACE.h"

C     == Routine Arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax,K
      INTEGER advection_scheme
      _RL  tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL  gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL  gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myIter, myThid
C     == Local variables ==
      INTEGER i,j
      _RL ab15,ab05
      _RL gTrtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C     Adams-Bashforth timestepping weights
      IF (myIter .EQ. 0) THEN
       ab15=1.0
       ab05=0.0
      ELSE
       ab15=1.5+abEps
       ab05=-(0.5+abEps)
      ENDIF

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

C-    Step forward tracer
      DO j=jMin,jMax
       DO i=iMin,iMax
         gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj)
     &            +deltaTtracer*gTracer(i,j,k,bi,bj)
       ENDDO
      ENDDO

      RETURN
      END
1	adcroft	1.7	C $Header: /u/gcmpack/models/MITgcmUV/model/src/timestep_tracer.F,v 1.6 2001/09/04 14:53:10 adcroft Exp $
2	adcroft	1.2	C $Name: $
3	adcroft	1.1
4			#include "CPP_OPTIONS.h"
5
6	adcroft	1.2	SUBROUTINE TIMESTEP_TRACER(
7	adcroft	1.4	I bi, bj, iMin, iMax, jMin, jMax, K,
8			I advection_scheme,
9			I tracer, gTracer,
10			O gTrNm1,
11	adcroft	1.2	I myIter, myThid )
12	adcroft	1.1	C /==========================================================\
13			C \| S/R TIMESTEP_TRACER \|
14			C \| o Step model tracer field forward in time \|
15			C \==========================================================/
16	adcroft	1.2	IMPLICIT NONE
17
18			C == Global variables ===
19	adcroft	1.1	#include "SIZE.h"
20			#include "EEPARAMS.h"
21			#include "PARAMS.h"
22	adcroft	1.4	#include "GAD.h"
23	jmc	1.5	#include "GRID.h"
24			#include "SURFACE.h"
25	adcroft	1.2
26	adcroft	1.1	C == Routine Arguments ==
27			INTEGER bi,bj,iMin,iMax,jMin,jMax,K
28	adcroft	1.4	INTEGER advection_scheme
29	adcroft	1.1	_RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
30			_RL gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
31			_RL gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
32			INTEGER myIter, myThid
33			C == Local variables ==
34			INTEGER i,j
35			_RL ab15,ab05
36	jmc	1.5	_RL gTrtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37	adcroft	1.1
38			C Adams-Bashforth timestepping weights
39	adcroft	1.3	IF (myIter .EQ. 0) THEN
40			ab15=1.0
41			ab05=0.0
42			ELSE
43	adcroft	1.4	ab15=1.5+abEps
44			ab05=-(0.5+abEps)
45	adcroft	1.3	ENDIF
46	adcroft	1.1
47	jmc	1.5	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
48
49			C- Step forward tracer
50			DO j=jMin,jMax
51			DO i=iMin,iMax
52	adcroft	1.7	gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj)
53			& +deltaTtracer*gTracer(i,j,k,bi,bj)
54	jmc	1.5	ENDDO
55			ENDDO
56	adcroft	1.1
57			RETURN
58			END