NOAA-EMC / WW3

WAVEWATCH III
Other
267 stars 545 forks source link

Problems with Excluded Points on grid boundaries in ww3 v6.07 #991

Open JMSnyder1980 opened 1 year ago

JMSnyder1980 commented 1 year ago

Discussed in https://github.com/NOAA-EMC/WW3/discussions/989

Originally posted by **JMSnyder1980** April 20, 2023 I have tried running ww3_grid with the following input files: [namelists.nml.txt](https://github.com/NOAA-EMC/WW3/files/11290834/namelists.nml.txt) [ww3_grid.nml.txt](https://github.com/NOAA-EMC/WW3/files/11290836/ww3_grid.nml.txt) [ww3_grid.inp.txt](https://github.com/NOAA-EMC/WW3/files/11290837/ww3_grid.inp.txt) With the following bathymetry files: [global_1p00_bathy.dat.txt](https://github.com/NOAA-EMC/WW3/files/11290842/global_1p00_bathy.dat.txt) [global_1p00_obst.dat.txt](https://github.com/NOAA-EMC/WW3/files/11290843/global_1p00_obst.dat.txt) When running ww3 v6.17 with the above files I get several excluded points especially around 180 degree longitude. This is a problem since in the ww3 output I get a gap around 180 that I do not want. I have never gotten this type of situation using ww3 3.14? Can you please look at my bathymetry files and input files to let me know how I can get rid of the excluded points. You can see these excluded points in the output log for ww3_grid as shown below: [test.log](https://github.com/NOAA-EMC/WW3/files/11290868/test.log)
EmilioEchevarria commented 1 year ago

Are you using the ww3_grid.inp or ww3_grid.nml? You only need one.

Also, how does your mask file look like? Looking at the output of your ww3_grid execution, it looks like you have 3s around the 180th meridian. 3 means "excluded point", is that what you want?

aliabdolali commented 1 year ago

At least in the grid.nml, I see your closure wrong. It should be SMPL. And the other user us right, do not put both nml and inp in your work directory. You need just one.

JMSnyder1980 commented 1 year ago

Are you using the ww3_grid.inp or ww3_grid.nml? You only need one.

Also, how does your mask file look like? Looking at the output of your ww3_grid execution, it looks like you have 3s around the 180th meridian. 3 means "excluded point", is that what you want? No it is not want I want. I am trying to eliminate those since when I display output I get a big gap at around 180. I do not want that. This gets eliminated when I using SMPL instead of the NONE but then the nested runs no longer work. Please refer to discussion #993 or issue #994 for this problem. How can I run the ww3 nests successfully without dealing with "excluded points" along the 180 boundary for the global run?

JMSnyder1980 commented 1 year ago

At least in the grid.nml, I see your closure wrong. It should be SMPL. And the other user us right, do not put both nml and inp in your work directory. You need just one.

i have put in SMPL but then the boundary conditions for nests do not pass properly. I wind up getting errors with sea points. See issue #994 or discussion #993. How can I resolve this issue?