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