【matlab程序】画海洋流场

clear;clc;
file =( '0227.nc');
lat=ncread(file,'latitude');
lon=ncread(file,'longitude');
u=ncread(file,'water_u');
v=ncread(file,'water_v');
x,y=meshgrid(lon,lat);
x=x'; y=y';
interval=4;
figure (1)
set(gcf,'color',1 1 1)
set(gcf,'position',150 300 500 500)
set(0,'defaultfigurecolor','w')
m_proj('Miller','lon',lon(1) lon(end),'lat',lat(1) lat(end));
m_gshhs_i('patch',0.6 0.6 0.6,'edgecolor','k');
m_grid('linest','-.','color','k');
hold on
m_quiver(x(1:interval:end,1:interval:end),y(1:interval:end,1:interval:end),u(1:interval:end,1:interval:end),v(1:interval:end,1:interval:end),0,'k');
title('海洋流场')
print('-djpeg','-r400','海洋流场.jpg');