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 etaHnm1 :: surface r-anomaly, etaH, at previous time level |
43 |
C dEtaHdt :: time derivative of total column height [r_unit/s = w unit] |
C dEtaHdt :: time derivative of total column height [r_unit/s = w unit] |
44 |
C PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step |
C PmEpR :: keep the fresh water input (=-EmPmR) of the previous time step |
45 |
COMMON /EXACT_ETA_LOCAL/ dEtaHdt, PmEpR |
COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR |
46 |
|
_RL etaHnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
47 |
_RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
_RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
48 |
_RS PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
49 |
#endif |
#endif |