def main():
# fetch shinnecock inlet test data
if not FORT14.is_file():
url = "https://www.dropbox.com/s/1wk91r67cacf132/"
url += "NetCDF_shinnecock_inlet.tar.bz2?dl=1"
g = urllib.request.urlopen(url)
tmpfile = tempfile.NamedTemporaryFile()
with open(tmpfile.name, 'b+w') as f:
f.write(g.read())
with tarfile.open(tmpfile.name, "r:bz2") as tar:
tar.extractall(PARENT / "data/NetCDF_Shinnecock_Inlet/")
# open mesh file
mesh = AdcircMesh.open(FORT14, crs=4326)
# init tidal forcing and setup requests
tidal_forcing = Tides()
tidal_forcing.use_constituent('M2')
tidal_forcing.use_constituent('N2')
tidal_forcing.use_constituent('S2')
tidal_forcing.use_constituent('K1')
tidal_forcing.use_constituent('O1')
mesh.add_forcing(tidal_forcing)
# set simulation dates
start_date = datetime(2015, 12, 14)
end_date = start_date + timedelta(days=5)
# instantiate AdcircRun object.
driver = AdcircRun(mesh, start_date, end_date)
# request outputs
driver.set_elevation_surface_output(sampling_rate=timedelta(minutes=30))
driver.set_velocity_surface_output(sampling_rate=timedelta(minutes=30))
# override the AdcircPy defaults so that the fort.15
# matches the original Shinnecock test case options
driver.timestep = 6.
driver.DRAMP = 2.
driver.TOUTGE = 3.8
driver.TOUTGV = 3.8
driver.smagorinsky = False
driver.horizontal_mixing_coefficient = 5.
driver.gwce_solution_scheme = 'semi-implicit-legacy'
# run parallel ADCIRC if binary is installed
if shutil.which('padcirc') is not None:
driver.run(PARENT / "outputs/example_1", overwrite=True)
# run serial ADCIRC if binary is installed
elif shutil.which('adcirc') is not None:
driver.run(PARENT / "outputs/example_1", overwrite=True, nproc=1)
# binaries are not installed, write to disk and exit
else:
msg = 'ADCIRC binaries were not found in PATH. ADCIRC will not run. '
msg += "Writing files to disk..."
warnings.warn(msg)
driver.write(PARENT / "outputs/example_1", overwrite=True)
tidal_forcing.use_constituent('N2')
tidal_forcing.use_constituent('S2')
tidal_forcing.use_constituent('K1')
tidal_forcing.use_constituent('O1')
mesh.add_forcing(tidal_forcing)
# set simulation dates
duration = timedelta(days=5)
start_date = datetime(2015, 12, 14)
end_date = start_date + duration
# instantiate driver object
driver = AdcircRun(mesh, start_date, end_date)
# request outputs
driver.set_elevation_surface_output(sampling_rate=timedelta(minutes=30))
driver.set_velocity_surface_output(sampling_rate=timedelta(minutes=30))
# override default options so the resulting `fort.15` matches the original Shinnecock test case options
driver.timestep = 6.0
driver.DRAMP = 2.0
driver.TOUTGE = 3.8
driver.TOUTGV = 3.8
driver.smagorinsky = False
driver.horizontal_mixing_coefficient = 5.0
driver.gwce_solution_scheme = 'semi-implicit-legacy'
if shutil.which('padcirc') is not None:
driver.run(OUTPUT_DIRECTORY, overwrite=True)
elif shutil.which('adcirc') is not None:
driver.run(OUTPUT_DIRECTORY, overwrite=True, nproc=1)
def main():
# fetch shinnecock inlet test data
if not FORT14.is_file():
url = "https://www.dropbox.com/s/1wk91r67cacf132/"
url += "NetCDF_shinnecock_inlet.tar.bz2?dl=1"
g = urllib.request.urlopen(url)
tmpfile = tempfile.NamedTemporaryFile()
with open(tmpfile.name, 'b+w') as f:
f.write(g.read())
with tarfile.open(tmpfile.name, "r:bz2") as tar:
tar.extractall(PARENT / "data/NetCDF_Shinnecock_Inlet/")
# open mesh file
mesh = AdcircMesh.open(FORT14, crs=4326)
# let's generate the tau0 factor
mesh.generate_tau0()
# let's also add a mannings to the domain (constant for this example)
mesh.mannings_n_at_sea_floor = np.full(mesh.values.shape, 0.025)
# init tidal forcing and setup requests
tidal_forcing = Tides()
tidal_forcing.use_constituent('M2')
tidal_forcing.use_constituent('N2')
tidal_forcing.use_constituent('S2')
tidal_forcing.use_constituent('K1')
tidal_forcing.use_constituent('O1')
mesh.add_forcing(tidal_forcing)
# set simulation dates
spinup_time = timedelta(days=2)
start_date = datetime(2015, 12, 14) + spinup_time
end_date = start_date + timedelta(days=3)
# instantiate AdcircRun object.
driver = AdcircRun(
mesh,
start_date,
end_date,
spinup_time,
)
# request outputs
driver.set_elevation_surface_output(sampling_rate=timedelta(minutes=30))
driver.set_velocity_surface_output(sampling_rate=timedelta(minutes=30))
# override defaults options
driver.timestep = 4.0
# run parallel ADCIRC if binary is installed
if shutil.which('padcirc') is not None:
driver.run(PARENT / "outputs/example_2", overwrite=True)
# run serial ADCIRC if binary is installed
elif shutil.which('adcirc') is not None:
driver.run(PARENT / "outputs/example_2", overwrite=True, nproc=1)
# binaries are not installed, write to disk and exit
else:
msg = 'ADCIRC binaries were not found in PATH. ADCIRC will not run. '
msg += "Writing files to disk..."
warnings.warn(msg)
driver.write(PARENT / "outputs/example_2", overwrite=True)