In [1]:
%%javascript
IPython.OutputArea.prototype._should_scroll = function(lines) {
    return false;
}
In [2]:
# Import packages
import numpy as np
import matplotlib.pyplot as pl
from flumyx import *
from flumyx.api import *
from collections import Counter
In [3]:
# Global configuration
trace_level = TL_WARNING # 0: None, 1: Fatal, 2: Error, 3: Warning, 4: Info, 5: Debug
# Verbose mode
#trace_level = TL_DEBUG

# Turbidites grain size templates distribution
nb_gr = 15
#              1   2   3   4   5   6   7   8   9  10  11  12  13  14  15
gr_temp = [ [  0,  1,  2,  3,  4,  8, 17, 25, 22,  8,  5,  3,  2,  0,  0 ], # Standard
            [  0,  0,  0,  0,  2,  3,  6, 13, 27, 26, 10,  8,  5,  0,  0 ], # Sandy
            [  0,  0,  0,  0,  0, 10, 13, 50, 22,  5,  0,  0,  0,  0,  0 ], # Fine grained
            [  2,  3,  3,  5, 10, 24, 29, 22,  2,  0,  0,  0,  0,  0,  0 ], # Muddy
            [  0,  0,  0,  0,  0,  0,  0,  0,  0,  8, 16, 28, 32, 16,  0 ], # Ultra Sandy
            [ 20, 30, 24, 17,  7,  2,  0,  0,  0,  0,  0,  0,  0,  0,  0 ], # Ultra Muddy
            [  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 50, 50,  0 ], # Ultra Cobble
            [  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,100,  0 ], # Single class 14
          ]
In [4]:
# Just a function which displays the Flumy version currently loaded
def printVersion(simu):
    old_tl = simu.getTracer().getTraceLevel()
    simu.getTracer().setTraceLevel(TL_DEBUG)
    simu.printVersion()
    simu.getTracer().setTraceLevel(old_tl)

# Just a function wich sets the required grain size distribution
def setGrainsizeTemplate(simu, index):
    if (index < 1 or index > 7):
        print("Error: invalid grain size template index")
        return
    for i in range(nb_gr):
        param = simu.getParameters()
        param.setParamIntIdx(KEY_SED_LOAD_VOL, i, gr_temp[index][i])
        print("Setting proportion (%) for grain size class ", i+1, " = ", gr_temp[index][i])
In [5]:
def runSimu(simu, factor=1, vargr=False, optim=False, proxy=False, tmpidx=0, rth=0.75) :
    
    print("Running simulation...")

    # Check input flags
    if (optim):
        if (proxy): print("proxy and optim cannot be True all together => using optim")
        proxy = False
        vargr = True
    if (proxy):
        optim = False
        vargr = True     
    
    # Print Flumy version
    printVersion(simu)

    # Reset the simulation (return at age 0)
    simu.getTracer().setTraceLevel(TL_WARNING)
    simu.clear()
    simu.getTracer().setTraceLevel(trace_level)
     
    # Configure the environement
    # Turbidite simulation
    simu.getParameters().setUsingTurbi()

    # Update parameters
    hmax = simu.getParameters().getParamDouble(KEY_CHNL_MAX_DEPTH)
    hmax = hmax * factor
    applyNexus(simu, hmax = hmax, adaptGrid=False)
    zul = 3 * hmax
    
    # Use variable grain size in PB/LAPs ?
    simu.getParameters().setParamBool(KEY_GR_PB, vargr)
    simu.getParameters().setParamBool(KEY_GR_PB_OPTIM, optim)
    simu.getParameters().setParamBool(KEY_GR_PB_PROXY, proxy)
    simu.getParameters().setParamDouble(KEY_GR_THRESHOLD, rth)

    # Set grain size template from pre-defined template index
    setGrainsizeTemplate(simu, tmpidx)

    # Simulation root seed
    seed = 665544
    
    # Launch the simulation up to zul
    print(simu.getParameters().getParamsDesc(all = True))
    runFlumy(simu, 0, seed, zul=zul, trace_level=trace_level)
    
    print("Simulation finished!")
In [6]:
# Create simulator and show Flumy version
simu = iSimulator.createSimulator()
printVersion(simu)
Flumy (8.000) / 2025-02-03
In [7]:
# Use original Variable Grain Size model with Standard template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=0)
showGrainDist(getGrainDist(simu, True), "Standard", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Standard")
Running simulation...
Flumy (8.000) / 2025-02-03
Error: invalid grain size template index
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 1 (1) [0, 100]
SED_LOAD_VOL_10 = 5 (5) [0, 100]
SED_LOAD_VOL_11 = 3 (3) [0, 100]
SED_LOAD_VOL_12 = 2 (2) [0, 100]
SED_LOAD_VOL_13 = 0 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 2 (2) [0, 100]
SED_LOAD_VOL_3 = 3 (3) [0, 100]
SED_LOAD_VOL_4 = 4 (4) [0, 100]
SED_LOAD_VOL_5 = 8 (8) [0, 100]
SED_LOAD_VOL_6 = 17 (17) [0, 100]
SED_LOAD_VOL_7 = 25 (25) [0, 100]
SED_LOAD_VOL_8 = 22 (22) [0, 100]
SED_LOAD_VOL_9 = 8 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [8]:
# Use original Variable Grain Size model with Sandy template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=1)
showGrainDist(getGrainDist(simu, True), "Sandy", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Sandy")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  0
Setting proportion (%) for grain size class  2  =  0
Setting proportion (%) for grain size class  3  =  0
Setting proportion (%) for grain size class  4  =  0
Setting proportion (%) for grain size class  5  =  2
Setting proportion (%) for grain size class  6  =  3
Setting proportion (%) for grain size class  7  =  6
Setting proportion (%) for grain size class  8  =  13
Setting proportion (%) for grain size class  9  =  27
Setting proportion (%) for grain size class  10  =  26
Setting proportion (%) for grain size class  11  =  10
Setting proportion (%) for grain size class  12  =  8
Setting proportion (%) for grain size class  13  =  5
Setting proportion (%) for grain size class  14  =  0
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 0 (1) [0, 100]
SED_LOAD_VOL_10 = 10 (5) [0, 100]
SED_LOAD_VOL_11 = 8 (3) [0, 100]
SED_LOAD_VOL_12 = 5 (2) [0, 100]
SED_LOAD_VOL_13 = 0 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 0 (2) [0, 100]
SED_LOAD_VOL_3 = 0 (3) [0, 100]
SED_LOAD_VOL_4 = 2 (4) [0, 100]
SED_LOAD_VOL_5 = 3 (8) [0, 100]
SED_LOAD_VOL_6 = 6 (17) [0, 100]
SED_LOAD_VOL_7 = 13 (25) [0, 100]
SED_LOAD_VOL_8 = 27 (22) [0, 100]
SED_LOAD_VOL_9 = 26 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [9]:
# Use original Variable Grain Size model with Fine Grained template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=2)
showGrainDist(getGrainDist(simu, True), "Fine Grained", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Fine Grained")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  0
Setting proportion (%) for grain size class  2  =  0
Setting proportion (%) for grain size class  3  =  0
Setting proportion (%) for grain size class  4  =  0
Setting proportion (%) for grain size class  5  =  0
Setting proportion (%) for grain size class  6  =  10
Setting proportion (%) for grain size class  7  =  13
Setting proportion (%) for grain size class  8  =  50
Setting proportion (%) for grain size class  9  =  22
Setting proportion (%) for grain size class  10  =  5
Setting proportion (%) for grain size class  11  =  0
Setting proportion (%) for grain size class  12  =  0
Setting proportion (%) for grain size class  13  =  0
Setting proportion (%) for grain size class  14  =  0
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 0 (1) [0, 100]
SED_LOAD_VOL_10 = 0 (5) [0, 100]
SED_LOAD_VOL_11 = 0 (3) [0, 100]
SED_LOAD_VOL_12 = 0 (2) [0, 100]
SED_LOAD_VOL_13 = 0 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 0 (2) [0, 100]
SED_LOAD_VOL_3 = 0 (3) [0, 100]
SED_LOAD_VOL_4 = 0 (4) [0, 100]
SED_LOAD_VOL_5 = 10 (8) [0, 100]
SED_LOAD_VOL_6 = 13 (17) [0, 100]
SED_LOAD_VOL_7 = 50 (25) [0, 100]
SED_LOAD_VOL_8 = 22 (22) [0, 100]
SED_LOAD_VOL_9 = 5 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [10]:
# Use original Variable Grain Size model with Muddy template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=3)
showGrainDist(getGrainDist(simu, True), "Muddy", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Muddy")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  2
Setting proportion (%) for grain size class  2  =  3
Setting proportion (%) for grain size class  3  =  3
Setting proportion (%) for grain size class  4  =  5
Setting proportion (%) for grain size class  5  =  10
Setting proportion (%) for grain size class  6  =  24
Setting proportion (%) for grain size class  7  =  29
Setting proportion (%) for grain size class  8  =  22
Setting proportion (%) for grain size class  9  =  2
Setting proportion (%) for grain size class  10  =  0
Setting proportion (%) for grain size class  11  =  0
Setting proportion (%) for grain size class  12  =  0
Setting proportion (%) for grain size class  13  =  0
Setting proportion (%) for grain size class  14  =  0
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 2 (0) [0, 100]
SED_LOAD_VOL_1 = 3 (1) [0, 100]
SED_LOAD_VOL_10 = 0 (5) [0, 100]
SED_LOAD_VOL_11 = 0 (3) [0, 100]
SED_LOAD_VOL_12 = 0 (2) [0, 100]
SED_LOAD_VOL_13 = 0 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 3 (2) [0, 100]
SED_LOAD_VOL_3 = 5 (3) [0, 100]
SED_LOAD_VOL_4 = 10 (4) [0, 100]
SED_LOAD_VOL_5 = 24 (8) [0, 100]
SED_LOAD_VOL_6 = 29 (17) [0, 100]
SED_LOAD_VOL_7 = 22 (25) [0, 100]
SED_LOAD_VOL_8 = 2 (22) [0, 100]
SED_LOAD_VOL_9 = 0 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [11]:
# Use original Variable Grain Size model with Ultra Sandy template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=4, rth=0.95)
showGrainDist(getGrainDist(simu, True), "Ultra Sandy", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Ultra Sandy")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  0
Setting proportion (%) for grain size class  2  =  0
Setting proportion (%) for grain size class  3  =  0
Setting proportion (%) for grain size class  4  =  0
Setting proportion (%) for grain size class  5  =  0
Setting proportion (%) for grain size class  6  =  0
Setting proportion (%) for grain size class  7  =  0
Setting proportion (%) for grain size class  8  =  0
Setting proportion (%) for grain size class  9  =  0
Setting proportion (%) for grain size class  10  =  8
Setting proportion (%) for grain size class  11  =  16
Setting proportion (%) for grain size class  12  =  28
Setting proportion (%) for grain size class  13  =  32
Setting proportion (%) for grain size class  14  =  16
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 0 (1) [0, 100]
SED_LOAD_VOL_10 = 16 (5) [0, 100]
SED_LOAD_VOL_11 = 28 (3) [0, 100]
SED_LOAD_VOL_12 = 32 (2) [0, 100]
SED_LOAD_VOL_13 = 16 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 0 (2) [0, 100]
SED_LOAD_VOL_3 = 0 (3) [0, 100]
SED_LOAD_VOL_4 = 0 (4) [0, 100]
SED_LOAD_VOL_5 = 0 (8) [0, 100]
SED_LOAD_VOL_6 = 0 (17) [0, 100]
SED_LOAD_VOL_7 = 0 (25) [0, 100]
SED_LOAD_VOL_8 = 0 (22) [0, 100]
SED_LOAD_VOL_9 = 8 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.95 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [12]:
# Use original Variable Grain Size model with Ultra Muddy template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=5)
showGrainDist(getGrainDist(simu, True), "Ultra Muddy", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Ultra Muddy")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  20
Setting proportion (%) for grain size class  2  =  30
Setting proportion (%) for grain size class  3  =  24
Setting proportion (%) for grain size class  4  =  17
Setting proportion (%) for grain size class  5  =  7
Setting proportion (%) for grain size class  6  =  2
Setting proportion (%) for grain size class  7  =  0
Setting proportion (%) for grain size class  8  =  0
Setting proportion (%) for grain size class  9  =  0
Setting proportion (%) for grain size class  10  =  0
Setting proportion (%) for grain size class  11  =  0
Setting proportion (%) for grain size class  12  =  0
Setting proportion (%) for grain size class  13  =  0
Setting proportion (%) for grain size class  14  =  0
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 20 (0) [0, 100]
SED_LOAD_VOL_1 = 30 (1) [0, 100]
SED_LOAD_VOL_10 = 0 (5) [0, 100]
SED_LOAD_VOL_11 = 0 (3) [0, 100]
SED_LOAD_VOL_12 = 0 (2) [0, 100]
SED_LOAD_VOL_13 = 0 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 24 (2) [0, 100]
SED_LOAD_VOL_3 = 17 (3) [0, 100]
SED_LOAD_VOL_4 = 7 (4) [0, 100]
SED_LOAD_VOL_5 = 2 (8) [0, 100]
SED_LOAD_VOL_6 = 0 (17) [0, 100]
SED_LOAD_VOL_7 = 0 (25) [0, 100]
SED_LOAD_VOL_8 = 0 (22) [0, 100]
SED_LOAD_VOL_9 = 0 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [13]:
# Use original Variable Grain Size model with Ultra Cobble template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=6)
showGrainDist(getGrainDist(simu, True), "Ultra Cobble", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Ultra Cobble")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  0
Setting proportion (%) for grain size class  2  =  0
Setting proportion (%) for grain size class  3  =  0
Setting proportion (%) for grain size class  4  =  0
Setting proportion (%) for grain size class  5  =  0
Setting proportion (%) for grain size class  6  =  0
Setting proportion (%) for grain size class  7  =  0
Setting proportion (%) for grain size class  8  =  0
Setting proportion (%) for grain size class  9  =  0
Setting proportion (%) for grain size class  10  =  0
Setting proportion (%) for grain size class  11  =  0
Setting proportion (%) for grain size class  12  =  0
Setting proportion (%) for grain size class  13  =  50
Setting proportion (%) for grain size class  14  =  50
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 0 (1) [0, 100]
SED_LOAD_VOL_10 = 0 (5) [0, 100]
SED_LOAD_VOL_11 = 0 (3) [0, 100]
SED_LOAD_VOL_12 = 50 (2) [0, 100]
SED_LOAD_VOL_13 = 50 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 0 (2) [0, 100]
SED_LOAD_VOL_3 = 0 (3) [0, 100]
SED_LOAD_VOL_4 = 0 (4) [0, 100]
SED_LOAD_VOL_5 = 0 (8) [0, 100]
SED_LOAD_VOL_6 = 0 (17) [0, 100]
SED_LOAD_VOL_7 = 0 (25) [0, 100]
SED_LOAD_VOL_8 = 0 (22) [0, 100]
SED_LOAD_VOL_9 = 0 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [14]:
# Use original Variable Grain Size model with Single Class 14 template
runSimu(simu, factor=1, vargr=True, optim=False, proxy=False, tmpidx=7)
showGrainDist(getGrainDist(simu, True), "Single class 14", 0.6)
fac,grain,age = simu.getBlock(1, zb=0)
showSection(grain = grain[50,:,:], size = 8, legend=1, title="Single class 14")
Running simulation...
[0] Flumy (8.000) / 2025-02-03
Setting proportion (%) for grain size class  1  =  0
Setting proportion (%) for grain size class  2  =  0
Setting proportion (%) for grain size class  3  =  0
Setting proportion (%) for grain size class  4  =  0
Setting proportion (%) for grain size class  5  =  0
Setting proportion (%) for grain size class  6  =  0
Setting proportion (%) for grain size class  7  =  0
Setting proportion (%) for grain size class  8  =  0
Setting proportion (%) for grain size class  9  =  0
Setting proportion (%) for grain size class  10  =  0
Setting proportion (%) for grain size class  11  =  0
Setting proportion (%) for grain size class  12  =  0
Setting proportion (%) for grain size class  13  =  0
Setting proportion (%) for grain size class  14  =  100
Setting proportion (%) for grain size class  15  =  0
SIM_TYPE = Turbidite (Fluvial)
AG_RAND_LAMBDA = 1 (1)
AV_LV_OB = 0 (0)
AV_NO_SPLAY = 0 (0)
COND_WELL_CLOSURE_ON = 0 (0)
COND_WELL_STOP_HONORED = 0 (0)
GR_PB = 1 (0)
GR_PB_OPTIM = 0 (0)
GR_PB_PROXY = 0 (0)
GR_QUICK_C0 = 0 (0)
MASS_BAL_ACTIVE = 0 (0)
AG_DRAP = 13 (13) [-2147483647, 2147483647]
AG_EP_DIST = 0 (0) [-2147483647, 2147483647]
AG_EP_FREQ = 0 (0) [-2147483647, 2147483647]
AG_EP_PERIOD = 24 (24) [1, 2147483647]
AG_EP_POISSON = 24 (24) [1, 2147483647]
AG_OB_DIST = 2 (2) [-2147483647, 2147483647]
AG_OB_FREQ = 2 (2) [-2147483647, 2147483647]
AG_OB_PERIOD = 24 (24) [1, 2147483647]
AG_OB_POISSON = 24 (24) [1, 2147483647]
AG_TYPE = 2 (2) [-2147483647, 2147483647]
AV_LOC_FREQ = 2 (2) [-2147483647, 2147483647]
AV_LOC_PERIOD = 264 (288) [1, 2147483647]
AV_LOC_POISSON = 264 (288) [1, 2147483647]
AV_REG_FREQ = 2 (2) [-2147483647, 2147483647]
AV_REG_PERIOD = 470 (470) [1, 2147483647]
AV_REG_POISSON = 470 (470) [1, 2147483647]
CHNL_SCALE_DIST = 0 (0) [-2147483647, 2147483647]
DOMAIN_NX = 251 (251) [20, 5000]
DOMAIN_NY = 201 (201) [20, 5000]
GR_PB_OPTIM_MAX = 13 (13) [3, 13]
GR_PB_OPTIM_MIN = 6 (6) [1, 10]
SED_LOAD_VOL_0 = 0 (0) [0, 100]
SED_LOAD_VOL_1 = 0 (1) [0, 100]
SED_LOAD_VOL_10 = 0 (5) [0, 100]
SED_LOAD_VOL_11 = 0 (3) [0, 100]
SED_LOAD_VOL_12 = 0 (2) [0, 100]
SED_LOAD_VOL_13 = 100 (0) [0, 100]
SED_LOAD_VOL_14 = 0 (0) [0, 100]
SED_LOAD_VOL_2 = 0 (2) [0, 100]
SED_LOAD_VOL_3 = 0 (3) [0, 100]
SED_LOAD_VOL_4 = 0 (4) [0, 100]
SED_LOAD_VOL_5 = 0 (8) [0, 100]
SED_LOAD_VOL_6 = 0 (17) [0, 100]
SED_LOAD_VOL_7 = 0 (25) [0, 100]
SED_LOAD_VOL_8 = 0 (22) [0, 100]
SED_LOAD_VOL_9 = 0 (8) [0, 100]
SIM_SEED = 165426111 (165426111) [1, 2147483647]
ZUL_TYPE = 0 (0) [-2147483647, 2147483647]
AG_EP_CONST = 2 (2) [-1e+30, 1e+30]
AG_EP_INIT_ELEV = 9999 (9999) [-1e+30, 1e+30]
AG_EP_LOGNORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EP_MAX = 2.6 (2.6) [-1e+30, 1e+30]
AG_EP_MIN = 1.4 (1.4) [-1e+30, 1e+30]
AG_EP_NORM_MEAN = 2 (2) [-1e+30, 1e+30]
AG_EP_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_EXP_DEC_THICK = 7993 (8107) [100, 1e+30]
AG_LV_WIDTH = 2 (2) [0.05, 1e+30]
AG_OB_CONST = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_LOGNORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_MAX = 2.6 (2.6) [0.03, 1e+30]
AG_OB_MIN = 1.4 (1.4) [0, 1e+30]
AG_OB_NORM_MEAN = 2 (2) [0, 1e+30]
AG_OB_NORM_STDEV = 0.6 (0.6) [0, 1e+30]
AG_OB_PEAT = 0 (0) [0, 99]
AG_OB_PEL_DRAP_INT = 0 (0) [0, 100]
AV_LOC_PROB1 = 0.5 (0.5) [0, 1]
AV_LOC_PROB2 = 0.9 (0.9) [0, 1]
CHNL_FLW_DIR = 90 (90) [-360, 360]
CHNL_MARGIN = 5 (5) [1, 100]
CHNL_MAX_DEPTH = 40 (40) [10, 500]
CHNL_SCALE_CONST = 100 (100) [80, 120]
CHNL_SCALE_LOGNORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_LOGNORM_STDEV = 20 (20) [0, 40]
CHNL_SCALE_MAX = 120 (120) [100, 150]
CHNL_SCALE_MIN = 80 (80) [50, 100]
CHNL_SCALE_NORM_MEAN = 100 (100) [90, 110]
CHNL_SCALE_NORM_STDEV = 20 (20) [0, 40]
CHNL_WAVELENGTH = 4000 (4000) [100, 20000]
CHNL_WIDTH = 800 (800) [20, 4000]
COND_WELL_CLOSURE = 1.4 (1.4) [0.1, 10]
COND_WELL_RELAX = 20 (20) [0, 100]
DOMAIN_DX = 200 (200) [1, 1e+30]
DOMAIN_DY = 200 (200) [1, 1e+30]
DOMAIN_GRID_DIR = 0 (0) [-360, 360]
DOMAIN_OX = 0 (0) [-1e+30, 1e+30]
DOMAIN_OY = 0 (0) [-1e+30, 1e+30]
DOMAIN_SLOPE = 0.005 (0.005) [0.0001, 0.5]
DOMAIN_ZREF = 0 (0) [-1e+30, 1e+30]
EROD_COEF = 4e-08 (4e-08) [0, 2e-07]
EROD_OFFSET = 0 (0) [-1e+30, 1e+30]
GR_PB_OPTIM_ALPHA = 0.25 (0.25) [0.05, 0.95]
GR_THRESHOLD = 0.75 (0.75) [0.05, 0.95]
MASS_BAL_COEF = 0.8 (0.8) [0, 1e+30]
SED_LOAD_MEAN = 0.1 (0.1) [0.001, 0.5]
Real Sinuosity Computation Interval         = 1
Number of Points for Curvatures Calculation = 3
Maximum Well Influence Distance Factor      = 15
Avulsions Conditioned by Wells              = Yes
Conditioning New Algorithm                  = Yes
Fill Abandoned Channels                     = Yes
Bank Erodibility Coefficient                = Not used
Using Old Filling Algorithm                 = No
Scour Factor                                = Not used
Friction Coefficient                        = Not used

Simulation finished!
No description has been provided for this image
No description has been provided for this image
In [ ]: