creating a fieldset with POP

[norafried]

I'm trying to create a fieldset with POP but I have problems with the B-Grid that POP has. I tried to follow the indexing example and adjust that to POP using from_pop (https://nbviewer.org/github/OceanParcels/parcels/blob/master/parcels/examples/documentation_indexing.ipynb) but I keep having issues with that. it always returns:

ValueError: On a B-grid, the dimensions of velocities should be the (top) corners of the grid cells, so the same for U, V and W

Is there an example for POP?

[erikvansebille]

@pdnooteboom might have some code to help you, as he has worked with pop data in his PhD work. Or you may have a look at the code of this paper at https://github.com/pdnooteboom/PO_res_error

[pdnooteboom]

Yes, you could check out specifically this part of the code, which uses a curvilinear grid with 1 degree horizontal resolution, or this part, which uses 0.1 degree horizontal resolution, tripolar grid with partial bottom cells.
I hope this helps.

[norafried]

Great, thanks. I think I found (one of) the errors. I can load the fieldset for one time step now. I have one file per time step containing all the variables. I tried to load them the same was as in the NEMO 3D tutorial. Not sure why that doesn't work. Is there any additional trick? Thanks for all your help!