Pyfvcom表面示例 - pyfvcom_surface_example

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt

from cmocean import cm
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.axes_grid1 import make_axes_locatable

from PyFVCOM.read import ncread

# Load the model output.
fvcom = 'sample.nc'

# Extract only the first 20 time steps.
dims = {
   'time': ':20'}

# And only these variables
varlist = ('lon', 'lat', 'nv', 'zeta', 'Times')

FVCOM = ncread(fvcom, vars=varlist, dims=dims, noisy=True)

得到：

File format: NETCDF4
Found nprocs 
Found partition 
Found x 
Found y 
Found lon (extracted [0:73683])
Found lat (extracted [0:73683])
Found xc 
Found yc 
Found lonc 
Found latc 
Found siglay 
Found siglev 
Found siglay_center 
Found siglev_center 
Found h_center 
Found h 
Found nv (extracted [0:<