pydeltasnow.core.deltasnow_snowpack_evolution

pydeltasnow.core.deltasnow_snowpack_evolution(Hobs, rho_max, rho_null, c_ov, k_ov, k, tau, eta_null, resolution)[source]

This is the main loop in the delta snow model. Should be called on HS chunks of consecutive nonzeros.

Parameters

  • Hobs (1D numpy.ndarray of floats) –

    Measured snow height. Needs to be in [m]. Must comply to the following constraints:

    • no nans

    • continuous entries (no missing dates)

  • rho_max (float, optional) – Maximum density of an individual snow layer produced by the deltasnow model in [kg/m3], rho_max needs to be positive. The default is 401.2588.

  • rho_null (float, optional) – Fresh snow density for a newly created layer [kg/m3], rho_null needs to be positive. The default is 81.19417.

  • c_ov (float, optional) – Overburden factor due to fresh snow [-], c_ov needs to be positive. The default is 0.0005104722.

  • k_ov (float, optional) – Defines the impact of the individual layer density on the compaction due to overburden [-], k_ov need to be in the range [0,1]. The default is 0.37856737.

  • k (float, optional) – Exponent of the exponential-law compaction [m3/kg], k needs to be positive. The default is 0.02993175.

  • tau (float, optional) – Uncertainty bound [m], tau needs to be positive. The default is 0.02362476.

  • eta_null (float, optional) – Effective compactive viscosity of snow for “zero-density” [Pa s]. The default is 8523356.

  • resolution (float) – Timedelta in hours between snow observations.

Raises

RuntimeError – If the snowpack evolution has somehow gone wrong.

Returns

SWE – Calculated SWE in [mm]. Same shape as Hobs.

Return type

1D numpy.ndarray of floats