model/src/init_predictor.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/init_predictor.F,v 1.3 1998/06/16 15:20:08 adcroft Exp $

#include "CPP_EEOPTIONS.h"

C     /==========================================================\
C     | S/R INIT_PREDICTOR                                       |
C     | o Copies Gu,Gv,Gt and Gs in U,V,theta and salt           |
C     |   This is to initialise the predictor-corrector time-    |
C     |   stepping scheme                                        |
C     \==========================================================/
      SUBROUTINE INIT_PREDICTOR( 
     I                       myThid )
      implicit none
! Common
#include "SIZE.h"
#include "DYNVARS.h"
C     == Routine Arguments ==
      INTEGER myThid
C     == Local variables ==
      INTEGER i,j,k,bi,bj

      DO bj=1,nSy
       DO bi=1,nSx
        DO k=1,Nz
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
           gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
           gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
           gTNm1(i,j,k,bi,bj)=theta(i,j,k,bi,bj)
           gSNm1(i,j,k,bi,bj)=salt(i,j,k,bi,bj)
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/init_predictor.F,v 1.3 1998/06/16 15:20:08 adcroft Exp $
2
3	#include "CPP_EEOPTIONS.h"
4
5	C /==========================================================\
6	C \| S/R INIT_PREDICTOR \|
7	C \| o Copies Gu,Gv,Gt and Gs in U,V,theta and salt \|
8	C \| This is to initialise the predictor-corrector time- \|
9	C \| stepping scheme \|
10	C \==========================================================/
11	SUBROUTINE INIT_PREDICTOR(
12	I myThid )
13	implicit none
14	! Common
15	#include "SIZE.h"
16	#include "DYNVARS.h"
17	C == Routine Arguments ==
18	INTEGER myThid
19	C == Local variables ==
20	INTEGER i,j,k,bi,bj
21
22	DO bj=1,nSy
23	DO bi=1,nSx
24	DO k=1,Nz
25	DO j=1-Oly,sNy+Oly
26	DO i=1-Olx,sNx+Olx
27	gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj)
28	gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj)
29	gTNm1(i,j,k,bi,bj)=theta(i,j,k,bi,bj)
30	gSNm1(i,j,k,bi,bj)=salt(i,j,k,bi,bj)
31	ENDDO
32	ENDDO
33	ENDDO
34	ENDDO
35	ENDDO
36
37	RETURN
38	END