ADCIRCWiki - User contributions [en]

Typical ADCIRC Parameter Selections

2019-06-05T13:00:58Z

JasonFleming: clarified SWAN columns

{| class="wikitable" border="1"
|-
!Region
!Mesh
!np
!GWCE Formulation
!Drag Law
!Drag Limit
!FFACTOR
!Advective Terms
!Steric
!Vert. Datum
!ESL
!ESL Node Count
!time step
!h0
!velmin
!tidal constituents
!convcr
!SWAN time step
!SWAN MXITNS
!SWAN NPNTS
|-
|Louisiana
|LA_v17a-WithUpperAtch_chk.grd
|1593521
|implicit
|Powell
|0.002
|0.0
|off
|0.228184
|navd
|0.05
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|Texas
|tx2008_r35h.grd
|3352598
|implicit
|Garratt
|0.002
|0.0
|off
|0.276300
|navd?
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-07
|1200
|20
|95
|-
|National
|hsofs.14
|1813443
|implicit or explicit
|Garratt
|0.0028
|0.0025
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|National (Drag Experiment)
|hsofs.14
|1813443
|explicit
|Garratt
|0.002
|0.0
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|MS/AL and FL Panhandle
|NGOM_RT_v18j_chk.grd
|2051346
|implicit
|Powell
|0.002
|0.0
|off
|0.230000
|navd
|0.02
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|North Carolina (low res)
|nc_inundation_v6d_rivers_msl.grd
|295328
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|0.5
|0.02
|0.02
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|North Carolina (high res)
|nc_inundation_v9.99a_w_rivers.grd
|624782
|explicit
|Garratt
|0.0028
|0.003
|on
|0.0
|msl
|none
|none
|0.5
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|Delmarva
|FEMA_R3_20110303_MSL.grd
|1875689
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|NY/NJ
|FEMA_R2_norivers_gcs_mNAVD.grd
|604790
|implicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-08
|1200
|10
|95
|-
|New England
|NAC2014_R01_ClosedRivers.grd
|3110470
|implicit
|Powell
|0.002
|0.0
|off
|0.10900
|msl
|none
|none
|?
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|n/a
|n/a
|n/a
|}

Typical ADCIRC Parameter Selections

2019-06-05T12:59:51Z

JasonFleming: fixed table formatting

{| class="wikitable" border="1"
|-
!Region
!Mesh
!np
!GWCE Formulation
!Drag Law
!Drag Limit
!FFACTOR
!Advective Terms
!Steric
!Vert. Datum
!ESL
!ESL Node Count
!time step
!h0
!velmin
!tidal constituents
!convcr
!time step
!MXITNS
!NPNTS
|-
|Louisiana
|LA_v17a-WithUpperAtch_chk.grd
|1593521
|implicit
|Powell
|0.002
|0.0
|off
|0.228184
|navd
|0.05
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|Texas
|tx2008_r35h.grd
|3352598
|implicit
|Garratt
|0.002
|0.0
|off
|0.276300
|navd?
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-07
|1200
|20
|95
|-
|National
|hsofs.14
|1813443
|implicit or explicit
|Garratt
|0.0028
|0.0025
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|National (Drag Experiment)
|hsofs.14
|1813443
|explicit
|Garratt
|0.002
|0.0
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|MS/AL and FL Panhandle
|NGOM_RT_v18j_chk.grd
|2051346
|implicit
|Powell
|0.002
|0.0
|off
|0.230000
|navd
|0.02
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|North Carolina (low res)
|nc_inundation_v6d_rivers_msl.grd
|295328
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|0.5
|0.02
|0.02
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|North Carolina (high res)
|nc_inundation_v9.99a_w_rivers.grd
|624782
|explicit
|Garratt
|0.0028
|0.003
|on
|0.0
|msl
|none
|none
|0.5
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|Delmarva
|FEMA_R3_20110303_MSL.grd
|1875689
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|NY/NJ
|FEMA_R2_norivers_gcs_mNAVD.grd
|604790
|implicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-08
|1200
|10
|95
|-
|New England
|NAC2014_R01_ClosedRivers.grd
|3110470
|implicit
|Powell
|0.002
|0.0
|off
|0.10900
|msl
|none
|none
|?
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|n/a
|n/a
|n/a
|}

Typical ADCIRC Parameter Selections

2019-06-05T12:56:00Z

JasonFleming: initial parameters table import

{| class="wikitable" border="1"
|-
! header 1
! header 2
! header 3
|-
| row 1, cell 1
| row 1, cell 2
| row 1, cell 3
|-
| row 2, cell 1
| row 2, cell 2
| row 2, cell 3
|}

{|
|Region
|Mesh
|np
|GWCE Formulation
|Drag Law
|Drag Limit
|FFACTOR
|Advective Terms
|Steric
|Vert. Datum
|ESL
|ESL Node Count
|time step
|h0
|velmin
|tidal constituents
|convcr
|time step
|MXITNS
|NPNTS
|-
|Louisiana
|LA_v17a-WithUpperAtch_chk.grd
|1593521
|implicit
|Powell
|0.002
|0.0
|off
|0.228184
|navd
|0.05
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|Texas
|tx2008_r35h.grd
|3352598
|implicit
|Garratt
|0.002
|0.0
|off
|0.276300
|navd?
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-07
|1200
|20
|95
|-
|National
|hsofs.14
|1813443
|implicit or explicit
|Garratt
|0.0028
|0.0025
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|National (Drag Experiment)
|hsofs.14
|1813443
|explicit
|Garratt
|0.002
|0.0
|off
|none
|msl
|none
|none
|2.0
|0.05
|0.05
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|MS/AL and FL Panhandle
|NGOM_RT_v18j_chk.grd
|2051346
|implicit
|Powell
|0.002
|0.0
|off
|0.230000
|navd
|0.02
|all
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|North Carolina (low res)
|nc_inundation_v6d_rivers_msl.grd
|295328
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|0.5
|0.02
|0.02
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|North Carolina (high res)
|nc_inundation_v9.99a_w_rivers.grd
|624782
|explicit
|Garratt
|0.0028
|0.003
|on
|0.0
|msl
|none
|none
|0.5
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-10
|1200
|20
|95
|-
|Delmarva
|FEMA_R3_20110303_MSL.grd
|1875689
|explicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|1200
|20
|95
|-
|NY/NJ
|FEMA_R2_norivers_gcs_mNAVD.grd
|604790
|implicit
|Garratt
|0.0035
|0.003
|on
|none
|msl
|none
|none
|1.0
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,q1
|1.00E-08
|1200
|10
|95
|-
|New England
|NAC2014_R01_ClosedRivers.grd
|3110470
|implicit
|Powell
|0.002
|0.0
|off
|0.10900
|msl
|none
|none
|?
|0.1
|0.01
|m2,s2,n2,k1,k2,o1,p1,q1
|1.00E-07
|n/a
|n/a
|n/a
|}

Examples and Test Cases

2019-06-05T12:51:40Z

JasonFleming:

Example 1: [https://adcirc.org/home/documentation/example-problems/quarter-annular-harbor-with-tidal-forcing-example/ Quarter Annular Harbor with Tidal Forcing (2D and 3D)]

Example 2: [https://adcirc.org/home/documentation/example-problems/shinnecock-inlet-ny-with-tidal-forcing-example Shinnecock Inlet with Tidal Forcing (2D)]

Example 3: [https://adcirc.org/home/documentation/example-problems/beaufort-inlet-nc-example/ Beaufort Inlet (2D)]

Example 4: [https://adcirc.org/home/documentation/example-problems/idealized-inlet-test-cases/ Idealized Inlet (2D)]

Example 5: [https://adcirc.org/home/documentation/example-problems/internal-overflow-boundaries-example/ Internal Overflow Boundary (2D)]

Example 6: [https://adcirc.org/home/documentation/example-problems/apes-wind-run-example/ APES Wind Run (NWS = 3) (2D)]

Example 7: [https://adcirc.org/home/documentation/example-problems/hurricane-isabel-example/ Hurricane Isabel Wind Run (NWS = 4) (2D)]

Example 8: [https://adcirc.org/home/documentation/example-problems/large-grid-tides-example/ Large Grid (~300,000 nodes) Tides Only]

Example 9: [https://adcirc.org/home/documentation/example-problems/katrina-run-2015-nws-20-example/ Hurricane Katrina (NWS=20)]

Example 10: [[Typical ADCIRC Parameter Selections]]

Examples and Test Cases

2019-06-05T12:49:01Z

JasonFleming: adding adcirc parameters table

Dynamic water level correction

2019-03-22T15:42:25Z

JasonFleming:

Improvements in forecasting and hindcasting of tropical cyclones have led to much more accurate meteorological data associated with these storms. These improvements, together with more detailed representations of the physics and greater model resolution, have improved storm surge model skill. This has led to a state in which storm surge forecast and hindcast errors are often dominated by what we refer to as unresolved drivers, i.e., physical processes that are not explicitly included in storm surge models. Common examples include baroclinic processes, major oceanic currents, precipitation, steric fluctuations and far-field atmospheric forcing. These phenomena typically have longer timescales than the storm surge itself and therefore often show up as biases in water level during the period leading up to, encompassing and following the surge event. [From T.G.Asher]

[include reference to Taylor's paper in references section at end]

AMS 2017 Presentation from Luettich, et al:

Representing Low Frequency, Spatially Varying Waterlevel Anomalies in Storm Surge Forecasting

https://ams.confex.com/ams/97Annual/videogateway.cgi/id/37411?recordingid=37411&uniqueid=Paper316033&entry_password=655934

[include link to Rick's powerpoint slides]
[include reference to Rick's powerpoint slides in references section at end]

In order to accommodate these unresolved drivers, a water level offset ingestion feature has been added to ADCIRC.

==Controlling Water Level Offset==

The water level offset feature is triggered by the presence of the &offsetControl namelist at the bottom of the fort.15 file. Here is an example of how this line is used:

&offsetControl offsetFileName="hsofs_offset_test_ft.dat" offsetMultiplier=0.3048, offsetRampStart=0.0, offsetRampEnd=259200.0 offsetRampReferenceTime="hotstart", offsetSkipSnaps=0 /

The code looks for the water level offset file in the same directory as the other fulldomain files, both in serial and in parallel. In parallel, each subdomain does domain decomposition of the offset data on the fly, so the only change to adcprep is the read/write of the new namelist to the subdomain fort.15 files.

The multiplier value allows us to easily change the entire dataset by single factor. It may be needed in a case where the data were provided in ft but adcirc requires meters, for example.

The offsetRampStart and offsetRampEnd parameter are the time in seconds for a linear ramp up period. The offsetRampReferenceTime parameter controls whether the ramp times are relative to the hotstart time or the coldstart time. Prior to the start of the ramp, the offset values are multiplied by 0.0, eliminating the offset until the ramp period starts.

The input parameter offsetSkipSnaps is used to skip the specified number of initial data sets in the offsets data file. It can also be used to insert blank snaps at the beginning of the run if it is negative.

==Water Level Offset Input File==

The water level offset input file starts with the following three lines:

# test file for 1ft offset on hsofs mesh during matthew
14400.0 # time increment between datasets -- read but ignored if constant offset
0.0 # default offset value

The first line is a comment line. The second line is the time interval in seconds between successive datasets. The third line is the default offset value; this allows us to create offset input files in a "sparse" format where only the non-default values are provided.

If there is only one dataset in the offset file, then the values are treated as temporally constant. If there is more than one dataset, ADCIRC linearly interpolates between them in time. If ADCIRC runs out of offset data, it continues with the last dataset as a set of constant values for the remainder of the run.

Because each offset dataset may have a different number of values, the datasets are separated by a hash tag ('#') in the 2nd column. When ADCIRC sees this, it knows it has reached the end of the dataset it is currently reading.

==Output==

If a water level offset input file is specified, then the time interpolated water level offset values will be written to a file called offset.63 on the same schedule as output to the fort.63. Currently, the offset.63 file can only be produced in ASCII format. This file can then be used as a diagnostic to check that the offset values have been applied as expected.

==Compatibility==

Although the offset is implemented as a pseudo barometric pressure, it is entirely independent of any/all meteorological forcing and meteorological data. The values do not show up in the fort.71 (barometric pressure stations) or fort.73 (fulldomain barometric pressure field) output files. The offset can be used with any meteorological forcing data, or with no meteorological forcing data at all.

==Example==

An ADCIRC model is coldstarted with a tide only run that uses a 10 day ramp, and runs for 15 days. How should the water level offset be configured?

In this case I suggest that you apply your bias correction after reaching full strength tidal forcing (i.e., after the 10 day tidal ramp). So to have your bias correction ramp completed at the end of the 15 day tidal spinup, I would start the ramp period at 14.5 days and end it at 15.0 days. That way you'll be able to clearly see the correction occurring when you look at the output data.

In summary, I suggest you set your ramp period as follows:

offsetRampStart =1252800.0, offsetRampEnd = 1296000.0

For the bias correction data file for this 15.0 day tidal spinup case, you'll only need a single data set, which will be treated as a constant value because there is only one dataset, so the format of the file will be for example

# bias correction values test for GFS 20161002 cycle 00Z
99999.9 # time increment in seconds; not used in this case
0.0 # default nodal value applied to any node not specified below
234 0.5 <-- node number and final bias correction value
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the dataset

After the spinup, in preparation for a forecast run, when you want to migrate from one set of bias correction values to another over a 6 hour time period, your offsetControl line would be something like this:

&offsetControl offsetFileName= "offset_migration_from_00Z_to_06Z.dat", offsetMultiplier = 1.0, offsetRampStart = 1252800.0, offsetRampEnd = 1296000.0, offsetRampReferenceTime = "coldstart", offsetSkipSnaps = 0 /

and the corresponding file would look like something like the following:

# bias correction values for GFS 20161002 from cycle 00Z to cycle 06Z
21600.0 # 6 hour time increment in seconds
0.0 # default nodal value applied to any node not specified below
234 0.5 <--- values match those at the end of the tidal spinup
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the first dataset
512 0.3 <-- new values for updated bias correction at 06Z, not necessarily defined at the same list of nodes
1001 -0.5
2346 0.74
4000 -0.1
# <-- hash tag in second column indicating the end of the second dataset

At the end of the 15.0 day tidal spinup in this example the bias correction at node 234 would be 0.5. Three hours later (after the hotstart) it would be 0.25. Six hours after the hotstart it would be zero.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 1001 would be 0.0. Three hours later (after the hotstart) it would be -0.25. Six hours after the hotstart it would be -0.5.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 4000 would be 0.1. Three hours later (after the hotstart) it would be 0.0. Six hours after the hotstart it would be -0.1.

==Frequently Asked Questions==

Question: When writing the &offsetControl namelist, should I write the parameters on a single line, comma separated?

Answer: I always put fortran namelist parameters on a single line and separate the assignments with commas, but other formatting may be possible depending on the compiler you are using. I've never had issues with the way I format these fortran namelist input lines, but I've never tried any other way so I'll defer to your compiler's documentation on that.

Question: Do I need to ramp the offset, i.e., is it wise to do so for some reason?

Answer: We haven't tried a sudden application of the water level offset to the solution so we don't yet know how that would behave. Our shortest ramp up period is 12 hours but you might get away with 6 hours or even less.

Question: Does the application of a dynamic water level offset surface affect the velocity solution (not just the elevation solution as it is primarily intended)?

Answer: Yes, in a physically realistic manner. As an example, when applying a water level offset to a back bay area with a nearby tidal inlet, the water level change in the back bay must be reflected in the velocity in the inlet in order to maintain continuity (conservation of mass). In the case of an initial application of a temporally steady water level offset, this velocity effect will be a one-time occurrence. On the other hand, for a time varying water level offset surface, the continuous adjustments to the water level will be continuously reflected in the velocity solution in physically realistic ways.

==References==

[Ricks 2017 AMS ppt]
[Taylor's manuscript]

Dynamic water level correction

2019-03-22T15:41:50Z

JasonFleming:

Improvements in forecasting and hindcasting of tropical cyclones have led to much more accurate meteorological data associated with these storms. These improvements, together with more detailed representations of the physics and greater model resolution, have improved storm surge model skill. This has led to a state in which storm surge forecast and hindcast errors are often dominated by what we refer to as unresolved drivers, i.e., physical processes that are not explicitly included in storm surge models. Common examples include baroclinic processes, major oceanic currents, precipitation, steric fluctuations and far-field atmospheric forcing. These phenomena typically have longer timescales than the storm surge itself and therefore often show up as biases in water level during the period leading up to, encompassing and following the surge event. [From T.G.Asher]

[include reference to Taylor's paper in references section at end]

AMS 2017 Presentation from Luettich, et al:

Representing Low Frequency, Spatially Varying Waterlevel Anomalies in Storm Surge Forecasting

https://ams.confex.com/ams/97Annual/videogateway.cgi/id/37411?recordingid=37411&uniqueid=Paper316033&entry_password=655934

[include link to Rick's powerpoint slides]
[include reference to Rick's powerpoint slides in references section at end]

In order to accommodate these unresolved drivers, a water level offset ingestion feature has been added to ADCIRC.

==Controlling Water Level Offset==

The water level offset feature is triggered by the presence of the &offsetControl namelist at the bottom of the fort.15 file. Here is an example of how this line is used:

#&offsetControl offsetFileName="hsofs_offset_test_ft.dat" offsetMultiplier=0.3048, offsetRampStart=0.0, offsetRampEnd=259200.0 offsetRampReferenceTime="hotstart", offsetSkipSnaps=0 /

The code looks for the water level offset file in the same directory as the other fulldomain files, both in serial and in parallel. In parallel, each subdomain does domain decomposition of the offset data on the fly, so the only change to adcprep is the read/write of the new namelist to the subdomain fort.15 files.

The multiplier value allows us to easily change the entire dataset by single factor. It may be needed in a case where the data were provided in ft but adcirc requires meters, for example.

The offsetRampStart and offsetRampEnd parameter are the time in seconds for a linear ramp up period. The offsetRampReferenceTime parameter controls whether the ramp times are relative to the hotstart time or the coldstart time. Prior to the start of the ramp, the offset values are multiplied by 0.0, eliminating the offset until the ramp period starts.

The input parameter offsetSkipSnaps is used to skip the specified number of initial data sets in the offsets data file. It can also be used to insert blank snaps at the beginning of the run if it is negative.

==Water Level Offset Input File==

The water level offset input file starts with the following three lines:

# test file for 1ft offset on hsofs mesh during matthew
14400.0 # time increment between datasets -- read but ignored if constant offset
0.0 # default offset value

The first line is a comment line. The second line is the time interval in seconds between successive datasets. The third line is the default offset value; this allows us to create offset input files in a "sparse" format where only the non-default values are provided.

If there is only one dataset in the offset file, then the values are treated as temporally constant. If there is more than one dataset, ADCIRC linearly interpolates between them in time. If ADCIRC runs out of offset data, it continues with the last dataset as a set of constant values for the remainder of the run.

Because each offset dataset may have a different number of values, the datasets are separated by a hash tag ('#') in the 2nd column. When ADCIRC sees this, it knows it has reached the end of the dataset it is currently reading.

==Output==

If a water level offset input file is specified, then the time interpolated water level offset values will be written to a file called offset.63 on the same schedule as output to the fort.63. Currently, the offset.63 file can only be produced in ASCII format. This file can then be used as a diagnostic to check that the offset values have been applied as expected.

==Compatibility==

Although the offset is implemented as a pseudo barometric pressure, it is entirely independent of any/all meteorological forcing and meteorological data. The values do not show up in the fort.71 (barometric pressure stations) or fort.73 (fulldomain barometric pressure field) output files. The offset can be used with any meteorological forcing data, or with no meteorological forcing data at all.

==Example==

An ADCIRC model is coldstarted with a tide only run that uses a 10 day ramp, and runs for 15 days. How should the water level offset be configured?

In this case I suggest that you apply your bias correction after reaching full strength tidal forcing (i.e., after the 10 day tidal ramp). So to have your bias correction ramp completed at the end of the 15 day tidal spinup, I would start the ramp period at 14.5 days and end it at 15.0 days. That way you'll be able to clearly see the correction occurring when you look at the output data.

In summary, I suggest you set your ramp period as follows:

offsetRampStart =1252800.0, offsetRampEnd = 1296000.0

For the bias correction data file for this 15.0 day tidal spinup case, you'll only need a single data set, which will be treated as a constant value because there is only one dataset, so the format of the file will be for example

# bias correction values test for GFS 20161002 cycle 00Z
99999.9 # time increment in seconds; not used in this case
0.0 # default nodal value applied to any node not specified below
234 0.5 <-- node number and final bias correction value
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the dataset

After the spinup, in preparation for a forecast run, when you want to migrate from one set of bias correction values to another over a 6 hour time period, your offsetControl line would be something like this:

&offsetControl offsetFileName= "offset_migration_from_00Z_to_06Z.dat", offsetMultiplier = 1.0, offsetRampStart = 1252800.0, offsetRampEnd = 1296000.0, offsetRampReferenceTime = "coldstart", offsetSkipSnaps = 0 /

and the corresponding file would look like something like the following:

# bias correction values for GFS 20161002 from cycle 00Z to cycle 06Z
21600.0 # 6 hour time increment in seconds
0.0 # default nodal value applied to any node not specified below
234 0.5 <--- values match those at the end of the tidal spinup
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the first dataset
512 0.3 <-- new values for updated bias correction at 06Z, not necessarily defined at the same list of nodes
1001 -0.5
2346 0.74
4000 -0.1
# <-- hash tag in second column indicating the end of the second dataset

At the end of the 15.0 day tidal spinup in this example the bias correction at node 234 would be 0.5. Three hours later (after the hotstart) it would be 0.25. Six hours after the hotstart it would be zero.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 1001 would be 0.0. Three hours later (after the hotstart) it would be -0.25. Six hours after the hotstart it would be -0.5.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 4000 would be 0.1. Three hours later (after the hotstart) it would be 0.0. Six hours after the hotstart it would be -0.1.

==Frequently Asked Questions==

Question: When writing the &offsetControl namelist, should I write the parameters on a single line, comma separated?

Answer: I always put fortran namelist parameters on a single line and separate the assignments with commas, but other formatting may be possible depending on the compiler you are using. I've never had issues with the way I format these fortran namelist input lines, but I've never tried any other way so I'll defer to your compiler's documentation on that.

Question: Do I need to ramp the offset, i.e., is it wise to do so for some reason?

Answer: We haven't tried a sudden application of the water level offset to the solution so we don't yet know how that would behave. Our shortest ramp up period is 12 hours but you might get away with 6 hours or even less.

Question: Does the application of a dynamic water level offset surface affect the velocity solution (not just the elevation solution as it is primarily intended)?

Answer: Yes, in a physically realistic manner. As an example, when applying a water level offset to a back bay area with a nearby tidal inlet, the water level change in the back bay must be reflected in the velocity in the inlet in order to maintain continuity (conservation of mass). In the case of an initial application of a temporally steady water level offset, this velocity effect will be a one-time occurrence. On the other hand, for a time varying water level offset surface, the continuous adjustments to the water level will be continuously reflected in the velocity solution in physically realistic ways.

==References==

[Ricks 2017 AMS ppt]
[Taylor's manuscript]

Dynamic water level correction

2019-03-22T15:39:22Z

JasonFleming:

Main Page

2019-03-22T15:36:53Z

JasonFleming:

Hey there. For editors and administrators, links to some general wiki information are provided toward the bottom of this page, while more specific editing help can be found at the "Help" link on the sidebar of this page. The following summary shows some of the major pages that need content. Blue links are pages that have been created (but still definitely need content). Red links are pages that haven't been created yet (clicking the link will let you create the page). Have fun, and welcome to the party!


{| style="border-style: solid; border-width:1px; width:100%; vertical-align:top; border-spacing:4px;"
|- style="vertical-align:top;
| style="width:50%" | 

{| width="100%" cellpadding="2" cellspacing="5" style="vertical-align:top; "
|-
! <div style="margin:0; background:#97bfde; font-size:120%; font-weight:bold; text-align:left; color:#13294B; padding:0.2em 0.4em;">Basics</div>
|-
|style="color:#000;"|
*[[adcirc:About|About ADCIRC]]
*[[Features]]
*[[Getting started]]
*[[FAQ]]
*[[Version information]]
|}
| style="width:50%" | 

{| width="100%" cellpadding="2" cellspacing="5" style="vertical-align:top; "
|-
! <div style="margin:0; background:#b9b7eb; font-size:120%; font-weight:bold; text-align:left; color:#000; padding:0.2em 0.4em;">User Guide</div>
|-
|style="color:#000;"|
*[[Compiling]]
*[[File formats]]
*[[Model setup]]
*[[Running]]
*[[Visualization]]
*[[Tools]]
|}
|- style="vertical-align:top;
|

{| width="100%" cellpadding="2" cellspacing="5" style="vertical-align:top; "
|-
! <div style="margin:0; background:#b7ebdb; font-size:120%; font-weight:bold; text-align:left; color:#000; padding:0.2em 0.4em;">How-To</div>
|-
|style="color:#000;"|
*[[Tutorials]]
*[[Examples and Test Cases]]
*[[Troubleshooting]]
*[[Grid Development and Editing]]
*[[Tips]]
|}
|

{| width="100%" cellpadding="2" cellspacing="5" style="vertical-align:top; "
|-
! <div style="margin:0; background:#7ea0ff; font-size:120%; font-weight:bold; text-align:left; color:#000; padding:0.2em 0.4em;">Technical</div>
|-
|style="color:#000;"|
*[[ADCIRC physics|Physics]]
*[[ADCIRC numerics|Numerics]]
*[[Code architecture]]
*[[Parallelization]]
*[[Coupling]]
|}
|}

== Contributor/Editor Resources ==
A quick way to learn Wiki syntax (the source code is called "wikitext") is to go to a Wikipedia page and click on the "Edit source" or "View source" link at the top of the page. I often use [https://en.wikipedia.org/wiki/Navier%E2%80%93Stokes_equations | the Navier Stokes page]. If something works on Wikipedia but not here, let Taylor know.
* [[Help:Style Guide|ADCIRC Wiki style guide]] (under development)
* [http://www.workingwithmediawiki.com/book/chapter3.html Nice list of key editing abilities]
* [https://meta.wikimedia.org/wiki/Help:Contents#For_editors General editing help]
* [https://meta.wikimedia.org/wiki/Help:Starting_a_new_page Starting a new page]
* [https://en.wikipedia.org/wiki/Wikipedia:Sandbox "Sandbox"] to play around with editing without affecting anything...this is on Wikipedia but should mostly work the same
* [https://meta.wikimedia.org/wiki/Help:Advanced_editing Advanced editing]
* [https://meta.wikimedia.org/wiki/Help:Reverting Reverting edits]
* [https://www.mediawiki.org/wiki/Help:Templates Templates help]
* Somewhat related: [https://en.wikipedia.org/wiki/Help:Directory Wikipedia help directory], [https://en.wikipedia.org/wiki/Template:Wikipedia_technical_help Wikipedia technical help list], [https://en.wikipedia.org/wiki/Help:Text_editor_support Wikipedia text editor help]

== Moderator/Admin Resources ==

* [https://www.mediawiki.org/wiki/Manual:Administrators Administrator manual]
* [https://www.mediawiki.org/wiki/Sysadmin_hub Sysadmin hub]
* [https://www.mediawiki.org/wiki/Manual:Configuration_settings Configuration settings]

== Wiki Management: Getting started ==
* [https://www.mediawiki.org/wiki/Special:MyLanguage/Manual:Configuration_settings Configuration settings list]
* [https://www.mediawiki.org/wiki/Special:MyLanguage/Manual:FAQ MediaWiki FAQ]
* [https://lists.wikimedia.org/mailman/listinfo/mediawiki-announce MediaWiki release mailing list]
* [https://www.mediawiki.org/wiki/Special:MyLanguage/Localisation#Translation_resources Localise MediaWiki for your language]
* [https://www.mediawiki.org/wiki/Special:MyLanguage/Manual:Combating_spam Learn how to combat spam on your wiki]
* [https://meta.wikimedia.org/wiki/Help:Contents User's Guide]

Dynamic water level correction

2019-03-22T15:28:38Z

JasonFleming: Created page with "AMS 2017 Presentation from Luettich, et al: Representing Low Frequency, Spatially Varying Waterlevel Anomalies in Storm Surge Forecasting https://ams.confex.com/ams/97Annual..."

AMS 2017 Presentation from Luettich, et al:

Representing Low Frequency, Spatially Varying Waterlevel Anomalies in Storm Surge Forecasting

https://ams.confex.com/ams/97Annual/videogateway.cgi/id/37411?recordingid=37411&uniqueid=Paper316033&entry_password=655934

[include link to Rick's powerpoint slides]
[include reference to Rick's powerpoint slides in references section at end]

[From T.G.Asher]: Improvements in forecasting and hindcasting of tropical cyclones have led to much more accurate meteorological data associated with these storms. These improvements, together with more detailed representations of the physics and greater model resolution, have improved storm surge model skill. This has led to a state in which storm surge forecast and hindcast errors are often dominated by what we refer to as unresolved drivers, i.e., physical processes that are not explicitly included in storm surge models. Common examples include baroclinic processes, major oceanic currents, precipitation, steric fluctuations and far-field atmospheric forcing. These phenomena typically have longer timescales than the storm surge itself and therefore often show up as biases in water level during the period leading up to, encompassing and following the surge event.

[include reference to Taylor's paper in references section at end]

In order to accommodate these unresolved drivers, a water level offset ingestion feature has been added to ADCIRC.

Controlling Water Level Offset

The water level offset feature is triggered by the presence of the &offsetControl namelist at the bottom of the fort.15 file. Here is an example of how this line is used:

&offsetControl offsetFileName="hsofs_offset_test_ft.dat" offsetMultiplier=0.3048, offsetRampStart=0.0, offsetRampEnd=259200.0 offsetRampReferenceTime="hotstart", offsetSkipSnaps=0 /

The code looks for the water level offset file in the same directory as the other fulldomain files, both in serial and in parallel. In parallel, each subdomain does domain decomposition of the offset data on the fly, so the only change to adcprep is the read/write of the new namelist to the subdomain fort.15 files.

The multiplier value allows us to easily change the entire dataset by single factor. It may be needed in a case where the data were provided in ft but adcirc requires meters, for example.

The offsetRampStart and offsetRampEnd parameter are the time in seconds for a linear ramp up period. The offsetRampReferenceTime parameter controls whether the ramp times are relative to the hotstart time or the coldstart time. Prior to the start of the ramp, the offset values are multiplied by 0.0, eliminating the offset until the ramp period starts.

The input parameter offsetSkipSnaps is used to skip the specified number of initial data sets in the offsets data file. It can also be used to insert blank snaps at the beginning of the run if it is negative.

Water Level Offset Input File

The water level offset input file starts with the following three lines:

# test file for 1ft offset on hsofs mesh during matthew
14400.0 # time increment between datasets -- read but ignored if constant offset
0.0 # default offset value

The first line is a comment line. The second line is the time interval in seconds between successive datasets. The third line is the default offset value; this allows us to create offset input files in a "sparse" format where only the non-default values are provided.

If there is only one dataset in the offset file, then the values are treated as temporally constant. If there is more than one dataset, ADCIRC linearly interpolates between them in time. If ADCIRC runs out of offset data, it continues with the last dataset as a set of constant values for the remainder of the run.

Because each offset dataset may have a different number of values, the datasets are separated by a hash tag ('#') in the 2nd column. When ADCIRC sees this, it knows it has reached the end of the dataset it is currently reading.

Output

If a water level offset input file is specified, then the time interpolated water level offset values will be written to a file called offset.63 on the same schedule as output to the fort.63. Currently, the offset.63 file can only be produced in ASCII format. This file can then be used as a diagnostic to check that the offset values have been applied as expected.

Compatibility

Although the offset is implemented as a pseudo barometric pressure, it is entirely independent of any/all meteorological forcing and meteorological data. The values do not show up in the fort.71 (barometric pressure stations) or fort.73 (fulldomain barometric pressure field) output files. The offset can be used with any meteorological forcing data, or with no meteorological forcing data at all.

Example

An ADCIRC model is coldstarted with a tide only run that uses a 10 day ramp, and runs for 15 days. How should the water level offset be configured?

In this case I suggest that you apply your bias correction after reaching full strength tidal forcing (i.e., after the 10 day tidal ramp). So to have your bias correction ramp completed at the end of the 15 day tidal spinup, I would start the ramp period at 14.5 days and end it at 15.0 days. That way you'll be able to clearly see the correction occurring when you look at the output data.

In summary, I suggest you set your ramp period as follows:

offsetRampStart =1252800.0, offsetRampEnd = 1296000.0

For the bias correction data file for this 15.0 day tidal spinup case, you'll only need a single data set, which will be treated as a constant value because there is only one dataset, so the format of the file will be for example

# bias correction values test for GFS 20161002 cycle 00Z
99999.9 # time increment in seconds; not used in this case
0.0 # default nodal value applied to any node not specified below
234 0.5 <-- node number and final bias correction value
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the dataset

After the spinup, in preparation for a forecast run, when you want to migrate from one set of bias correction values to another over a 6 hour time period, your offsetControl line would be something like this:

&offsetControl offsetFileName= "offset_migration_from_00Z_to_06Z.dat", offsetMultiplier = 1.0, offsetRampStart = 1252800.0, offsetRampEnd = 1296000.0, offsetRampReferenceTime = "coldstart", offsetSkipSnaps = 0 /

and the corresponding file would look like something like the following:

# bias correction values for GFS 20161002 from cycle 00Z to cycle 06Z
21600.0 # 6 hour time increment in seconds
0.0 # default nodal value applied to any node not specified below
234 0.5 <--- values match those at the end of the tidal spinup
1223 0.6
1789 -0.2
4000 0.1
# <-- hash tag in second column indicating the end of the first dataset
512 0.3 <-- new values for updated bias correction at 06Z, not necessarily defined at the same list of nodes
1001 -0.5
2346 0.74
4000 -0.1
# <-- hash tag in second column indicating the end of the second dataset

At the end of the 15.0 day tidal spinup in this example the bias correction at node 234 would be 0.5. Three hours later (after the hotstart) it would be 0.25. Six hours after the hotstart it would be zero.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 1001 would be 0.0. Three hours later (after the hotstart) it would be -0.25. Six hours after the hotstart it would be -0.5.

At the end of the 15.0 day tidal spinup in this example the bias correction at node 4000 would be 0.1. Three hours later (after the hotstart) it would be 0.0. Six hours after the hotstart it would be -0.1.

Frequently Asked Questions

Question: When writing the &offsetControl namelist, should I write the parameters on a single line, comma separated?

Answer: I always put fortran namelist parameters on a single line and separate the assignments with commas, but other formatting may be possible depending on the compiler you are using. I've never had issues with the way I format these fortran namelist input lines, but I've never tried any other way so I'll defer to your compiler's documentation on that.

Question: Do I need to ramp the offset, i.e., is it wise to do so for some reason?

Answer: We haven't tried a sudden application of the water level offset to the solution so we don't yet know how that would behave. Our shortest ramp up period is 12 hours but you might get away with 6 hours or even less.

Question: Does the application of a dynamic water level offset surface affect the velocity solution (not just the elevation solution as it is primarily intended)?

Answer: Yes, in a physically realistic manner. As an example, when applying a water level offset to a back bay area with a nearby tidal inlet, the water level change in the back bay must be reflected in the velocity in the inlet in order to maintain continuity (conservation of mass). In the case of an initial application of a temporally steady water level offset, this velocity effect will be a one-time occurrence. On the other hand, for a time varying water level offset surface, the continuous adjustments to the water level will be continuously reflected in the velocity solution in physically realistic ways.

References

[Ricks AMS ppt]
[Taylor's manuscript]