updatedzc

UpdateGeopotentialHeightOnCGrid

__call__(zs, ut, vt, gz, ws, dt)

Parameters:

Name	Type	Description	Default
dp_ref		layer thickness in Pa	required
zs	FloatFieldIJ	surface height in m	required
ut	FloatField	horizontal wind (TODO: covariant or contravariant?)	required
vt	FloatField	horizontal wind (TODO: covariant or contravariant?)	required
gz	FloatField	geopotential height on model interfaces	required
ws	FloatFieldIJ	surface vertical wind implied by horizontal motion over topography	required
dt	Float	timestep over which to evolve the geopotential height, in seconds	required

update_dz_c(dp_ref, zs, area, ut, vt, gz, gz_x, gz_y, ws, *, dt)

Step dz forward on c-grid Ensures gz is monotonically increasing in z at the end Args: dp_ref: zs: area: ut: vt: gz: gz_x: gz with corners copied to perform derivatives in x-direction gz_y: gz with corners copied to perform derivatives in y-direction ws: lagrangian (parcel-following) surface vertical wind implied by lowest-level gz change note that a parcel moving horizontally across terrain will be moving in the vertical (eqn 5.5 in documentation) dt:

xy_flux(gz_x, gz_y, xfx, yfx)

Compute first-order upwind fluxes of gz in x and y directions.

Parameters:

Name	Type	Description	Default
gz_x		gz with corners copied to perform derivatives in x-direction	required
gz_y		gz with corners copied to perform derivatives in y-direction	required
xfx	out	contravariant c-grid u-wind interpolated to layer interfaces, including metric terms to make it a "volume flux"	required
yfx	out	contravariant c-grid v-wind interpolated to layer interfaces	required

Returns:

Name	Type	Description
fx		first-order upwind x-flux of gz
fy		first-order upwind y-flux of gz