Navigation: Home › Das Modell › Modellentwicklung › Entwicklungsdetails
Entwicklungsdetails ab WaSiM 6.4
| 2005 | 2006 | 2007 | 2008 | 2009 | 2010 | 2011 | 2012 |
| ! | = Bugfix (Patch), Fehlerkorrektur |
| + | = neue Fähigkeiten |
| - | = herausgenommene Features |
| o | = sonstige Änderungen |
 | 2012 |
| WaSiM 9.0.8, 01.05.2012 | |
| ==> Official Release | |
| ! | checking and correcting all error codes according to the new documentation |
| ! | calculation of ETP for open water: Radiation was not correctly applied (there was some radiation spared for soil heat flux, which is now used for evaporation too) |
| ! | DynamicPhenology_4: altitude correction is now working for non-dynamic sample days (was not working at all for DP4 previously) |
| WaSiM 9.0.7, 30.04.2012 | |
| + | Implemented implicit solution schema for heat transfer model: Now this scheme is delivering exactly the same results as the explicit solution. Thus allowing to take advantage of an implicit solution. There are almost no numerical instabilities when using much longer time steps. Tests revealed that for daily time step the sub time step for explicit solutions should not be larger than 1800 seconds, whereas for implicit solutions the time step can still be as large as a day without any substantial loss of accuracy. For hourly resolution the situation is similar: Explicit solution should not use larger time steps than 1200 seconds whereas for implicit solutions 3600 seconds are still ok. However, some artificial effects will be seen in both versions because of the larger temporal temperature gradients within a day. So shorter time steps are recommended as down to 900s for both methods when running in an hourly time step. For even shorter input data resolution even shorter sub time steps may apply. |
| + | sub-timestep-optimizer (based on the mesh Fourier number) for heat transfer implemented (and commented, because it is often useless when melting energy and it's fast change over time must be considered, so the old approach of the explicit solution is kept active). |
| WaSiM 9.0.6, 10.03.2012 | |
| ! | fixed a severe implementation error from version 9.0.4: Since the order of initialization for snow model and soil model was changed, there were some settings for snow model related entities which were set do their default values in unsatzonmodel (so the soil model acted as if no snow model was active). This is fixed now by changing the initialization order back to the original order and extending the snow model by a call to the same zone identification method as in the soil model, thus making sure that the snow model uses the same zones for glaciers as are used in the soil model. |
| ! | making sure that the glacier zones used in the soil model were taken from the correct grid: if the traditional ice_firn grid AND the new glacierized cells (fraction of cells with glaciers) were read in both, until now the soil model used the glacierized cells grid, no matter if the static or dynamic glacier model was used. Thus, if both grids were different (what is the case most probably) the resulting runoff could be wrong at all. This is fixed now: the soil model will also look for some spezial grids only available when using the dynamic glacier model prior to decide which glacier grid to use. |
| WaSiM 9.0.5, 28.02.2012 | |
| ! | When using Wendling Evaporation algorithm (Method 3) and having no radiation as input data (modelled by sunshine duration instead), the model crashed because internally the DEM was not retrieved for radiation modelling. |
| ! | corrected interpolation of HAMON coefficients (so far only a shift of 15 days was possible) |
| WaSiM 9.0.4, 27.02.2012 | |
| ! | thorough debugging of the 32-bit version: Chrashes when using the lake model and abstractions, and inflows could occur due to some wrong memory allocations when parsing the control file. Should be fixed now. |
| ! | glacier statistics and ELS handling in snow model: Now this will be controlled by the same zones as used in the unsaturated zone model and in the routing model, since the runoff of the glaciers must be available in those sub-models in the expected zonal structure. Thus, other zones may be used in the snow model for other statistics like snow storage etc. (all not [only] glacier model related statistics and all non-glacier-runoff linked variables even from the glacier model). There is no user action (i.e. control file change) required The addresses will be taken automatically from the unsatzonmodel. For this purpose the order of initialization of snow model and unsatzonmodel was changed internally. Now the snowmodel is initialized after the unsatzonmodel, since otherwise the addresses of the zone-related objects (taken form unsatzonmodel) were not yet available. |
| WaSiM 9.0.3, 21.02.2012 | |
| + | optimizing the performance of the heat transfer model: Now the sub time step for heat transfer and water transport are handled separately. If the sub time step for heat transfer is half or less of the water transfer sub time step, the heat transfer sub time step is made an integer factor of the water transport sub time step (i.e. if heat transfer sub time step would be 48 seconds and water transfer sub time steps would be 600 seconds, the integer factor would be 13 (600/48 = 12.5 ? rounded up to 13), so the effective heat transfer time step is now 600/13 = 46.15 seconds. This allows the heat transfer to be de-coupled from the sub time steps of the water transport. Now, in this example in every water transport sub time step, the heat transfer is completely computed 13 times, thus reducing the computational time for water transport by factor 13, which affects strongly the entire run time of the model (experimental results: speed up is approx. 2 times). |
| WaSiM 9.0.2, 18.02.2012 | |
| !/o | changing the parsing mechanism for routing descriptions from pure old C to C++ streams and strings. Thereby simplifying the code considerably and removing some implicit errors with copying strings. |
| WaSiM 9.0.1, 16.02.2012 | |
| + | a brand new sub-model is available: the 1D vertical heat transfer. It is coupled to the solution of the Richards equation. So it will be available in the Richards-version of WaSiM only. For this first version only heat diffusion is regarded (no advection by infiltrating water etc.). Also, no heat transfer trough snow packs, no lateral heat diffusion, no heat advection by lateral groundwater flows etc. are implemented. However, the temperatures in the soil column can be modelled now - and thus the effects of frozen soil on runoff generation, infiltration etc. There is one disadvantage: The computational times will multiply depending on the configuration (number and thickness of soil layers, especially thickness of the uppermost layers down to approx. 10 cm) by a factor of up to 10 to 20 compared to model runs without heat model. For daily time steps, the increase may be even larger, since internally the time step may be reduced to some seconds in order to keep the numerical solution schema stable. |
| 2011 |
| WaSiM 8.10.4, 14.10.2011 | |
| ! | check in unsatzonmodel if root depth is less than soil depth |
| ! | error fixed in precipitation correction when using regions for interpolation (no wind dependent correction was done for snow) |
| ! | fixed an error in week day dependent abstraction rules (on Mondays, no abstractions were taken from the reservoirs at all) |
| WaSiM 8.10.3, 20.09.2011 | |
| + | date dependent abstraction rules implemented |
| + | time dependent abstraction rules implemented |
| + | writing standard grids rightly after reading (in order to get the changes made by nearest neighbor-fill algorithm) |
| ! | corrected error in GLMB initialization when reading grids from disk |
| o | minor changes to error message when reading row and col parameters |
| o | OpenMP: in debug-mode, the output when a balance is wrong in a cell must be put into an OpenMP-critical section |
| WaSiM 8.10.2, 19.06.2011 | |
| ! | Make sure that no negative volumes and areas can occur. This also prevents the total runoff from -1.#IND values from glacierized areas with zero area. |
| o | writing balance errors into the log file: make this section critical for OMP (only important for debugging the application |
| WaSiM 8.10.1, 25.05.2011 | |
| ! | Expression parser results were not returned correctly to WaSiM values resulting in SDISP=0 etc. for all time steps. Error is fixed now. |
| WaSiM 8.10.0, 16.04.2011 | |
| ==> Official Release | |
| + | expression parser implemented into the silting-up model... More » |
| ! | fixed a bug when searching the control file for some optional parameters which could cause WaSiM to loop infinitely |
| + | extension of meteo interpolation method 12. Now three new parameters can be used optionally (as a group):
|
| WaSiM 8.9.5, 19.02.2011 | |
| + | added 8 new optional parameters to the control file for [SiltingUpModel]. The parameters are used for the following equations... More » |
| + | added a method in section [SiltingUpModel] in the control file.
|
| WaSiM 8.9.4, 13.02.2011 | |
| + | added a parameter to the permafrost model: k_sat_reduced will be taken instead of 1e-12 (taken in former versions as constant value) to define the infiltration capacity into frozen soil layers. |
| WaSiM 8.9.3, 09.02.2011 | |
| ! | some minor bug fixed for dynamic glaciers: firn modelling used n-1 layers instead of all n layers and SSTO and SLIQ were not set to 0 when making old snow to firn. |
| WaSiM 8.9.2, 08.02.2011 | |
| + | adding an optional parameter to the 5 regression parameters in elevation dependent regression method 10 for meteo interpolation. This is a fix lapse rate valid for the uppermost range. Applicable for regions where we have no information on, e.g. precipitation in high elevations. So we can set a constant lapse rate to be able to reduce precipitation in elevations of e.g. higher than 2500m. The new parameter set would look like this:
[precipitation] 10 # method $inpath//precip.txt AdditionalColumns=0 800 1400 200 1 300 -0.002 # regression parameters for precipitation: # lower inversion, upper inversion, tolerance, overlap (0/1), # max cluster range, lapse rate fix for upper regression... |
| WaSiM 8.9.1, 06.02.2011 | |
| ! | Irrigation method 4 (ETR/ETP > x): when transpiration was taken from interception and/or snow only, the soil transpiration was zero and thus ETR/ETP was always < x. This is fixed now: ETR/ETP will be 1.0 in these cases. |
| ! | Permafrost: transpiration extraction not possible from frozen soil now (below thaw depth) |
| ! | After thaw depth increases, the soil will get it's old k_sat value instead of keeping the reduced value of 1e-12 m/s forever. |
| ! | Now infiltration into first soil layer is controlled by the reduced k_sat value (of 1e-12) if thaw depth is within this layer. Before, the first soil layer could be saturated even if permafrost was in layer 1. This could have strong effects on soil water when the first soil layer was considerably thicker than a few mm. |
| 2010 |
| WaSiM 8.9.0, 15.09.2010 | |
| + | Now standard grids can be configured to create output statistics as well as to be written as any other grid (like variable grids) at any requested interval. This is useful for grids which will usually be created and changed internally by default. Now they can be read as standard grids (or even created and given a default value). More » |
| + | adding RSE, ALbedo, LAI[layer], VCF[layer], RSC[layer], RSI[layer], Z0[layer] and ROOT[layer] to the detailed output in the otherdata file for the control cell. |
| WaSiM 8.8.3, 24.08.2010 | |
| ! | Error in initialization of lakes: Lakes with a certain shape stretching from the upper right hand corner to the lower left hand corner with the upper left cell being the most left cell of this lake at all were initialized wrongly: this cell did not count as lake internally. |
| WaSiM 8.8.2, 06.08.2010 | |
| ! | Output of “otherdata” statistics: there was some rubbish beyond the last data column (depending on how many soil layers there were). These columns are removed now. |
| ! | Glacierized cells caused an error in the topmodel version when using no lake grid. |
| ! | Glacierized and unglacierized cells were swapped for balances. |
| WaSiM 8.8.1, 09.07.2010 | |
| + | Implementing a soil stability index into WaSiM-Topmodel version. Documentation on demand (will be part of the new documentation - coming soon). |
| WaSiM 8.8.0, 16.06.2010 | |
| ==> Official Release | |
| + | Sub-time step for the sub-model group "soil model, surface routing, groundwater" implemented in the Richards version. Thus, these three sub-models can run with a shorter time step (as the soil model and surface model already did) in a loop, e.g. in areas where the gradients are too high for the common model time step; for instance, when running the model in 1 m spatial resolution with strong local gradients of up to 0.5 in a moor or marsh region with high hydraulic conductivities in an hourly temporal resolution. Internally, soil model and surface routing model will keep using their dynamically estimated sub-time step, which has now an upper limit set by the new sub-time step of the mentioned model group; now it can not be longer than that sub-time step, but can well be shorter. The required model setting in the control file is a single parameter (the number of sub-time steps) in a new section [soil_surface_groundwater_substeps]. Values to start with are 1 (default), 2 (half of the common time step), 3 etc. Please set this value carefully. Too high values can decrease the model performance dramatically, since unsatzonmodel and surface routing are called each time! |
| ! | made warnings "no station in search radius" thread safe (module could crash when two or more messages were printed at the same time from different processors). |
| o | routing module: changed the WARNING that no oberved runoff is used for calculating the efficiency criterion into an INFO. |
| WaSiM 8.7.0, 25.05.2010 | |
| + | Dynamic glacier model implemented. This is a major extension to the existing glacier model:
|
| + | new interpolation method 12 (Thiessen + lapse rate). |
| ! | air pressure calculated for actual elevation for radiation extinction instead of set to 970 as fix value (only small effects on results). |
| ! | Topmodel version crashed because it could not find the MPI-regions (which are now created by default as in the Richards-Version). |
| ! | Topmodel version crashed when cancel the model run with option 2 because output stream was deleted before the last message was printed to the stream. This is fixed now. |
| WaSiM 8.6.3, 10.04.2010 | |
| ==> Official Release | |
| ! | bug fix for initialization of DynamicPhenology_2: When starting the model under cold start conditions, DP2 switches to DP1 and causes an access violation when trying to access a data structure not valid for DP2 but DP3 only. Now the data structure is accessed for DP3 only, so DP2 will work properly even when using DP1 for the first model year. |
| WaSiM 8.6.2, 25.03.2010 | |
| o | Changing diagnostic levels for some messages from WARNING to INFO in order to allow FEWS a smoother run. |
| + | Added an optional output statistics for thaw depth. To enable this statistics, after the last parameter (readgrids) in section [SurfaceRoutingModel] of the control file a common parameter set for statistical files has to be added (i.e. a filename or a variable containing a filename and a 4-digit code specifying the output like: $outpath//sfstsr//$grid//.//$code//$year $hour_mean # statistics for surface storage in mm per sub catchment) |
| WaSiM 8.6.1, 06.03.2010 | |
| ! | Initialization of shallow lakes from read in grids corrected for such cells which were fallen dry in the previous simulation. |
| ! | Handling of complete or incomplete filenames changed when using "/" instead of "\" in Windows environments. Now absolute filenames beginning with "\\" or "//" are also allowed (server-path pattern). |
| ! | When using WaSiM for interpolation only, now the land-use grid is no longer required to run the model at all. |
| + | Added a new command line parameter -diaglevel:< level >. < level > can be one of the following: DEBUG/Debug/debug: all output except some very special display specific output is also sent to the XML-Diag file INFO/Info/info: only error messages, warnings and infos are sent to the - XML-Diag file (information on start time, finish time, finished model day and some more) WARNING/Warning/warning: only error messages and warnings are sent to the XML-Diag file ERROR/Error/error: only recoverable errors are sent to the XML-Diag file. Actually no error is recoverable, so this setting will have no other effect than FATAL. FATAL/Fatal/fatal: only fatal error messages will be sent to the XML-Diag file. |
| + | Added a new command line parameter -diagfile: < filename >. This parameter gives the filename of the output XML-diagnostics. It is optional. If omitted, the default filename will be WaSiM_diag.xml |
| o | Some internal changes to the stream handling were made in order to allow suppressing of messages for the XML file (also the message headers like INFO:, WARNING:, a.o. are suppressed now). |
| ! | When reading the observed runoff file, the analysis of the columns is done on the first valid data row (which is the first row starting with a number). All other rows are expected to have exactly the same number of values like this first row. When this is not true, WaSiM may shift dates and values. This is still a potential problem, but now at least a misreading of a floating point value into an integer date (year, month, day or hour) field will be recognized as an error. WaSiM will stop with an appropriate error message. However, if the file contains integer values only, then no error will occur until the end of the file is read in. The best case will be that the shift in values is not resulting in the exact number of values for the last file row; so the described error may still occur. |
| WaSiM 8.5.6, 11.01.2010 | |
| ! | When a Pond was overrunning, QD was wrongly assumed to be in mm instead of m/s. This error led to extremely high runoff values for direct flow. It is corrected now. |
| ! | The storage file was appended to the old file and not re-created when no routing module was active. This is corrected now: When the unsaturated zone model is active only, then the storage file is created there, otherwise it is appended. |
| o | internal optimization: strstreams replaced with stringstreams in order to avoid deprecated functions |
| o | changed the internal stream class to make it compatibel with ANSI C++ (and thus operable under Linux) |
| ! | The snow and rain fractions of precipitation are no input state grids, so they are no longer required when starting with initial conditions. |
| 2009 |
| WaSiM 8.5.5, 16.12.2009 | |
| + | Implementation of an additional diagnostics output as XML file. This file with the fix name WaSiM_diag.xml will be created for each model run (and it will overwrite existing files with the same name). It contains the same information as the console output but in XML format and with additional attributes for the severity level: 0 = FATAL (non recoverable errors) 1 = ERROR (recoverable errors) 2 = WARNING 3 = INFO (most output as written on the console) 4 = DEBUG (currently not used) |
| WaSiM 8.5.4, 06.12.2009 | |
| ! | Corrected detailed output "otherdata" in unsatzonmodel: exactly 16 columns shold be written and not as many columns as there are soil layers. This is corrected now. |
| + | Surface routing: Now it is possible to rout surface runoff across sub-basin borders (but only if the other cell has a valid zone code too; so no routing into nodata-cells is possible). Sub-basin borders may change by rising water tables in very flat regions. Also this configuration allows some “special effects“ when running the model for comparison with some special locations: each location may then be represented by a single-cell sub-basin. However, no channel routing should be used for such applications with single cell locations, since a single cell basin should not be set as a river cell - surface routing would not work then for this cell. |
| WaSiM 8.5.3, 29.11.2009 | |
| ! | Because of limited numerical precision the comparison between the hydraulic head and the upper limit of an aquifer could result in confined conditions for the uppermost aquifer. This caused some saturated cells with high gradients to the neighboring cells to generate direct runoff when the aquifer thickness was set to 0 within the permafrost extension. Now, the aquifer is assumed to be confined only when the hydraulic head exceeds the elevation of the upper limit of the aquifer by at least 0.01m. |
| WaSiM 8.5.2, 22.11.2009 | |
| + | Groundwater model extended by an active layer algorithm for deep permafrost: lateral flows in the thawed soil layer may be calculated with the groundwater model now. The transmissivity for the aquifer is limited to the thickness of the active layer. This will NOT work for shallow permafrost but only for deep permafrost. In multi-aquifer models, the lower aquifers are not affected by this change (the deeper aquifers are assumed to be unfrozen all the time). |
| WaSiM 8.5.1, 17.11.2009 | |
| ! | Content of irrigation grid was always 0 because the grid was written after setting the irrigation amount to 0 in unsatzonmodel again. Now it is written directly in the irrigation model after calculating the irrigation. |
| WaSiM 8.5.0, 29.09.2009 | |
| ==> Official Release | |
| + | New method for demand driven irrigation implemented: method 4 does not use critical suction values for irrigation start and stop but rather the ratio of ETR/ETP... More » |
| + | New module [RegionalSuperposition] implemented: short description of the new module... More » Note that the use of the region grid in its old way is no longer possible! |
| + | Optimization of model performance: - unsatzon-model: implemented look-up table for h(theta) and k(theta) search as binary search (speed up will be up to 20%!) - exposition dependent radiation correction and shading: now shading algorithm works like an active ray tracer (formerly as a passive shading algorithm). Speed up may be up to 10% for very large grid dimensions (e.g. 1 Million cells). |
| + | Permafrost algorithms completely implemented and tested. |
| + | Some preparations have been done for integrating WaSiM in FEWS-FOEN, the Flood Early Warning System used at the Swiss Federal Office of Environment: - Now elevation dependent regressions are directly calculated without further need of using the helptools regr or regress. - Now write grids can be written in ESRI ASCII format. - Now default input and output directories can be defined. |
| ! | (27.09.2009) An internal temporary variable in the snow model was declared as integer but should have been float. This lead to some sudden but short peaks in the snow melt. This is fixed now. |
| ! | (25.09.2009) Thiessen interpolation used wrong coordinates when an anisotropy was defined, resulting in incorrect interpolation results. Now the Thiessen method uses unrotated coordinates again. This error was already some years old, obviously nobody used this method or if, then only for a single station. In the latter case is no difference in using the rotated or the unrotated coordinate system(). |
| ! | So far the runoff dynamics in surface routing was limited by the amount of water stored in the active cell (however, water from other cells could flow into each cell, but the outflow was limited). Now this limitation is overcome. But there is another point to pay attention: since WaSiM does not use a hydrodynamic flow model (but kinematic wave instead), the parameter for setting the maximum allowed sub-time step has a second meaning now. After the time given with this parameter, all flow velocities are newly calculated. Thus, it is possible to pour the entire grid cell content within the period of the major sub-time step (the maximum allowed one). All following sub-time steps will then operate with the new velocities and thus with the new water content until the next major sub-time step is reached. Consequently, some calibrating work is needed for this parameter now (as well as for the Manning roughness parameter for each land-use type). |
| WaSiM 8.4.3, 19.08.2009 | |
| o | Now MPI-regions are defined without taking into account the river cells. This ensures a fast gathering of grids and stacks. The disadvantage is that re-infiltration from rivers into lakes can not be used in the MPI-Version. However, this may not be dramatically since usually the effects of re-infiltrating water are small. |
| ! | Evaporation when using HAMON: wind speed is not needed as input any longer (was never used, but searched for). |
| ! | Groundwater when using multiple layers: problems are fixed for the case that the uppermost layer falls dry and the lower layer also becomes unconfined. Until now a leakage was calculated from the uppermost to the second aquifer. This was because of the pressure head in the uppermost layer which was limited to not fall below the bottom of the layer. When the head in the lower aquifer decreases below the boundary between the two layers, there was calculated a positive head difference, resulting in a leakage (even if there was no water at all...). This is fixed. Now the groundwater model has a zero balance even when running with multiple layers. |
| WaSiM 8.4.2, 09.08.2009 | |
| ! | Now explicitly looking for „temperatur_14“ and „air_humidity_14“ as these meteorological variables are used in Haude-evaporation Schema. This also implies looking for „temperature_14“ in the snow model as well when no other temperature was found. |
| ! | When using lakes but not explicite the lake model (by using a lake grid), then the model ran into an access violation when the lake were overflowing. This is fixed now. |
| + | Now the model components surface routing and channel routing will be skipped if the model is run in MPI-mode and the actual instance is NOT of MPI-rank 0. This means that only the instance with rank 0 will run these submodels. Thus, MPI-implementations may run with all submodels now (which has still to be tested with irrigation, surface routing, and silting up). |
| WaSiM 8.4.1, 10.07.2009 | |
| + | (10.07.2009) Time step splitting for reservoir/lake handling implemented. This is especially usefull for small reservoirs in case of floods when a long time step would cause strong nonlinearities in volume-outflow-relations (e.g. the reservoir could be emptied in one time step and filled again in the next one or never reach the maximum filling at all, even if the flood is serious. Changes in the control file are... More » |
| + | (05.07.2009) Permafrost Thawing algorithms implemented (not yet the core in snowmodel, but all the input/output stuff as well as the setting of the saturated hydraulic conductivity within the unsaturated zone model). Changes in the control file are... More » |
| ! | (26.06.2009) bugfix in handling of POND: when MAXPOND is 0, then this spot is NOT handled as a POND at all. Hence infiltration excess may discharge as direct runoff whereas it was set to 0 before and stored into the soil somehow. |
| WaSiM 8.4.0, 24.04.2009 | |
| + | (24.04.2009) Now storage capacity s0 (drainable porosity) of the groundwater model is calculated in the unsatzonmodel in order to have a better match between both modules. However, the s0-value is calculated for the soil layer where the groundwater table is actually within, so different theta-values in other layers may still lead to oscillations, but with a far lesser aplitude than before. |
| + | A new parameter for non-persistent default grids is added: the last parameter says for how many intervals the new default grid should be valid until the original default is valid again (this is important when reading a new constant head boundary with rising groundwater tables, because the s0-value as described above may not be as expected. When making the boundary condition grid valid for 2 or 3 intervals, the oscillation will be damped sufficiently. |
| + | (23.04.2009) The option persistent=0|1 for standardgrids was introduced (in combination with option periodicity). This means, that a non-persistent standard grid will be valid for a single time step only until it is replaced by the old standard grid read in at model start... More » |
| + | (20.04.2009) MPI in snowmodel was implemented (still except for the glaciermodel!) |
| + | Implemented MPI-compatibility in irrigation module: since MPI is not suiteable here, the inputgrids must be gathered and (after this submodel) be scattered to all the instances. |
| + | The reduce function for statistics was changed to a function using a custom MPI::Op-style function. |
| + | (19.04.2009) MPI was extended up to evaporation (potential, snow, real). The longwave outgoing radiation calculation was changed: the cloudiness at night, which is at day estimated by the sunshine duration, was formerly estimated by taking the average value of the sunshine duration for all intervals of bright day for an entire subbasin. Now, it is the average value of each cell only. Thus, we are avoiding synchronization overhead between MPI-instances. And anyhow: using a subbasin mean value is a somewhat artificial limitation, so it is even more to the basic concepts of gridded modelling now. |
| + | (16.04.2009) Now MPIGatherFloatGrid and MPIGatherDoubleGrid is available. |
| + | (13.04.2009) First experimental MPI version of WaSiM available. Note that Message Passing Interface (MPI) is needed to use WaSiM in computer clusters and supercomputers. As a first step MPI was implemented for the interpolation of meteorological data, except for writing the result grids -> meaning that it cannot be used for real modelling until at least the results of the interpolation are gathered in rank 0 instance. |
| WaSiM 8.3.2, 17.03.2009 | |
| ! | When reading grids in the [unsatzon_model] with readcode 2, then not only the groundwater table below surface will be read in but also the POND grid (if ponds are enabled). This will lead to a proper initialization of lakes. Because lakes are in very close connection to groundwater. This feature is rather a bug fix then a new feature. All other grids are initialized depending on the groundwater table (i.e. soil water content, hydraulic heads, etc.). |
| WaSiM 8.3.1, 11.03.2009 | |
| ! | When initializing the unsaturated zone, a cell with a lake code was always set to saturated, even if the lake was not fully extended. So the margins of the lake which were empty and dry had a saturated soil which led to a wave of groundwater exfiltration into the lake after model start and to a suddenly much higher mean groundwater table for the subbasin than in the old model run. Now only those cells which are belong to a lake and have a positive POND volume will be initialized in this way. When old grids are read in, this should be the state of the old grids anyway. |
| WaSiM 8.3.0, 02.03.2009 | |
| ! | Balancing of GW: Until now the last (lowest) grid row and the last (rightmost) grid column were not regarded in calculations of fluxes for tracer mixing and balance checks. |
| ! | Balancing of GW: Storage change was wrongly subtracted from balance of fluxes (added instead of subtracted), so the balances were double the total flux. |
| + | Grid for groundwater balance has been implemented -> It must be added in the control file for each aquifer after the flow grid together with the writecode for the grid, which should be 10 or 20 or 13 or 23 in order to have the sum over the model period or each year, respectively. |
| + | Internal precision in groundwater model has been increased from single to double in order to allow much smaller errors in mass balances when small head changes occur only; e.g. an actual head change of 0.1mm at a total head of 2000m a.s.l. would require more than 7 significant digits for keeping precise results. |
| + | New balance check of groundwater implemented: The change in storage is calculated based on the change in hydraulic heads and the storage coefficient. Additionally, the fluxes into and out of the element are balanced as div(grad(h)). The difference of both entities is taken as the calculation error. The error should be near zero, unless there are constant head boundary conditions. Constant flux boundary conditions will also be considered. The averaged annual sum of those balance will serve as a water balance correction term when using boundary conditions in the groundwater model. Attention: This method cannot be applied when an external groundwater model should be coupled. |
| + | Periodical re-reading of standard grids during the model run has been implemented. For setting up this feature, the paragraph [standard_grids] may contain optional parameters for each grid: after the name, the identification and the fillcode, three more parameters may follow a keyword "periodicity": the number of periods to wait until a new grid is read in, the period (H=hour, D=day, M=month, Y=year) and the number of sub intervals of this period. "periodicity 2 D 12" means that each other day (second day, i.e.) at hour 12 a new grid is re-read. The grids must have a filename matching the original filename of the standard grid extended by the date-time-scheme. In the above example, the grid name may be m500.bh1, so the next standard grids must have the names: m500.bh1_1984010212 m500.bh1_1984010412 m500.bh1_1984010612 m500.bh1_1984010612 etc. This feature will be interesting for reading known changing boundary conditions for the groundwater model as well as for changing landuse in long term simulations. The only prerequisite is the creation of all the needed grids in advance to the model run. If a grid is missing, the model run will stop with an exception (no results are written!). |
| WaSiM 8.2.11, 19.02.2009 | |
| ! | Macropore infiltration regarded for INFEX output and balancing of PONDs was incorrect since version 8.2.6. |
| ! | QD from snowmelt was added two times under certain conditions (since 8.2.1). |
| ! | Balancing of PONDs was incorrect since 8.2.1, even in 8.2.10. |
| WaSiM 8.2.10, 15.02.2009 | |
| ! | Unsatzonmodel: Now the POND grid is correctly filled (there was a bug introduced in version 8.2.6). |
| ! | GW-Model: The intermediate head for GW-level was written to a wrong grid. As a result of this the iteration took too long. Now corrected! |
| ! | GW-Model: When using boundary conditions (constant head in *.bh*-Grid), then the groundwater level was not correctly written out (except when initializing the unsatzonmodel with readgrids=2, so the real groundwater surface was taken for initialization which should have matched the boundary condition, too). Now, the boundary condition will be regarded correctly. Attention: There is still a known problem! When using certain van-Genuchten parameters, then the groundwater level may strongly oscillate when using such boundary conditions. The work on this issue is still in progress. |
| WaSiM 8.2.9, 10.01.2009 | |
| ==> Official Release | |
| ! | Altitude dependent shift in phenology: fix initialized with -0.025 -0.025 0.025 0.025 for the old style version of the land-use table. Now it is changed to 0.025 0.025 -0.025 -0.025 (so the phenophases will start later in higher elevations and end earlier instead of start earlier and end later). |
| ! | Calculation of snowmelt: the melt amount was defined internally as an integer value, so all melt amounts were rounded to integer values (which was not obvious from the statistics because there we have an average value for a subbasin). Now the melt amount is a floating point value again (bug existed since parallelizing WaSiM with OpenMP in Version 7.5). |
| + | WaSiM-Helptools for Linux compiled for 64-bit systems -> provided for download according to the latest documentation (1999, which is still valid). |
| 2008 |
| WaSiM 8.2.8, 09.12.2008 | |
| + | New grid write codes in lowest significant digit: 6 = at each end of a month, 7 = at the end of a year. |
| + | New naming scheme for all periodically written grids: |
| ! | DynamicPhenology 1, 2, 3, and 4: Theta-Scaling default was 1.0, which was not compatible with the new algorithm. Now, Theta-Scaling defaults to 0.0. When running the model without soilmodel (or with topmodel), then Theta-Scaling of 1 may have resulted in too large forcing units (in fact, a value of 1 means a doubling of forcing units). |
| WaSiM 8.2.7, 05.12.2008 | |
| ! | DP4: Start day for vegetation period was calculated correctly but internally never used for evaluating entities with vegetation period dependent average or sum grids. --> corrected. Attention: When using DynamicPhenology_2 or DynamicPhenology_4, be sure to start the model run on or better at least one day before DP2/4_t0_dorm, e.g. on October 1st, when DP4_t0_dorm is October, 2nd. Otherwise WaSiM would use DP1 also for the following model year (because dormation with its chilling units could not be calculated correctly). |
| WaSiM 8.2.6, 02.12.2008 | |
| ! | Bugfix in initialization of Dynamic Phenology 4: MTSTART-stack was only created if at least one other landuse had DP2. |
| ! | Time of starting leave fall was set to the default value for following model years (not the first, but the next ones). |
| ! | Bugfix in macropore infiltration; too much infiltration was calculated due to wrong comparison between infiltration amount in mm and capacity in m/s. |
| ! | Bugfix in Qdir-Statistics: INFEX counted twice; now corrected to count only once. |
| ! | Silting-up/unsatzon/surfacerouting: when no silting up was active but surface routing was, the surface routing storage was not added to the precipitation input. |
| WaSiM 8.2.5, 30.11.2008 | |
| + | Surface runoff routing: kelsqd and kelsqi from unsatzonmodell are now used. This was not entirely implemented yet (only a default value of 3 for kelsqi ond nothing for kelsqd, meaning QD was only transformed by the translation effect, once a river was reached but not by SLR-effects). |
| + | Unsatzonmodel: if surface routing is active, then all surface runoff and interflow is now NOT routed into the SLR-cascade but rather statistically accumulated "as is". The transforming takes place in the surface routing model only. Thus, the balances may be checked much easier. Also, the difference in interflow statistics between unsatzonmodel and surface routing will give a hint on the effect of the SLR-cascade. |
| ! | DynamicPhenology4: actual Julian day was internally compared with sample point ID for start leave fall instead of with the corresponding Julian day to this sample point. This lead to unpredictable behaviour in certain cases. |
| WaSiM 8.2.4, 27.11.2008 | |
| ! | Bugfix in generating lookup-tables for van-Genuchten-parameterization when using Par_m (default: m=1-1/Par_n --> no problems there). This error may lead to -1.#IND-results in the soil. |
| WaSiM 8.2.3, 23.11.2008 | |
| ! | Some bugs in setting default parameters for parameters not given in land-use table corrected (value for last sample point was not read in for some variables). |
| ! | Modul silting-up: bug in equation for calculating grain size distribution (*+ instead +) corrected. |
| ! | Modul surfrout: for tau-calculation: slope was used as sqrt(S), should be S --> only. corrected. |
| ! | Modul unsatzon: tracer balance for surface runoff: QDSU now considered. |
| ! | Modul unsatzon: QDGrid: now INFEX considered correctly (was set AFTER setting QDgrid, so it was always zero before). |
| WaSiM 8.2.2, 12.11.2008 | |
| ! | Output of infex and qinf in unsatzonmodell corrected. |
| WaSiM 8.2.1, 09.11.2008 | |
| ! | Outputs of infex and qinf in [unsatzon_model] have been corrected; now |
| ! | Bug in DP4: when no cell has DP4, then an exception is thrown when accessing the moving average grid (now it is switched off when no DP4 is present). |
| + | Silting Up: Algorithm only works for cells which have a valid TillageDays-entry in the land-use table. However, the submodel will be called for all cells. |
| WaSiM 8.2.0, 04.11.2008 | |
| + | Implementation of DynamicPhenology4 (DP4): This method allows to estimate the start of the leave fall. Required new parameters in the land-use table are... More » |
| WaSiM 8.1.2, 26.10.2008 | |
| - | Conditional compilation removed for surface routing and silting-up module [#ifedef #endif - block definitions of _UNIBW_ removed]. |
| ! | Changed measure of INFX from [m/s] into [mm/timestep]. |
| ! | Calculation of weights for multiple flow paths corrected. |
| + | Calculations of flow directions and aspects have been changed. Now all 8 neighboring cells are regarded to calculate the main and minor flow directions. |
| + | New method parameter for flowpath calculation: 1 = multiple flowpaths without regarding water content, 2 = single flowpaths without regarding water content, 3 = multiple flowpaths with regarding water content, 4 = single flowpaths with regarding water content |
| + | New algotithms for silting-up module have been implemented (taking into account stones, other infiltration rates, a.s.o); GrainSizeDist must contain values for clay, silt, sand and stones1, stones2, stones3, stones4 (the last ranging from 2-6.3mm, 6.3-20mm, 20-63mm and 63-200mm resp.) |
| + | Soil table in the control file has been extended by the van Genuchten parameter m. Par_m can be defined for each soil type (analogous to Par_n); default value: Par_m = 1 - 1/Par_n (if omitted). |
| + | A minimum flow sheet depth for regarding crop roughness for surface routing can be specified. |
| WaSiM 8.1.1, 04.10.2008 | |
| + | Implementation of "Silting Up Modul" has been completed (Alpha-Version). |
| + | Now INFEX grid and statistics are written out with the amount of non-infiltrating water from the soil model. [Before: INFEX was always zero since Green-Ampt approach was switched off.] |
| + | Now filenames for layered entities (e.g. interception evaporation or potential and real evaporation) can be created with relative pathnames. [Before: The dots in the pathnames confused the algorithm to add the layer-ID to the name.] |
| WaSiM 8.1.0, 19.09.2008 | |
| + | Implementation of surface routing with dynamically calculated roughness: either use conserving (e.g. meadows) or non-conserving (e.g. conventional corn crops) land-use types or other types (e.g. forest) to parameterize k-modelling... More » |
| ! | Radioactive decay for tracer concentrations for groundwater recharge within the unsaturated zone will be only done after mixing with groundwater for unsaturated zone layers (before: The decayed and then mixed concentrations were transferred to the groundwater model leading to a faster decay (appr. 30-50% faster). |
| ! | Clipping drainage runoff to the maximum allowed runoff from that layer (drainable water as for interflow). |
| WaSiM 8.0.11, 02.08.2008 | |
| + | New writegrid codes: 1xxx and 2xxx are now available to indicate that an additional grid should be written with balances for the vegetation period only. The dates for start and stop of the vegetation period are individually calculated for each cell by using the landuse parameters and lokal climate... More » |
| WaSiM 8.0.10, 30.07.2008 | |
| ==> Official Release | |
| ! | During initialization no water is taken from reservoirs. Now the start of a new model run can be smoothly connected to an older run. |
| ! | No initial routing is done when readgrids is > 0, since the file "storage_richards.ftz" already contains all the channel content. |
| WaSiM 8.0.9, 29.07.2008 | |
| ! | Fixed initialization of ETR for lake evaporation (so far incorrect when using dynamic time steps). |
| + | Waterlevel-volume-relation (W-V-function) for lakes has been tweaked: the density of sampling points was increased from 5cm(0..2m), 20cm (2..12m) and 50cm (>12m) to 1cm (0..4m), 2cm (4..14m) and 5cm (>14m), respectively. Thus, the errors in the balance for lakes will be minimized. |
| + | Basics for vegetation period based statistics implemented: using new landuse dependent parameters like StartVegetationPeriod/StopVegetationPeriod and ResetVegetationPeriod enables the model to calculate the starting point, the endpoint, and the duration of the vegetation period. |
| WaSiM 8.0.8, 19.07.2008 | |
| + | Irrigation amount will be available as grid output now. The name and the writegrid-code must be given as two new parameters in the unsatzon model (not in irrigation!). So this affects all users and all control files. |
| ! | Simply add the suggested new lines for the file name and the write code and add the variables if symbols are used. |
| + | Evaporation from lakes: LAI and VCF as well as RSC for lakes are no longer needed. rsc is now set to 0 internally and the algorithms for calculating the effective resistancies are aware of possibly zero values. |
| + | Also for LAI and VCF values of 0 are now handled specially: no divisions by zero are possible now. |
| WaSiM 8.0.7, 17.07.2008 | |
| ! | When running without interception model, an access violation could occur in etpetr (because it tried to use interception grids. |
| ! | When using sunshine duration with values not equal to 0 for a dark period of the day, radiation exposition correction could deliver invalid corrected temperatures. |
| ! | Interception with one vegetation layer did not use potential interception evaporation until now but rather used the old interception model. Now, when using one layer of vegetation only and using method 2, the potential interception evaporation will be used instead (correctly). |
| WaSiM 8.0.6, 24.06.2008 | |
| ! | Correction of subcatchment-mapping in routing module (for lake modelling) when subcatchments were not enumerated from 1 to n in steps of 1. |
| ! | Same correction in tracer mixing for a special case of foreland flooding. |
| WaSiM 8.0.5, 15.06.2008 | |
| ! | Write a lot of more information for error 152 to be able to find misconfigured landuse or soil types faster. |
| WaSiM 8.0.4, 04.06.2008 | |
| ! | Balance for EVAR corrected (was calculated based on cells with EVAR>0 only, so the average was valid for those cells only, but not for all cells, and thus much too high). |
| ! | Correction for Linux when reading in the abstraction rules for reservoirs and rivers. |
| WaSiM 8.0.3, 24.05.2008 | |
| + | WaSiM Topmodel version is now also able to use the new lake model. For lake cells, no Topmodel approach is applied (like for glacier cells), so also no groundwater flow balance will be taken into account for lakes in this model version. Note, that the pond grid has NO connection to the soil, so ALL precipitation will be stored in the pond grid and ALL evaporation (only ETP, no EVAR/EVAP here!) will be taken from the POND grid, too. |
| WaSiM 8.0.2, 17.05.2008 | |
| + | Check for connection of ZL(modus=intern_with_rule) against AL (modus=intern_with_rule) (not only against AL(modus=intern). |
| + | Adding AL amount to the target subbasin for modus=intern_with_rule (like for modus=intern). |
| WaSiM 8.0.0, 01.05.2008 | |
| History of lake model and time variant abstractions: | |
| 24.04.-01.05.2008: Test and documentation, compilation of release and linux versions, regression tests against the older versions. | |
| ! | (24.04.2008): Now irrigation code 0 will not reduce tracer concentrations any more. |
| ! | (20.04.2008): Corrections in Alpha 8.0.0 (abstractions and irrigations could lead to negative reservoir content or at least to balance errors) |
| + | (19.04.2008): Taking into account targetCapacity for internal abstractions. Variable abstraction rules implemented for reservoirs and abstractions (still untested). |
| + | (16.04.2008): ReadAbstractionRule implemented for Abstractions (not yet for reservoirs). |
| + | (15.04.2008): Routing-model: added a check for correctly connecting internal abstractions with internal inflows. Each internal abstraction must be connected to an internal inflow in another sub-basin of the model domain. |
| + | Setting targed data (sub-catchment area and relative area as well as sub-catchment number) in abstraction class in order to be able to limit the abstraction when the target channels capacity is reached. |
| ! | Deleting semaphore file after reading grids (to gain some time). |
| ! | Now irrigation amount is subtracted from the last reservoir in a routing description. |
| + | (13.04.2008): Implemented abstraction as read in from file. |
| ! | For using observed discharges: When discharges are read in for routing (either for using observed runoff for routing input and/or for calculating r-square), the input file is checked in two ways now: |
| 1. Only existing columns will be used (if there are more sub-catchments than columns, then the column mappings for such sub-catchments with mapped columns which does not exist in the input file are ignored. | |
| 2. If a "measured" value is -9999, this value is ignored, too. | |
| + | (06.04.2008): Routing: AL_class implemented (refactoring of AL, which was handled within routingclass before); |
| Now AL-handling in routingclass before SP handling. So the AL amount is taken from the storage in the same interval; | |
| SP_class: take abstraction_sum amount into account for balance and error calculations and write it to the result file, too. | |
| + | (04.04.2008): Corrections in routing model (interpolation in W-V-Table), some bugfixes in routing for SP elements. |
| ! | Balance for SP elements corrected by irrigation amount (was not correctly taken into account for balance output only). |
| + | (24.03.2008): Initializing PONDs and Vakt-Values from each other, depending on readgrids, including tracer-concentrations. |
| + | (21.03.2008): Partly implementing the initialization of POND grid from routing. The Grid handling is implemented into abflussmodellklasse, routing_model and SP_class (which is newly derived from SP_struct) and handles now all reservoir- or lake-specific functions as methods. |
| + | (18.03.2008): Lake Evapotranspiration depending on POND content implemented in evapotranspiration model (getETP). |
| + | (10.03.2008): Initialization for etpetr, unsatzonmodel and gwflowmodel implemented. |
| Other changes/bugfixes: | |
| ! | (04.03.2008): Inflow files and observed runoff were not processed if they were write protected. This was changed, both file types may be write protected now. |
| ! | (04.03.2008): Preventing the model from a crash when no observed runoff was found but the parameter "use observed runoff" was set (method 2 in routing model). |
| WaSiM 7.10.3, 02.03.2008 | |
| ! | Another bugfix in phenology: z0-values of 0 still caused -1.#IND-values. Now, such values are limited to >=0.01m. |
| + | Buffering input of observed runoff data and external inflows similar to the input of meteorological data. |
| WaSiM 7.10.2, 28.02.2008 | |
| ! | Bugfix for dynamic phenology 3: When entering a new year, the initial forcing sum is set to 0 for interpolation of phenology now. Thus values of -1.#IND are avoided in ETP (caused by negative values for some phenology parameters due to interpolation between the last forcing value from the last year and the (larger) forcing value in the new year, e.g. for z0, which then was input into a log()-function, resulting in an error for negative arguments). |
| WaSiM 7.10.1, 16.01.2008 | |
| ! | Reading parameters of landuse and soil tables between "{" and "}" only. |
| ! | Corrected initialization of baseflow when groundwater model is disabled. Without this correction extremely high discharge values can be calculated during initialization and also for the first model time step. |
| ! | Coupling external groundwater table: If the external groundwater table is below the soil column, the groundwater table is internally set to the bottom of the soil column (prior it was tried to set the groundwater table to even lower layers, resulting in access violations or other errors). |
| + | Error handling for reading grids: multiple tests when grid cannot be opened. |
| + | Grid reading: If nodata values in input grids are smaller than -9999, they are translated into -9999. This may be important for external coupling, when the other model generates nodata values of e.g. -1e+038 |
| WaSiM 7.10.0, 06.01.2008 | |
| ! | Tracer handling in [unsaturated zone]: some bugs in balances for tracers were fixed... More » |
| + | Implementing two new modules: silting up and advanced surface routing as empty frames. It is implemented all the handling of fluxes for water and tracers between the modules. Solely the algorithms of the new modules are still missing... More » |
| 2007 |
| WaSiM 7.9.11, 13.12.2007 | |
| ==> Official Release | |
| ! | When generating QB without groundwater model and when using a dynamic timestep at the same time, the baseflow was generated without consideration of the splitfactor. |
| WaSiM 7.9.10, 08.12.2007 | |
| ! | If there were additional regions in the south of the last zone in the zone grid, then the region check failed. |
| WaSiM 7.9.9, 28.11.2007 | |
| ! | Bugfix in output of R-square at the end of the routing output file: When using a selection of subbasins AND discontinuous subbasin codes (e.g. 12,17,28 without subbasins 1 to 11 and 13 to 16 a.s.o) then the R-square values were complete rubbish until now... |
| WaSiM 7.9.8, 21.11.2007 | |
| + | Two new parameters implemented into the phenological method "DynamicPhenology_2": DP2_Offset_1 (defaults to -3.4 if not present) and DP2_Offset_2 (defaults to 10.4 if not present). |
| WaSiM 7.9.7, 18.11.2007 | |
| + | Now the POND grid can be considered in the [groundwater_flow] model in order to calculate the hydraulic head and lake growth by exfiltrating water. |
| WaSiM 7.9.6, 09.11.2007 | |
| ! | Fixing a bug for statistics output: When using a set of subcatchments for output only (codes 5 or 6) and at the same time the subcatchments in the zone grid are not numbered from 1..n without missing value, then the output was for invalid zones or even complete rubbish. |
| ! | Extended error handling when using a river link code without corresponding routing description. |
| WaSiM 7.9.5, 05.11.2007 | |
| ! | GW-Reinfiltration was handled with wrong riverflow_table-indices. Now, the indirection for river-links works correctly --> caused segmentation faults under LINUX. |
| ! | Regression with multiple regions: Now zonegrid masks the region grid too in order to avoid errors due to distant regions which have a positive weight in the weighting grid but no region code. |
| WaSiM 7.9.4, 28.10.2007 | |
| ! | Korrekte Ausgabe der Grids gwh1..n, gwx1..n und gwy1..n; diese Ausgaben wurden im Multithread-Modus unterdrückt. |
| ! | Nochmalige Korrektur der ksat_topmodel/ksat-Eingabe (war evtl. ein Linker-Problem). |
| ! | Korrektur für Multithread-Modus bei [unsatzon_model] Methode 3 (dynamischer Zeitschritt). Hier traten Fehler auf, weil eine Variable nicht threadsafe war. |
| WaSiM 7.9.3, 23.10.2007 | |
| ! | Korrektur an GWIN/GWEX-Berechnung, wenn die Bilanzfehler durch numerische Ungenauigkeiten in einem Iterationsschrit im Richardsmodell mithilfe von GWIN/GWEX kompensiert werden. |
| + | Vorbereitung eines detaillierteren Loggings mit Hilfe eines in der Steuerdatei festlegbaren Parameters - der aber bisher noch nicht eingelesen wird - und auch Logeinträge werden noch nicht erzeugt. |
| WaSiM 7.9.2, 18.10.2007 | |
| o | Ausgabe ETR[0] als Grid und als Statistik: enthält nun die totale reale Evapotranspiration als Summe aus ETR[1..n] + EI[1..n] + EVAR + EVRS |
| WaSiM 7.9.1, 16.10.2007 | |
| ! | Korrektur bei Verwendung von Regionen: Steuerdatei wurde geschlossen und anschliessend wurde der Abschnitt für die externe Kopplung nicht gefunden –> Steuerdatei wird jetzt wieder geöffnet vor der Suche nach den Kopplungsparametern. |
| WaSiM 7.9.0, 30.09.2007 | |
| + | dynamischer Zeitschritt in Richardsmodell implementiert: neue Methode für Unsatzonmodell: 3=dynamische Zeitschrittsteuerung |
| zusätzliche Parameter: | |
| o | minimum dynamic time step in sec |
| o | result statistics of the number of substeps |
| o | result grids containing the number of substeps |
| ! | ksat_topmodel wird nun nicht mehr irrtümlich gelesen, wenn ksat gesucht wird. |
| WaSiM 7.8.3, 20.09.2007 | |
| + | Stressfaktor für dynamische Phänologie: fehlt er oder ist er <= 0, dann wird ohne Stressbeeinflussung gerechnet. |
| ! | Temperatur für dynamische Phänologie wird nicht mehr in Expokorr und Verdunstungsmodell verändert (dynamische Phänologie wird VOR Expokorr und ETP berechnet -> bei Stundenschritten also im Intervall 24 des Tages) |
| WaSiM 7.8.2, 10.09.2007 | |
| ! | Korrektur an ETR-Entnahme aus gesättigten Bodenschichten, bei negativen Differenzen zwischen aktuellem und Grenz-Theta |
| WaSiM 7.8.1, 23.07.2007 | |
| + | Bodenfeuchteabhängigkeit für DynamicPhenology implementiert |
| ! | Korrekturen an der Berechnung des GW-Standes, um Schwingungen im GW-Stand und in hhyd, qu und theta zu vermeiden; Entnahme von Verdunstungswasser aus oberster gesättigter Schicht entgegen der Entnahmedichte, wenn Wurzeln ins GW ragen. |
| WaSiM 7.8.0, 30.06.2007 | |
| + | Dynamic Phenology 3 implemented: läuft mit Topmodel- und Richardsversion |
| WaSiM 7.7.1, 22.05.2007 | |
| + | neue Interpolationsmethode 9 zum Einlesen von Grids anstatt zu interpolieren: die Tabellendatei enthält Dateipfade und -namen der einzulesenden Grids. |
| WaSiM 7.7.0, 06.05.2007 | |
| + | Sequential model 24b for dynamic phenology implemented |
| + | variable grids must now have a fourth parameter indicating the default value. This was important for reading e.g. chilling units at the model startup. |
| WaSiM 7.6.0, 22.04.2007 | |
| + | dynamische Phänologie, es wird an einem definierten Tag (in der Steuerdatei pro Landnutzung festlegbar) der TSTART-Wert auf -1 gesetzt und die Forcingunits auf 0. Somit kann man auch längere Jahre simulieren. |
| WaSiM 7.5.1a, 10.04.2007 (inoffizielle Zwischenversion) | |
| + | dynamic phenology implementation started; should compile and run with usual control files dynamic phenology is not yet complete: |
| - TSTART and FORC grids must be created when required | |
| - Forcing units must be calculated | |
| - actual parameters for TSTART must be set when the conditions are meet | |
| + | implementation of dynamic phenology: reading the landuse table is implemented |
| todo: creating FORC and TSTART-grids and doing the forcing calculation as well as setting the actual parameter in the LUtable when beiing in a TStart...TStart+Delta time slice | |
| WaSiM 7.5.1, 06.03.2007 | |
| ! | in Precicorr war Array entf mit 1 Element zuwenig reserviert, gab Fehler bei nur 1 P-Station |
| 2006 |
| WaSiM 7.5.0, 18.11.2006 | |
| ==> offizieller Release | |
| + | kleinere Korrekturen für die Ausgabe der Nachrichten an den Bildschirm (wurden im Paralelbetrieb bei den sections zerhackt) |
| + | Modapuzr weiter parallelisiert (GW-Model, Routing etc. als sections definiert) |
| + | Unsatzon is now threadsafe (there were some errors when updating riverflow_table-content) |
| - | gwflow will probably never be faster using multiple threads than it is now, so the omp-directives are commented there (in gauss_seidel_iteration only) |
| + | OpenMP-Threading für Unsatzon arbeitet jetzt, jedoch sind in der Rechnung noch einige winzige Abweichungen im Resultat zu sehen --> irgendwo wird noch eine gemeinsame Variable überschrieben... |
| WaSiM 7.2.7, 08.11.2006 | |
| ==> offizieller Release | |
| ! | Korrekturen am Externalbalance-Grid sowie an SB1-Statistik |
| WaSiM 7.2.6, 08.11.2006 | |
| ! | Probleme mit externer Kopplung korrigiert (GW-Stand über Gelände verursachte Fehler) |
| + | Unsatzon: alle ZoneFZT_Q*-Zugriffe nach ausserhalb der Schleife verlegt, leichte Unterschiede in den Ergebnissen aufgrund Speicherung der QD und QI-Werte neu als mm /*timestep) und nicht mehr als m/s, Der Mittelwet wird sozusagen auf bereits umgerechneten Werten berechnet. Früher wurde Qd und Qi erst in m/s für das TG berechnet und dann * timestep genommen, um auf mm zu kommen. Ausserdem sind kleinere arithmetische Unterschiede zwischen GCC und VC-Ergebnissen feststellbar, was aber schon immer so war. |
| + | Unsatzon: threadsafe except for FZTzone_Q*, which still must be changed -> qd[0]..qd[z]calculation is not independent on the grid cell, so we either introduce a new stack or we move the qd-calculation outside the main parallel loop. |
| WaSiM 7.2.5, 30.10.2006 (Seitenzweig, interne Version) | |
| Unsatzonmodel: | |
| + | 9 Funktionen threadsave gemacht (es fehlen noch: new_concentrations; infiltration; wertzufuegen; zoneFZT_Q*-Zugriffe); alle in ::run aufgerufenen Funktionen (innerhalb Zeilenschleife) threadsafe gemacht |
| ! | kleinere Bugs bzgl. Irrigation-Bilanz (DELTAGW-Handling) korrigiert |
| ! | Initialisierung SB bei irrigation korrigiert |
| ! | im Schneemodell einen Bug gefixt, der durch Parallelisierung reingekommen ist (Schneeanteil lokal definiert und als Referenz-Parameter zurückgegeben aus Calculate_Melt-..) |
| WaSiM 7.2.4, 27.10.2006 (Seitenzweig, interne Version) | |
| ! | Korrigierte Kopplung mit Externem GW (Theta bei steigendem GW korrekt aufgefüllt) |
| ! | Werte in Storagedatei z.T. Null, weil Format %lf bei exponentialwerten von z.B. 1.5e-7 0 rausschreibt |
| ! | SB_2 Grid Inhalt jetzt korrekt als absoluter Wassergehalt der Bodensäule |
| + | OpenMP-Implementierung für Meteoint, ETP, Expokorr, Interzeption, Schneemodell und einige Teile Unsatzonmodell (dort 6 Funktionen von 13 bereits ok, die restlichen 7 müssen noch threadsave gemacht werden) |
| + | OpenMP in Meteointerpol und ExpoKorr drin |
| WaSiM 7.2.3, 21.03.2006 | |
| + | Ausgabe der Schneeverdunstung |
| WaSiM 7.2.2, 18.03.2006 | |
| ! | Korrigierte Version, nun bilanzrein |
| + | Detailausgaben ETR, EI, EVAR, EVAP |
| + | Interzeption für mehrere Vegetations-Layer implementiert |
| o | Anpassungen für Sleep-Funktion (neu: _sleep in VC++ 2005) |
| WaSiM 7.1, 11.03.2006 | |
| ! | 7.0 -> 7.1 Header in wasim.h geändert |
| WaSiM 7.0, 11.03.2006 | |
| + | detaillierte Ausgabe für ETR aus Bodenmodell |
| + | Berücksichtigung von Makroporen |
| + | separate Parametrisierung von Interzeptionsverdunstung und Evaporation |
| + | diverse zusätzliche Ausgaben |
| 2005 |
| WaSiM 6.8.0, 20.12.2005 | |
| ! | Bugfixes für Dicken[0] in unsatzonmodel |
| + | Multi-Layer-Lannutzungstabelle implementiert |
| + | Aufnehmen der Parameter für Makroporenfluss in die Steuerdatei und Lesen derselben |
| + | k-recession auch für neue Bodenparameterisierung ermöglicht |
| WaSiM 6.7.0, 18.11.2005 | |
| + | variable Bodenartentabelle implementiert (nicht nur intern, nun auch über Steuerdatei nutzbar) |
| WaSiM 6.6.2, 04.11.2005 (interne Version) | |
| + | Bodentabelle wird intern mit Detailparametern für alle Schichten verwendet |
| WaSiM 6.6.1, 22.10.2005 | |
| + | in LU-Table und Verdunstungsmodell wurde anaerobe Reduktion der Verdunstung frei parametrisierbar gemacht (TReduWet und LimitReduWet werden nun in der Steuerdatei gelesen und in unsatzonmodel::evapotranspiration berücksichtigt) |
| WaSiM 6.5.1, 20.09.2005 | |
| + | koppelbare Version |
| o | LU-Tabelle als gekapselte Klasse -> dies als Vorbereitung auf die anstehende Erweiterung, um die Module davon frei zu halten. |
| WaSiM 6.4, 01.07.2005 | |
| ==> offizielles Release | |
| ! | Ablageort Storagefile bei Readgrids = 1 in Unsatzonmodel (war neuer Bug in Version 6.1) |
| ! | Infiltration/Exfiltration zwischen Grundwasser und Fliess-Gewässer in Unsatzonmodel korrigiert (nur aufgetreten, wenn die Teileinzugsgebiete nicht fortlaufend nummeriert waren) |
| ! | k-Wert-Berechnung verbessert im Richardsmodell |
| + | interne van-Genuchten-Tabelle intern mit mehr Stützstellen |
| + | Dicke erste Bodenschicht neu 0.05 anstatt 0.1 m |
| + | mehr interne Buffergrössen für seltenere Plattenzugriffe |
| + | Entnahme Irrigationwasser aus Reservoirs (wenn vorhanden) |
| + | Neue Initialisierung GW-Level (Steuerfile-Anpassung!) |
| + | Schneemodell: variable Albedowerte, Modellierung Schneedeckenalter, Anteil Schnee/Regen, variabler Strahlungsanteil für Ice_firn_model --> In Steuerdatei neue Parameter beachten |
| + | Unsatzonmodel: neue Ausgabegrössen als Grids (2 Bodenwassergehalte sowie Wurzeltiefe) sowie Ausgabe Wurzeltiefe als Statistikfile, Initialisierung GW-Level inSteuerdatei festglegt |
| o | zulässigen Iterationsfehler verkleinert in Richards-Ansatz (von 1e-3 zu 1e-4) |
WaSiM-News
01.05.2012
New Model Version
WaSiM 8.10.0 -> Mehr »