42 |
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] |
43 |
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 |
44 |
COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR |
COMMON /EXACT_ETA_LOCAL/ etaHnm1, dEtaHdt, PmEpR |
45 |
_RL etaHnm1(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
_RL etaHnm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
46 |
_RL dEtaHdt(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
_RL dEtaHdt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
47 |
_RL PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL PmEpR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
48 |
#endif |
#endif |
49 |
|
|
74 |
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 |
75 |
C rStarFacW :: same but for West face |
C rStarFacW :: same but for West face |
76 |
C rStarFacS :: same but for South face |
C rStarFacS :: same but for South face |
77 |
|
C pStarFacK :: rStarFacC**atm_kappa (for atmosphere in p* coords) |
78 |
C rStarFacNm1C, etc. :: prior values |
C rStarFacNm1C, etc. :: prior values |
79 |
C rStarExpC :: column expansion factor = h^n+1/h^n , Centered |
C rStarExpC :: column expansion factor = h^n+1/h^n , Centered |
80 |
C rStarExpW :: column expansion factor = h^n+1/h^n , Western face |
C rStarExpW :: column expansion factor = h^n+1/h^n , Western face |
83 |
C rStarDhWDt:: relative time derivative of h_West_face (u.point) |
C rStarDhWDt:: relative time derivative of h_West_face (u.point) |
84 |
C rStarDhSDt:: relative time derivative of h_South_face (v.point) |
C rStarDhSDt:: relative time derivative of h_South_face (v.point) |
85 |
COMMON /RSTAR_CHANGE/ |
COMMON /RSTAR_CHANGE/ |
86 |
& rStarFacC, rStarFacW, rStarFacS, |
& rStarFacC, rStarFacW, rStarFacS, pStarFacK, |
87 |
& rStarFacNm1C, rStarFacNm1W, rStarFacNm1S, |
& rStarFacNm1C, rStarFacNm1W, rStarFacNm1S, |
88 |
& rStarExpC, rStarExpW, rStarExpS, |
& rStarExpC, rStarExpW, rStarExpS, |
89 |
& rStarDhCDt,rStarDhWDt,rStarDhSDt |
& rStarDhCDt,rStarDhWDt,rStarDhSDt |
90 |
_RL rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
91 |
_RL rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
92 |
_RL rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
93 |
|
_RL pStarFacK (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
94 |
_RL rStarFacNm1C (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacNm1C (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
95 |
_RL rStarFacNm1W (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacNm1W (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
96 |
_RL rStarFacNm1S (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL rStarFacNm1S (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |