Coordinate conversion using run_telfile.py
The problem:
In TELEMAC v8p4, coordinate conversion using run_telfile.py doesn't work properly converting from lat-long.
(e.g. run_telfile.py alter --proj EPSG:4326/EPSG:27700 <LL.slf> <OSGB36.slf>)
This is because the pyproj transformer assumes default lat-long coordinates are ordered latitude, longitude (i.e. y,x)
The fix:
In pretel/manip_telfile.py, edit the following line
meshx, meshy = pyproj.transform(in_proj, out_proj, meshx, meshy)
Add the option always_xy=True:
meshx, meshy = pyproj.transform(in_proj, out_proj, meshx, meshy, always_xy=True)