pydda.cost_functions.calculate_mass_continuity

pydda.cost_functions.calculate_mass_continuity(u, v, w, z, dx, dy, dz, coeff=1500.0, anel=1)[source]

Calculates the mass continuity cost function by taking the divergence of the wind field.

All arrays in the given lists must have the same dimensions and represent the same spatial coordinates.

Parameters
u: Float array

Float array with u component of wind field

v: Float array

Float array with v component of wind field

w: Float array

Float array with w component of wind field

dx: float

Grid spacing in x direction.

dy: float

Grid spacing in y direction.

dz: float

Grid spacing in z direction.

z: Float array (1D)

1D Float array with heights of grid

coeff: float

Constant controlling contribution of mass continuity to cost function

anel: int

= 1 use anelastic approximation, 0=don’t

Returns
J: float

value of mass continuity cost function