model/src/correction_step.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/correction_step.F,v 1.4 1998/06/16 15:20:08 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 "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.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 hxFac,hyFac,rRhoNil

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

C     Rotate theta/gT/gTnm1
      DO j=jMin,jMax
       DO i=iMin,iMax
        theta(i,j,k,bi,bj)=gTNm1(i,j,k,bi,bj)
        gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)
       ENDDO
      ENDDO

C     Rotate salt/gS/gSnm1
      DO j=jMin,jMax
       DO i=iMin,iMax
        salt(i,j,k,bi,bj)=gSNm1(i,j,k,bi,bj)
        gSNm1(i,j,k,bi,bj)=gS(i,j,k,bi,bj)
       ENDDO
      ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/correction_step.F,v 1.4 1998/06/16 15:20:08 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 "EEPARAMS.h"
18	#include "PARAMS.h"
19	#include "GRID.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 hxFac,hyFac,rRhoNil
30
31	C On/off scaling paramters
32	hxFac = pfFacMom
33	hyFac = pfFacMom
34
35	rRhoNil=1. / rhonil
36
37	C Step forward zonal velocity
38	DO j=jMin,jMax
39	DO i=iMin,iMax
40	uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
41	& -deltaTmomhxFacrRhonil *pSurfX(i,j)
42	& )*_maskW(i,j,k,bi,bj)
43	gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
44	ENDDO
45	ENDDO
46
47	C Step forward meridional velocity
48	DO j=jMin,jMax
49	DO i=iMin,iMax
50	vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
51	& -deltaTmomhyFacrRhonil *pSurfY(i,j)
52	& )*_maskS(i,j,k,bi,bj)
53	gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
54	ENDDO
55	ENDDO
56
57	C Rotate theta/gT/gTnm1
58	DO j=jMin,jMax
59	DO i=iMin,iMax
60	theta(i,j,k,bi,bj)=gTNm1(i,j,k,bi,bj)
61	gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)
62	ENDDO
63	ENDDO
64
65	C Rotate salt/gS/gSnm1
66	DO j=jMin,jMax
67	DO i=iMin,iMax
68	salt(i,j,k,bi,bj)=gSNm1(i,j,k,bi,bj)
69	gSNm1(i,j,k,bi,bj)=gS(i,j,k,bi,bj)
70	ENDDO
71	ENDDO
72
73	RETURN
74	END