Run global UnWW3 to provide boundary to Atlantic set up #9

janahaddad · 2024-05-01T15:21:00Z

https://github.com/noaa-ocs-modeling/SurgeTeamCoordination/issues/431

...continue documenting this task here

janahaddad · 2024-05-01T15:27:04Z

To do:

explore WW3-ERA5 boundary forcing option

yunfangsun · 2024-05-07T17:42:00Z

https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form
The available variables for wave:

janahaddad · 2024-05-09T14:12:43Z

See meeting notes here: oceanmodeling/ufs-weather-model#83

janahaddad · 2024-06-17T16:08:21Z

@yunfangsun you've run this and generated a boundary condition, correct? The next step is to verify that the BC you created is correct?

Can you add some info here:

summary of steps to generate the BC file from global run
where is the run directory and output and generated BC located
what are next steps

yunfangsun · 2024-06-24T21:34:40Z

Hi @aliabdolali ,

I have used your global wave model configuration and run a simulation, when I generate the whole filed by using ww3_ounf.inp

vi ww3_ounf.inp

$ -------------------------------------------------------------------- $
$ WAVEWATCH III Grid output post-processing                            $
$--------------------------------------------------------------------- $
$ Time, time increment and number of outputs
$
$ 10800=Every 3 hours. Use a large number to go all the way to end.
$  20220901 000000 10800 999
$ 3600=Every 1 hour. Use a large number to go all the way to end.
  20220915 000000 3600 720
$
$--------------------------------------------------------------------- $
$ Output request flags identifying fields as in ww3_shel input and
$ section 2.4 fo the manual.
$
$ N for namelist
  N
$
HS FP
$  DPT HS FP DIR SPR
$
$--------------------------------------------------------------------- $
$ Output type 4 [3,4] (version netCDF)
$        and   variable type 4 [2 = SHORT, 4 = REAL]
$ Output type 0 1 2 [0,1,2,3,4,5]   (swell partition)
$ variables T [T] or not [F] in the same file
$
  3 4
  0 1 2
  F
$
$ -------------------------------------------------------------------- $
$ File prefix
$ number of characters in date
$ IX, IY range
$
ww3.
6
1 655113  1 1
$       nodeCount = 4996111 ;
$       elementCount = 2470094 ;
$
$ For each field and time a new file is generated with the file name
$ ww3.date_xxx.nc , where date is a conventional time idicator with S3
$ characters,
$ and xxx is a field identifier.
$
$
$ -------------------------------------------------------------------- $
$ End of input file                                                    $

It only generate one time step of results

 ncdump -v time ww3.202209_hs.nc
netcdf ww3.202209_hs {
dimensions:
        level = 1 ;
        node = 655113 ;
        element = 1290999 ;
        time = UNLIMITED ; // (45 currently)
        noel = 3 ;
variables:
        float longitude(node) ;
                longitude:units = "degree_east" ;
                longitude:long_name = "longitude" ;
                longitude:standard_name = "longitude" ;
                longitude:valid_min = -180.f ;
                longitude:valid_max = 360.f ;
                longitude:axis = "X" ;
        float latitude(node) ;
                latitude:units = "degree_north" ;
                latitude:long_name = "latitude" ;
                latitude:standard_name = "latitude" ;
                latitude:valid_min = -90.f ;
                latitude:valid_max = 90.f ;
                latitude:axis = "Y" ;
        double time(time) ;
                time:long_name = "time" ;
                time:standard_name = "time" ;
                time:units = "days since 1990-01-01T00:00:00" ;
                time:axis = "T" ;
                time:calendar = "standard" ;
        int tri(element, noel) ;
        short MAPSTA(node) ;
                MAPSTA:long_name = "status map" ;
                MAPSTA:standard_name = "status map" ;
                MAPSTA:units = "1" ;
                MAPSTA:valid_min = -32 ;
                MAPSTA:valid_max = 32 ;
        float hs(time, node) ;
                hs:long_name = "significant height of wind and swell waves" ;
                hs:standard_name = "sea_surface_wave_significant_height" ;
                hs:globwave_name = "significant_wave_height" ;
                hs:units = "m" ;
                hs:_FillValue = 9.96921e+36f ;
                hs:scale_factor = 1.f ;
                hs:add_offset = 0.f ;
                hs:valid_min = 0.f ;
                hs:valid_max = 64.f ;

// global attributes:
                :WAVEWATCH_III_version_number = "7.14" ;
                :WAVEWATCH_III_switches = "PDLIB SCOTCH NOGRB DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 DB1 MLIM FLD1 TR0 BS0 WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 O14 O15 IC0 IS0 REF0" ;
                :SIN4\ namelist\ parameter\ BETAMAX = 1.33f ;
                :product_name = "ww3.202209_hs.nc" ;
                :area = "GLOBAL" ;
                :latitude_resolution = "n/a" ;
                :longitude_resolution = "n/a" ;
                :southernmost_latitude = "-78." ;
                :northernmost_latitude = "90." ;
                :westernmost_longitude = "-180." ;
                :easternmost_longitude = "180." ;
                :minimum_altitude = "-12000 m" ;
                :maximum_altitude = "9000 m" ;
                :altitude_resolution = "n/a" ;
                :start_date = "2022-09-15 00:00:00" ;
data:

 time = 11945, 11945.0416666667, 11945.0833333333, 11945.125,
    11945.1666666667, 11945.2083333333, 11945.25, 11945.2916666667,
    11945.3333333333, 11945.375, 11945.4166666667, 11945.4583333333, 11945.5,
    11945.5416666667, 11945.5833333333, 11945.625, 11945.6666666667,
    11945.7083333333, 11945.75, 11945.7916666667, 11945.8333333333,
    11945.875, 11945.9166666667, 11945.9583333333, 11946, 11946.0416666667,
    11946.0833333333, 11946.125, 11946.1666666667, 11946.2083333333,
    11946.25, 11946.2916666667, 11946.3333333333, 11946.375,
    11946.4166666667, 11946.4583333333, 11946.5, 11946.5416666667,
    11946.5833333333, 11946.625, 11946.6666666667, 11946.7083333333,
    11946.75, 11946.7916666667, 11946.8333333333 ;

and the ww3_ounp.inp could give me the correct outputs, could you please suggest me what I should modify in ww3_ounf.inp?

Thank you!

aliabdolali · 2024-06-24T22:31:06Z

from what you pasted here, I see 45 time steps in the netcdf file.
time = UNLIMITED ; // (45 currently)

you can see those time steps if you try: ncdump -tv time ww3.202209_hs.nc

yunfangsun · 2024-06-27T14:01:08Z

ncdump -tv time ww3.202209_hs.nc
netcdf ww3.202209_hs {
dimensions:
	level = 1 ;
	node = 655113 ;
	element = 1290999 ;
	time = UNLIMITED ; // (45 currently)
	noel = 3 ;
variables:
	float longitude(node) ;
		longitude:units = "degree_east" ;
		longitude:long_name = "longitude" ;
		longitude:standard_name = "longitude" ;
		longitude:valid_min = -180.f ;
		longitude:valid_max = 360.f ;
		longitude:axis = "X" ;
	float latitude(node) ;
		latitude:units = "degree_north" ;
		latitude:long_name = "latitude" ;
		latitude:standard_name = "latitude" ;
		latitude:valid_min = -90.f ;
		latitude:valid_max = 90.f ;
		latitude:axis = "Y" ;
	double time(time) ;
		time:long_name = "time" ;
		time:standard_name = "time" ;
		time:units = "days since 1990-01-01T00:00:00" ;
		time:axis = "T" ;
		time:calendar = "standard" ;
	int tri(element, noel) ;
	short MAPSTA(node) ;
		MAPSTA:long_name = "status map" ;
		MAPSTA:standard_name = "status map" ;
		MAPSTA:units = "1" ;
		MAPSTA:valid_min = -32 ;
		MAPSTA:valid_max = 32 ;
	float hs(time, node) ;
		hs:long_name = "significant height of wind and swell waves" ;
		hs:standard_name = "sea_surface_wave_significant_height" ;
		hs:globwave_name = "significant_wave_height" ;
		hs:units = "m" ;
		hs:_FillValue = 9.96921e+36f ;
		hs:scale_factor = 1.f ;
		hs:add_offset = 0.f ;
		hs:valid_min = 0.f ;
		hs:valid_max = 64.f ;

// global attributes:
		:WAVEWATCH_III_version_number = "7.14" ;
		:WAVEWATCH_III_switches = "PDLIB SCOTCH NOGRB DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 DB1 MLIM FLD1 TR0 BS0 WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 O14 O15 IC0 IS0 REF0" ;
		:SIN4\ namelist\ parameter\ BETAMAX = 1.33f ;
		:product_name = "ww3.202209_hs.nc" ;
		:area = "GLOBAL" ;
		:latitude_resolution = "n/a" ;
		:longitude_resolution = "n/a" ;
		:southernmost_latitude = "-78." ;
		:northernmost_latitude = "90." ;
		:westernmost_longitude = "-180." ;
		:easternmost_longitude = "180." ;
		:minimum_altitude = "-12000 m" ;
		:maximum_altitude = "9000 m" ;
		:altitude_resolution = "n/a" ;
		:start_date = "2022-09-15 00:00:00" ;
data:

 time = "2022-09-15", "2022-09-15 01", "2022-09-15 02", "2022-09-15 03",
    "2022-09-15 04", "2022-09-15 05", "2022-09-15 06", "2022-09-15 07",
    "2022-09-15 08", "2022-09-15 09", "2022-09-15 10", "2022-09-15 11",
    "2022-09-15 12", "2022-09-15 13", "2022-09-15 14", "2022-09-15 15",
    "2022-09-15 16", "2022-09-15 17", "2022-09-15 18", "2022-09-15 19",
    "2022-09-15 20", "2022-09-15 21", "2022-09-15 22", "2022-09-15 23",
    "2022-09-16", "2022-09-16 01", "2022-09-16 02", "2022-09-16 03",
    "2022-09-16 04", "2022-09-16 05", "2022-09-16 06", "2022-09-16 07",
    "2022-09-16 08", "2022-09-16 09", "2022-09-16 10", "2022-09-16 11",
    "2022-09-16 12", "2022-09-16 13", "2022-09-16 14", "2022-09-16 15",
    "2022-09-16 16", "2022-09-16 17", "2022-09-16 18", "2022-09-16 19",
    "2022-09-16 20" ;
}

but```
ncdump -tv time ww3.202209_spec.nc
netcdf ww3.202209_spec {
dimensions:
time = UNLIMITED ; // (384 currently)
station = 325 ;
string40 = 40 ;
frequency = 50 ;
direction = 36 ;
variables:
double time(time) ;
time:long_name = "julian day (UT)" ;
time:standard_name = "time" ;
time:units = "days since 1990-01-01 00:00:00" ;
time:conventions = "Relative julian days with decimal part (as parts of the day)" ;
time:axis = "T" ;
time:calendar = "standard" ;
int station(station) ;
station:long_name = "station id" ;
station:_FillValue = -2147483647 ;
station:axis = "X" ;
int string40(string40) ;
string40:long_name = "station_name number of characters" ;
string40:_FillValue = -2147483647 ;
string40:axis = "W" ;
char station_name(station, string40) ;
station_name:long_name = "station name" ;
station_name:content = "XW" ;
station_name:associates = "station string40" ;
float longitude(time, station) ;
longitude:long_name = "longitude" ;
longitude:standard_name = "longitude" ;
longitude:globwave_name = "longitude" ;
longitude:units = "degree_east" ;
longitude:scale_factor = 1.f ;
longitude:add_offset = 0.f ;
longitude:valid_min = -180.f ;
longitude:valid_max = 360.f ;
longitude:_FillValue = 9.96921e+36f ;
longitude:content = "TX" ;
longitude:associates = "time station" ;
float latitude(time, station) ;
latitude:long_name = "latitude" ;
latitude:standard_name = "latitude" ;
latitude:globwave_name = "latitude" ;
latitude:units = "degree_north" ;
latitude:scale_factor = 1.f ;
latitude:add_offset = 0.f ;
latitude:valid_min = -90.f ;
latitude:valid_max = 180.f ;
latitude:_FillValue = 9.96921e+36f ;
latitude:content = "TX" ;
latitude:associates = "time station" ;
float frequency(frequency) ;
frequency:long_name = "frequency of center band" ;
frequency:standard_name = "sea_surface_wave_frequency" ;
frequency:globwave_name = "frequency" ;
frequency:units = "s-1" ;
frequency:scale_factor = 1.f ;
frequency:add_offset = 0.f ;
frequency:valid_min = 0.f ;
frequency:valid_max = 10.f ;
frequency:_FillValue = 9.96921e+36f ;
frequency:axis = "Y" ;
float frequency1(frequency) ;
frequency1:long_name = "frequency of lower band" ;
frequency1:standard_name = "frequency_of_lower_band" ;
frequency1:globwave_name = "frequency_lower_band" ;
frequency1:units = "s-1" ;
frequency1:scale_factor = 1.f ;
frequency1:add_offset = 0.f ;
frequency1:valid_min = 0.f ;
frequency1:valid_max = 10.f ;
frequency1:_FillValue = 9.96921e+36f ;
frequency1:content = "Y" ;
frequency1:associates = "frequency" ;
float frequency2(frequency) ;
frequency2:long_name = "frequency of upper band" ;
frequency2:standard_name = "frequency_of_upper_band" ;
frequency2:globwave_name = "frequency_upper_band" ;
frequency2:units = "s-1" ;
frequency2:scale_factor = 1.f ;
frequency2:add_offset = 0.f ;
frequency2:valid_min = 0.f ;
frequency2:valid_max = 10.f ;
frequency2:_FillValue = 9.96921e+36f ;
frequency2:content = "Y" ;
frequency2:associates = "frequency" ;
float direction(direction) ;
direction:long_name = "sea surface wave to direction" ;
direction:standard_name = "sea_surface_wave_to_direction" ;
direction:globwave_name = "direction" ;
direction:units = "degree" ;
direction:scale_factor = 1.f ;
direction:add_offset = 0.f ;
direction:valid_min = 0.f ;
direction:valid_max = 360.f ;
direction:_FillValue = 9.96921e+36f ;
direction:axis = "Z" ;
float efth(time, station, frequency, direction) ;
efth:long_name = "sea surface wave directional variance spectral density" ;
efth:standard_name = "sea_surface_wave_directional_variance_spectral_density" ;
efth:globwave_name = "directional_variance_spectral_density" ;
efth:units = "m2 s rad-1" ;
efth:scale_factor = 1.f ;
efth:add_offset = 0.f ;
efth:valid_min = 0.f ;
efth:valid_max = 1.e+20f ;
efth:_FillValue = 9.96921e+36f ;
efth:content = "TXYZ" ;
efth:associates = "time station frequency direction" ;
float dpt(time, station) ;
dpt:long_name = "depth" ;
dpt:standard_name = "depth" ;
dpt:globwave_name = "depth" ;
dpt:units = "m" ;
dpt:scale_factor = 1.f ;
dpt:add_offset = 0.f ;
dpt:valid_min = -100.f ;
dpt:valid_max = 10000.f ;
dpt:_FillValue = 9.96921e+36f ;
dpt:content = "TX" ;
dpt:associates = "time station" ;
float wnd(time, station) ;
wnd:long_name = "wind speed at 10m" ;
wnd:standard_name = "wind_speed" ;
wnd:globwave_name = "wind_speed" ;
wnd:units = "m s-1" ;
wnd:scale_factor = 1.f ;
wnd:add_offset = 0.f ;
wnd:valid_min = 0.f ;
wnd:valid_max = 100.f ;
wnd:_FillValue = 9.96921e+36f ;
wnd:content = "TX" ;
wnd:associates = "time station" ;
float wnddir(time, station) ;
wnddir:long_name = "wind direction" ;
wnddir:standard_name = "wind_from_direction" ;
wnddir:globwave_name = "wind_from_direction" ;
wnddir:units = "degree" ;
wnddir:scale_factor = 1.f ;
wnddir:add_offset = 0.f ;
wnddir:valid_min = 0.f ;
wnddir:valid_max = 360.f ;
wnddir:_FillValue = 9.96921e+36f ;
wnddir:content = "TX" ;
wnddir:associates = "time station" ;
float cur(time, station) ;
cur:long_name = "sea water speed" ;
cur:standard_name = "sea_water_speed" ;
cur:globwave_name = "sea_water_speed" ;
cur:units = "m s-1" ;
cur:scale_factor = 1.f ;
cur:add_offset = 0.f ;
cur:valid_min = 0.f ;
cur:valid_max = 100.f ;
cur:_FillValue = 9.96921e+36f ;
cur:content = "TX" ;
cur:associates = "time station" ;
float curdir(time, station) ;
curdir:long_name = "direction from of sea water velocity" ;
curdir:standard_name = "direction_of_sea_water_velocity" ;
curdir:globwave_name = "direction_of_sea_water_velocity" ;
curdir:units = "degree" ;
curdir:scale_factor = 1.f ;
curdir:add_offset = 0.f ;
curdir:valid_min = 0.f ;
curdir:valid_max = 360.f ;
curdir:_FillValue = 9.96921e+36f ;
curdir:content = "TX" ;
curdir:associates = "time station" ;

// global attributes:
:product_name = "ww3.202209_spec.nc" ;
:area = "GLOBAL" ;
:data_type = "OCO spectra 2D" ;
:format_version = "1.1" ;
:southernmost_latitude = "n/a" ;
:northernmost_latitude = "n/a" ;
:latitude_resolution = "n/a" ;
:westernmost_longitude = "n/a" ;
:easternmost_longitude = "n/a" ;
:longitude_resolution = "n/a" ;
:minimum_altitude = "n/a" ;
:maximum_altitude = "n/a" ;
:altitude_resolution = "n/a" ;
:start_date = "2022-09-15 00:00:00" ;
:stop_date = "2022-10-01 00:00:00" ;
:field_type = "hourly" ;
data:

time = "2022-09-15", "2022-09-15 01", "2022-09-15 02", "2022-09-15 03",
"2022-09-15 04", "2022-09-15 05", "2022-09-15 06", "2022-09-15 07",
"2022-09-15 08", "2022-09-15 09", "2022-09-15 10", "2022-09-15 11",
"2022-09-15 12", "2022-09-15 13", "2022-09-15 14", "2022-09-15 15",
"2022-09-15 16", "2022-09-15 17", "2022-09-15 18", "2022-09-15 19",
"2022-09-15 20", "2022-09-15 21", "2022-09-15 22", "2022-09-15 23",
"2022-09-16", "2022-09-16 01", "2022-09-16 02", "2022-09-16 03",
"2022-09-16 04", "2022-09-16 05", "2022-09-16 06", "2022-09-16 07",
"2022-09-16 08", "2022-09-16 09", "2022-09-16 10", "2022-09-16 11",
"2022-09-16 12", "2022-09-16 13", "2022-09-16 14", "2022-09-16 15",
"2022-09-16 16", "2022-09-16 17", "2022-09-16 18", "2022-09-16 19",
"2022-09-16 20", "2022-09-16 21", "2022-09-16 22", "2022-09-16 23",
"2022-09-17", "2022-09-17 01", "2022-09-17 02", "2022-09-17 03",
"2022-09-17 04", "2022-09-17 05", "2022-09-17 06", "2022-09-17 07",
"2022-09-17 08", "2022-09-17 09", "2022-09-17 10", "2022-09-17 11",
"2022-09-17 12", "2022-09-17 13", "2022-09-17 14", "2022-09-17 15",
"2022-09-17 16", "2022-09-17 17", "2022-09-17 18", "2022-09-17 19",
"2022-09-17 20", "2022-09-17 21", "2022-09-17 22", "2022-09-17 23",
"2022-09-18", "2022-09-18 01", "2022-09-18 02", "2022-09-18 03",
"2022-09-18 04", "2022-09-18 05", "2022-09-18 06", "2022-09-18 07",
"2022-09-18 08", "2022-09-18 09", "2022-09-18 10", "2022-09-18 11",
"2022-09-18 12", "2022-09-18 13", "2022-09-18 14", "2022-09-18 15",
"2022-09-18 16", "2022-09-18 17", "2022-09-18 18", "2022-09-18 19",
"2022-09-18 20", "2022-09-18 21", "2022-09-18 22", "2022-09-18 23",
"2022-09-19", "2022-09-19 01", "2022-09-19 02", "2022-09-19 03",
"2022-09-19 04", "2022-09-19 05", "2022-09-19 06", "2022-09-19 07",
"2022-09-19 08", "2022-09-19 09", "2022-09-19 10", "2022-09-19 11",
"2022-09-19 12", "2022-09-19 13", "2022-09-19 14", "2022-09-19 15",
"2022-09-19 16", "2022-09-19 17", "2022-09-19 18", "2022-09-19 19",
"2022-09-19 20", "2022-09-19 21", "2022-09-19 22", "2022-09-19 23",
"2022-09-20", "2022-09-20 01", "2022-09-20 02", "2022-09-20 03",
"2022-09-20 04", "2022-09-20 05", "2022-09-20 06", "2022-09-20 07",
"2022-09-20 08", "2022-09-20 09", "2022-09-20 10", "2022-09-20 11",
"2022-09-20 12", "2022-09-20 13", "2022-09-20 14", "2022-09-20 15",
"2022-09-20 16", "2022-09-20 17", "2022-09-20 18", "2022-09-20 19",
"2022-09-20 20", "2022-09-20 21", "2022-09-20 22", "2022-09-20 23",
"2022-09-21", "2022-09-21 01", "2022-09-21 02", "2022-09-21 03",
"2022-09-21 04", "2022-09-21 05", "2022-09-21 06", "2022-09-21 07",
"2022-09-21 08", "2022-09-21 09", "2022-09-21 10", "2022-09-21 11",
"2022-09-21 12", "2022-09-21 13", "2022-09-21 14", "2022-09-21 15",
"2022-09-21 16", "2022-09-21 17", "2022-09-21 18", "2022-09-21 19",
"2022-09-21 20", "2022-09-21 21", "2022-09-21 22", "2022-09-21 23",
"2022-09-22", "2022-09-22 01", "2022-09-22 02", "2022-09-22 03",
"2022-09-22 04", "2022-09-22 05", "2022-09-22 06", "2022-09-22 07",
"2022-09-22 08", "2022-09-22 09", "2022-09-22 10", "2022-09-22 11",
"2022-09-22 12", "2022-09-22 13", "2022-09-22 14", "2022-09-22 15",
"2022-09-22 16", "2022-09-22 17", "2022-09-22 18", "2022-09-22 19",
"2022-09-22 20", "2022-09-22 21", "2022-09-22 22", "2022-09-22 23",
"2022-09-23", "2022-09-23 01", "2022-09-23 02", "2022-09-23 03",
"2022-09-23 04", "2022-09-23 05", "2022-09-23 06", "2022-09-23 07",
"2022-09-23 08", "2022-09-23 09", "2022-09-23 10", "2022-09-23 11",
"2022-09-23 12", "2022-09-23 13", "2022-09-23 14", "2022-09-23 15",
"2022-09-23 16", "2022-09-23 17", "2022-09-23 18", "2022-09-23 19",
"2022-09-23 20", "2022-09-23 21", "2022-09-23 22", "2022-09-23 23",
"2022-09-24", "2022-09-24 01", "2022-09-24 02", "2022-09-24 03",
"2022-09-24 04", "2022-09-24 05", "2022-09-24 06", "2022-09-24 07",
"2022-09-24 08", "2022-09-24 09", "2022-09-24 10", "2022-09-24 11",
"2022-09-24 12", "2022-09-24 13", "2022-09-24 14", "2022-09-24 15",
"2022-09-24 16", "2022-09-24 17", "2022-09-24 18", "2022-09-24 19",
"2022-09-24 20", "2022-09-24 21", "2022-09-24 22", "2022-09-24 23",
"2022-09-25", "2022-09-25 01", "2022-09-25 02", "2022-09-25 03",
"2022-09-25 04", "2022-09-25 05", "2022-09-25 06", "2022-09-25 07",
"2022-09-25 08", "2022-09-25 09", "2022-09-25 10", "2022-09-25 11",
"2022-09-25 12", "2022-09-25 13", "2022-09-25 14", "2022-09-25 15",
"2022-09-25 16", "2022-09-25 17", "2022-09-25 18", "2022-09-25 19",
"2022-09-25 20", "2022-09-25 21", "2022-09-25 22", "2022-09-25 23",
"2022-09-26", "2022-09-26 01", "2022-09-26 02", "2022-09-26 03",
"2022-09-26 04", "2022-09-26 05", "2022-09-26 06", "2022-09-26 07",
"2022-09-26 08", "2022-09-26 09", "2022-09-26 10", "2022-09-26 11",
"2022-09-26 12", "2022-09-26 13", "2022-09-26 14", "2022-09-26 15",
"2022-09-26 16", "2022-09-26 17", "2022-09-26 18", "2022-09-26 19",
"2022-09-26 20", "2022-09-26 21", "2022-09-26 22", "2022-09-26 23",
"2022-09-27", "2022-09-27 01", "2022-09-27 02", "2022-09-27 03",
"2022-09-27 04", "2022-09-27 05", "2022-09-27 06", "2022-09-27 07",
"2022-09-27 08", "2022-09-27 09", "2022-09-27 10", "2022-09-27 11",
"2022-09-27 12", "2022-09-27 13", "2022-09-27 14", "2022-09-27 15",
"2022-09-27 16", "2022-09-27 17", "2022-09-27 18", "2022-09-27 19",
"2022-09-27 20", "2022-09-27 21", "2022-09-27 22", "2022-09-27 23",
"2022-09-28", "2022-09-28 01", "2022-09-28 02", "2022-09-28 03",
"2022-09-28 04", "2022-09-28 05", "2022-09-28 06", "2022-09-28 07",
"2022-09-28 08", "2022-09-28 09", "2022-09-28 10", "2022-09-28 11",
"2022-09-28 12", "2022-09-28 13", "2022-09-28 14", "2022-09-28 15",
"2022-09-28 16", "2022-09-28 17", "2022-09-28 18", "2022-09-28 19",
"2022-09-28 20", "2022-09-28 21", "2022-09-28 22", "2022-09-28 23",
"2022-09-29", "2022-09-29 01", "2022-09-29 02", "2022-09-29 03",
"2022-09-29 04", "2022-09-29 05", "2022-09-29 06", "2022-09-29 07",
"2022-09-29 08", "2022-09-29 09", "2022-09-29 10", "2022-09-29 11",
"2022-09-29 12", "2022-09-29 13", "2022-09-29 14", "2022-09-29 15",
"2022-09-29 16", "2022-09-29 17", "2022-09-29 18", "2022-09-29 19",
"2022-09-29 20", "2022-09-29 21", "2022-09-29 22", "2022-09-29 23",
"2022-09-30", "2022-09-30 01", "2022-09-30 02", "2022-09-30 03",
"2022-09-30 04", "2022-09-30 05", "2022-09-30 06", "2022-09-30 07",
"2022-09-30 08", "2022-09-30 09", "2022-09-30 10", "2022-09-30 11",
"2022-09-30 12", "2022-09-30 13", "2022-09-30 14", "2022-09-30 15",
"2022-09-30 16", "2022-09-30 17", "2022-09-30 18", "2022-09-30 19",
"2022-09-30 20", "2022-09-30 21", "2022-09-30 22", "2022-09-30 23" ;
}


The simulation finished the whole of September, but the HS postprocessing could only show 45 hours, and the spectra could show all time steps. I think there must be something wrong in my ww3_ounf.inp.

aliabdolali · 2024-06-27T15:35:43Z

please attach ww3_shel.inp(nml) and ww3_ounf.inp(nml) here

yunfangsun · 2024-06-27T16:03:52Z

ww3.shel.nml

! -------------------------------------------------------------------- !
! WAVEWATCH III ww3_shel.nml - single-grid model                       !
! -------------------------------------------------------------------- !


! -------------------------------------------------------------------- !
! Define top-level model parameters via DOMAIN_NML namelist
! -------------------------------------------------------------------- !
&DOMAIN_NML
  DOMAIN%IOSTYP  = 0
  DOMAIN%START   = '20220901 000000'
  DOMAIN%STOP    = '20221015 000000'
/

! -------------------------------------------------------------------- !
! Define each forcing via the INPUT_NML namelist
! -------------------------------------------------------------------- !
&INPUT_NML
  INPUT%FORCING%ICE_CONC      = 'F'
  INPUT%FORCING%CURRENTS      = 'F'
  INPUT%FORCING%WINDS         = 'T'
/

! -------------------------------------------------------------------- !
! Define the output types point parameters via OUTPUT_TYPE_NML namelist
! -------------------------------------------------------------------- !
&OUTPUT_TYPE_NML
  TYPE%FIELD%LIST          = 'WND HS FP DP PTP PDIR CUR ICE'
  TYPE%POINT%FILE          = 'ww3_points.list'
/

! -------------------------------------------------------------------- !
! Define output dates via OUTPUT_DATE_NML namelist
! -------------------------------------------------------------------- !
&OUTPUT_DATE_NML
  DATE%FIELD          = '20220901 000000' '3600' '20221015 000000'
  DATE%POINT          = '20220901 000000' '3600' '20221015 000000'
!  DATE%RESTART        = '20210401 000000' '86400' '20210402 000000'
/

! -------------------------------------------------------------------- !
! Define homogeneous input via HOMOG_COUNT_NML and HOMOG_INPUT_NML namelist
! -------------------------------------------------------------------- !
&HOMOG_COUNT_NML

&HOMOG_INPUT_NML
/

! -------------------------------------------------------------------- !
! WAVEWATCH III - end of namelist                                      !
! -------------------------------------------------------------------- !

ww3_ounf.inp

$ -------------------------------------------------------------------- $
$ WAVEWATCH III Grid output post-processing                            $
$--------------------------------------------------------------------- $
$ Time, time increment and number of outputs
$
$ 10800=Every 3 hours. Use a large number to go all the way to end.
$  20220901 000000 10800 999
$ 3600=Every 1 hour. Use a large number to go all the way to end.
  20220915 000000 3600 720
$
$--------------------------------------------------------------------- $
$ Output request flags identifying fields as in ww3_shel input and
$ section 2.4 fo the manual.
$
$ N for namelist
  N
$
HS FP
$  DPT HS FP DIR SPR
$
$--------------------------------------------------------------------- $
$ Output type 4 [3,4] (version netCDF)
$        and   variable type 4 [2 = SHORT, 4 = REAL]
$ Output type 0 1 2 [0,1,2,3,4,5]   (swell partition)
$ variables T [T] or not [F] in the same file
$
  3 4
  0 1 2
  F
$
$ -------------------------------------------------------------------- $
$ File prefix
$ number of characters in date
$ IX, IY range
$
ww3.
6
1 655113  1 1
$       nodeCount = 4996111 ;
$       elementCount = 2470094 ;
$
$ For each field and time a new file is generated with the file name
$ ww3.date_xxx.nc , where date is a conventional time idicator with S3
$ characters,
$ and xxx is a field identifier.
$
$
$ -------------------------------------------------------------------- $
$ End of input file                                                    $
$ -------------------------------------------------------------------- $

aliabdolali · 2024-06-27T16:09:10Z

change $ 20220901 000000 10800 999 to $ 20220901 000000 3600 9999 in ounf.nml

aliabdolali · 2024-06-27T16:24:21Z

in case you did not fix it, send me an invite and we can take a look at your log.ww3 and see what's wrong there

yunfangsun · 2024-06-27T17:04:29Z

Hi @aliabdolali
Thank you! Now I can produce the whole fields with time steps. I thought the program didn't read the line starting with $.

aliabdolali · 2024-06-27T18:32:05Z

It shouldn't, it is either a bug or some hidden character. That is why we switched to nml at ERDC

yunfangsun · 2024-07-01T15:22:49Z

For this week, I will try to set up to use the spectra from the outputs of the global wave model as open boundary conditions to drive the Atlantic domain with the coarse-resolution SCHISM mesh.

aliabdolali · 2024-07-01T15:28:54Z

For this week, I will try to set up to use the spectra from the outputs of the global wave model as open boundary conditions to drive the Atlantic domain with the coarse-resolution SCHISM mesh.

awesome, keep up posted on your progress or if you faced an issue.

gseroka · 2024-09-30T21:27:36Z

@janahaddad thanks for the tag up this morning and suggestion to post here for UFS Coastal Applications Team (CAT) project. Several UFS CAT testers (a couple at Colorado School of Mines, and a couple at Univ of South Florida) are interested in getting WW3 outputs for boundary conditions/forcing for their testing of SCHISM and FVCOM.

The location is Mid-Atlantic Bight (NY/NJ Harbor).

The time periods are:
(June spinup) Jul 1 - Sep 30, 2021: includes hurricanes (Elsa, Henri, Ida)
(Dec '21 spinup) Jan 1 - Mar 31, 2022: includes Nor’easter (Jan 14-17; Jan 28-29)

This is an initial comment to track this issue and start the discussion. Thank you.

gseroka · 2024-10-09T19:37:45Z

@janahaddad thanks for the tag up this morning and suggestion to post here for UFS Coastal Applications Team (CAT) project. Several UFS CAT testers (a couple at Colorado School of Mines, and a couple at Univ of South Florida) are interested in getting WW3 outputs for boundary conditions/forcing for their testing of SCHISM and FVCOM.

The location is Mid-Atlantic Bight (NY/NJ Harbor).

The time periods are: (June spinup) Jul 1 - Sep 30, 2021: includes hurricanes (Elsa, Henri, Ida) (Dec '21 spinup) Jan 1 - Mar 31, 2022: includes Nor’easter (Jan 14-17; Jan 28-29)

This is an initial comment to track this issue and start the discussion. Thank you.

@yunfangsun just following up on the above. Let me know if you were able to find out anything about the NOAA operational WW3 output (issues in wave spectra data but only prior to 2020 or 2016?) from the HPSS tapes. Thanks!

yunfangsun · 2024-10-10T02:46:56Z

Hi @gseroka ,

I have just rechecked, the global GFS wave spectra product from EMC is available from March 21, 2021, to June 26, 2022, version 16.2 from June 27, 2022, to Nov 29, 2022, version 16.3 from Nov 30, 2022, to now.

The spectra bin number is 50, direction bin number is 36.

The locations of the buoy around the Mid-Atlantic Bight are as follows:

The spectra bin is 50, the direction bin is 36. The global wave model is updated every 6 hours, with 384-hour results, we need to pick up the first 6 hours of each run to interpolate the open boundary conditions.

And before March 21, 2021, there are not enough buoy outputs from the global GFS wave model outputs.

I have downloaded the wave spectra from HPSS on Hera for Jun 1 - Sep 30, 2021, and Dec 1, 2021 - Mar 31, 2022.

And you can find some samples of the spectra on Hercules in the folder /work2/noaa/nosofs/yunfangs/hpss/2021/gfs.20210601/06/wave/station

gseroka · 2024-10-10T11:39:11Z

Hi @gseroka ,

I have just rechecked, the global GFS wave spectra product from EMC is available from March 21, 2021, to June 26, 2022, version 16.2 from June 27, 2022, to Nov 29, 2022, version 16.3 from Nov 30, 2022, to now.

The spectra bin number is 50, direction bin number is 36.

The locations of the buoy around the Mid-Atlantic Bight are as follows:

The spectra bin is 50, the direction bin is 36. The global wave model is updated every 6 hours, with 384-hour results, we need to pick up the first 6 hours of each run to interpolate the open boundary conditions.

And before March 21, 2021, there are not enough buoy outputs from the global GFS wave model outputs.

I have downloaded the wave spectra from HPSS on Hera for Jun 1 - Sep 30, 2021, and Dec 1, 2021 - Mar 31, 2022.

And you can find some samples of the spectra on Hercules in the folder /work2/noaa/nosofs/yunfangs/hpss/2021/gfs.20210601/06/wave/station

@yunfangsun thank you so much, highly appreciated. I will take a look and I will probably have some questions. Thank you again!

gseroka · 2024-10-10T17:14:01Z

@yunfangsun I added a short meeting tomorrow morning with you and Ayumi to discuss this. Let me know if this time doesn't work, thanks!

yunfangsun · 2024-10-10T17:32:18Z

Hi @gseroka ,

Yes, it works for me, should @janahaddad also be involved?

gseroka · 2024-10-10T17:34:11Z

Hi @gseroka ,

Yes, it works for me, should @janahaddad also be involved?

Sure, I'll add her! (as optional)

janahaddad added enhancement New feature or request and removed enhancement New feature or request labels May 1, 2024

janahaddad added this to ufs-coastal project May 1, 2024

janahaddad assigned yunfangsun May 1, 2024

janahaddad moved this to Todo in ufs-coastal project May 20, 2024

janahaddad mentioned this issue Jun 10, 2024

SCHISM+WW3 H. Ian testing-suite: Phase 3 oceanmodeling/ufs-weather-model#97

Open

4 tasks

janahaddad moved this from Todo to In Progress in ufs-coastal project Jul 1, 2024

janahaddad moved this from In Progress to Backlog in ufs-coastal project Sep 23, 2024

janahaddad mentioned this issue Oct 11, 2024

UFS-CAT Round 4 support: Create wave BCs for NY Harbor SCHISM and FVCOM domains #14

Closed

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Run global UnWW3 to provide boundary to Atlantic set up #9

Run global UnWW3 to provide boundary to Atlantic set up #9

janahaddad commented May 1, 2024

janahaddad commented May 1, 2024

yunfangsun commented May 7, 2024

janahaddad commented May 9, 2024

janahaddad commented Jun 17, 2024

yunfangsun commented Jun 24, 2024

aliabdolali commented Jun 24, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jul 1, 2024

aliabdolali commented Jul 1, 2024

gseroka commented Sep 30, 2024

gseroka commented Oct 9, 2024

yunfangsun commented Oct 10, 2024

gseroka commented Oct 10, 2024

gseroka commented Oct 10, 2024

yunfangsun commented Oct 10, 2024

gseroka commented Oct 10, 2024

Run global UnWW3 to provide boundary to Atlantic set up #9

Run global UnWW3 to provide boundary to Atlantic set up #9

Comments

janahaddad commented May 1, 2024

janahaddad commented May 1, 2024

yunfangsun commented May 7, 2024

janahaddad commented May 9, 2024

janahaddad commented Jun 17, 2024

yunfangsun commented Jun 24, 2024

aliabdolali commented Jun 24, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jun 27, 2024

aliabdolali commented Jun 27, 2024

yunfangsun commented Jul 1, 2024

aliabdolali commented Jul 1, 2024

gseroka commented Sep 30, 2024

gseroka commented Oct 9, 2024

yunfangsun commented Oct 10, 2024

gseroka commented Oct 10, 2024

gseroka commented Oct 10, 2024

yunfangsun commented Oct 10, 2024

gseroka commented Oct 10, 2024