Skip to content

dyn_core

AcousticDynamics

Fortran name is dyn_core Performs the Lagrangian acoustic dynamics described by Lin 2004

__init__(comm, stencil_factory, quantity_factory, grid_data, damping_coefficients, grid_type, nested, stretched_grid, config, phis, wsd, state, checkpointer=None)

Parameters:

Name Type Description Default
comm Communicator

object for tile or cubed-sphere inter-process communication

required
stencil_factory StencilFactory

creates stencils

required
quantity_factory QuantityFactory

creates quantities

required
grid_data GridData

metric terms defining the grid

required
damping_coefficients DampingCoefficients

damping configuration

required
grid_type

grid geometry used

required
nested

if the grid contains a nested, high-res region

required
stretched_grid

if the grid is stretched so tile faces cover different areas

required
config AcousticDynamicsConfig

configuration settings

required
pfull

atmospheric Eulerian grid reference pressure (Pa)

required
phis FloatFieldIJ

surface geopotential height

required
checkpointer Optional[Checkpointer]

if given, used to perform operations on model data at specific points in model execution, such as testing against reference data

None

gz_from_surface_height_and_thicknesses(zs, delz, gz)

Parameters:

Name Type Description Default
zs in

surface height

required
delz in

layer thickness

required
gz out

geopotential height

required

interface_pressure_from_toa_pressure_and_thickness(delp, pem, ptop)

Parameters:

Name Type Description Default
delp in

pressure thickness of atmospheric layer

required
pem out

interface pressure

required
ptop in

pressure at top of atmosphere

required

p_grad_c_stencil(rdxc, rdyc, uc, vc, delpc, pkc, gz, dt2)

Update C-grid winds from the backwards-in-time pressure gradient force

When this is run the C-grid winds have almost been completely updated by computing the momentum equation terms, but the pressure gradient force term has not yet been applied. This stencil completes the equation and Arakawa C-grid winds have been advected half a timestep upon completing this stencil.

Parameters:

Name Type Description Default
rdxc in
required
rdyc in
required
uc inout

x-velocity on the C-grid, has been updated due to advection but not yet due to pressure gradient force

required
vc inout

y-velocity on the C-grid, has been updated due to advection but not yet due to pressure gradient force

required
delpc in

vertical delta in pressure

required
pkc in

pressure if non-hydrostatic, (edge pressure)**(moist kappa) if hydrostatic

required
gz in

height of the model grid cells (m)

required
dt2 in

half a model timestep (for C-grid update) in seconds

required

zero_data(mfxd, mfyd, cxd, cyd, heat_source, diss_estd, first_timestep)

Parameters:

Name Type Description Default
mfxd out

mass flux in x direction

required
mfyd out

mass flux in y direction

required
cxd out

courant number in x direction

required
cyd out

courant number in y direction

required
heat_source out

heat accumulated from diffusion of kinetic energy gets applied at the end of all acoustic steps

required
diss_estd out

how much energy is dissipated, is mainly captured to send to the stochastic physics (in contrast to heat_source)

required
first_timestep in

is this the first acoustic timestep?

required