model/src/correction_step.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/correction_step.F,v 1.1 1998/06/01 20:36:13 adcroft Exp $

#include "CPP_EEOPTIONS.h"

C     /==========================================================\
C     | S/R CORRECTION_STEP                                      |
C     | o Corrects the horizontal flow fields with the surface   |
C     |   pressure gradient.                                     |
C     \==========================================================/
      SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,
     I                       K, pSurfX, pSurfY,
     I                       myThid )
      implicit none
! Common
#include "SIZE.h"
#include "DYNVARS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "EEPARAMS.h"
#include "CG2D.h"
C     == Routine Arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      INTEGER K
      _RL  pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL  pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER myThid
C     == Local variables ==
      INTEGER i,j
      _RL ab15,ab05,hxFac,hyFac,rRhoNil

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

C     On/off scaling paramters
      hxFac = pfFacMom
      hyFac = pfFacMom

      rRhoNil=1. / rhonil

C     Step forward zonal velocity
      DO j=jMin,jMax
       DO i=iMin,iMax
        uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
     &       -deltaTmom*hxFac*rRhonil *pSurfX(i,j)
     &                    )*_maskW(i,j,k,bi,bj)
        gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
       ENDDO
      ENDDO

C     Step forward meridional velocity
      DO j=jMin,jMax
       DO i=iMin,iMax
        vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
     &       -deltaTmom*hyFac*rRhonil *pSurfY(i,j)
     &                    )*_maskS(i,j,k,bi,bj)
        gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
       ENDDO
      ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/correction_step.F,v 1.1 1998/06/01 20:36:13 adcroft Exp $
2
3	#include "CPP_EEOPTIONS.h"
4
5	C /==========================================================\
6	C \| S/R CORRECTION_STEP \|
7	C \| o Corrects the horizontal flow fields with the surface \|
8	C \| pressure gradient. \|
9	C \==========================================================/
10	SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,
11	I K, pSurfX, pSurfY,
12	I myThid )
13	implicit none
14	! Common
15	#include "SIZE.h"
16	#include "DYNVARS.h"
17	#include "PARAMS.h"
18	#include "GRID.h"
19	#include "EEPARAMS.h"
20	#include "CG2D.h"
21	C == Routine Arguments ==
22	INTEGER bi,bj,iMin,iMax,jMin,jMax
23	INTEGER K
24	_RL pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
25	_RL pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26	INTEGER myThid
27	C == Local variables ==
28	INTEGER i,j
29	_RL ab15,ab05,hxFac,hyFac,rRhoNil
30
31	C Adams-Bashforth timestepping weights
32	ab15=1.5+abeps
33	ab05=-0.5-abeps
34
35	C On/off scaling paramters
36	hxFac = pfFacMom
37	hyFac = pfFacMom
38
39	rRhoNil=1. / rhonil
40
41	C Step forward zonal velocity
42	DO j=jMin,jMax
43	DO i=iMin,iMax
44	uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
45	& -deltaTmomhxFacrRhonil *pSurfX(i,j)
46	& )*_maskW(i,j,k,bi,bj)
47	gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
48	ENDDO
49	ENDDO
50
51	C Step forward meridional velocity
52	DO j=jMin,jMax
53	DO i=iMin,iMax
54	vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
55	& -deltaTmomhyFacrRhonil *pSurfY(i,j)
56	& )*_maskS(i,j,k,bi,bj)
57	gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
58	ENDDO
59	ENDDO
60
61	RETURN
62	END