model/src/timestep.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/timestep.F,v 1.12 1998/06/10 16:05:39 adcroft 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 "EEPARAMS.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
      IF (tempStepping) THEN
       DO j=jMin,jMax
        DO i=iMin,iMax
         gTNm1(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))
        ENDDO
       ENDDO
      ENDIF

C     Step forward salt
      IF (saltStepping) THEN
       DO j=jMin,jMax
        DO i=iMin,iMax
         gSNm1(i,j,k,bi,bj)=salt(i,j,k,bi,bj)
     &    +deltaTtracer*(ab15*gS(i,j,k,bi,bj)+ab05*gSNm1(i,j,k,bi,bj))
        ENDDO
       ENDDO
      ENDIF

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