Few check on trackeddy_V1-0_global_eddies.nc

check_track_eddy_dataset.py

from netCDF4 import Dataset
from matplotlib import pyplot as plt

with Dataset('trackeddy_V1-0_global_eddies.nc') as h:
# with Dataset('http://dapds00.nci.org.au/thredds/dodsC/ks32/CLEX_Data/TrackEddy_eddies/v1-0/netcdf/trackeddy_V1-0_global_eddies.nc') as h:
    # Variable
    time = h.variables["time"][:]
    neddy = h.variables["neddy"][:]
    longitude = h.variables["center_mass_loc_x"][:]
    latitude = h.variables["center_mass_loc_y"][:]
    # Order
    # I didn't understand how data are ordering, i wonder there is a reason
    # that explain why there are not sorting by time or position
    print(time[:10])
    print(neddy[:10])
    print(longitude[:10])
    print(latitude[:10])
    
    # Detection over land
    fig =plt.figure()
    ax = fig.add_subplot(111)
    ax.set_title('Detection over italian land')
    ax.set_aspect('equal'), ax.grid()
    m = (longitude < 15) * (longitude > 10) * (latitude >41) * (latitude < 44)
    ax.plot(longitude[m], latitude[m], '.', markersize=.5)
    
    # 0 longitude missing detection
    longitude = (longitude + 180) % 360 -180
    m = (longitude < 10) * (longitude > -10) * (latitude >-40) * (latitude < -30)
    fig =plt.figure()
    ax = fig.add_subplot(111)
    ax.set_title('No detection at longitude 0°')
    ax.set_aspect('equal'), ax.grid()
    ax.plot(longitude[m], latitude[m], '.', markersize=.5)

    #
    plt.show()