sbrus89 / ww3_utils Goto Github PK

Hi!

I have been checking the file plot_grid.py in order to recreate the plots resulting from the joint use of OceanMesh2D (which provides the general mesh file) and AlphaBetaLab (which provides the obstruction files).

I have seen that plot_grid.py works nicely for the "bot" and "mask" files. However, the two "obstruction" files that anyone can get from AlphaBetaLab have a very different format that cannot be read by np.loadtxt(). For instance, they typically start with things like:

$WAVEWATCH III SHADOW OBSTRUCTIONS
18015
$ ilon ilat of the cell. lon: -0.0, lat: -43.0
1   95
$ sizes of the cell in km
    40.628       44.286       49.408       56.200       64.946       66.929       64.335       61.661       58.470       56.605       60.148       63.041       65.602       68.436       60.299       52.583       46.650       42.290       42.290       46.650       52.583       60.299       68.436       65.602       63.041       60.148       56.605       58.470       61.661       64.335       66.929       64.946       56.200       49.408       44.286       40.628
$ mean alpha: 0.9371637038137896
$ mean beta: 0.9448303374500736
$alpha by ik, ith

In addition, they have many other lines and parts along the file that are not directly readable, like those of the alpha and beta ik, ith.

Then, the question is, are you using a different tool or format to represent the obstructions file with plot_grid.p?

Thank you so much!