Alaskan Winter Storm with Ice - Revision history

Wpringle at 23:43, 8 June 2020

2020-06-08T23:43:52Z

← Older revision		Revision as of 19:43, 8 June 2020
Line 1:		Line 1:
	This example tests ADCIRC version 55 (and beyond). It tests the simulation of the storm tides in a regional Alaska domain under astronomical and atmospheric forcing in November 2011 during a strong winter storm in the presence of sea (affecting the surface wind drag)<ref name=Brian>Joyce, B.R., Pringle, W.J., Wirasaet, D., Westerink, J.J., Van der Westhuysen, A.J., Grumbine, R., Feyen, J., 2019. High resolution modeling of western Alaskan tides and storm surge under varying sea ice conditions. Ocean Model. 141, 101421. doi:10.1016/j.ocemod.2019.101421</ref>. The results of interest are the global elevations, velocities and meteorology. The test finishes in about 5 minutes in serial ADCIRC for two weeks of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_alaska_ice-2d GitHub test suite].		This example tests ADCIRC version 55 (and beyond). It tests the simulation of the storm tides in a regional Alaska domain under astronomical and atmospheric forcing in November 2011 during a strong winter storm in the presence of sea ice (affecting the surface wind drag)<ref name=Brian>Joyce, B.R., Pringle, W.J., Wirasaet, D., Westerink, J.J., Van der Westhuysen, A.J., Grumbine, R., Feyen, J., 2019. High resolution modeling of western Alaskan tides and storm surge under varying sea ice conditions. Ocean Model. 141, 101421. doi:10.1016/j.ocemod.2019.101421</ref>. The results of interest are the global elevations, velocities and meteorology. The test finishes in about 5 minutes in serial ADCIRC for two weeks of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_alaska_ice-2d GitHub test suite].

	== Mesh ==		== Mesh ==

Wpringle: /* Options/Features Tested */

2020-06-08T22:46:33Z

Options/Features Tested

← Older revision		Revision as of 18:46, 8 June 2020
Line 16:		Line 16:
	*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).		*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).
	*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
	*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: ~~Change~~ "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).		*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Set "WarnElev", the warning elevation level, to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).
	*[[Fort.15_file_format#Namelists\|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true.		*[[Fort.15_file_format#Namelists\|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true.

	== References ==		== References ==
	<references />		<references />

Wpringle at 22:45, 8 June 2020

2020-06-08T22:45:49Z

← Older revision		Revision as of 18:45, 8 June 2020
Line 17:		Line 17:
	*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
	*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Change "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).		*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Change "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).
	*[[Fort.15_file_format#Namelists\|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true		*[[Fort.15_file_format#Namelists\|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true.

	== References ==		== References ==
	<references />		<references />

Wpringle: /* Options/Features Tested */

2020-06-08T22:45:38Z

Options/Features Tested

← Older revision		Revision as of 18:45, 8 June 2020
Line 17:		Line 17:
	*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
	*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Change "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).		*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Change "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).
			*[[Fort.15_file_format#Namelists\|&metControl namelist]]: Set "rhoAir", to 1.29193 (density of air at 0 deg C for 1013 mbar); set "WindDragLimit" equal to 0.0025; set "invertedBarometerOnElevationBoundary" to true (in Alaska extremely large-scale low pressure systems persist and cross over the open boundaries, so it is important to have the inverted barometer condition along the elevation specified boundary); set "outputWindDrag" to true

	== References ==		== References ==
	<references />		<references />

Wpringle: /* Options/Features Tested */

2020-06-08T22:40:20Z

Options/Features Tested

← Older revision		Revision as of 18:40, 8 June 2020
Line 8:		Line 8:
	*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).		*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).
	*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.		*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.
	*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient<ref name=Brian></ref><\.		*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient<ref name=Brian></ref>.
	*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.		*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.
	*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.		*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.

Wpringle at 22:40, 8 June 2020

2020-06-08T22:40:11Z

← Older revision		Revision as of 18:40, 8 June 2020
Line 1:		Line 1:
	This example tests ADCIRC version 55 (and beyond). It tests the simulation of the storm tides in a regional Alaska domain under astronomical and atmospheric forcing in November 2011 during a strong winter storm in the presence of sea (affecting the surface wind drag)<ref name=~~Foreman~~>Joyce, B.R., Pringle, W.J., Wirasaet, D., Westerink, J.J., Van der Westhuysen, A.J., Grumbine, R., Feyen, J., 2019. High resolution modeling of western Alaskan tides and storm surge under varying sea ice conditions. Ocean Model. 141, 101421. doi:10.1016/j.ocemod.2019.101421</ref>. The results of interest are the global elevations, velocities and meteorology. The test finishes in about 5 minutes in serial ADCIRC for two weeks of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_alaska_ice-2d GitHub test suite].		This example tests ADCIRC version 55 (and beyond). It tests the simulation of the storm tides in a regional Alaska domain under astronomical and atmospheric forcing in November 2011 during a strong winter storm in the presence of sea (affecting the surface wind drag)<ref name=Brian>Joyce, B.R., Pringle, W.J., Wirasaet, D., Westerink, J.J., Van der Westhuysen, A.J., Grumbine, R., Feyen, J., 2019. High resolution modeling of western Alaskan tides and storm surge under varying sea ice conditions. Ocean Model. 141, 101421. doi:10.1016/j.ocemod.2019.101421</ref>. The results of interest are the global elevations, velocities and meteorology. The test finishes in about 5 minutes in serial ADCIRC for two weeks of simulation. Find the test at the [https://github.com/adcirc/adcirc-cg-testsuite/tree/v55/adcirc/adcirc_alaska_ice-2d GitHub test suite].

	== Mesh ==		== Mesh ==
Line 8:		Line 8:
	*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).		*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).
	*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.		*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.
	*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient.		*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient<ref name=Brian></ref><\.
	*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.		*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.
	*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.		*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.

Wpringle: /* Options/Features Tested */

2020-06-08T22:39:17Z

Options/Features Tested

← Older revision		Revision as of 18:39, 8 June 2020
Line 16:		Line 16:
	*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).		*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).
	*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
	*[[~~fort~~.~~13_file~~#~~Internal_Tide_Energy_Conversion~~\|~~internal_tide_friction~~]]: ~~Namelist to change~~ the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).		*[[Fort.15_file_format#Namelists\|&WarnElevControl namelist]]: Change "WarnElev", the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).

	== References ==		== References ==
	<references />		<references />

Wpringle: /* Options/Features Tested */

2020-06-08T22:37:53Z

Options/Features Tested

← Older revision		Revision as of 18:37, 8 June 2020
Line 5:		Line 5:

	== Options/Features Tested ==		== Options/Features Tested ==
	*<code>[[ICS]]</code> = -~~22: Uses the Mercator~~ projection ~~with a coordinate rotation to remove the pole singularity (need to provide a [[fort.rotm]])~~.		*<code>[[ICS]]</code> = 20: Equal-Area cylindrical projection.
	*<code>[[IM]]</code> = ~~513113~~: Uses the ~~fully~~ implicit scheme for the gravity wave term (computational time step is 12 minutes).		*<code>[[IM]]</code> = 513111: Uses the implicit scheme for the linear component of the gravity wave term (computational time step is 4 minutes).
	*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 10 tidal constituents.		*<code>[[NTIP]]</code> = 2: Equilibrium tide + self-attraction and loading tide (read from a [[fort.24 file]]) forcing for 8 tidal constituents.
	*<code>[[NWS]]</code> = ~~-14~~: Reads from GRIB2 files that specify the global atmospheric forcing (6-hourly ~~CFS~~ reanalysis data) ~~in addition to OWI ASCII files that specify the 3~~-~~hourly atmospheric forcing in~~ the ~~Hurricane Katrina landfall region~~.		*<code>[[NWS]]</code> = 14014: Reads from GRIB2 files that specify the global atmospheric forcing and sea-ice concentration (6-hourly CFSv2 reanalysis data). Sea-ice concentration affects the wind drag coefficient.
	*<code>[[WTIMINC]]</code> = 21600, ~~10800~~: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the ~~OWI met~~ data (3 hours) ~~in seconds~~.		*<code>[[WTIMINC]]</code> = 21600, 21600: First value gives the temporal interval of the GRIB2 met data (6 hours), second value gives the temporal interval of the GRIB2 ice data (6 hours) - these should always be the same.
	*<code>[[A00, B00, C00]]</code> = 0.5, 0.5, 0: Ensures that the ~~fully~~ implicit scheme is stable with a large time step.		*<code>[[A00, B00, C00]]</code> = 0.4, 0.4, 0.2: Ensures that the implicit scheme is stable with a fairly large time step.
	*<code>[[ESLM]]</code> = -0.2: Enables the Smagorinsky turbulence closure with a coefficient of 0.2.		*<code>[[ESLM]]</code> = -0.2: Enables the Smagorinsky turbulence closure with a coefficient of 0.2.
	*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]].		*<code>[[NOUTGE]]</code> = 5: Outputs the global elevations into a netCDF4 [[fort.63 file]].
Line 16:		Line 16:
	*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).		*<code>[[NOUTGM]]</code> = 5: Outputs the global meteorology into a netCDF4 [[fort.73 file]] (pressure) and a netCDF4 [[fort.74 file]] (velocity).
	*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Spatially varying linear wave drag [[fort.13 file]] attribute accounting for energy conversion due to internal tide generation in the deep ocean.
	*[[fort.13_file#~~Quadratic_Friction_coefficient~~\|~~quadratic_friction_coefficient~~]]: ~~Spatially varying quadratic bottom friction~~ [~~[fort.13 file]~~] ~~attribute~~.		*[[fort.13_file#Internal_Tide_Energy_Conversion\|internal_tide_friction]]: Namelist to change the warning elevation level to 30-m (elevations reach beyond 20-m [default] but remain below 30-m).

	== References ==		== References ==
	<references />		<references />

Wpringle: /* Mesh */

2020-06-08T22:31:50Z

Mesh

← Older revision		Revision as of 18:31, 8 June 2020
Line 2:		Line 2:

	== Mesh ==		== Mesh ==
	The mesh ~~is a coarse representation of Alaska (Gulf of Alaska, Bering Sea, and Chukchi Sea) with a minimum resolution of 5 km, comprised of 15,876 vertices and 27,757 triangular elements~~ generated using the OceanMesh2D Alaska [https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/Examples/Example_8_AK.m Example_8_AK.m].		The mesh was generated using the OceanMesh2D Alaska [https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/Examples/Example_8_AK.m Example_8_AK.m]. The domain encompasses the Gulf of Alaska, Bering Sea, and Chukchi Sea with a minimum resolution of 5 km, comprised of 15,876 vertices and 27,757 triangular elements.

	== Options/Features Tested ==		== Options/Features Tested ==

Wpringle: /* Mesh */

2020-06-08T22:30:50Z

Mesh

← Older revision		Revision as of 18:30, 8 June 2020
Line 2:		Line 2:

	== Mesh ==		== Mesh ==
	The mesh is a coarse representation of ~~the spherical Earth~~ with minimum resolution of ~~approximately 50~~ km, comprised of 27,~~330~~ vertices and 50,~~859~~ triangular elements.		The mesh is a coarse representation of Alaska (Gulf of Alaska, Bering Sea, and Chukchi Sea) with a minimum resolution of 5 km, comprised of 15,876 vertices and 27,757 triangular elements generated using the OceanMesh2D Alaska [https://github.com/CHLNDDEV/OceanMesh2D/blob/Projection/Examples/Example_8_AK.m Example_8_AK.m].

	== Options/Features Tested ==		== Options/Features Tested ==