- 积分
- 632
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
clc;
close all;clear;
mmap_path='C:\Users\youlin\Documents\MATLAB\matlab\m_map\';
%nc_path = 'D:\dt_global_twosat_madt_h_20000315_20140106.nc';
nc_path = 'H:\AVISO\AGC\two_sat_merged\h\2000\dt_global_twosat_madt_h_20000315_20140106.nc';
lon=[121 124];
lat=[24 27];
data=double(ncread(nc_path,'adt'));
lon=double(ncread(nc_path,'lon'));
lat=double(ncread(nc_path,'lat'));
time=double(ncread(nc_path,'time'));
data=data';
pdata=data(435:480,450:540);
%pdata(data>150)=nan;
pdata=pdata*100
[lon lat]=meshgrid(lon,lat);
figure;hold on;warning off;
m_proj('Equidistant cylindrical','lon',[121 124],'lat',[24 27]);
m_grid('linestyle','none','tickdir','out','linewidth',3);
[c h]=m_contour(lat,lon,pdata,5);
shading flat,colorbar
m_gshhs_i('patch',[.5 .5 .5]);
出现error
Error in m_contour (line 37)
data(i)=NaN;
Error in adt (line 29)
[c h]=m_contour(lat,lon,pdata,5);
哪里需要修改,谢谢!
|
|