卫星高度计.zip_卫星高度_卫星高度计_波浪_波浪场_风场

clc;clear; % addpath(genpath('F:\MATLAB学习资料\m_map1.4\m_map.')); lon1=ncread('201612__300km.nc','longitude'); lat1=ncread('201612__300km.nc','latitude'); lonmin=119; lonmax=125; in1=find(lon1>=lonmin & lon1<=lonmax); latmin=30; latmax=36; in2=find(lat1>=latmin & lat1<=latmax); lon_space1 = lonmin:2:lonmax; lat_space1 = latmin:2:latmax; [lat1,lon1]=meshgrid(lat_space1,lon_space1); lon_space = lonmin:1:lonmax; lat_space = latmin:1:latmax; [lat,lon]=meshgrid(lat_space,lon_space); Y=2016; figure x0=0.04;y0=0.05;wd=0.23;ht=0.28; xgap=0.00;ygap=0.03; for mon=1:12 if(mon<=4) xp=x0+(wd+xgap)*(mon-1);yp=y0+(ht+ygap)*2; elseif(mon<=8 && mon>4) xp=x0+(wd+xgap)*(mon-5);yp=y0+(ht+ygap)*1; else xp=x0+(wd+xgap)*(mon-9);yp=y0; end subplot('position',[xp,yp,wd,ht]); nc=[num2str(Y), num2str(mon,'%02d'),'__300km.nc']; mean_wsp=ncread(nc,'mean_wsp'); wsp=squeeze(mean_wsp); h=wsp(in1,in2); hh=interp2(lat1,lon1,h,lat,lon) m_proj('mercator','lon',[119 124],'lat',[31 36]); hold on; m_pcolor(lon,lat,hh);shading interp;colorbar; hold on colormap(jet); caxis([3 10]); colorbar('delete'); % set(get(colorbar,'ylabel'),'String','Wind Speed (m s^-^1)',... % 'FontSize',8,'FontWeight','bold'); hold on m_grid('box','on','YaxisLocation','left','YLim',[31 36],'YTick',31:1:36,... 'XaxisLocation','bottom','XLim',[119 124],'XTick',119:1:124,'fontsize',7); m_gshhs_i('patch',[.7 .7 .7]); % m_coast('patch',[.7 .7 .7]); grid on; % m_text(-120.5,35.0, 'Southern California','FontSize',12,'fontweight','bold','color','k'); % hold off; % title(strcat('Altimate',num2str(Y),num2str(m,'%02d')), 'FontSize',14,'fontweight','bold'); % print('-dtiff','-r300',[ 'Altimate_', num2str(Y), num2str(m,'%02d'), '.tif']); % close all end % figure % m_proj('mercator','lon',[-125.5 -116.5],'lat',[29.2 36.5]); % hold on; % color_list = linspace(0,5,10); % [h_color,h_handle] = m_contourf(lon,lat,hs,color_list); % set(h_handle,'LineWidth',1.5) % colorbar('fontsize',14,'fontweight','bold'); % caxis([color_list(1) color_list(end)]); % colormap('jet'); % clabel(h_color,h_handle,'FontSize',15); % m_grid('box','on','Fontsize',14,'fontweight','bold','tickdir','in','XaxisLocation','bottom','YaxisLocation','left'); % m_gshhs_i('patch',[.7 .7 .7]); % % m_coast('patch',[.7 .7 .7]); % grid on; % m_text(-120.5,35.0, 'Southern California','FontSize',12,'fontweight','bold','color','k'); % hold off; % title(strcat('Altimate',num2str(Y),num2str(m,'%02d')), 'FontSize',14,'fontweight','bold'); % print('-dtiff','-r300',[ 'Altimate_', num2str(Y), num2str(m,'%02d'), '.tif']); % close all