296 |
C \ev |
C \ev |
297 |
CEOP |
CEOP |
298 |
|
|
|
#include "PACKAGES_CONFIG.h" |
|
|
|
|
299 |
C Macros that override/modify standard definitions |
C Macros that override/modify standard definitions |
300 |
#include "GRID_MACROS.h" |
#include "GRID_MACROS.h" |
301 |
|
|
320 |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
C Note: The code needs terms like MIN(hFac,hFac(I+1)) |
321 |
C On some platforms it may be better to precompute |
C On some platforms it may be better to precompute |
322 |
C hFacW, hFacE, ... here than do MIN on the fly. |
C hFacW, hFacE, ... here than do MIN on the fly. |
323 |
C gravitySign - indicates whether gravity points in the opposite |
C gravitySign - indicates the direction of gravity relative to R direction |
324 |
C direction of R or not. |
C (= -1 for R=Z (Z increases upward, -gravity direction ) |
325 |
C ( = +1 for R=Z (gravity points downward in Z) |
C (= +1 for R=P (P increases downward, +gravity direction) |
|
C ( = -1 for R=P (gravity points upward in P) |
|
326 |
C rkFac - Vertical coordinate to vertical index orientation. |
C rkFac - Vertical coordinate to vertical index orientation. |
327 |
C ( -1 same orientation, 1 opposite orientation ) |
C ( -1 same orientation, 1 opposite orientation ) |
|
C ( vertical coord == m -> rkFac = 1 ) |
|
|
C ( vertical coord == Pa -> rkFac = -1 ) |
|
328 |
C maskH - cell Center full-column mask (= 2D mask) |
C maskH - cell Center full-column mask (= 2D mask) |
329 |
C maskC - cell Center land mask |
C maskC - cell Center land mask |
330 |
C maskW - West face land mask |
C maskW - West face land mask |
368 |
INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
INTEGER klowC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
369 |
|
|
370 |
COMMON /GRID_R/ |
COMMON /GRID_R/ |
371 |
|
& cosfacU,cosfacV,sqcosfacU,sqcosfacV, |
372 |
|
& gravitySign,rkFac, recip_rkFac, |
373 |
& dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU, |
& dxC,dxF,dxG,dxV,dyC,dyF,dyG,dyU, |
374 |
& R_low,Ro_surf,HFacC,HFacW,HFacS, |
& R_low,Ro_surf,HFacC,HFacW,HFacS, |
375 |
& recip_dxC,recip_dxF,recip_dxG,recip_dxV, |
& recip_dxC,recip_dxF,recip_dxG,recip_dxV, |
381 |
& maskH, maskC,maskW,maskS, |
& maskH, maskC,maskW,maskS, |
382 |
& recip_rA,recip_rAw,recip_rAs,recip_rAz, |
& recip_rA,recip_rAw,recip_rAs,recip_rAz, |
383 |
& tanPhiAtU, tanPhiAtV, |
& tanPhiAtU, tanPhiAtV, |
|
& cosfacU,cosfacV,sqcosfacU,sqcosfacV, |
|
384 |
& drC,drF,recip_drC,recip_drF,rC,rF, |
& drC,drF,recip_drC,recip_drF,rC,rF, |
385 |
& gravitySign, |
& xC0, yC0 |
386 |
& rkFac, recip_rkFac, xC0, yC0 |
_RL cosfacU(1-Oly:sNy+Oly,nSx,nSy) |
387 |
|
_RL cosfacV(1-Oly:sNy+Oly,nSx,nSy) |
388 |
|
_RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy) |
389 |
|
_RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy) |
390 |
|
_RL gravitySign |
391 |
|
_RL rkFac |
392 |
|
_RL recip_rkFac |
393 |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
394 |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
395 |
_RS dxG (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS dxG (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
433 |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
_RS maskS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:Nr,nSx,nSy) |
434 |
_RS tanPhiAtU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS tanPhiAtU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
435 |
_RS tanPhiAtV (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS tanPhiAtV (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
_RL cosfacU(1-Oly:sNy+Oly,nSx,nSy) |
|
|
_RL cosfacV(1-Oly:sNy+Oly,nSx,nSy) |
|
|
_RL sqcosfacU(1-Oly:sNy+Oly,nSx,nSy) |
|
|
_RL sqcosfacV(1-Oly:sNy+Oly,nSx,nSy) |
|
436 |
_RS drC (1:Nr) |
_RS drC (1:Nr) |
437 |
_RS drF (1:Nr) |
_RS drF (1:Nr) |
438 |
_RS recip_drC (1:Nr) |
_RS recip_drC (1:Nr) |
440 |
_RS saFac (1:Nr) |
_RS saFac (1:Nr) |
441 |
_RS rC (1:Nr) |
_RS rC (1:Nr) |
442 |
_RS rF (1:Nr+1) |
_RS rF (1:Nr+1) |
|
_RL gravitySign |
|
|
_RS rkFac |
|
|
_RS recip_rkFac |
|
443 |
_RS xC0 |
_RS xC0 |
444 |
_RS yC0 |
_RS yC0 |
445 |
|
|