39 |
INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
40 |
|
|
41 |
#ifdef EXACT_CONSERV |
#ifdef EXACT_CONSERV |
42 |
C hDivFlow :: Div. Barotropic Flow at current time [transport unit m3/s] |
C dEtaHdt :: time derivative of total column height [r_unit/s = w unit] |
43 |
COMMON /EXACT_ETA_LOCAL/ hDivFlow |
C PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step |
44 |
_RL hDivFlow(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
COMMON /EXACT_ETA_LOCAL/ dEtaHdt, PmEpR |
45 |
|
_RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
46 |
|
_RS PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
47 |
#endif |
#endif |
48 |
|
|
49 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
52 |
C center (Tracer point) |
C center (Tracer point) |
53 |
C hFac_surfW :: idem, West interface (U point) |
C hFac_surfW :: idem, West interface (U point) |
54 |
C hFac_surfS :: idem, South interface (V point) |
C hFac_surfS :: idem, South interface (V point) |
|
C PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step |
|
55 |
COMMON /SURF_CHANGE/ |
COMMON /SURF_CHANGE/ |
56 |
& hFac_surfC, hFac_surfW, hFac_surfS, |
& hFac_surfC, hFac_surfW, hFac_surfS |
|
& PmEpR |
|
57 |
_RS hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS hFac_surfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
58 |
_RS hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS hFac_surfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
59 |
_RS hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS hFac_surfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
_RS PmEpR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
60 |
|
|
61 |
C-- COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate |
C-- COMMON /RSTAR_CHANGE/ transient variables used with r* Coordinate |
62 |
C rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H |
C rStarFacC :: = dr/dr* = ratio of r-thickness / r*-thickness = h^n / H |