GSU not precise anough to find nearest point for active boundaries in high resolution grid #1516
Description
Describe the bug
The GSU is not able to find the right nearst point for active boundaries when grid is at a resolution of 30meters.
This due to the EPS value which is not small enough
To Reproduce
providing a nest.ww3 and running ww3_shel, the program will crash wit the error message :
*** WAVEWATCH III ERROR IN W3IOBC : POINT 'X' 'Y' NOT ACTIVE SEA POINT (1) ***
It happens when calling W3GFCL which is itself called by W3GFPT in w3iobcmd.F90 :
IF ( W3GFPT( GSU, XBPI(I), YBPI(I), IX, IY, DCIN=0.1, DEBUG=.FALSE. ) ) THEN
Expected behavior
It should be able to find the right cell. To make it works, we need to specify EPS=1.0E-8 instead of the default value which is 1.0E-6
IF ( W3GFPT( GSU, XBPI(I), YBPI(I), IX, IY, EPS=1.0E-8, DCIN=0.1, DEBUG=.FALSE. ) ) THEN