windrot_geo2proj.gr3 missing for surface forcing

[dkobashi]

Hi,

I created surface forcing files (sflux) using pyschism's era5 routine, but had an error when running SCHISM saying windrot_geo2proj.gr3 is missing.
I checked the manual and it looks like windrot_geo2proj.gr3 is needed to run SCHISM as I use epsg:4326 that is based on lon/lat. However, pyschism era5 script did not create the file.

Does pyschism have a function to create windrot_geo2proj.gr3? or does anyone have a working example for Python to create the gr3 file?
It appears that ERA5().write() does not have an option to create windrot_geo2proj.gr3.

Thanks in advance.

-DJ

[dkobashi]

Maybe I was a bit confused.

My coordinate is epsg:4326 for everything including grid. So I assume I don't need windrot_geo2proj.gr3 much like the fact that I did not create hgrid.ll either (just copy hgrid.gr3 to hgrid.ll). From what I checked for Test_ICM_ChesBay test case, the format of windrot_geo2proj.gr3 is pretty much the same as hgrid.gr3 excluding boundary conditions. So I just copy hgrid.gr3 to windrot_geo2proj.gr3 and remove the boundary condition part just like hgrid.ll.

Will see if this works.

Bests,

-DJ

[pmav99]

you need to replace the Z values of the hgrid.gr3 (i.e. the depth) with zeroes, too. There is an example here: https://github.com/schism-dev/pyschism/blob/b8f4bdde9b7bec0b0225a6f9fcbb8673d1d805ce/examples/Gr3/gen_gr3_constant.py