model/src/timestep.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/timestep.F,v 1.9 1998/06/08 21:43:02 cnh Exp $

#include "CPP_OPTIONS.h"

C     /==========================================================\
C     | S/R TIMESTEP                                             |
C     | o Step model fields forward in time                      |
C     \==========================================================/
      SUBROUTINE TIMESTEP( bi, bj, iMin, iMax, jMin, jMax,
     I                     K,
     I                     myThid )
      implicit none
! Common
#include "SIZE.h"
#include "DYNVARS.h"
#include "PARAMS.h"
#include "GRID.h"
C     == Routine Arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      INTEGER K
      INTEGER myThid
C     == Local variables ==
      INTEGER i,j
      _RL ab15,ab05

C     Adams-Bashforth timestepping weights
      ab15=1.5+abeps
      ab05=-0.5-abeps

C     Step forward zonal velocity (store in Gu)
       DO j=jMin,jMax
        DO i=iMin,iMax
         gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
     &    +deltaTmom*(ab15*gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj)
#ifdef ALLOW_CD
     &    +guCD(i,j,k,bi,bj)
#endif
     &        )*_maskW(i,j,k,bi,bj)
        ENDDO
       ENDDO
C     Step forward meridional velocity (store in Gv)
       DO j=jMin,jMax
        DO i=iMin,iMax
         gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
     &    +deltaTmom*(ab15*gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj)
#ifdef ALLOW_CD
     &    +gvCD(i,j,k,bi,bj)
#endif
     &        )*_maskS(i,j,k,bi,bj)
        ENDDO
       ENDDO
C     Step forward temperature
       DO j=jMin,jMax
        DO i=iMin,iMax
         theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj)
     &    +deltaTtracer*(ab15*gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj))
         gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)
        ENDDO
       ENDDO

      RETURN
      END
1	adcroft	1.10	C $Header: /u/gcmpack/models/MITgcmUV/model/src/timestep.F,v 1.9 1998/06/08 21:43:02 cnh Exp $
2	cnh	1.1
3	adcroft	1.10	#include "CPP_OPTIONS.h"
4	cnh	1.1
5			C /==========================================================\
6			C \| S/R TIMESTEP \|
7			C \| o Step model fields forward in time \|
8			C \==========================================================/
9	adcroft	1.4	SUBROUTINE TIMESTEP( bi, bj, iMin, iMax, jMin, jMax,
10	adcroft	1.7	I K,
11	adcroft	1.4	I myThid )
12	cnh	1.1	implicit none
13			! Common
14			#include "SIZE.h"
15			#include "DYNVARS.h"
16			#include "PARAMS.h"
17			#include "GRID.h"
18			C == Routine Arguments ==
19			INTEGER bi,bj,iMin,iMax,jMin,jMax
20	adcroft	1.4	INTEGER K
21	cnh	1.1	INTEGER myThid
22			C == Local variables ==
23	adcroft	1.4	INTEGER i,j
24	adcroft	1.7	_RL ab15,ab05
25	cnh	1.1
26			C Adams-Bashforth timestepping weights
27			ab15=1.5+abeps
28			ab05=-0.5-abeps
29
30	adcroft	1.7	C Step forward zonal velocity (store in Gu)
31	cnh	1.1	DO j=jMin,jMax
32			DO i=iMin,iMax
33	adcroft	1.7	gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
34	cnh	1.1	& +deltaTmom(ab15gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj)
35	cnh	1.9	#ifdef ALLOW_CD
36			& +guCD(i,j,k,bi,bj)
37			#endif
38	adcroft	1.7	& )*_maskW(i,j,k,bi,bj)
39	cnh	1.1	ENDDO
40			ENDDO
41	adcroft	1.7	C Step forward meridional velocity (store in Gv)
42	cnh	1.1	DO j=jMin,jMax
43			DO i=iMin,iMax
44	adcroft	1.7	gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
45	cnh	1.1	& +deltaTmom(ab15gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj)
46	cnh	1.9	#ifdef ALLOW_CD
47			& +gvCD(i,j,k,bi,bj)
48			#endif
49	adcroft	1.7	& )*_maskS(i,j,k,bi,bj)
50	cnh	1.1	ENDDO
51			ENDDO
52			C Step forward temperature
53			DO j=jMin,jMax
54			DO i=iMin,iMax
55	adcroft	1.8	theta(i,j,k,bi,bj)=theta(i,j,k,bi,bj)
56	cnh	1.1	& +deltaTtracer(ab15gT(i,j,k,bi,bj)+ab05*gTNm1(i,j,k,bi,bj))
57	adcroft	1.8	gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)
58	cnh	1.1	ENDDO
59			ENDDO
60
61			RETURN
62			END