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

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

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

revision 1.1 by adcroft, Wed Sep 19 02:43:27 2001 UTC revision 1.2 by cnh, Wed Sep 26 18:09:13 2001 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    CBOP
7    C     !ROUTINE: ADAMS_BASHFORTH2
8    C     !INTERFACE:
9        SUBROUTINE ADAMS_BASHFORTH2(        SUBROUTINE ADAMS_BASHFORTH2(
10       I                     bi, bj, K,       I                     bi, bj, K,
11       U                     gTracer, gTrNm1,       U                     gTracer, gTrNm1,
12       I                     myIter, myThid )       I                     myIter, myThid )
13  C     /==========================================================\  C     !DESCRIPTION: \bv
14  C     | S/R ADAMS_BASHFORTH2                                     |  C     *==========================================================*
15  C     | o Extrapolate tendancies forward in time using           |  C     | S/R ADAMS_BASHFORTH2                                      
16  C     |   quasi-second order Adams-Bashforth method.             |  C     | o Extrapolate tendancies forward in time using            
17  C     \==========================================================/  C     |   quasi-second order Adams-Bashforth method.              
18        IMPLICIT NONE  C     *==========================================================*
19    C     \ev
20    
21    C     !USES:
22          IMPLICIT NONE
23  C     == Global variables ===  C     == Global variables ===
24  #include "SIZE.h"  #include "SIZE.h"
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 22  C     == Global variables === Line 28  C     == Global variables ===
28  #include "GRID.h"  #include "GRID.h"
29  #include "SURFACE.h"  #include "SURFACE.h"
30    
31    C     !INPUT/OUTPUT PARAMETERS:
32  C     == Routine Arguments ==  C     == Routine Arguments ==
33    C     bi,bj,K :: Tile and level indices
34    C     gTracer :: Tendency at current time  ( generally units of quantity/sec )
35    C     gTrNm1  :: Tendency at previous time ( generally units of quantity/sec )
36    C     myIter  :: Current time step number
37    C     myThid  :: Thread number of this thread
38        INTEGER bi,bj,K        INTEGER bi,bj,K
39        _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)
40        _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)
41        INTEGER myIter, myThid        INTEGER myIter, myThid
42    
43    C     !LOCAL VARIABLES:
44  C     == Local variables ==  C     == Local variables ==
45    C     i,j        :: Loop counters
46    C     ab15, ab05 :: Adams bashforth extrapolation weights.
47        INTEGER i,j        INTEGER i,j
48        _RL ab15,ab05        _RL ab15,ab05
49        _RL gTrtmp        _RL gTrtmp
50    CEOP
51    
52  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
53        IF (myIter .EQ. 0) THEN        IF (myIter .EQ. 0) THEN
# Line 41  C     Adams-Bashforth timestepping weigh Line 58  C     Adams-Bashforth timestepping weigh
58         ab05=-(0.5+abEps)         ab05=-(0.5+abEps)
59        ENDIF        ENDIF
60    
61  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----
62    
63  C-    Compute effective G-term with Adams-Bashforth weights:  C-    Compute effective G-term with Adams-Bashforth weights:
64        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22