riem_solver_c

NonhydrostaticVerticalSolverCGrid

Fortran subroutine Riem_Solver_C

Semi-implicit solver for pressure, vertical velocity, and dz (not a Riemann solver)

accounts for: Vertically-propagating sound wave and straining terms vertical non-hydrostatic pressure gradient force change in layer interface heights due to straining/compression by sound waves

__call__(dt2, cappa, ptop, hs, ws, ptc, q_con, delpc, gz, pef, w3)

Solves for the nonhydrostatic terms for vertical velocity (w) and non-hydrostatic pressure perturbation after C-grid winds advect and heights are updated.

Parameters:

Name	Type	Description	Default
dt2	in	acoustic timestep in seconds	required
cappa	in	???	required
ptop	in	pressure at top of atmosphere	required
hs	in	surface height in m	required
ws	in	vertical velocity of the lowest level	required
ptc	in	potential temperature	required
q_con	in	total condensate mixing ratio	required
delpc	in	vertical delta in pressure	required
gz	inout	geopotential height	required
pef	out	full hydrostatic pressure	required
w3	in	vertical velocity	required

finalize(pe2, pem, hs, dz, pef, gz, ptop)

Enforce vertical boundary conditions.

The top of atmosphere pressure is constant. At bottom of domain, the height should be equal to hs, the surface elevation.

Parameters:

Name	Type	Description	Default
pe2	in	nonhydrostatic perturbation pressure defined on interfaces	required
pem	in	total hydrostatic pressure defined on interface (including condensate)	required
hs	in	surface elevation	required
dz	in		required
pef	out		required
gz	out		required

precompute(delpc, cappa, w3, w, gz, dm, q_con, pem, dz, gm, pm, ptop)

Parameters:

Name	Type	Description	Default
delpc	in		required
cappa	in		required
w3	in		required
w	out		required
gz	in		required
dm	out	delta mass, mass of gridcell per unit area	required
q_con	in		required
pem	out	total hydrostatic pressure defined on interface (including condensate)	required
dz	out		required
gm	out	gamma parameter, Cp/Cv, used to compute pressure gradient force using potential temperature and ideal gas law	required
pm	out	hydrostatic cell mean pressure, derivation in documentation (Chapter 4? 7?) TODO: identify chapter reference, will be sent by Lucas	required