Skip to content

Commit 0753743

Browse files
SierdSierd
authored
Merge branch 'add_profiler' into main
2 parents d91a74a + 0d54f7d commit 0753743

File tree

3 files changed

+370
-0
lines changed

3 files changed

+370
-0
lines changed

‎.devcontainer/Dockerfile‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,7 @@
11
FROM continuumio/miniconda3
22

33
# Set the working directory
4+
45
# WORKDIR /workspace
56

67
# Install Python
@@ -10,6 +11,7 @@ RUN conda install python -y
1011
# COPY . .
1112

1213
# Install the necessary packages
14+
1315
# RUN pip install .
1416

1517
# Clean up

‎aeolis/examples/vanWesten2024/blowout/blowout_short.txt‎

Lines changed: 155 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,155 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
%% AeoLiS model configuration %%
3+
%% Date: 2020-04-22 12:26:11 %%
4+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5+
6+
%% -------------------- [Timing] ----------------------------- %%
7+
dt = 3600.000 % [s] Timeinterval between timesteps
8+
tstart = 36979200 % [s] Starting time of the simulation (03-04-2013)
9+
tstop = 37000000 % 66674520 % [s] End time of the simulation
10+
restart = None % [s] Interval for which to write restart files
11+
refdate = 2012-01-01 00:00 % [-] Reference datetime in netCDF output
12+
13+
%% -------------------- [Grid files *.grd] ------------------- %%
14+
xgrid_file = x_NPZK.grd % Filename of ASCII file with x-coordinates of grid cells
15+
ygrid_file = y_NPZK.grd % Filename of ASCII file with y-coordinates of grid cells
16+
bed_file = z_NPZK.grd % Filename of ASCII file with bed level heights of grid cells
17+
ne_file = zne_NPZK.grd % Filename of ASCII file with non-erodible layer
18+
veg_file = veg_NPZK.grd % Filename of ASCII file with initial vegetation density
19+
wave_mask = wave_mask_NPZK.grd % Filename of ASCII file with mask for wave height
20+
tide_mask = tide_mask_NPZK.grd % Filename of ASCII file with mask for tidal elevation
21+
vver_mask = vver_mask_NPZK.grd % Filename of ASCII file with mask for vertical growth (V_ver)
22+
23+
%% -------------------- [Time series *.txt] ------------------ %%
24+
tide_file = tide_hvh_2012_2022_ref20120101.txt % Filename of ASCII file with time series of water levels
25+
wave_file = waves_euro_2012_2022_ref20120101.txt % Filename of ASCII file with time series of wave heights (and periods)
26+
wind_file = wind_NPZK_2012_2023.txt % Filename of ASCII file with time series of wind velocity and direction
27+
28+
%% -------------------- [Output] ----------------------------- %%
29+
output_times = 604800.000 % [s] Timeinterval [s] between *.nc output file
30+
output_vars = zb zs zne ustar ustarn ustars ustars_avg ustarn_avg u us un SWL TWL rhoveg moist_avg % [-] List of output variables for in *.nc output file
31+
visualization = T
32+
33+
%% -------------------- [Flags Processes] -------------------- %%
34+
process_wind = T % [T/F] Enable the process of wind
35+
process_shear = T % [T/F] Enable the process of wind shear
36+
process_tide = T % [T/F] Enable the process of tides
37+
process_wave = T % [T/F] Enable the process of waves
38+
process_runup = T % [T/F] Enable the process of wave runup
39+
process_moist = T % [T/F] Enable the process of moist
40+
process_mixtoplayer = T % [T/F] Enable the process of mixing
41+
process_threshold = T % [T/F] Enable the process of threshold
42+
process_transport = T % [T/F] Enable the process of transport
43+
process_bedupdate = T % [T/F] Enable the process of bed updating
44+
process_avalanche = T % [T/F] Enable the process of avalanching
45+
process_separation = F % [T/F] Enable the including of separation bubble
46+
process_vegetation = T % [T/F] Enable the process of vegetation
47+
process_wet_bed_reset = T
48+
49+
%% -------------------- [Flags Threshold] -------------------- %%
50+
th_grainsize = T % [T/F] Enable wind velocity threshold based on grainsize
51+
th_bedslope = F % [T/F] Enable wind velocity threshold based on bedslope
52+
th_moisture = T % [T/F] Enable wind velocity threshold based on moisture
53+
th_sheltering = F % [T/F] Enable wind velocity threshold based on roughness
54+
th_nelayer = T % [T/F] Enable wind velocity threshold based on a non-erodible layer
55+
56+
%% -------------------- [General physics] -------------------- %%
57+
g = 9.81 % [m/s^2] Gravitational constant
58+
v = 0.000015 % [m^2/s] Air viscosity
59+
rhoa = 1.225 % [kg/m^3] Air density
60+
rhog = 2650.0 % [kg/m^3] Grain density
61+
rhow = 1025.0 % [kg/m^3] Water density
62+
porosity = 0.4 % [-] Sediment porosity
63+
cpair = 0.0010035 % [MJ/kg/oC] Specific heat capacity air
64+
65+
%% -------------------- [Sediment] --------------------------- %%
66+
nlayers = 1 % [-] Number of bed layers
67+
layer_thickness = 0.05000 % [m] Thickness of bed layers
68+
nfractions = 1 % [-] Number of sediment fractions
69+
grain_dist = 1.0 % [-] Initial distribution of sediment fractions
70+
grain_size = 0.00025 % [m] Average grain size of each sediment fraction
71+
72+
%% -------------------- [Wind and shear] --------------------- %%
73+
wind_convention = nautical % [-] Convention used for the wind direction in the input files (cartesian or nautical)
74+
alfa = 0 % [deg] Real-world grid cell orientation wrt the North (clockwise)
75+
k = 0.003 % [m] Bed roughness
76+
z = 10. % [m] Measurement height of wind velocity
77+
kappa = 0.41 % [-] Von Kármán constant
78+
h = None % [m] Representative height of saltation layer
79+
L = 100. % [m] Typical length scale of dune feature (perturbation)
80+
l = 1. % [m] Inner layer height (perturbation)
81+
m = 1. % [-] Factor to account for difference between average and maximum shear stress
82+
83+
%% -------------------- [Transport] -------------------------- %%
84+
bi = 0.050000 % [-] Bed interaction factor
85+
method_transport = bagnold % [-] Name of method to compute equilibrium sediment transport rate
86+
method_grainspeed = duran_full % [-] Name of method to assume/compute grainspeed
87+
Aa = 0.085 % [-] Constant in formulation for wind velocity threshold based on grain size
88+
Cb = 1.5 % [-] Constant in bagnold formulation for equilibrium sediment concentration
89+
Ck = 2.78 % [-] Constant in kawamura formulation for equilibrium sediment concentration
90+
Cl = 6.7 % [-] Constant in lettau formulation for equilibrium sediment concentration
91+
Cdk = 5. % [-] Constant in DK formulation for equilibrium sediment concentration
92+
93+
%% -------------------- [Solver] ----------------------------- %%
94+
T = 1. % [s] Adaptation time scale in advection equation
95+
solver = steadystate % [-] Numerical solver of advection scheme
96+
CFL = 1. % [-] CFL number to determine time step in explicit scheme
97+
accfac = 1. % [-] Numerical acceleration factor
98+
scheme = euler_backward % [-] Name of numerical scheme (euler_forward, euler_backward or crank_nicolson)
99+
max_error = 0.000001 % [-] Maximum error at which to quit iterative solution in implicit numerical schemes
100+
max_iter = 1000 % [-] Maximum number of iterations at which to quit iterative solution in implicit numerical schemes
101+
% max_bedlevel_change = 0.01 % [m] Maximum bedlevel change after one timestep. Next timestep dt will be modified (use 999. if not used)
102+
103+
%% -------------------- [Boundary conditions] ---------------- %%
104+
boundary_onshore = constant % [-] Name of onshore boundary conditions (flux, constant, uniform, gradient)
105+
boundary_lateral = circular % [-] Name of lateral boundary conditions (circular, constant ==noflux)
106+
boundary_offshore = constant % [-] Name of offshore boundary conditions (flux, constant, uniform, gradient)
107+
offshore_flux = 0. % [-] Factor to determine offshore boundary flux as a function of Q0 (= 1 for saturated flux , = 0 for noflux)
108+
constant_offshore_flux = 0. % [kg/m/s] Constant input flux at offshore boundary
109+
onshore_flux = 0. % [-] Factor to determine onshore boundary flux as a function of Q0 (= 1 for saturated flux , = 0 for noflux)
110+
constant_onshore_flux = 0. % [kg/m/s] Constant input flux at offshore boundary
111+
lateral_flux = 0. % [-] Factor to determine lateral boundary flux as a function of Q0 (= 1 for saturated flux , = 0 for noflux)
112+
sedimentinput = 0. % [-] Constant boundary sediment influx (only used in solve_pieter)
113+
114+
%% -------------------- [Rotating shear-grid] ---------------- %%
115+
dx = 2.0 % [m] Size of dimension of the computational grid that is generated for a rotating shear computation in x-direction
116+
dy = 2.0 % [m] Size of dimension of the computational grid that is generated for a rotating shear computation in x-direction
117+
buffer_width = 5.0 % [m] Width of the bufferzone around the rotational grid for wind perturbation
118+
119+
%% -------------------- [Vegetation] ------------------------- %%
120+
sigma = 4.2 % [-] Ratio between basal area and frontal area of roughness elements
121+
beta = 130.000000 % [-] Ratio between drag coefficient of roughness elements and bare surface
122+
gamma_vegshear = 16. % [-] Roughness factor for the shear stress reduction by vegetation
123+
avg_time = 86400. % [s] Indication of the time period over which the bed level change is averaged for vegetation growth
124+
dzb_interval = 86400. % [s] Interval used for calcuation of vegetation growth
125+
hveg_max = 1. % [m] Max height of vegetation
126+
dzb_opt = 0. % [m/year] Sediment burial for optimal growth
127+
V_ver = 3. % [m/year] Vertical growth
128+
V_lat = 0. % [m/year] Lateral growth
129+
germinate = 0. % [1/year] Possibility of germination per year
130+
lateral = 0. % [1/year] Posibility of lateral expension per year
131+
veg_gamma = 0.05 % [-] Constant on influence of sediment burial
132+
veg_sigma = 0. % [-] Sigma in gaussian distrubtion of vegetation cover filter
133+
vegshear_type = raupach % [-] Choose the Raupach grid based solver (1D or 2D) or the Okin approach (1D only)
134+
veg_min_elevation = -10
135+
136+
%% -------------------- [Seperation] ------------------------- %%
137+
c_b = 0.2 % [-] Slope at the leeside of the separation bubble # c = 0.2 according to Durán 2010 (Sauermann 2001: c = 0.25 for 14 degrees)
138+
mu_b = 20. % [deg] Minimum required slope for the start of flow separation
139+
sep_filter_iterations = 0 % [-] Number of filter iterations on the separation bubble (0 = no filtering)
140+
zsep_y_filter = F % [T/F] Boolean for turning on/off the filtering of the separation bubble height in y-direction
141+
142+
%% -------------------- [Soil moisture] ---------------------- %%
143+
Tdry = 5400.000 % [s] Adaptation time scale for soil drying
144+
eps = 0.2 % [m] Minimum water depth to consider a cell "flooded"
145+
method_moist = belly_johnson % [-] Name of method to compute wind velocity threshold based on soil moisture content
146+
147+
%% -------------------- [Waves] ------------------------------ %%
148+
Tbedreset = 86400. % [s] Adaptation time for resetting morphology in the Swash-zone
149+
xi = 0.3 % [-] Surf similarity parameter
150+
facDOD = 0.5 % [-] Ratio between depth of disturbance and local wave height
151+
152+
%% -------------------- [Avalanching] ------------------------ %%
153+
theta_dyn = 44. % [degrees] Initial Dynamic angle of repose, critical dynamic slope for avalanching
154+
theta_stat = 45. % [degrees] Initial Static angle of repose, critical static slope for avalanching
155+
max_iter_ava = 100 % [-] Maximum number of iterations at which to quit iterative solution in avalanching calculation

0 commit comments

Comments
 (0)