/[MITgcm]/MITgcm/model/src/correction_step.F
ViewVC logotype

Annotation of /MITgcm/model/src/correction_step.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.1 - (hide annotations) (download)
Mon Jun 1 20:36:13 1998 UTC (26 years ago) by adcroft
Branch: MAIN 
Swapped the time-stepping algorithm around (just a little bit).
We now officially use the predictor-corrector terminology.
We make the prediction step at the end of the dynamics() section
and store the result in the gUNm1, gVNm1, gTNm1 arrays.
The "tricky" part is that at the beginning of the dynamics section,
where the "correction" is made, theses arrays must be initialised
at the beginning of any run. A new routine init_predictor() does this.
This is "all" in preparation for implicit diffusion. Let's hope
it's enough...
1 adcroft 1.1 C $Header$
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     & -deltaTmom*hxFac*rRhonil *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     & -deltaTmom*hyFac*rRhonil *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     C Step forward temperature velocity
62     DO j=jMin,jMax
63     DO i=iMin,iMax
64     theta(i,j,k,bi,bj)=gTNm1(i,j,k,bi,bj)
65     gTNm1(i,j,k,bi,bj)=gT(i,j,k,bi,bj)
66     ENDDO
67     ENDDO
68    
69     RETURN
70     END
  ViewVC Help
Powered by ViewVC 1.1.22