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

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

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


Revision 1.15 - (show annotations) (download)
Fri May 20 16:34:19 2011 UTC (13 years ago) by jmc
Branch: MAIN
Changes since 1.14: +15 -8 lines
leave tracer unchanged outside OB interior region: This has no effect
 on the solution but just to prevent unrealistic tracer value ouside OB.

1 C $Header: /u/gcmpack/MITgcm/model/src/timestep_tracer.F,v 1.14 2009/06/26 23:10:09 jahn Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: TIMESTEP_TRACER
9 C !INTERFACE:
10 SUBROUTINE TIMESTEP_TRACER(
11 I bi, bj, iMin, iMax, jMin, jMax, k,
12 I advection_scheme, deltaTloc,
13 I tracer, gTracer,
14 I myIter, myThid )
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | S/R TIMESTEP_TRACER
18 C | o Step model tracer field forward in time
19 C *==========================================================*
20 C \ev
21
22 C !USES:
23 IMPLICIT NONE
24 C == Global variables ===
25 #include "SIZE.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28 #include "GRID.h"
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C == Routine Arguments ==
32 INTEGER bi,bj,iMin,iMax,jMin,jMax,k
33 INTEGER advection_scheme
34 _RL deltaTloc
35 _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
36 _RL gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
37 INTEGER myIter, myThid
38
39 C !LOCAL VARIABLES:
40 C == Local variables ==
41 INTEGER i,j
42 CEOP
43
44 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
45
46 C- Step tracer forward in time and store provisional value in gTracer array
47 c DO j=jMin,jMax
48 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)
52 & +deltaTloc*gTracer(i,j,k,bi,bj)
53 #ifdef ALLOW_OBCS
54 & *maskInC(i,j,bi,bj)
55 #endif
56 ENDDO
57 ENDDO
58
59 RETURN
60 END

  ViewVC Help
Powered by ViewVC 1.1.22