utils.solver.flux_divergence

utils.solver.flux_divergence(face_fluxes, dz)

Compute finite-volume tendency from face fluxes and layer thicknesses.

Parameters

face_fluxes : (array - like, shape(n_layers + 1))

Fluxes at cell faces.

dz : (array - like, shape(n_layers))

Layer thicknesses.

Returns

tendency : (ndarray, shape(n_layers))

Per-layer tendency -(flux_out - flux_in) / dz.

Raises

: ValueError

If face_fluxes does not have length n_layers + 1.