Navigation: Home ›  Produkte ›  WaSiM-Richards

WaSiM (Richards)

Important: version 10.02.03 and higher provide netCDF support. That requires a few additional libraries (DLLs) to be installed on your system. Please download and install either the lib-package for 32-bit or the lib-package for 64-bit before using WaSiM on Windows and/or download and install the latest netCDF package from the internet. On Linux, you need to install a netCDF package as well (see development details for more information)

Actual release version (may still be unstable -> it is recommended to always use the latest patch, if one is provided!)

Version 10.07.02(2023-11-12) More »
ice wedge degradation and many patches
  • new functionality for ice wedge degradation: ExcessIceFactor is initialised as additional water equivalent and hgeo and theta are initialized according to the excess ice content. When running unsatzonmodel, actual theta_res and theta_sat are set for each layer in an ice wedge. All functions of the richards-solution (in unsatzon and soiltable) were extended by handling different thetasat and thetares values for ice wedges (e.g. if the normal theta_sat is 0.4 and thetares is 0.1 but with excess ice it is 0.65 and 0.05 resp., a return theta of 0.25 from the h_calc function (which is exactly 50% of useable water in the soil with the original theta values) is changed into 0.35 (which is also 50% of usable water in the soil with the excess ice - here, including the ice of course). After each model run, the changed DHM is written back to the input grid directory of the original dhm with the name extension _out. This could be used for coupling several runs (rename it first and remove the _out-extension!) IceWedge degradation can be used by
  • adding a new control file section like this:

    1 # do_ArcVeg; 0=no, 1=yes
    $outpath//$ExcessIceStack # stack, factor that expresses the excess ice in that soil layer
    $Writestack # Writecode for this stack

  • adding two standard grids, containing the ice wedges upper and lower depth for each cell (-9999 if no ice wedge is present)

    $set $IceWedgeUpperDepth = $grid//.iwdgup
    $set $IceWedgeLowerDepth = $grid//.iwdglo
    xx # number of standard grids
    $inpath_grid//$IceWedgeUpperDepth IceWedgeUpperDepth fillcode = 0
    $inpath_grid//$IceWedgeLowerDepth IceWedgeLowerDepth fillcode = 0

  • fixed routing with -mode:routing for lakes with abstraction rule "HydraulicConnection". No routing took place at all because there usually was no connection (both water levels were 0 but channel bottom level was above). To make this work, the routing description must contain two other values: A0 and W0 (A0 = initial area of the lake in m2, W0 is initial water table height in m a.s.l.). Example:

    ...and SP 1 ( file = $outpath//spv_01_sub09.//$code//$year , V0 = 4.69261e+10, A0 = 473.25E6, W0 = 395.107, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 )

    0 and W0 must follow V0 in this order. A0 and W0 must be provided for all lakes that are connected hydraulically.
    V0, A0 and W0 can/must be taken from the last entry in the lake-statistics file of the previous run manually, since -mode:routing doesn't read these Values from the richards-storage file.
    • and several bug fixes (see details of model development)
-> older versions More »

Patched versions (stable versions only are listed here)

Version 10.05.05 (2020-11-15) Mehr »
  • new (only Richards-version): implementing changes into the main branch of WaSiM that were made by Matvey Debolskiy during his PhD thesis to better model the frozen depth in a permafrost soil (only matters when heat transfer is used)
  • bug fix: interpolation method 9 did not handle memory for the filename read from the gridlist correctly (was allocated too late)
  • and several bug fixes (see details of model development)
-> ältere Versionen Mehr »


Beispiel-Datensatz (Inputdaten für WaSiM) Mehr »
mit Grids und Zeitreihen für ein Beobachtungsjahr
back top home



Neuer Release verfügbar:
WaSiM 10.07.02 -> Mehr »