Bug fixes and changes to SSIIM 2 for Windows

27. February 2000:
Added to the map graphics: Number of cells in the vertical direction and roughness/depth ratio.

29. February 2000:
Fixed bug in the discharge editor for removal of inflow surfaces. Added F2 option for rapid removal of inflow surfaces.

9. March 2000:
Fixed a bug in reading the unstruc file for large values of the grid coordinates. Added short-cut keys for viewing text, map and the grid editor.

16. March 2000:
Fixed several bugs in the sediment computation module, for the upstream boundary cells.

21. March 2000:
Fixed bugs in the specification of the water discharge and the time-dependent input of water and sediment discharge.

22. March 2000:
Fixed bugs in the specification of the water discharge, in the dialog box.
Added zero-gradient boundary condition for outflow - only for steady water discharge.

27. March 2000:
Fixed two bugs in the boundary conditions for the sediment calculation.
Fixed a bug in reading the water levels on the D data set in the timei file.

3. April 2000:
Fixed bugs in the boundary conditions for the sediment calculation.

26. April 2000:
Fixed bugs for sediment continuity control algorithm, for inflow concentrations given in the timei file.
Fixed bugs in GUI: number of cells in vertical direction+ bed roughness.
Added user-defined thickness of top sediment layer: On F 106 data set (in meters) Default: Maximum sediment size.
Added sediment continuity control for algorithm preventing blocking of downstream outlet.
Corrected algorithm for sediment concentration limits for one sediment size
Added F 105 data set: integer for iterations between updates of the water surface for transient computations. Default: 10.

6. September 2000:
Fixed bug in sediment transport module - boundary condition.

26. October 2000:
Fixed bug for computation of turbulent kinetic energy, when the flow is supercritical.

11. December 2000:
Fixed a bug for reading the F 106 data set
Fixed a bug for computing the total sediment concentration when only having one sediment size
Fixed a bug in the grid editor, when giving the number of grid lines for a block
Added the Froude number as a parameter for the Map graphics

19. February 2001:
Fixed a bug for hexahedral cells with one or more surface of zero area.
Fixed a bug for dynamic generation of grids for more than one block.
Fixed a bug for reading the RESULT file when doing transient computations.
Added beddll.dll DLL file.
Added new algorithms for sediment computations with dynamic grid and a hydrostatic/2D water flow computation (not tested).
Added the possibility of fixing two points instead of one when updating the water surface.
Added an algorithm for wetting/drying, when starting with a partially dry geometry.

13. March 2001:
Implemented an improved algorithm for initializing the variables in newly wetted cells.
Fixed two bugs in the contour map graphics: for roughness/depth and Froude number.

20. March 2001:
Made a change to the Rhie/Chow interpolation close to the wall. The new algorithm is used if F 113 1 is present in the control file. If not, the old algorithm is used.

10. April 2001:
Made a file called 'habitat', which is written every time the result file is written. This file contains the availability function for making a habitat preference function from the HSI index.

11. April 2001:
Made a new pressure-gradient algorithm for cases where there are very shallow areas inside the grid. This is invoked by using F 114 9 or 10 (two versions).

14. April 2001:
Made sand slide algorithm for multiple grain sizes. Not tested.

16. April 2001:
Fixed a bug for the contour map graphics for cells that change form.
Started implementing the vegdll.dll file, for computing effects of vegetation on the water flow.
Moved allocation of arrays for contour map generation, fixing a bug where the program crashes when it shows the contour map on PC's with little RAM.
Implemented depth-averaged velocity and vegetation in contour map graphics.

17. April 2001:
Changed the method to compute the maximum depth for generation of unstructured grid, fixing a bug for reading result and bedres files at the same time.

18. July 2001:
Fixed a bug for calling a function in the vegdll.dll.

25. August 2001:
Fixed a bug for reading time series using only the D data in the timei file.
Fixed a bug for reading the bedrough file.

26. January 2002:
Fixed a bug that makes the program crash if the graphics are viewed right after startup.
Added the tsc2dll.dll file.
Fixed some bugs for computing sediment transport without bed changes (F 68 2). Added alternative pressure gradient algorithms for triangular cells
Added locally denser grid with quadrilateral cells (not tested)
Fixed a bug in the sediment transport algorithm, for computing the critical shear stress of a particle.

25. March 2002:
Fixed bugs in the contour map graphics.

30. May 2002:
Fixed an error in the setting of a switch for the compiler, for the multi-thread library. The error prevented the interres file to be written properly. Unfortunately, the changed switch also caused the program to run about 8 % slower.

4. November 2002:
Improved data structure written to the bedres file.

6. November 2002:
Added the following data sets in the 'control' file:
F 138: Mininum value of the turbulent kinetic energy (default 1.0e-20)
F 139: Minimum value of u+ in the wall laws for the velocity and turbulent kinetic energy (default 1.0).
F 140: Flag for correction of epsilon values in the wall laws, if u+ goes below the value on the F 139 data set. Used if the integer is 1 or larger. Default 0.
F 141: Maximium value of epsilon. Default 100.0.
Improved data structure written to the bedres file. Also note that using the F 99 data set may make it impossible to read both the result and the bedres file for restarting the computations.

12. November 2002:
Fixed a bug in the computation of d90 in the sediment transport algorithm.

14. December 2002:
Added error messages to be written to the boogie file and the window if the program tries to allocate more memory than what is available on the computer.

19. February 2003:
Added scales in the Map graphics and the Grid editor.

5. March 2003:
Fixed a bug in the computation of the bed roughness based on the sediment grain size distribution.

7. April 2003:
Fixed one bug for the update of the water surface in the sediment computations. Fixed several bugs for the F 143 option.

8. April 2003:
Fixed a bug for the update of the water surface. Added the possibility to move/scale the graphics in smaller increments, by using ALT+arrows or ALT+PgUp/PgDown.

9. May 2003:
Added interpolation in the scale dialog box for the contour map graphics.

20. May 2003:
Fixed a bug in the function calling the bed form height computation in the BEDDLL dll.

17. June 2003:
Implemented F 143 2, using an average pressure used for the water surface update, instead of the linear interpolation algorithm, which may have a bug.

23. July 2003:
Added scales for the profile graphics.

30. July 2003:
Fixed bugs in the F 143 option.
Added F 148 option, giving the maximum local slope of the water surface.

31. July 2003:
Added algorithm to use zero-gradient outflow boundary conditions also for the time-dependent computations. Previously, this only worked for the steady computation. The zero-gradient boundary condition is set in the dialog box of the discharge editor, in the "update" check box.

5. August 2003:
Seem to have fixed bugs in the algorithm for extrapolating values to cells that are wetted. There is still a bug in the F 143 data set option, giving continuity errors.

7. August 2003:
Fixed the bug in the F 143 data set option, related to update of water surface for multiple blocks. Implemented possibilities to write Tecplot files, but this is not tested yet.

8. August 2003:
Fixed a bug in writing the 2d Tecplot file. Changed the F 48 7 option, so that a 2d Techplot compatible file is written instead of the old Tecplot file. This is used in combination with the P 10 data set, and multiple files are written while the program is computing.

11. August 2003:
Added modules writing 2D and 3D Tecplot files in combination with the P 10 data set. If the F 48 data set is 7 or 8, the 2D Tecplot file will be written. If the F 48 data set is 8 or 9, the 3D (FE) Tecplot file will be written.

13. August 2003:
Fixed a bug in writing the Tecplot file for the F 48 8 option.

12. September 2003:
Fixed a bug in the inflow/outflow boundary condition for the sediment concentration. The bug was in the program from the 5th August -03, and caused the sediment flux at the boundary to increase by a factor equal to the number of grid cells in the vertical direction.

2. October 2003:
Added code for writing more than 999 files in time-dependent computations. The limit is now 999 999 files.

4. December 2003:
Removed two bugs: One was the F 141 data set, which gave non-convergence for epsilon for high velocities. The other bug was in the map graphics for display of the bed roughness for non-sediment transport computations.

7. February 2004:
Added warnings written to the boogie file if memory allocation errors occurred.

11. February 2004:
Added two new algorithms for water surface changes:
F 36 3: the initial water surface is moved up/down equally in all cells according to downstream changes in water surface.
F 36 4: the water surface elevation is computed solving a partial differential equation for the local slope as a function of the pressure gradients.
Note that the algorithms are not tested.

12. February 2004:
Changed the slope limiter in the F 36 4 algorithm from 1:100 to 1:30.
Included the relaxation factor from the G 6 data set in the F 36 4 algorithm, instead of the hard-coded value of 0.1.

13. February 2004:
Added the F 41 and F 42 data set to the control file. The F 41 data set gives the sediment thickness in meters. The F 42 data set gives the level of non-erodable material. These data sets can be used instead of the koomin file.

1. March 2004:
Implemented an improved algorithm for computing changes in the grain size distribution at the bed. Also improved continuity defect controls.

24. March 2004:
Added residuals to the Map graphics.
Changed the Tecplot 2D file, so that the variables are user-defined in the G 24 data set.
Started making a kind of Multigrid algorithm, which was not finished.

25. March 2004:
Fixed a bug in the number of iterations for transient sediment compuations.

26. March 2004:
Fixed bugs in writing the Tecplot 2D file.
Changed the structure of the Tecplot 3D file to use user-defined variables in the G 24 data set.

2. April 2004:
Fixed a bug in computation of residuals for the velocity when the SOU scheme was used.
Added the option F 64 13. This is the same as F 64 11, except that the data on the F 94 data set is used to give different cell generation criteria for new and old cells. The first float on the F 94 data set is the minimum depth for creation of new cells, while the second is the minimum depth for old cells. This is an attempt to improve the stability of the wetting/drying algorithm.

5. April 2004:
Fixed a bug in the F 64 data set, for the 13 option.

6. April 2004:
Fixed a bug in writing the Tecplot files.

19. April 2004:
Fixed a bug in writing the Tecplot 2D file, related to the water depth/water surface/bed level/bed movements/sediment thickness.
Did some changes in the algorithms for the F 36 4 option. Note that these are not finished yet.

20. April 2004:
Added an error check for the F 2 data set in the control file.

21. April 2004:
Fixed a bug in writing the Tecplot 3D file, related to the water depth/water surface/bed level/bed movements/sediment thickness.

13. May 2004:
Added the sfdifdll DLL file. This file is used in a new algorithm for computing the location of the free surface. The DLL is used if F 116 1 is present in the control file. The algorithm is also coded in the main program, and invoked if the option F 36 4 is used.
Added a part in the F 36 4 algorithm where the downstream water level is horizontal, instead of only being moved at one point. Then the third integer on the G 6 data set is used to specify the outflow group from the discharge editor.

19. May 2004:
Fixed some bugs in writing the bed level, water level and sediment thickness to the Tecplot files.

22. May 2004:
Fixed a bug in the SFDIFDLL.DLL file.
Fixed a bug for moving the water surface over the whole downstream cross-section for F 116 > 0.
Added the option F 116 2, giving better stability and lower convergence speed for the new free surface algorithm.

24. May 2004:
Added the functionality of the F 109 option to SSIIM 2.

29. May 2004:
Increased the maximum value for the F 21 data set from 2 to 20.

2. June 2004:
Added a new criteria for wetting/drying : F 64 111.

11. June 2004:
Added a surface smoothing option: F 154, reading an integer and a float. The smoothing is invoked if the integer is 1. The float is a relaxation of the smoothing, with 0.0 being no smoothing, and 1.0 full smoothing.

19. June 2004:
Fixed a bug in the cell shape change algorithm, which sometimes caused a cell to collapse to zero volume.

21. June 2004:
Fixed another bug in the cell shape change algorithm.
Changed the F 64 13 option, so that the smallest number on the F 94 data set is first. Also combined the wetting/drying modification with the one-cell criteria of F 64 11.

28. June 2004:
Fixed two bugs in the grid generation algorithm. One was related to the number of connection points between multiple blocks. The other was related to the limits of local scour and the koomin file.
Added the option of pausing the program before/after some of the algorithms in the program. This is done from the menu. This option does not work on the speed-optimized version of the program on the web.

30. July 2004:
Changed the algorithm for computing the location of the inflow/outflow as the grid changes.
Added a new algorithm for the SIMPLE corrections: F 157 2. The old algorithm is now invoked by using F 157 1 instea of F 7 P.

2. August 2004:
Improved the F 157 2 algorithm.
Fixed a bug related to the loops and brackets in the Simple Correction algorithms. Compiler bug.

10. September 2004:
Modified the algorithms for generation of the unstructured grid and the changes in cell shapes.
16. September 2004:
Made a new version of the Rhie and Chow algorithm, to improve velocities in triangular cells. This is invoked by using F 113 7 in the control file.

21. September 2004:
Modified the Rhie and Chow algorithms for F 113 7. Also added a pressure-gradient limiter. The new version without the limiter is invoked by using F 113 6. The purpose of the algorithms is to reduce the problem of unphysical directions of the velocity vectors in triangular cells.

15. October 2004:
Fixed some smaller bugs in the F 37 1 option. Added the F 37 2 option, where the bed boundary condition for the sediments are given as a pick-up rate flux instead of a concentration.

26. October 2004:
Fixed some bugs for computation of the thickness of the upper sediment layer, given a finite thickness of the movable bed.

9. November 2004:
Fixed more bugs for the computation of the thickness of the active sediment layer.

29. November 2004:
Fixed a bug for the F 64 13 option.
Fixed a bug for the F 37 4 option.

3. December 2004:
Changed a depth criteria for the F 64 13 algorithm.

8. December 2004:
Changed the depth criteria for the F 64 13 algorithm. The value can be given on the F 162 data set. The default is zero.

11. January 2004:
Fixed a bug for the bed changes at the edges of the boundary during wetting/drying.
Fixed a problem for determining the number of grid cells in the vertical direction, close to the maximum values.
Fixed a problem for distribution of bed movement from the cell center to the cell corners if there is limited sediment available or limited space for deposition.
Added the F 163 data set: Algorithm for limited space/sediment in bed movement algorithm:
F 163 0 (default) Algorithm only used for limited amount of sediments
F 163 1 Algorithm also used for limited space for deposits
F 163 2 Algorithm not used
Added the F 164 data set: solver options:
F 164 0 (default) solver in streamwise lines
F 164 1 solver in both streamwise and lateral lines
F 164 2 solver only along lateral lines
F 164 3 solver in streamwise lines, but starting downstream.
Added the F 165 data set: numbering the cells from the left bank instead of the right bank if the integer on the data set is unity.

13. January 2005
Fixed an allocation error for an array in the pressure-regression algorithm for the free water surface.
Added the F 166 data set: If F 166 1 is included, the grid will be regenerated also after each water surface uptdate, even though bed changes are also computed and the grid is updated after each bed change.

14. January 2005:
Added extra debug print-out for free surface algorithm if F 1 E and F 36 3 options are used.

15. January 2005:
Fixed a bug in the SOU algorithm for the water flow computation.

19. January 2005:
Fixed a bug for the F 36 3 option.
Included an error message for when the reference cell dries up when computing the free water surface

21. January 2005:
Fixed a bug in the F 37 2 algorithms

26. January 2005:
Fixed a bug in the setting of the inflow boundary conditions for the turbulence

18. February 2005:
Added new variables to the BEDDLL file: indexes + d50/d_mean.

23. February 2005:
Fixed a bug in setting the epsilon values at the inflow boundary.

1. March 2005:
Fixed a bug for the F 106 data set, which was not used for computations with only one sediment size.

9. March 2005:
Found a bug in the wall laws for k, for the side walls. Not fixed yet.

13. April 2005:
Fixed the bug in the side wall for k.
Added the F 169 data set, for giving in hiding/exposure parameters to the beddll.dll file.
Added a new multi-grid solver for the pressure-correction equation. It is invoked by using the F 170 data set, where the parameter is an integer. The integer gives the number of levels in the multi-grid algorithm. The algorithm is not completely tested yet.

27. May 2005:
Added the turbdll.dll file.
Added the possibility to use different eddy-viscosity in horizontal and vertical direction.
Added modification of the eddy-viscosity in the vertical direction from the low-Reynolds number turbulence models by Rodi and Mansour (F 24 8) and Nagano and Tagawa (F 24 9).

31. May 2005:
Added the F 127 data set. Default values: F 172 0 2.0 2.33. If the first integer is 1, then an algorithm will be used to increase the shear stress as a function of the bed form/roughness height in relation to the water depth. The first floating point is the maximum value for the average depth divided by the depth above the bed form/roughness. The second float is the factor the ratio is raised to for giving the increase in the shear stress. Note that this algorithm is not tested at all yet.

1. June 2005:
Fixed a bug in the F 24 8 and 9 algorithms.
Added an option 3 in the F 172 data set.

2. June 2005:
Fixed bugs in the F 24 8 and 9 algorithm.

6. June 2005:
Fixed a bug in the F 24 8 and 9 algorithm.

25. April 2006:
Fixed several bugs in the algorithm interpolating results between changing grids.
Added a multi-grid algorithm for the pressure-correction equation: Add F 168 8 and 10 instead of 0 for the pressure-correction equation on the K 5 data set.
Added the F 182 data set: reduction of critical shear stress on sloping beds. An integer is read. If it is 1, then this will be the same as F 7 B. If the integer is 2, then the algorithm will only be applied to cells that has dry neighbours. Which will typically be the side banks of a channel.
Modified the F 36 4 algorithm.
Added a new free surface algorithm: F 36 8. The algorithm is a bit similar to F 36 4, but it uses a 9-point computation molecule instead of the usual 5-point.
Added a low-depth Reynolds number model: F 24 10. If the water depth is very shallow, the eddy-viscosity is computed by a zero-equation model instead of the k-epsilon model.
Added the option P 11 4. Then the total number of cells, surfaces and points in the grid is written to the boogie file for each iteration.
Added specification of the water surface in the Grid Editor. Three points are given graphically, and the program interpolates the surface between these.

25. September 2006:
Added outblocking of region from the G 13 data set. The format is the same as SSIIM 1, but the first and the two last integers are not used. They still have to be given, though.
Added F 191 data set, specifying the accuracy when gluing blocks together. Default 0.00001 meters. The number can be set higher for larger-scale geometries.
Added a new algorithm for interpolation of values between changing grids. It is invoked on the second parameter on the F 147 data set.
Modified some parts of the algorithm for interpolation of values between changing grids, when the old or the new 2D cell array has only one cell.
Added algorithm for increasing water depth in shallow areas: Last integer on the F 159 data set has to be set to 2.
Collected array allocation in fewer functions.

27. September 2006:
Fixed a bug in the generation of nested grids.
Added the nestdll.dll file.
Fixed a bug in the F 36 3 option when not using the timei file.

29. September 2006:
Added option of viewing nested blocks or non-nested blocks.
Added F 192 data set, giving the ratio of number of vertical cells in the nested grid to the number of grid cells in the coarse grid.

30. September 2006:
Fixed a bug related to activating nested grid for computations.
Fixed an error message for allocation of nested surfaces.

2. October 2006:
Fixed a bug in the metric for the nested block.

3. October 2006:
Modified the interpolation to the nested block.

5. October 2006:
Fixed bugs in the G 25 data set algorithms, and in the F 188 data set.

31. October 2006:
Fixed a bug in the F 90 4 algorithm, for the computation of the bed form height and the critical bed shear stress. The reduced shear stress was used as input for the functions. Now the total shear stress is used.
Added calls to the new beddll

22. January 2007:
Tried to fix a bug in the writing of 3D Tecplot files when having multiple blocks.

23. January 2007:
Fixed a bug in the writing of 3D Tecplot files.
Fixed a bug causing program crashes for non-sediment computations.

30. January 2007:
Coded a new sand slide algorithm. Invoked if the first integer is 200 on the F 56 data set. Not tested.
Added a new wall law: F 15 8, a three-layer wall law for smooth and rough boundaries. Not tested.

16. February 2007:
Added the slide2dll file. The dll is invoked by giving a non-zero integer on the F 198 data set.

19. February 2007:
Changed the limiter for the F 15 9 option from -0.06 to -0.1.

20. February 2007:
Extended the numbering system for multiple Tecplot files.

21. February 2007:
Added legend for contour map and vectors.

13. March 2007:
Added density and turbulent kinetic energy as parameters to the slide2dll dll.

26. March 2007:
Added the ParaView.vtk file to be written from the user interface or the main program (F 48 10).
Added Rodi's wall laws for k for normal rough wall laws (F 15 11).

27. March 2007:
Fixed a bug in the ParaView.vtk file, so that also velocity vectors can be shown (Filter-Glyph)

3. April 2007:
Added a fix to the zero-gradient outflow boundary condition.
Added logarithmic velocity profile for the inflow. This is invoked by adding F 202 1 to the control file.
10. April 2007:
Fixed a bug in the upstream and right profile for the update of the water surface, option F 36 2.
Removed the F 140 data set, as it gave incorrect bed shear stress.
13. April 2007:
Changed the F 202 data set, for the logarithmic inflow, so that two parameters are read. First an integer and then a float. The float is the roughness of the velocity profile.

2. May 2007:
Fixed a bug in the algorithm for the inflow/outflow boundary when the grid changes.
Fixed several bugs in the F 36 8 and 9 algorithms.

3. May 2007:
Fixed a bug in the F 56 200 algorithm.

14. May 2007:
Fixed a bug in the F 56 algorithms.

28. August 2007:
Added the F 209 option, where a scaling depth is given for the generation of the unstructured grid. Typically, the maximum depth is given.
Changed the algorithm for the surface interpolation from given points in the grid editor, so that multiple consecutive triangles are used. Not tested nor finished yet.

30. August 2007:
Fixed a bug where the program crashed when using the F 112 data set and not using the F 36 or the F 37 data sets.

4. September 2007:
Fixed a bug for the F 209 data set.

5. September 2007:
Fixed a bug in the graphics when using the F 112 data set.

11. September 2007:
Increased the maximum number of grid cells in the vertical direction to 500.
Added some k-omega extensions (F 24 16), but not tested these.

12. September 2007:
Fixed and allocation bug for using the F 115 data set in a steady-state computation.
Increased the default vertical number of levels to 99 for the discharge editor.
Fixed some bugs in the interpolation of vegetation parameters.
Added level in vegetation map graphics.

13. September 2007:
Fixed more bugs in the vegetation interpolations.

15. September 2007:
Fixed more bugs in the vegetation interpolations, and added some error messages for the vegdata file.

18. September 2007:
Fixed a bug in combination of running free water surface computations without sediments and using the F 168 data set.

20. September 2007:
Fixed a problem in using timei parameters after end of specified time.

24. September 2007:
Fixed a bug in reading the bedrough file.

25. September 2007:
Fixed a bug in reading the M data set in the control file.

29. September 2007:
Added the source2dll file.

1. October 2007:
Made a parallel version. Use the F 206 data set with an integer to specify how many processors you want to use. Default: F 206 0 (sequential version).

10. October 2007:
Fixed an error in reading vegdata files without T data sets.
Fixed an error reading the "fracres" file.

25. October 2007:
Added circle segmenet to the grid editor for the parallel version of SSIIM 2.

31. October 2007:
Fixed some bugs in the parallel algorithms of SSIIM 2, including the multi-grid algorithm and the pressure source.

6. December 2007:
Added a new interpolation option for the interres file: F 48 5. This gives only the coordinates and the water surface velocity in the interres file.

12. December 2007:
Fixed an interpolation bug in the wetting/drying algorithms, which caused the program to crash.

17. January 2008:
Changed the scaling of the graphics from the geodata file.

19. January 2008:
Modified the viewing of the geodata points in the Grid Editor, so that the user can vary the number of points.

24. January 2008:
Modified the k-omega model (F 24 16), and updated the RNG model (F 24 1).

15. February 2009:
Fixed a bug in the inflow/outflow boundary condition when using dynamically changing grids.

6. March 2008:
Improved the interpolation algorithm between grid changes to reduce the sediment continuity defect.
Made a new parallel algorithm for the SOU scheme.

7. March 2008:
Fixed a bug in the new SOU scheme for the parallel algorithm.

3. April 2008:
Increased the number of surfaces for an inflow/outflow group from 2000 to 20 000.

11. April 2008:
Fixed a bug for steady state lowering of the water surface.

23. April 2008:
Added the I option to the F 2 data set. This runs through the bed sediment routine one time, initializing the roughness.

2. May 2008:
Improved the speed to write unstructured Tecplot and Paraview files.

5. May 2008:
Fixed a bug where the P 10 parameter in the control file didn't work when a time-dependent solution converged in a time step.
Introduced the F 222 data set, reading an integer. If it is zero, an algorithm used to lower the outflow cross-section is not used. The default value of the integer is 1.

7. May 2008:
Added two new F 48 options: 12 and 14. 14 produces an interres file with the depht*velocity value for each point. 12 produces a 2D paraview file of the cross-section, with the avereage velocity as the variable. The files are written whenever the "result" file is written.

13. May 2008:
Fixed a bug in the computation of the residuals for the OpenMP algorithms.

31. May 2008:
Fixed some bugs in the OpenMP algorithms of the program.

2. June 2008:
Fixed a bug in the non-parallel algorithm of the program, introduced on the 31st of May.

3. May 2008:
Fixed a bug in the non-parallel algorihtm of the program. Implemented F 164 27, which uses a more consistent but slower parallel solver.

29. July 2008:
Fixed two bugs in the 64 bits version of SSIIM 2: One about printing the correct amount of allocated memory to the boogie file and the other was an error in the Grid Editor, when selecting points with large index numbers.

2. October 2008:
Improved the F 36 8 free surface algorithm.
Improved interpolation of variables between grid changes.
Improved algorithm to write 3D Tecplot files when using multiple blocks.
Added grid index numbers when choosing "legend" for the grid editor and the discharge editor.

7. October 2008:
Fixed a bug for time-dependent inflow from multiple sediment sources (timei file).
Fixed a bug for reading the "fracres" file in initially dry cells.
Fixed a bug for viewing bed changes with the F 68 2 option.
Added printout to user interface and boogie file for sediment computations using F 68 2.
Fixed a bug for the P 10 data set used together with the F 68 2 option.
Added the "Nagel" file for post-processing: A point cloud of values. Time-depenent writing of the file using F 48 11 together with P 10 in the control file.

20. October 2008:
Fixed a bug in the interpolation of variables between changing grids, causing instabilities.
Note: the pressure-limiter on the F 136 data set may give instabilities.
Added the Nagel file, giving x,y,z coordinates, velocities, k,epsilon, eddy-viscosity, pressure and sediment concentration (total) for all grid cells. The file is used for post-processing.

17. November 2008:
Added option to add 5 geodata points between two existing geodata points. Choose this option in the menu and click on the two existing geodata points.
Added new algorithm for F 78 data set when multiple sediment sizes are used.

18. November 2008:
Added option to view index and level of geodata points.

1. December 2008:
Implemented the F 78 data set for mulitple sediment fractions.

9. December 2008:
Fixed a bug in the bed levels when writing the 2d Tecplot file.

10. December 2008:
Did some more fixes related to the 2d Tecplot file.
Increased the maximum number of grid cells in the vertical direction from 99 to 999 for the 64 bits version.

12. December 2008:
Fixed a bug with the sediment fractions when erosion of the whole active layer occur.

7. January 2009:
Fixed a bug in the bed interpolation algorithm.
Improved algorithms for grid generation: F 64 3 and F 64 4.
Added P 16 data set. Four intergers are read, defining an area in the 2D structured grid seen in plan view. The contour map graphics will only show this area. The purpose of the data set is to reduce the time used for graphics updates with large grids.

15. January 2009:
Fixed a bug in the graphics scaling for 2D depth-averaged grids.
Added new algorithms for free water surface: F 36 14,15,16 and two smoothing algorithms: F 178 3.
Fixed a bug when applying the new free surface algorithms with the multi-grid solver.
Added time in the profile graphics, and also fixed the velocity scale. This is shown when "Legend" is chosen in the menu.
Fixed two bugs when writing pressure and some other variables to the Paraview files.

3. March 2009:
Fixed a bug in the free surface algorithm for the F 36 2 option. The bug would only be noticed when the pressure deviates from the hydrostatic distribution in the vertical direction.

6. March 2009:
Fixed a bug when using small values on the F 209 data set in combniation with the F 64 options 8, 11 or 13.
Improved the accuracy of the free surface solver for the F 36 options 4 and 7.

25. March 2009:
Fixed a bug in the interpolation of values between changing grids.

3. April 2009:
Improved the use of the bedrough file.

11. May 2009:
Simplified the conversion of SSIIM 1 files to SSIIM 2 files, by enabeling the koosurf file to be read into the grid editor.
Added option for stabilizing velocities in triangular cells: F 235 10.
Added the option of specifying the discharges in the control file instead of the unstruc file. This is done on the F 237 data set. The data set takes one integer and one float. The integer is the discharge group, and the float is the discharge. For example: F 237 1 0.02 means a discharge of 0.02 m3/s for discharge group 1. The discarge groups are numbered from 1 as the first group.

15. May 2009:
Fixed a bug in the F 113 7 option.

15. June 2009:
Inreased the maximum number of G 13 data sets from 9 to 499.

31. August 2009:
Improved the accuracy of the computation of the average area across the center of a cell. Most improvements for cells on steep bed or surface slopes with a high ratio of hight/(width or length). The improvement leads to better stability in these cases.

7. September 2009:
Fixed a bug in the F 56 200 algorithm.
Fixed bugs in the nested grid algorithms.

8. September 2009:
Fixed a bug for the parallel algorithms when not specifying the F 206 data set.

9. October 2009:
Fixed a bug in the algorithm for rapidly moving water surfaces (F 36 1 and 15), that was introduced in June 2009.
Added more variables to the graphics: secondary current angle, bed slope and potential sediment concentration.

27. October 2009:
Added secondary current angle (G 24 C) and bed slope (G 24 B) to the Tecplot and ParaView files that are written from the main program (Menu: Write result file").

25. November 2009:
Made changes in the sfdifdll.dll file and the free surface algorithms (F 36 4/7/8/9/10).

4. December 2009:
Added extra print-out in case of errors in the geometry.
Modified the F 36 7-10 algorithm.
Added F 233 6 option, using depth-averaged pressures for the F 36 2 and 7-10 algorithms.
Modified the F 36 2 option, reducing the water elevation changes.

20. January 2010:
Modified the zero-gradient downstream boundry condition.
Modified the F 36 7 option.
Added the F 187 option, with an integer above 1, as a limiter for the free water surface at the downstream boundary.

25. January 2010:
Modified the sediment inflow from the M data set, so that it can also be used together with a timei file.

26. January 2010:
Fixed some bugs in the F 187 option for the free surface.
Added a test for the first integer on the S data set in the control file.

27. January 2010:
Modified the sediment inflow from the M data set, for multiple fractions, so that it can also be used together with a timei file.
Added a test to check if the time-dependent sediment computation is started without specifying a time step.

6. February 2010:
Fixed a bug in the computation of the changes in the thickness of the sediment layers if only one size of sediments is used and the thickness of the sediment layer is small.

8. February 2010:
Improved the STL file written for the bed.

11. February 2010:
Improved the STL file written for the bed.

15. February 2010:
Fixed a bug for using the "fracres" file with F 37 2.
Fixed a bug for generating the ParaView file with the variables depth, sediment thickness, bed elevation or water level.
Increased the number of decimals in the number for the water volume, written to the boogie file.
Made the F 249 1 data set, allowing for negative sediment thicknesses.
Made the F 250 data set. This reads an integer pointing to the inflowing discharge group. The effect is the same as the third integer on the G 6 data set. Works only for F 36 7.

8. March 2010:
Added the option 2 on the F 202 data set. This will give an inflow velocity distribution with lower velocities and shear stresses in the shallow areas.

17. March 2010:
Fixed a bug in the water surface computation, for F 36 2, when using two reference points on the G 6 data set.

22. March 2010:
Modified the reading of the F 2 and F 7 data sets in the control file, so that the program stops reading the data set when a non-printable charachter is read. Note that this also applies to tabs, so there should not be tabs between the data set identifying integer and the letters coming afterwards.

24. March 2010:
Fixed a bug for writing the Tecplot2D file from the main user interface when the variable is water depth, sediment thickness, bed level changes, water level and bed level.

25. March 2010:
Fixed a bug in the algorithm invoked by the third integer on the F 246 data set, for SSIIM versions that are compiled without DLL's. The bug is only in the 64 bits Windows version and the Unix versions.

26. March 2010:
Added a new free surface algorithm: F 36 17. It is a bit similar to the F 36 15 algorithm.
Added logarithmic inflow with time-dependent water discharge. Note that this is not tested yet.
Modified the structure of the boundary conditions for the water discharge.

30. March 2010:
Rewrote the algorithm for inflow boundary conditions for the sediments, fixing a bug that were present if no sediments were specified in the timei file or on the M data sets in the control file.
Added a new solver algorithm: F 164 29.

7. April 2010:
Modified the F 78 1 option for 1 sediment size, so that it should give less continuity defects. Added a new explicit F 78 function, invoked with the integer 3, for example F 78 3 0.6. The new algorithm should be less accurate than the old function, but give better sediment continuity. The F 78 3 option should work for both one and multiple sediment sizes. Note that the F 78 1 function only works for one sediment size. The F 78 3 option is not tested in cobmination with the F 37 1 option.
Fixed a bug in the inflow boundary condition for the water flow computation.
Rewrote the outflow boundary condition function for the sediments. The new function fixes some continuity problems when using the M data set in the control file for specifying the sediment inflow in combination with the F 37 1 option. The new function includes an algorithm for zero gradient boundary, which was not in the old function. This new part is not tested, though.

9. April 2010:
Fixed three bugs for the computation of the sediment inflow. Two bugs were in the F 37 1 option, and one was for the logarithmic velocity distribution (F 202).

16. April 2010:
¨ Updated the F 36 17 option, including an addition of the F 254 data set, giving a weighting factor between fluxes and velocities in the algorithm.

17. April 2010:
Added a new solver algorithm: F 164 30.
Added a new bed change algorithm: F 255 1.

19. April 2010:
Fixed a bug in the F 164 30 option.

20. April 2010:
Added a new data set: F 256, giving an iteger for iterations for erosion in the bed changes. Default 100.

23. April 2010:
Added two new options to the F 222 data set. Option 2 forces the outflow cross-section to remain at the same bed elevation. Option 3 forces the four corners in the cells where the outflow is to remain at the same bed elevation.

26. April 2010:
Fixed a bug in reading the F 256 data set. Changed the print-out to the boogie file so that it will only print OpenMP times if the F 1 data set is P. Also, the program will only print out warning messages for erosion under movable bed if the F 1 data set is S.

27. April 2010:
Fixed two bugs for reading the F 255 and F 256 data set on the IBM AIX version of SSIIM 2.
Added a new solver: F 164 31. This only works for the 1st order upwind scheme.

28. April 2010:
Added 2nd order upwind solver for F 164 31.

29. April 2010:
Fixed a small bug in the upwind solver for F 164 31.
Added a new solver: F 164 32.

12. May 2010:
Added a new option for the third integer on the F 246 data set: -3. This will reduce the possibilities of grid splitting.
Added writing the "kookant" file when writing the unstruc file. This will produce the coordinates of the edges of the 2D domain for the first block. The file can be renamed to koordina and used when reading old results, avoiding strange lines around the boundary.
Added an algorithm raising water level at the inflow cross-section, preventing it from drying up. The algorithm is invoked by F 187 -1.
Fixed the G 50 option for the solver for the pressure-correction equation (9th parameter).
Added a new F 235 option: 14. This will increase the depth of triangular cells to keep them more stable.
Added print-out to the boogie file of the 2D indexes of the reference cell(s) on the G 6 data set.
Added the 2D indexes of the inflow/outflow areas to the map graphics when the cell numbers are shown.

20. May 2010:
Changed some of the structure in the function for the bed changes.
Changed the test for errors reading the F 200 data set, so that an integer above zero must first be read.
Added more information at the end of the boogie file regarding problems related to no surface updates taking place.
Added error message in the user interface window when a reference cell dries up.
Fixed a bug in the F 187 -1 algorithm.
Added a test for the correctness of the two first integers on the G 62 data set. The test also works for the G 6 data set, if the third integer is non-zero.

22. May 2010:
Added a F 84 4 function. The bed concentration is then a function of a constant times the dimensionless shear stress, T* in van Rijn's formula. The constant is the first parameter on the F 83 data set (default 0.053). The option is only meant to be used for debugging purposes.

25. May 2010:
Added an advice in the boogie file for errors in G 6/G 62 data sets in the control file.

27. May 2010:
Added a new algorithm for the grid generation: 2nd integer on the F 159 data set set to 11.

28. May 2010:
Added a new algorithm for the grid generation: 2nd integer on the F 159 data set set to 12. Fixed a bug in the F 164 31 algorithm.

31. May 2010:
Moved the F 256 integer to the F 242 data set. Then the same parameter is used on the same data set for both SSIIM 1 and SSIIM 2. Note that the default value is 1000 for SSIIM 1 and 100 for SSIIM 2.

1. June 2010:
Changed the solver structure, to make it more consistent for large F 164 values. Changed the bugtest for correctness of coordinates in the koomin file, solving a problem with large x and y values.

2. June 2010:
Fixed a bug for the combination of F 206 and F 235 data sets.

5. June 2010:
Implemented a change in the F 222 data set so that F 222 2 applies to both inflow and outflow sections, not only outflow.

7. June 2010:
Fixed a bug when reading more than 199 W 6 data sets.
Modified the parallel version of the coefficient generation equation for the SIMPLE method.
Hard-coded the use of a sequential solver for the sediments instead of a parallel solver when using the F 164 31/32 options.

8. June 2010:
Implemented a more consistent parallel solver for the sediment concentration, for F 164 31.

11. June 2010:
Improved parallel consistency for the F 235 option.
Added the possibility to make a 2D ParaView file, similar to the 2D Tecplot file, if using F 48 11 in the control file.

14. June 2010:
Fixed a bug in the F 164 31 function.

22. June 2010:
Fixed a small bug in a warning messages written to the boogie file, regarding the discharges.
Added a new solver option: F 164 33, which is a slight modificaiton of the F 164 31 algorithm to reduce the computational time. Note that the new option is not tested yet.

24. August 2010:
Changed the wetting/drying algorithms at startup, so that it is possible to get new wetted areas if the depth in the koordina file allowes it (F 112 1). Also, the G 13 data set now has to be present in the control file after the grid is generated, during the wetting/drying. In earlier versions, the G 13 data set was not needed after the first unstruc file had been made.

25. August 2010:
Fixed a bug when using F 206 0 together with the F 235 data set with integers between 8 and 23.

15. October 2010:
Changed the sediment continuity tests so that the values are given in kg instead of volume units.
Changed the SSIIM menu, so that the variables are chosen from two sections.
Changed the SSIIM Map graphics, so that grid sides without wall laws are drawn in blue instead of red. This would typically be for inflow/outflow areas. This is similar to SSIIM 1. Added two DLLs: flocdll, which includes a formula for flocculation, where the sediments can be moved between the size groups. The vdmldll includes new algorithms for variable bed sediment density and multiple sediment layers. None of the new DLLs are tested yet.

23. October 2010:
Added reading/writing of the con2res file, which is equivalent to the conres file in SSIIM 1. The file is written/read from the menu when the F 37 data set in the control file is above zero. It can also be read using the F letter on the F 2 data set in the control file, or written using and O (letter) on the F 2 data set.
Added the G 26 data set, controlling which blocks are used in the free surface algorithm. An integer is first read, giving the number of further integers on the data set. Then this number of integer are read. Typically, one integer for each block. If the integer is 1, the free water surface is computed. If it is 0, the free water surface is not computed for this block.
Added the F 261 data set. An integer is read. If it is 1, a coupling from the flow field of the nested grid to the coarse grid is computed.

25. October 2010:
Fixed a bug and improved the speed of the nested grid bed interpolation algorithms.

13. November 2010:
Changed the resolution of the grid editor, so that the graphics is faster when using large grids.
Added bed sediment fraction to the graphics. (N is used in G 24).
Added a print-out to the boogie file about the discharges and boundary condition for each discharge group.
Modified the logarithmic inflow boundary condition, so that it can also be used for time-dependent flow.
Simplified the code for the boundary condition for the water flow.
Fixed many bugs for the vdmldll algorithms, but there are still more left.
Added a check for the S data sets in the control file, that the grain sizes are given as the coarsest first, starting with integer 1.

15. November 2010:
Fixed some bugs in the sediment continuity tests for the F 37 3 option.

20. November 2010:
Changed parameters in the call to the FlocDLL.

27. November 2010:
Fixed bugs in reading the fracres file.

30. November 2010:
Fixed a bug in the algorithm for more than two bed layers.
Added an integer sent to the vdmldll function to compute the bed sediment fraction/porosity. This is given on the F 264 data set.

1. December 2010:
Fixed a bug in reading the fracres file for more than 2 bed layers.
Fixed a number of bugs relating to more than 2 sediment layers.
Changed the map graphics to easier see variables for more than 2 sediment layers.

2. December 2010:
Fixed a bug for the inflow boundary condition combined with the F 37 3 option.

9. December 2010:
Fixed bugs in the graphcis for multiple sediment layers.
Added an error message for not using F 64 11 when using the F 37 data set.

10. December 2010:
Modified parameters sent to the vdmldll file.

14. December 2010:
Changed the reading of the bedres file, so that the grid is only regenerated if the F 112 1 option is used. Also removed the exit statement if trying to read a result file that is not there.
Fixed a bug in the F 159 data set when writing/reading the bedres file.

15. December 2010:
Modified the call to the vegdll function.

18. December 2010:
Added a new bed smoothing algorithm: F 144 6 1.0.

21. December 2010:
Added more options to the menu, for different bed sediment distributions and multiple layers.
Modified the writing of Tecplot and ParaView files for the bed sediments.
Added a new initialization of the bed sediment fraction when using F 264 -1.

1. January 2011:
Added print-out of a 2D ParaView file from the menu.

Modified some of the parameters in the 2D ParaView file so that also values in dry areas can be shown. The parameters are: bedform height, bed roughness, bed sediment fraction, bed solid fraction and bed sediment sieve, water depth, bed movements, water level, bed level and sediment thickness.

11. January 2011:
Fixed a bug in the bed concentration algorithm for more than 2 bed layers.
Fixed a bug in the sediment concentration for the F 37 3 option.
Added a geotechnical-based sand slide algorithm: F 266. One integer and two floats are read. The algorithm is used if the integer is 1. The first float is the angle of repose for the material. The second float is the cohesion of the material (Pa). Not tested.

17. January 2011:
Added a test for the use of the F 255 data set in the control file in combination with the F 205 data set and the F 266 data set.

19. January 2011:
Fixed a bug in specifying attraction functions in the Grid Editor.
Fixed a bug and added warning message for the Roughness Editor.
Added varying number of cells in the Roughness Editor.
Added print-out of a bedrough.t file from the menu.
Added the possibility to view the geodata points as lines in the Map view and the Roughness Editor.
Added the possibility to break up lines from the geodata points, by giving a z value smaller than -1.0e6.
Added squares for the Map graphics, as an alternative to the contours.

20. January 2011:
Changed the F 261 data set, so that an integer and a float is read. The float is a relaxation factor for the nested grid interpolation.

22. January 2011:
Fixed a bug in the interpolation from the coarse to the nested grid for the bed cells.
Made new options for the transfer of variables from the nested to the coarse grid:
F 261 10 0.0: transfer of only bed elevation changes.
F 261 11 0.1: new algorithm for transfer of velocities.

25. January 2011:
Fixed a bug in the F 37 1 option for the nested grid.
Added a new F 261 option: 18.
Modified the print-out from the menu to the Tecplot file, so that the block viewed in the Map Graphis is written.
Fixed a bug in the water continuity when interpolating fluxes from the coarse to the nested grid.

26. January 2011:
Fixed a bug in writing the Tecplot 2D file from the menu.

28. January 2011:
Improved the writing of the Tecplot 2D file from the main program, adding a time stamp and enabling writing of the nested grid.
Added the F 268 option, which is an index for which block is to be written to the 2D Tecplot file.

7. February 2011:
Fixed a bug in the use of the cyclic file.
Added the possibility to use an "inflow" file, similar to the "innflow" file for SSIIM2. The "inflow" file can only be used when the inflow discharge is in group 1 and there is no inflow from other discharge groups.
Error in the User's Manual: The description of the F 109 data set states that the tan of the angle of repose is specified. This is actually the inverse of the tan of the angle of repose. (This can also be found when comparing the code in the beddll file with the original formula by Brooks)
Added new options for the F 182 data set: 5 is the empirical function by Prof. Dey. 6 is the same function, but only applied to the border cells.
Fixed a bug in the Brooks formula (F 7 B or F 182 1/2).

10. February 2011:
Fixed a bug in the Dey formula (F 182 5/6), causing non-symmetric erosion.

11. February 2011:
Fixed a bug in writing the 3D Tecplot file from the main program. Added an index to the 3D Tecplot file written from the main program that is 0 for the coarse grid and 1 for the nested grid.

12. February 2011:
Improved the error check for the N data sets in the control file.

14. February 2011:
Fixed a bug for writing the Tecplot 3D file from the main program for the nested grid.
Improved the bed interpolation algorithm from the nested to the coarse grid.

16. February 2011:
Fixed a bug for writing the bed levels and bed level changes to the 2D Tecplot file from the program menu.

17. February 2011:
Fixed a bug for testing of the existence of the "fracres" file.
Fixed a bug in the Dey formula (F 182 5/6).

23. February 2011:
Fixed a bug in the computations of the slope angles in the Dey formula.
Removed the functions called with the F 240 data set.
Fixed a bug for the upstream momentum in the nested grid, when using the F 206 data set.
Renumbered and improved the output from the F 1 P option.
Parallelized the interpolation of the boundary surfaces. Parallelization invoked when using F 269 0.
Added the Lane (1955) formula for the reduction in critical shear stress. Invoked when using F 182 7/8.

26. February 2011:
Fixed the bug with the F 36 15 free surface option and the F 206 data set.
Added a new sand slide algorithm: F 270 201 a b. The a and b are tan(angle of repose). The angle a is for uphill flow and b is for downhill flow. The data set can be used instead of the F 56 data set.
Fixed a bug for the F 81 data set.

2. March 2011:
Fixed a bug where the concentrations in the bed cells were set to zero before the time-dependent computation of the sediment transport (F 2 S).

3. March 2011:
Improved the interpolation of sediment concentration between grid changes.

4. March 2011:
Improved the sediment continuity control for the boundary condition with time-varying water discharges.

5. March 2011:
Added the alternative Pk algorithm of Kato and Launder (1993): F 272 1.
Added wall laws with pressure gradiens: F 273 1. Will only work for F 15 0 (default).

7. March 2011:
Fixed bugs in the F 273 1 algorithm.
Added new algorithm F 273 4: formula by Nakayama and Koyama (1983).

9. March 2011:
Fixed a bug in the F 273 4 formula.

10. March 2011:
Fixed a bug in the F 264 -1 formula.

26. March 2011:
Fixed a bug in the graphics for the surface points in the Grid Editor.

28. March 2011:
Fixed a bug in the F 264 -1 formula.

9. April 2011:
Added the F 275 data set, which reads two floats: the two constants in the epsilon equation of the k-epsilon model, which are by default 1.44 and 1.92.
Added the option of using spatially varying values of the angles of repose. This is done by adding a F 274 1 data set in the control file. The values are given in a file called 'bedangle'. This file has one set of four parameters for each bed cell. The four parameters are read on one line. But first two integers are read on the same line, which indicate which cell it is. The four parameters are: angle of repose for the critical shear stress for upstream slopes, angle of repose for critical shear stress for downstream slopes, limiter for reduction of critical shear stress, and angle of repose for sand slide algorithm. The first three parameters can also be given values on the F 109 data set, for all the cells. The variation in angle of repose for the sand slide algorithm is only coded for the F 56 200 option. Note that all angles in the bedangle file are given in degrees.

12. April 2011:
Fixed a bug in the F 269 0 parallelization algorithm for the interpolation of the discharge locations.
Rewrote both the interpolation algorithm and its parallelization.
Added a test to see if dry cells exist in the grid when writing the unstruc file. If so, a warning message is written to the boogie file and the written unstruc file is given the extension .dry.

18. April 2011:
Fixed a bug in the interpolation algorithm for the discharge locations that occurred when using UTM coordinates in the grid.

28. April 2011:
Fixed a bug in the F 37 2 algorithm.
Fixed a bug for computing bed level changes for the F 68 2 option.

3. May 2011:
Added another continuity control for the sediments.
Added code that sets the highest level in the discharge group equal to the maximum number of grid cells in the vertical direction. This is now default. The code is not used if F 276 0 is given in the control file.

4. May 2011: (v. 101)
Changed the code for the nested grid to allow for a second level of nested grids (not tested).
Fixed a sediment continuity bug for the F 37 2 option.
Modified the walldll calls (F 120) so that SSIIM 2 uses the same functions as is used in the SSIIM 1 DLL.

9. May 2011: (v. 102)
Added a new wall law option, which is invoked by F 15 19 in the control file. This gives smooth walls on the sides and rough on the bed.
Did some slight improvements of the nested grid algorithms.
Fixed a bug in calling the walldll DLL.

10. May 2011: (v. 103)
Fixed a bug in calling the walldll DLL.

11. May 2011: (v. 104)
Added variations of the F 36 7 algorithm. This is invoked with the F 278 data set. An integer is read. Options that seemed to work was 6 and 7. Option 6 is for sub-critical flow, while option 7 has some formulas for supercritical flow.

31. May 2011: (v. 105)
Fixed several bugs in the F 37 3 algorithms. Note that this option will by default use the density function in the vdmldll. To use an algorithm in the main SSIIM program, use the F 264 -1 option.
3. June 2011: (v. 106)
Added the F 279 option. An integer is read. If it is 1 then the bed shear stress will be computed from the maximum of the k value in the bed cell and the k value in the cell above the bed cell. If the integer is 2, then the bed shear stress is computed from the maximum of the bed cell and its neighbours (max 3) if the bed cell is at the boundary of the geometry.

4. June 2011: (v. 107)
Fixed two bugs in the F 261 option, with the nested grids.
Added two options: F 261 16, which interpolates from the fine grid to the coarser grid. F 261 17 interpolateds only from the 2nd level nested grid to the 1st level nested grid.
Added parallelization of the generation of a nested interpolation array. This can be cancelled by using an integer 0 on the second parameter of the F 269 data set. Added this second integer to the F 269 data set. The first integer is a similar option for another nested grid algorithm for interpolation of the boundary conditions between the nested blocks.
Added the F 280 option, giving different options in how the fluxes on the nested grid border are computed. The default is to use an average of the nested and the coarse grid velocities. F 280 1 uses an upwind scheme, with priority to the coarse grid.

7. June 2011: (v. 108)
Added relaxation factors to the F 261 16 and 17 options.

8. June 2011: (v. 108)
Added x and y coordinates for the cell together with distance to bed in the turbdll DLL. The parameters can be used to compute the distance to the wall in a modified turbulence model.

9. June 2011: (v. 109)
Modified the call to the turbdll.dll.
Modified the bed level interpolation between the nested blocks.

18. June 2011: (v. 110)
Added new options for the F 246 data set:
If the first integer is 1 and the water is coming into the cell, or the integer is 2, a limiter will be invoked on the water depth in each cell. Different types of limiters can be used, depending on what the third integer is.
If the second integer is 1, the ap term in the equation is increased if the Froude number is between 0.9 and 1.1. Around supercritical flow the above ap term sometimes can be close to zero, causing a crash. If the second integer is 2, the ap term will never be below 1.0. If the second integer is 3, the ap term will be increased if the Froude number is between 0.5 and 1.5.
The third integer decides the magnitude and conditions of the depth limiter invoked by the first integer. If the integer is above zero, the limiting value is decided by a formula in the sfdifdll.dll file, which can be coded by the user. Negative values indicate different alternative hard-coded formulas, not in the DLL. The following options for the limiter are:
-1: the depth if the Froude number is above 1
-2: the depth * 1.1 if the Froude number is above 1
-3: the depth * the Froude number if the Froude number is above 1
-4: the depth if the Froude number is above 0.9
-5: same as option -3, except the Froude number is taken as the maximum of the cell and its neigbours
-6: the depth if the Froude number is above 0.5
The float is a limiter on the allowable surface slope used in the determination of the anb coefficients for the F 36 7 algorithm.
Added a new option for the F 90 data set: 6. Instead of reducing the bed shear stress due to the dune shear stress partition, it reduces the difference between the bed shear stress and the critical shear stress for the particles.

21. June 2011: (v. 111)
Fixed a bug in the computation of computed time. The bug was in the free surface function.
Added the computed time as a print-out to the boogie file if F 1 D is specified in the control file.
Added a stabilizing function for the F 36 7 algorithm, increasing the ap coefficient when the Froude nubmer is above 1. The function is used if F 283 1 is given in the control file. Note that this option can be used together with the F 246 data set.

24. June 2011: (v. 112)
Fixed a bug in writing Tecplot/Paraview files from more than two sediment layers.
Added completely smooth wall laws for all sides: F 15 5.

1. July 2011: (v. 112)
Fixed a bug in the setting of the pressure for the G 62 data set.
Replaced the F 36 10 algorithm with a new algorithm that is a mixture of F 36 2 and F 36 7.

13. July 2011: (v. 113)
Added more data sets to the timei file: The H data set only reads the water levels upstream and downstream. The Q data set only reads the discharges. The C data set only reads the sediment concentrations. The H data sets must be the first of these three in the file.
Added error messages for the timei file.
Added more functions to the beddll file, with new parameters. The different functions are invoked by the F 285 data set.
Added the "inspace" file, which gives a spatially varying upstream sediment concentration.
Added time and discharge to the Tecplot 3D file.
Added 2nd order time discretization, F 208. (not tested)
Fixed initialization of the grain size distribution of more than 2 bed layers for the F 37 3 algorithm.

15. July 2011: (v. 113)
Fixed two bugs: Reading the I data set in the timei file and using the beddll file when not using a variable bed solid fraction.

19. July 2011: (v. 114)
Fixed a bug in the F 36 10 algorithm, related to the number of grid cells in the lateral direction.

25. July 2011: (v. 115)
Added two parameters to the beddll function computeBedConcentration1.
Added a test for reasonable concentration values from the inspace file.

3. August 2011: (v. 115)
Fixed a bug in the F 36 10 algorithm.

5. August 2011: (v. 115)
Added an algorithm correcting the residual for the pressure if the inflow is not equal to the outflow. The algorithm is invoked automatically.
Added an algorithm where the outflow discharge is increased/decreased as a function of the changes in the water levels in the reservoir. The algorithm is invoked with F 287 2 in the control file. Note that this is not tested yet.

30. August 2011: (v. 116)
Moved the call to the flocdll function, making the algorithm explicit. This fixed a bug and also enabled the use of the flocdll for 1 sediment size, and for both options F 37 2 and F 37 3.

23. September 2011: (v. 117)
Added two options for computing the free surface in the nested grid. F 238 3 0.0 computes the surface from the coarse cells. F 238 4 0.0 scales the pressure in the nested grid so that the average nested surface corresponds to the average coarse grid surface.

26. September 2011: (v. 118)
Fixed a bug in the F 261 18 option.

17. October 2011: (v. 119)
Fixed a bug when reading long timei files.

1. November 2011: (v. 121)
Removed the option of writing porosity files.
Added the possibility to write fracres.new files from the menu.
Added the F 290 option: two intergers are read. If the first one is above zero, the vertical eddy-viscosity will be increased by the integer multiplied with the dune height/depth.
Added the F 291 option: Used to try to prevent instabilities in areas that are drying up. Three integers are read. The default is F 291 0 0 0, meaning none of the algorithms are used. The first integer invokes a relaxation for the velocities after the solver if it is 1. The second integer invokes an ordinary relaxation procedure for the velocities it is 1, and only in 2D cells. If the second integer is 2, the relaxation is also used for 3D cells that are partially dry. If the integer is 3, the relaxation is also used on the neighbour cells. The third integer is a relaxation for the velocity-correction equation for cells that are partially dry.

7. November 2011:
New User's Manual. The Chapters 5.12 and 5.16 are recommended.
(v. 123)
Added F 48 16, which in combination with the interpol file gives profiles of concentrations, similar to what F 48 12 gives for velocities.
Added F 48 17, giving a file with the x and y coordinates of the banks of the geometry.

14. November 2011: (v. 127)
Modified the depth-averaging call in the F 233 algorithm.
Added an inner time step: F 292. An integer and a float are read. If the integer is 1, the float is used as an inner time step to stabilize the solution.
Added more options for the surface smoothing algorithm: F 154. Option 3 will not smooth dry areas, whereas option 1 will. Options 2 and 4 are the same as 1 and 3 respectively, except the boundary is also smoothed. Option 5 adds a diffusive term to the surface solver instead of direct smoothing of the surface. The float is then used as a diffusive factor.
Added print-out of the residuals to a file called 'sigurd'. This file will be written if F 293 1 is given in the control file.

22. November 2011: (v. 128)
Fixed a bug in the computation of time-dependent inflow of sediment concentration.
Fixed a bug in writing the paraview file with option F 48 16.
Added code for reading a time-dependent inspace file. An F 294 data set has to be read, with an integer giving the size of the array for the max number of time steps. Then a C option is used instead of B in the inspace file, where the time step in seconds is read before the other numbers.

24. November 2011: (v. 129)
Modified the reading of the inspace file. Now, the time step is given on a T data set, and the C data set has the same information as the B data set. Example:
A 1
T 0.0
C 1 0.0 0.0 0.0 0.0 0.0 0.001 0.001 0.001 0.002 0.003
C 2 0.0 0.0 0.0 0.0 0.0 0.001 0.001 0.001 0.002 0.003
T 100.0
C 1 0.0 0.0 0.0 0.0 0.0 0.001 0.001 0.001 0.002 0.003
C 2 0.0 0.0 0.0 0.0 0.0 0.002 0.003 0.002 0.002 0.001

26. November 2011: (v. 132)
Fixed a bug in the boundary condition for the inspace file.

5. December 2011: (v. 133)
Added an option 2 for the F 281 data set. This is similar to option 1, except the increase in eddy-viscosity is done for all equations, not only the velocity equations. Also, the increase is done for both the vertical and the horizontal directions.

6. December 2011: (v. 134)
Added the option F 287 3. This is the same as F 287 2, but the outflow water discharge will not go below what is specified in the timei file.

7. December 2011: (v. 135)
Fixed bug in the Tecplot file, for showing the grain size distribution at the bed.

14. December 2011: (v. 136)
Added error messages in the Grid Editor for making grids with larger horizontal dimensions than what is allocated on the G 1 data set.
Added error message to the boogie file for changing the G 1 data set after the unstruc file was created. This will only work for new unstruc files.
Added a warning message in a dialog box for writing the unstruc file when having dry cells in the geometry.
Restructured the reading of the geodata points.
Added an option "K" for the F 2 data set in the control file. This invokes an algorithm where the geodata points are used to compute the bed roughness. An average distance between the bed level and the geodata points inside each 2D cell is computed. This distance is multiplied with a factor and set equal to the roughness. The default factor is 1.0, but other values can be given on the F 295 data set. A diffusion equation is solved to interpolate values in cells with no geodata points. A bedrough.new file is written afterwards.

19. December 2011: (v. 137)
Fixed a bug in the F 2 K algorithm.

13. January 2012: (v. 138)
Added the F 297 data set, which reads an integer between 2 and 4. The integer is the version of the unstruc file that is written. 2 is default. 3 is a new version where the G 1 data set numbers can increase. Also, it contain information about the dry areas, making it not necessary to use a koordina file simultaneously. Note that version 3 is not tested much yet.

19. January 2012: (v. 138)
Added boundary lines in the Map graphics when choosing black colour.

28. January 2012: (v. 138)
Fixed one bug in the F 37 3 function, and added bugtraps for bed fractions above unity and invalid sum of bed grain size fractions. There is still a continuity bug in the F 37 3 function.
Changed sediment concentration formats written to the tecplot files, to enable values under 1 ppm.

30. January 2012: (v. 138)
Fixed a bug in the continuity control for the sediment transport in the F 37 3 function.

4. February 2012: (v.139)
Increased the number of geodata points it is possible to add manually in the Grid Editor from 3 to 5000. If more than 5000 are needed, save the geodata file and restart the program.

7. February 2012: (v. 139)
Fixed a bug for writing the geodata file.

9. February 2012: (v. 139)
Added code that increased the two last integers on the F 65 data set to minimum 2. This will prevent problems with too low values.

24. February 2012: (v. 140)
Fixed two bugs for writing the Froude number and the water depth to the 2D Tecplot file, from the main program. Not the user interface.

29. February 2012. (v. 141)
Fixed a bug in the smooth wall law for the F 15 19 option.
Removed the relaxation for the wall law for epsilon in the wall cells.

8. March 2012. (v. 142)
Improved the OpenFOAM files written from the program.
Improved the F 64 0 algorithm, so that multiple similar points are avoided.

28. March 2012: (v. 142)
Improved F 64 2 and 11 algorithms so that wedge cells are allowed in OpenFOAM files if F 299 1 is used in the control file. Improved content of k, p_rgh, epsilon, alpha1 and nut files to avoid error messages in OpenFOAM. Note that wedge cells do not seem to work for a 3D grid in OpenFOAM.

11. April 2012: (v. 143)
Added more sediment sizes: up to 50.
Added two new grid generation functions: F 64 14 and 15.
Added a new Tecplot function: F 48 18, which writes both the coarse and the nested grid.
Simplified the functions for generation of the a_nb coefficients.

16. April 2012: (v. 143)
Added an error message if the fracres or the N data sets in the control file are not used.
Fixed a bug for more than 19 sediment fractions.
Initiated the bed grain size distribution with the N data sets after allocation of arrays.

17. April 2012: (v. 143)
Fixed a bug introduced 16. April, for computing only water flow.
Added new discretization schemes according to Waterson and Deconinck (1995). These are explict source terms added to the implicit first-order upwind scheme. The numbers are used on the K 6 data set:
10: QUICK
11: Cubic upwind
12: SMART
13: H-QUICK
14: van Leer
15: Superbee
16: Minmod
Note that the implementation is only working for F 206 0, and it is not tested.

20. April 2012: (v. 143)
Fixed a bug in the combination of F 206 0 and F 292 1.

23. April 2012: (v. 143)
Added the K 6 disretization function for F 206 larger than 0. Note that this is not parallelized, so it may run slow.

28. April 2012: (v. 144)
Fixed a bug in the discretization schemes (K 6 10-16).

3. May 2012: (v. 144)
Added thicker and grey lines for the bed velocity vectors when plotted with the surface vectors.
Scaled the velocity vectors with the geometry for the map graphics.

Fixed two bugs in the high-order schemes invoked by 10-17 values on the K 6 data set.

5. May 2012: (v. 144)
Fixed a small bug in the SOU scheme (K 6 1).

11. May 2012: (v. 144)
Fixed a bug in the F 78 data set. The parameter on the data set was set to a fixed value of 0.6, which could not be changed. Now, the value on the data set is used.

26. May 2012: (v. 145)
Increased the maximum number of points in the inflow file to 499x49 (lateral and vertical).

2. June 2012: (v. 145)
Added the F 287 6 option. Two floating point numbers are read on the same line right after the 6 integer. The two floats are by default 0.5 and 2.0. The computed outflow will then be within the limits of 0.5 times the discharge in the timei file and 2.0 times the disharge in the timei file. If F 287 5 is used, no floats are read and the default values are used.

6. June 2012: (v. 145)
Added the F 301 option, reading a floating point number. The number is the criteria for program crash, related to the maximum residual. Default is 1.0e10.
Fixed a problem in a bugtrap for checking the grain size distribution in the top layer.

11. June 2012: (v. 146)
Changed the default unstruc file version to 3 (F 297).
Added an exit dialog box. To end the program use the menu option File->exit.

14. July 2012: (v. 147)
Fixed a bug for the F 238 3 0.0 option (free surface of nested grid) when some part of the grid has dried up. An alternative algorithm is coded as F 238 5 0.0.
Separated the index for which variables are written to files and which are shown in the graphics.

14. August 2012: (v. 148)
Upgraded the Windows compiler to version 12. This lead to a reduction in the file size.
Fixed a bug in the error messages for the inspace file when using multiple sediment sizes.

15. August 2012: (v. 148)
Changed the code so that the OpenMP library libiomp5mt.lib is linked statically into the program, instead of needing the libiomp5md.dll as is the default with the new verson of the Intel compiler.
Added output of xyz coordinates for location of negative concentrations from the inspace file.

18. August 2012: (v. 148)
Modified the print-out at the end of the boogie file for F 1 D and sediment computation to make it more understandable.
Simplified the computation of the bed grain size distribution for F 37 2.
Changed the print-out of the timei file to the boogie file to only be done if there is an error in the boogie file.

23. August 2012: (v. 148)
Added the F 304 2 option. This adds limiters for the F 37 2 algorithm, giving better stability and less continuity problems, especially when the sediment thickness is small.

27. August 2012: (v. 149)
Fixed a bug in the F 37 2 algorithm, which was introduced 18 August 2012. Please do not use version 148 with the F 37 2 option.

17. October 2012: (v. 151)
Fixed a problem with the libiomp5md.dll for the 64 bits version.
Added the F 306 1 option in the control file. With this option, the sediment grain size distribution in the bedres file will be used instead of the fracres or the N data sets when starting a new computation.
Allowed the bedres file to be read even if the F 37 data set is not specified in the control file.
Initialized the grain size distribution with only the coarsest size if no fracres/bedres/N-data sets are used.
Increased the maximum number of geodata points to 500 million.
Modified the print-out of error messages from the computation of the changes in the grain size distribution.
Fixed a bug in the nested grid interpolation that occurred if some of the nested grid was outside the coarse grid.
Changed the walldll function for SSIIM2. The dll function is now different from the one used in SSIIM 1.
Added a new function for the wall laws, invoked by F 305 1. Only F 15 values of 0, 5 and 8 are implemented in the new function.
Parallelized the bed interpolation of vertical values from the geodata file to the grid.

19. October 2012: (v. 152)
Added four new F 48 options: 88, 89, 91 and 92. They will do the same as the options 8, 9, 11 and 12, except the result and bedres files will not be written.

22. October 2012: (v. 152)
Fixed a bug for estimating the water discharge for times longer than what is given in the timei file.
Modified the squares graphics.

25. October 2012: (v. 152)
Modified the roughness editor, so all given roughness values can be seen with crosses.
Added a second float on the F 47 data set, which gives a distance which outide a faster and less accurate method is used for choosing points for interpolation. Typical values could be the same size as the grid cell. The use of this parameter may give a faster bed interpolation.
Added print-out of the grid of the last block when the program crashes during grid generation due to too small memory allocations. Further instructions for restart is in the boogie file.
Changed the logic to avoid two grid generations at startup.
Added the location of the G 6 data points in the Discharge Editor.

14. November 2012: (v. 152)
Fixed a bug in writing the bedres file.
Fixed a bug in reading the unstruc file.
Added the F 78 10 option, which uses the average diameter for the computation, instead of using each size fraction.
Changed the F 78 option for 1 sediment size, so that the 2nd parameter on the F 78 data set is used, instead of a hard-coded value of 0.6.
Modified two tests for the G 6 data set.

20. November 2012: (v. 153)
Changed the code for the G 6/62 data sets, fixing some bugs.
Added print-out to the boogie file of the 2D indexes resulting from the G 6/62 data sets.
Added a termination of the programme if not sufficient memory was allocated on the F 65 data set for the multi-grid solver.

23. November 2012: (v. 154)
Fixed bugs in the printing of the paraview 2D VTK file.
Added the F 48 19 option: printing of one 2D paraview file for each block in the grid. This is similar to the F 48 18 option for the Tecplot files.

28. November 2012: (v. 155)
Fixed a bug in reading the koomin file. This bug was introduced 24. October 2012.
Made an error message for the problem of different sediment thickness in the fracres/bedres file compared with the grid/koomin file. Made two algorithms to overcome the problem. F 310 1 adjusts the sediment layer thicknesses according to the koomin file values. F 310 2 adjusts the values from the koomin file according to the thicknesses of the sediment layers.

18. February 2013: (v. 157)
Added MPI connections. Fixed a bug in the generation of the koosurf files for the MPI implementation F 2 UA.

21. February 2013: (v. 157)
Fixed bugs in the MPI version of SSIIM 2, when a time step is used and for the final print-out of results.

28. February 2013: (v. 158)
Fixed a number of bugs in the MPI version nof SSIIM 2, for the connection between the blocks. Fixed bugs in the paraview files written from the MPI version. Added the F 315 1 data set (not tested), which causes the upstream and downstream cross-section of each block not to move laterally, which they would have done similar to the banks. This only works for the MPI version.

14. March 2013: (v. 158)
Added the F 316 data set, which gives different options of removing areas of high pressure at the downstream side of MPI blocks. An integer is read. Option 8 seems to give the best results.
Added the F 317 data set. Two integers are read. The first integer is the block number. The second integer is the 2D cell index in the lateral direction at upstream cross-section for a reference cell for the pressure. The parameter is used in the F 316 data set to specify a reference cell for the pressure algorithm.
Example: F 317 3 148. If the cross-section has 200 cells, cell no. 148 is used in block 3 to be reference for the pressure. In the F 316 8 algorithm, the pressure in this cell will be used in all the cells with too high pressure.
Removed the line written to the boogie file with the cell number of the cell with the highest residual and its information.

20. March 2013: (v. 158)
Changed the default relaxation parameter for the free surface from 1.0 to 0.1 (G 6 data set).
Disabled error control checks for the G 6/62 data sets for the MPI version. The reference cell for the MPI version is given on the F 317 data set, not the G 6 data set.

3. April 2013: (v. 159)
Fixed a bug for computing turbulent diffusion for the sediment concentration after update of the grid from the free water surface computation.
Wrote more information to the boogie file when the solution diverges, and also a 2D ParaView file with graphics from the crash.

5. April 2013: (v. 159)
Improved the speed for writing 3D ParaView files via the F 48 data set. Also changed the format of the variables to an exponential form in the file.

8. April. 2013: (v. 159)
Fixed a bug that caused the program to crash on some occations.

8. May 2013: (v. 160)
Rewrote many of the free surface algorithms (F 36). The new algorithms are invoked by F 323 1 in the control file. (Default F 323 0). The new function includes a new algorithm for the free surface, which is chosen by F 36 9, in combination with F 323 1. The new algorithm is based on the Shallow water equations. The F 322 data set reads two integers. If the first one is 1, then the time term is included in the continuity equation of the algorithm. If the second one is 1, then the momentum term is included excplicitly in the Shallow water equation. Default: F 322 0 0. Note that even though the new free surface algorithm is based on the Shallow water equations, it only computes changes in the free surface. The Navier-Stokes equations are still used to compute the velocities.
Added the F 286 2 option, that stabilizes the free surface computations by reducing the pressure-corrections.

16. May 2013: (v. 160)
Found a bug in the F 287 data set: The outflow flux was not changed. Fixed the bug and added a new option: F 287 7. The option does the same as the 3 option, but an algorithm is added where the flux caused by the water surface elevaton changes is included in the SIMPLE corrections.
Fixed a bug in the longitudinal profile graphics, occurring when using a large numer of cells in the grid.

9. July 2013: (v. 161)
Made an algorithm where the values on the F 94 data set are computed as 0.05 % of the average cell length for the first parameter and 5 % of the average cell length for the second parameter. The algorithm is used when the first float on the F 94 data set is negative. This is now the default. The values are written to the boogie file if F 1 D is used in the control file.
Added two more F 235 options: 15 and 16. The option 15 makes the cells at the boundary more horizontal. Option 16 makes extremely thin cells more deep.

7. August 2013: (v. 161)
Fixed a bug in the coefficient generation for the power-law scheme.

14. August 2013: (v. 161)
Modified the print-out of result files.
Added the F 54 data set, which is similar to what is used in SSIIM 1. The default is 2e-3, though.

30. August 2013: (v. 161)
Fixed allocation bugs for generating the grid without using an F 37 1/2 data set in the control file.
Fixed bug in map graphcs when not using F 37 1/2 in the control file.
Print of error message to boogie file when the integers on the F 314 data set is above 1.

7. September 2013: (v. 161)
Fixed another bug that caused the program to crash if the F 37 data set was not used.

10. September 2013: (v. 162)
Fixed a bug in the F 329 data set.

13. September 2013: (v. 162)
Fixed a bug in the arrays of the surfaces in the grid, emerging when computing the water flow with free surface and not sediment transport.
Fixed bugs in the print-out of result files, paraview files etc.

12. December 2013: (v. 163)
Changed the default grid algorithm from F 64 0 to F 64 11.
Changed a number of algorithms in the soil slide functions (F 236).

13. February 2014: (v. 164)
Fixed a bug in the Map graphics, which showed lines between the geodata points.

17. February 2014: (v. 164)
Fixed a bug in the Discharge Editor, which showed lines between the geodata points.
Added a velocity damping function in shallow areas. Invoked by F 328 1 and F 341 1.
Added an algorithm for more changes in the shapes of cells. It is invoked by using the number 6 on the 5th integer of the F 159 data set.
Added the G 22 data set, to be used when modelling bridges. Four integers and three floats are read. The first four integers denote the cells with the bridge, in the horizontal plane. The first float is the lowest elevation of the bridge deck, the second float is the highest elevation. The last float is a drag/solidfraction coefficient. Up to 19 G 22 data sets can be used.

6. March 2014: (v. 165)
Fixed a bug when regenerating the grid and there is an F 64 data set in the control file. The bug was related to changing the default value of the F 64 data set from 0 to 11, which took place on the 10th of December 2013. This means the bug is in versions 163 and 164.

10. May 2014: (v. 168)
Implemented the Engelund-Hansen formula with F 84 6 in the control file. This is a different code than F 10 E, which is not correctly implemented. F 84 6 does not use the Shields function to decide if there is erosion or not, according to the Engelund-Hansen formula.

13. June 2014: (v. 169)
Version 168 did not call the DLLs, which is fixed in the 32 bits 169 version.

16. June 2014: (v. 170)
Fixed a bug for the time term in the k equation when running only on one core. The bug lead to instabilities.

24. July 2014: (v. 171)
Fixed bugs in writing the ParaView and Tecplot files, particularly for bed parameters, for example grain size fractions, and also particulaly for the Linux version.

30. August 2014: (v. 172)
Fixed a bug for interpolating sediment concentrations when the grid changes.
New feature in the grid editor: generate a straight line with NoMovePoints between the two last made NoMovePoints. This is the Generate->NoMovePLine option in the menu.
Added an option to write the average sediment concentration at the outflow section to the timeo file. This is done by specifying an "a" on the O data set in the timei file. This is thought to be useful when simulating a river with several dams. Then each section between the dam can be simulated at the time, and the concentrations from the upstream reach can be given as input to the downstream reach using the timei/timeo files.
Added a new option for the F 234 data set, with an integer 3. The function increases or decreases the outflow discharge according to the velocity gradient at the boundary and the computed water level in comparison to the specified level.
Added a new option: F 347. If an integer 1 is read, then the free surface water level is rised at the location of a sand slide, with a similar magnitude as the bed level rises due to the slide.
Modified the F 64 112 option, so that the removal of the shallow cells are done for cells smaller than the 2nd parameter on the F 64 data set. Added a convergence criteria for the F 36 9 option: If the average water surface elevation change is under 0.1 %, then the water surface and grid is not updated.
Added more checks for errors in the geodata file.

29. October 2014: (v. 173)
Fixed a bug in the F 292 data set in the control file, which only emerged if the F 206 data set was not in the file, or the integer on the F 206 data set was set to zero.

20. Mach 2015: (v. 175)
Fixed a bug in writing the 2D paraview file.
Fixed a bug in writing inflow sediment concentrations with the a option in the timei file.

2. April 2015: (v. 176)
Added the F 356 data set: reads an integer which is the interval for plotting geodata points. This interval is also used when writing the geodata.new file. The purpose is to reduce the size of very large geodata files.
Added an option in the Grid Editor to define a triangular region for which the geodata points are deleted. The deletion is done from the Generate option in the menu. The deletion also affects the points in the geodata.new file that can be written from the program.
Recoded the writing of the depth-averagaed velocities to the 2D ParaView file.

24. September 2015: (v. 180)
Removed the control for decreasing sediment sizes on the S data sets if the F 260 data set is used.
Added the F 365 data set. If an integer equal to 1 is read (F 365 1), then the points in the interpol file will show up in the map graphics, if they have been read.

24. January 2016: (v. 182)
Used a new bedll with more functions.
Fixed a bug in the F 102 algorithm, changing the cell shapes.

15. April 2016: (v. 184)
Fixed a bug in reading the timei file.

20. May 2016: (v. 186)
Implemented a change in the G 6 data set. If the third integer is non-zero, indicating a discharge group, the first integer can be set to zero for transient sediment computations. The program will find a reference cell.
Implemented the alldata file. This is a combination of the unstruc, result, bedres and con2res files. The file can be used for restarts. Reading the alldata file instead of the unstruc file on the F 2 data set is done by using a B instead of a U.

10. August 2016: (v. 188)
Implemented the Schmidt number (F 12) for the sediment computations.

12. December 2023 (v. 550)
Fixed bugs in the GridEditor.

4. June 2024 (v. 586)
Fixed a bug in the GridEditor that caused a crash if the geodata file is read before a new grid block is made.