In [1]:

    

from shapely.geometry.polygon import Polygon
from shapely.geometry import MultiPolygon
import pyugrid
import json


    

In [2]:

    

#url = 'http://www.smast.umassd.edu:8080/thredds/dodsC/FVCOM/NECOFS/Forecasts/NECOFS_GOM3_FORECAST.nc'
url = 'http://www.smast.umassd.edu:8080/thredds/dodsC/FVCOM/NECOFS/Forecasts/NECOFS_FVCOM_OCEAN_MASSBAY_FORECAST.nc'


    

In [3]:

    

ug = pyugrid.UGrid.from_ncfile(url)


    

In [ ]:

    

lon = ug.nodes[:,0]
lat = ug.nodes[:,1]
nv = ug.faces[:]


    

In [ ]:

    

mp = MultiPolygon([Polygon(zip(lon[element],lat[element])) for element in nv])


    

In [ ]:

    

with open('ugrid.json','w') as f:
    f.write(json.dumps(mp))