/[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.1 by adcroft, Fri Feb 2 21:35:19 2001 UTC revision 1.2 by adcroft, Tue May 29 14:01:37 2001 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6          SUBROUTINE TIMESTEP_TRACER(
7         I                     bi, bj, iMin, iMax, jMin, jMax, K, tauAB,
8         U                     tracer, gTracer, gTrNm1,
9         I                     myIter, myThid )
10  C     /==========================================================\  C     /==========================================================\
11  C     | S/R TIMESTEP_TRACER                                      |  C     | S/R TIMESTEP_TRACER                                      |
12  C     | o Step model tracer field forward in time                |  C     | o Step model tracer field forward in time                |
13  C     \==========================================================/  C     \==========================================================/
14        SUBROUTINE TIMESTEP_TRACER(        IMPLICIT NONE
15       I                     bi, bj, iMin, iMax, jMin, jMax, K,  
16       U                     tracer, gTracer, gTrNm1,  C     == Global variables ===
      I                     myIter, myThid )  
       implicit none  
 ! Common  
17  #include "SIZE.h"  #include "SIZE.h"
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PARAMS.h"  #include "PARAMS.h"
20    
21  C     == Routine Arguments ==  C     == Routine Arguments ==
22    C     tauAB - Adams-Bashforth timestepping weight: 0=forward ; 1/2=Adams-Bashf.
23        INTEGER bi,bj,iMin,iMax,jMin,jMax,K        INTEGER bi,bj,iMin,iMax,jMin,jMax,K
24          _RL  tauAB
25        _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)
26        _RL  gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
27        _RL  gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
# Line 27  C     == Local variables == Line 31  C     == Local variables ==
31        _RL ab15,ab05        _RL ab15,ab05
32    
33  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
34  Caja  IF (myIter .EQ. 0) THEN        ab15 = 1.+ tauAB
35  Caja   ab15=1.0        ab05 = -tauAB
 Caja   ab05=0.0  
 Caja  ELSE  
        ab15=1.5+abeps  
        ab05=-0.5-abeps  
 Caja  ENDIF  
36    
37  C     Step forward temperature  C     Step forward temperature
38        DO j=jMin,jMax        DO j=jMin,jMax
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22