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!
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!
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!
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!
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!
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!
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!
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!
In [ ]: