identifier = -1
print_output('Simulation identifier : ' + str(identifier))
ramptime = 250000
t_start = -ramptime # Simulation start time relative to tidal_forcing
t_end = ramptime + t_start # Simulation duration in sec
Dt = 100. # Crank Nicolson timestep
t_export = 1000.0 # Export time if necessary
wd_alpha = 0.5 # Wetting and drying
#mu_manning = 0.02 # Bed friction
lat_coriolis = 51 # Coriolis calculation parameters
CG_2d = FunctionSpace(mesh2d, 'CG', 1)
lagoon_input = input_barrages.input_barrage("inputs/LagoonSpecs.dat")
bathymetry_2d, h_viscosity, mu_manning = tools.thetis_support_scripts.initialise_fields(
mesh2d, inputdir, outputdir, identifier, manning=True)
coriolis_2d = tools.thetis_support_scripts.coriolis(mesh2d, lat_coriolis)
with timed_stage('initialisation'):
# --- create solver ---
solver_obj = solver2d.FlowSolver2d(mesh2d, bathymetry_2d)
options = solver_obj.options
options.cfl_2d = 1.0
options.use_nonlinear_equations = True
options.simulation_export_time = t_export
options.simulation_end_time = ramptime
options.coriolis_frequency = coriolis_2d
options.output_directory = outputdir
# ctrl_ID = int(sys.argv[2])
#
#
#
# if ctrl_ID == 0 : operation = None
# elif ctrl_ID == 1 : operation = 'cardiff_ebb'
# elif ctrl_ID == 2 : operation = 'cardiff_two-way'
identifier = 0
operation = None
bathymetry_2d, h_viscosity, elev_init, uv_init, mu_manning = tools.thetis_support_scripts.initialise_fields(
mesh2d, inputdir, outputdir, identifier, manning=True)
coriolis_2d = tools.thetis_support_scripts.coriolis(mesh2d, lat_coriolis)
lagoon_swansea_input = input_barrages.input_barrage(
"../swansea_2018_copy/inputs/LagoonSpecs.dat")
lagoon_cardiff_input = input_barrages.input_barrage(
"inputs/LagoonSpecs.dat") # inputs/LagoonSpecs.dat
lagoon_status = pickle.load(
open(inputdir + "/barrage_status_" + str(identifier) + ".p",
"rb")) #should add one?
if operation != None:
adaptive_control = pickle.load(
open(inputdir + "/" + operation + ".p", "rb"))
holding_timing = pickle.load(
open(inputdir + "/holding_" + operation + ".p", "rb"))
else:
adaptive_control = None