Add ability to subset a .static.nc file
Closed this issue · 1 comments
MiCurry commented
No reason why one shouldn't be able to. It'll also help fix problems with the current static interpolations on convex regions or channel regions.
We should be able to process any mesh and use its sphere radius for marking its boundary points in LimitedArea.mark_boundary; however, it appears that mark_boundary assumes that it is working upon a unit sphere.
For the sphere radius of static files, mark_boundary will fail when it tries to take arccosine of the angle between the current cell and its edges, which will cause k to not be defined.
The following lines in LimitedArea.mark_boundary have no effect:
temp = np.linalg.norm(pta)
cross = pta / temp