53 |
LOGICAL momDissip_In_AB |
LOGICAL momDissip_In_AB |
54 |
LOGICAL momStartAB |
LOGICAL momStartAB |
55 |
INTEGER i,j |
INTEGER i,j |
|
_RL ab15,ab05 |
|
56 |
_RL phxFac,phyFac, psFac |
_RL phxFac,phyFac, psFac |
57 |
_RL gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
58 |
_RL gVtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL gVtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
70 |
#endif |
#endif |
71 |
CEOP |
CEOP |
72 |
|
|
73 |
C Adams-Bashforth timestepping weights |
C Start AB with low-order timestepping weights |
74 |
momStartAB = nIter0.EQ.0 |
momStartAB = nIter0.EQ.0 |
|
IF (myIter .EQ. 0) THEN |
|
|
ab15=1.0 |
|
|
ab05=0.0 |
|
|
ELSE |
|
|
ab15=1.5+abeps |
|
|
ab05=-0.5-abeps |
|
|
ENDIF |
|
75 |
|
|
76 |
C-- explicit part of the surface potential gradient is added in this S/R |
C-- explicit part of the surface potential gradient is added in this S/R |
77 |
psFac = pfFacMom*(1. _d 0 - implicSurfPress) |
psFac = pfFacMom*(1. _d 0 - implicSurfPress) |