This testcase is based on:
https://svn.oss.deltares.nl/repos/delft3d/branches/research/Deltares/20151209_34238_slib3D/examples/19_slib3d_McKay

The information below belongs to the original McKay testcase:
============================================================

Simulatie met Gensig-model met aanpassingen op verzoek van Sitoni/Talmon (valsnelheid zand wordt beinvloed door slibconcentratie (fallsand).:
Recent overleg tussen Rob En Arno hebben nog wat aanvullende aanpassingen in de subroutine BNGHAM en FALLSAND opgeleverd.

RVW: Q_H model; waterstandsrand hard: alfa = 0.0

bin_cvisc = 0.252

Initiele conditie met schuine waterstand en put.
2 fracties: zandconcentratie (963, was 1000 g/l), let op dit moet stof 1 zijn!!!
en uniforme slibconcentratie (fines) ( 193, was 180 g/l), let op dit moet stof lsed zijn (stof 2)


Bij de uitstroomrand ontstaat er een zandhoop, zonder toevoeging van een put bij de uitstroomrand.

Fallvelocity mud = 0.0 (in code);

Csoil = 1193.

Power in fallsand veranadert:

nf_floc=2.15
power     = 2./(3.-nf_floc)

Deze waarde staat in de invoer, maar ook hard in de code
om de parameter niet door te hoeven geven naar de BNGHAM-subroutine.

   FlcNf    = 2.15
   FlcK2    = 1.
   FlcK3    = 5.
   FlcK4    =  .0070
   FlcD0    = 10.e-06
   FlcD     = 0.44
   FlWsmin  = 0.3e-3
   EroSchel = false
   EroSns   = 0.
   EroScv   = 3.e-9
   EroScu   = 400.
   EroSd50  = 20.e-6
   cons_mud = true
   cons_kk  = 2.0e-14
   cons_ksig= 2.0e8
   cons_eta = 1.0e5
   bingham     = true
   bin_cvisco  = 0.252
   bin_cnvisco = 2.35
   bin_cyield  = 30.0
   bin_abingh  = 1000
   ero_intfc   = true
   conc_intfc  = 25.0
   
   In de sourcecode wordt irobed op 1 gezet en de valsnelheid van de fines op 0.0.
   
   Let op consolidatie uitgeschakeld in deze .sed file door cons_kk op nul te zetten.
   
   Let op bij plotten settling velocity. Moet de grootheid op het generieke sigma rooster zijn, de andere
   is afkomstig uit FALLVE i.p.v. uit SETTOT. Valsnelheid moet het gedrag van dudz vertonen bij de bodem.
   
   Voor deze condities ontstaan er golven in de oplossing. Daarom simulatie herhaald met andere formulering voor
   de Binghmstress. Fractale dimensie aangepast van 2.15 naar 1.5 in subroutine Bngham en bin_cyield van 37.8 naar 30.0 in kk3_2frac.sed..
   
   Simulatie gedaan op verzoek van Arno Talmon zonder turbulentuiemodel, dan verdwijnen slingeringen.
   Constant. Vicoww =10-3 en Dicoww=10-4.