oceanmodeling / WW3

WAVEWATCH III
Other
0 stars 0 forks source link

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

Open janahaddad opened 7 months ago

janahaddad commented 7 months ago

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

...continue documenting this task here

janahaddad commented 7 months ago

To do:

yunfangsun commented 6 months ago

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

Screenshot 2024-05-07 at 1 39 46 PM
janahaddad commented 6 months ago

See meeting notes here: https://github.com/oceanmodeling/ufs-coastal/discussions/83

janahaddad commented 5 months ago

@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:

yunfangsun commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago
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 commented 5 months ago

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

yunfangsun commented 5 months ago

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 commented 5 months ago

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

aliabdolali commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago

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

yunfangsun commented 5 months ago

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 commented 5 months ago

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 commented 1 month ago

@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 commented 1 month ago

@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 commented 1 month ago

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:

Screenshot 2024-10-09 at 10 31 25 PM

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 commented 1 month ago

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 commented 1 month ago

@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 commented 1 month ago

Hi @gseroka ,

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

gseroka commented 1 month ago

Hi @gseroka ,

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

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