pydda.cost_functions.calculate_mass_continuity

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

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