/[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.3 by edhill, Fri Oct 31 20:35:32 2003 UTC revision 1.6 by jmc, Fri Apr 15 14:11:01 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
 #include "PACKAGES_CONFIG.h"  
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CBOP  CBOP
7  C     !ROUTINE: ADAMS_BASHFORTH2  C     !ROUTINE: ADAMS_BASHFORTH2
8  C     !INTERFACE:  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     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
14  C     *==========================================================*  C     *==========================================================*
15  C     | S/R ADAMS_BASHFORTH2                                        C     | S/R ADAMS_BASHFORTH2                                      
16  C     | o Extrapolate tendancies forward in time using              C     | o Extrapolate tendencies forward in time using            
17  C     |   quasi-second order Adams-Bashforth method.                C     |   quasi-second order Adams-Bashforth method.              
18  C     *==========================================================*  C     *==========================================================*
19  C     \ev  C     \ev
# Line 25  C     == Global variables === Line 24  C     == Global variables ===
24  #include "SIZE.h"  #include "SIZE.h"
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
 #include "GAD.h"  
 #include "GRID.h"  
 #include "SURFACE.h"  
27    
28  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
29  C     == Routine Arguments ==  C     == Routine Arguments ==
30  C     bi,bj,K :: Tile and level indices  C     bi,bj,k :: Tile and level indices
31  C     gTracer :: Tendency at current time  ( generally units of quantity/sec )  C     gTracer :: Tendency at current time  ( generally units of quantity/sec )
32  C     gTrNm1  :: Tendency at previous time ( generally units of quantity/sec )  C     gTrNm1  :: Tendency at previous time ( generally units of quantity/sec )
33  C     myIter  :: Current time step number  C     myIter  :: Current time step number
34  C     myThid  :: Thread number of this thread  C     myThid  :: Thread number of this thread
35        INTEGER bi,bj,K        INTEGER bi,bj,k
36        _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)
37        _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)
38        INTEGER myIter, myThid        INTEGER myIter, myThid
# Line 51  C     ab15, ab05 :: Adams bashforth extr Line 47  C     ab15, ab05 :: Adams bashforth extr
47  CEOP  CEOP
48    
49  C     Adams-Bashforth timestepping weights  C     Adams-Bashforth timestepping weights
50        IF (myIter .EQ. 0) THEN        IF ( myIter.EQ.0 ) THEN
51         ab15=1.0         ab15=1.0
52         ab05=0.0         ab05=0.0
53        ELSE        ELSE

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6

  ViewVC Help
Powered by ViewVC 1.1.22