cuill

@yangz888 Hi Zizang, You can find the revised script here (on Frontera): /work2/07933/cuivims/ICOGS3D/Forecast/PostP/NewIO/gen_geojson.py The command line to run it is the same as the old one. The only difference is...

@yangz888 Thanks for the updates.

@erin-nagel @saeed-moghimi-noaa @yangz888 The python script Zizang is using is the same as VIMS' version. So there are no changes to these files. Maybe check if the file size is...

@BahramKhazaei-NOAA I tested with Pacific hgrid. Both ranges generated the same bctides. Would you mind sharing your hgrid?

@BahramKhazaei-NOAA It looks like the conversion from [0, 360) to [-180, 180) is wrong. It should be: `lon[lon>180] -= 360 `

@BahramKhazaei-NOAA We mainly use this Fortran script to generate vqs vgrid for STOFS3D: https://github.com/schism-dev/schism/blob/master/src/Utility/Pre-Processing/STOFS-3D-Atl-shadow-VIMS/Pre_processing/Vgrid/gen_vqs.f90 pylib has a script to generate this type of vgrid: https://github.com/wzhengui/pylibs/blob/master/pyScripts/gen_vqs.py

The downloaded data was saved to a temporary folder created by `tempfile.TemporaryDirectory()`. After generating sflux files, this folder will be deleted. Downloading ERA5 has been slow in the last couple...

This is something wrong with the boundary information. Is it possible to share the hgrid.gr3 file?

The hgrid is missing information of exterior land boundary/island for land boundaries: ``` 16 = number of land boundaries (including islands) 1743 = Total number of land boundary nodes 753...

@calquigs This line `npts, ibtype = map(int, buf.readline().split()[:2])` takes header information for each land or island boundary, for example, **line 848638** in your hgrid.gr3, and assigns the first two strings...