Navigation: Home ›  Das Modell ›  Modellentwicklung ›  Meilensteine

Meilensteine der Modellentwicklung

1994 - 1996 Modellversion I mit erweitertem Sättigungsflächenansatz nach TOPMODEL.
Weiterhin: Bereitstellung von diversen WaSiM-Helptools für das Pre- und Postprocessing (z.B. Tanalys)
1997 - 1998 Modellversion II mit Bodenwasserdynamik in geschichteten Böden nach RICHARDS.
Weiterhin: Implementierung eines 2.5 dimensionalen Grundwassermodells mit Exfiltration und Reinfiltration; Berücksichtigung von Bewässerungsmanagement sowie Tracer- und Salztransport
2000 Erweiterung des Schneemodells
2001 Prototyp einer grafischen Modelloberfläche (siehe auch WaSiM-GUI)
2002 Erweiterung der WaSiM-Ausgabeoptionen
2005 - 2006
  • WaSiM-Adapter für das Flood Early Warning System FEWS der Firma Deltares (NL);
  • Erweiterung des Gletschermodells (dynamische Berechnung der Gletscherausdehnung);
  • Ausdehnung der Zeitschrittsteuerung auf das Grundwassermodell
  • Erweiterungen im Verschlämmungsmodul (u.a. Einführung expression parser);
  • Erweiterte Modellausgaben für Variablen der Vegetationsentwicklung;
  • Ausbau der Methoden zur Meteodaten-Interpolation (z.B. multiple lapse rates);
  • Erweiterte Regelwerke zur Abfluss- und Speicherregulierung (erhöhte zeitliche Flexibilität)
  • Implementierung 1D Vertical Heat Transfer Model, welches an das Richards-Bodenmodell gekoppelt ist (einschliesslich Auftauen und Gefrieren der Bodenschichten). Es wird Wärmeleitung udn Advektion berücksichtigt
  • Erweiterung im GLetschermodell: Berücksichtigung von schuttbedeckten Gletschern durch ein Debris-Grid
  • Neuer Abschnitt in der Steuerdatei für [special_output]: Individuell kann die Ausgabe fast aller intern verfügbaren Grössen für einzelne Zellen und/oder Bodenschichten zusammengestellt werden
  • WaSiM arbeitet nun intern ausschliesslich mit double-precision Werten, wodurch das Modell bis zu 20% schneller und z.T. genauer wird (z.B. im Grundwasser und Bodenmodell)
  • Erweiterung der MPI Funktionalität in einigen anderen Modulen (Gletschermodell, Grundwassermodell), um das Modell so bereit zu machen für den Einsatz auf Supercomputern (kann produktiv genutzt werden, einige Module arbeiten anch wie vor nur auf einer der Instanzen, z.B. surface routing und GW-Modell)
  • Im Schneemodell wurden 8 neue Methoden aufgenommen, die alle möglichen Kombinationen der neuen Optionen (auch mit alten Methoden) darstellen: wind driven lateral snow redistribution, gravitational snow redistribution (kleine Lawinen) und energy balance approach
  • WaSiM kann nun auch Stacks im ASCII-Format lesen und schreiben. Somit können alle Inputs als ASCII-Dateien zur Verfügung gestellt werden und auch alle Ergebnisse können als ASCII-Dateien erhalten werden (leichteres Postprocessing, aber auch längere Lese- und Schreibzeiten!)
  • Ein neuer Betriebsmoidus, welcher speziell für Spin-Up-Läufe für das Temperaturmodell nützlich ist, aber auch generell verwendet werden kann, wurde eingeführt: Der Hydrotop/HRU-Modus (HRU = Hydrologic Response Units oder Hydrotope). Weil dieser Modus vor allem für den SPinUp verwendet wird, kann die HRU-Generierung über den neuen, optionalen Abschnitt [SpinUp] gesteuert werden. HRUs können automatisch als eindeutige Kombinationen beliebiger Grids erzeugt werden (z.B. Höhenstufen, Landnutzungen, Bodenarten, Gefällebereiche usw.
  • Zur Unterstützung des SpinUp-Modus wurde ein Pre-SpinUp Modus für die Abschätzung der mittleren Jahrestemperatur Bodens implementiert. Vor dem normalen SpinUp-Lauf läuft lediglich die Temperatur-Interpolation (nur für die HRU-Masterzellen). Anhand des Ergebnisses wird die mittlere Temperatur der Permafrostoberfläche geschätzt und davon ausgehend die Temperatur des Bodens bis zur Untergrenze des Bodenprofils. Wenn kein Permafrost vorhanden ist, wird auf diese Weise die normale Bodentempertatur geschätzt (ausgehend von der mittleren Jahrestemperatur und dem geothermischen Wärmefluss).
  • Reading GIF Images (specifically: Radar precipitation measurements) as precipitation input, regardless of their grid domain and cell size. Values are then interpolated using IDW method to the model domain and model resolution
  • Soil depth and groundwater layer thickness are always identical: aquifer thickness is taken from soil tabel. ALso, easy scaling of soil thickness (and aquifer thickness) is possible per coild type by using a new optional parameter in the soil table
  • MPI memory optimization: until now, each instance used to allocate the complete memory that was required for a single instance. Now, each instance only allocates the stripe it is working on (for most grids and stacks - some grids are allocated completely in the master instance only and stripe-wise in the other ranks). The effect is dramatic: WaSiM can now run in much more instances per nodes (on big workstations and super computers) even for extremely large model domains with several GByte of allocated RAM for a single instance.
  • Some important bugs where fixed (interpolation and precipitation correction related and glacier runoff routing related bugs - see details in releasenotes).
  • Surface routing works with a linear set of equations and Gauss-Seidel solution algorithm now. This is much faster and lakes or other backwater playec have much more stable surfaces now. However, sionce this algorithm is a Finite Differences approach, some artifacts may show up (like water flowing sometimes over dry dams when using too les iteration steps)
  • Layered snow model introduced. This model component actually only works with heat transfer model switched on, since the snow layers are handled as soil and water/ice at the same time. So the temperature is quite important for all computations of heat and water fluxes
  • A new abstraction rule for reservoirs allows the physically based routing between hydraulically connected lakes (abstraction rule of type HydraulicConnection). Thus, fluxes can accure in both directions. Sample application: three peripheral Jura Lakes in Switzerland.
  • External coupling reactivated. Works now also with MPI. In addition to the classic functionallity, the semaphorfile optionally contains a command that can be interpreted by WaSiM or an external application (continue, cancel, save, full stop).
  • reading and writing of netCDF files implemented. netCDF files must follow the classic data format with one default group. Compression is supported for both writing and readin. Grid time series can be stored in a single, compressed netCDF file.
  • Re-infiltration is now available in MPI version by re-ordering river link cells in a separate grid row and processing all cells of a single link in a single MPI rank from highest to lowest cell, thus following the natural flow paths. This feature is also used in non-MPI-versions, so the re-infiltration will show slightly different results if compared to older versions. The main advantage of the re-ordering of river link cells and processing them separately from the normal cells is the fact that this can be used in MPI without too big performance losses (and that it is possible now at all to use re-infiltration in MPI).
  • Extended netCDF handling implemented: until now, the variable to read must have had the dimensions in order var(time,x,y), whereas the y variable was considered top-down. New behavior: var(time,x,y) and var(time,y,x) are possible. x-dimension is still always considered to be read from west to east, while y variable can be from south to north or vice versa. var(time,y,x) means, that the grid was written in rows to the netCDF file, starting either from the top or from the bottom of the grid. WaSiM will read the y-dimension variable and analyze the first two entries. If the second entry is larger than the first entry, the rows are considered to start at the lower left hand corner, while when the second y value is smaller than the first y value, the rows are considered to start in the upper left hand corner.
  • Extension of the simple interception model by a physically based snow canopy interception model (implemented by Kristian Förster, Uni Hannover, Germany) and by a canopy energy balance model for a snow filled canopy interception storage (implemented by Matthias Kopp, TU Munich, Germany). Both modules allow a much better snow pack modelling in forested areas, where the ground snow cover varies often quite remarkably from open field conditions.
  • Data assimilation as new optional sub model available. Grids can be read in at arbitrary times, e.g. snow water equivalent grids that are available on a weekly basis but only in winter months. SWE values can be split into liquid and solid components based on the present SWE at the actual or neighboring cells. Also, reservoir water table values can be read in to adjust for accumulated errors in inflow/outflow computation, thus bringing the model back to realistic lake water tables.
  • sub time steps of down to one minute can explicitely be used by adding a fifth time header column "mm" (for minutes) in each input file and in each output (graphlines and all time series processing tools are compatible with this change)
  • time variabel melt parameters for all T-index and T-u-index melt methods can be defined in the snow model section
  • command line option -mode:routing invokes the "routing-only" mode of WaSiM
back top home



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