/[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.8 by jmc, Mon Nov 7 20:20:31 2011 UTC revision 1.9 by jmc, Tue Feb 19 13:42:19 2013 UTC
# Line 7  CBOP Line 7  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, kSize,
11       U                     gTracer, gTrNm1,       U                     gTracer, gTrNm1,
12       O                     AB_gTr,       O                     AB_gTr,
13       I                     startAB, myIter, myThid )       I                     startAB, myIter, myThid )
# Line 29  C     == Global variables === Line 29  C     == Global variables ===
29  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
30  C     == Routine Arguments ==  C     == Routine Arguments ==
31  C     bi,bj,k :: Tile and level indices  C     bi,bj,k :: Tile and level indices
32    C     kSize   :: 3rd dimension of tracer and tendency arrays
33  C     gTracer ::  in: Tendency at current time  ( generally units of quantity/sec )  C     gTracer ::  in: Tendency at current time  ( generally units of quantity/sec )
34  C             :: out: Extrapolated Tendency at current time ( same units )  C             :: out: Extrapolated Tendency at current time ( same units )
35  C     gTrNm1  ::  in: Tendency at previous time             ( same units )  C     gTrNm1  ::  in: Tendency at previous time             ( same units )
# Line 37  C     AB_gTr  :: Adams-Bashforth tendenc Line 38  C     AB_gTr  :: Adams-Bashforth tendenc
38  C     startAB :: number of previous time level available to start/restart AB  C     startAB :: number of previous time level available to start/restart AB
39  C     myIter  :: Current time step number  C     myIter  :: Current time step number
40  C     myThid  :: my Thread Id. number  C     myThid  :: my Thread Id. number
41        INTEGER bi,bj,k        INTEGER bi, bj, k, kSize
42        _RL  gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)
43        _RL  gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)
44        _RL  AB_gTr (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  AB_gTr (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45        INTEGER startAB        INTEGER startAB
46        INTEGER myIter, myThid        INTEGER myIter, myThid
# Line 62  C     Adams-Bashforth extrapolation fact Line 63  C     Adams-Bashforth extrapolation fact
63  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
64    
65  C-    Compute effective G-term with Adams-Bashforth weights:  C-    Compute effective G-term with Adams-Bashforth weights:
66        DO j=1-Oly,sNy+Oly        DO j=1-OLy,sNy+OLy
67         DO i=1-Olx,sNx+Olx         DO i=1-OLx,sNx+OLx
68          AB_gTr(i,j) = abFac          AB_gTr(i,j) = abFac
69       &              *( gTracer(i,j,k,bi,bj)- gTrNm1(i,j,k,bi,bj) )       &              *( gTracer(i,j,k,bi,bj)- gTrNm1(i,j,k,bi,bj) )
70          gTrNm1(i,j,k,bi,bj) = gTracer(i,j,k,bi,bj)          gTrNm1(i,j,k,bi,bj) = gTracer(i,j,k,bi,bj)
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22