- 积分
- 237
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-10-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
matlab 利用SODA数据绘制南中国海SCS四季流场图
download SODA 2015年数据,读取其中的u,v方向的current矢量,在选取目标海域绘制流场图
这里选取了[4 7 10 1]月作为代表月,展示spring,summer,autumn,winter 流场
数据太大,需要自己下载哦
- clc; clear; close;
- path = '/xx]);
- lat = ncread([path,'soda2015.nc'],'latitude');
- lon = ncread([path,'soda2015.nc'],'longitude');
- u = ncread([path,'soda2015.nc'],'u');
- v = ncread([path,'soda2015.nc'],'v');
- u_surface = squeeze(u(:,:,1,:));
- v_surface = squeeze(v(:,:,1,:));
- u_surface(u_surface == -1.000000020040877e+20) = nan;
- v_surface(v_surface == -1.000000020040877e+20) = nan;
- U_surface = sqrt(u_surface.^2 + v_surface.^2);
- lonmin = 98;
- lonmax = 130;
- latmin = 0;
- latmax = 30;
- ix = find(lon>lonmin & lon<lonmax);
- iy = find(lat>latmin & lat<latmax);
- tlon = lon(ix);
- tlat = lat(iy);
- u_sf = u_surface(ix,iy,:);
- v_sf = v_surface(ix,iy,:);
- U_sf = U_surface(ix,iy,:);
- %U_sf(U_sf = inf) = nan;
- u_sf(u_sf < -1.0000000e+2) = nan;
- v_sf(v_sf < -1.0000000e+2) = nan;
- %U_sf(U_sf < 0) = nan;
- U_sf(U_sf > 100) = nan;
- U_sp = U_sf(:,:,4);U_sm = U_sf(:,:,7);
- U_fl = U_sf(:,:,10);U_wt = U_sf(:,:,1);
- [xx,yy] = meshgrid(tlon,tlat);
- size = 1;
- %% [spring summer fail winter 1 4 7 10]
- close
- figure
- subplot(221)
- contourf(xx,yy,U_sp',66,'LineColor','none'); colorbar;caxis([0 1]);
- mycolor=zeros(3,202);
- mycolor(1,:)=[0:0.01:1 zeros(1,101)+1];
- mycolor(2,:)=[0:0.01:1 1:-0.01:0];
- mycolor(3,:)=[zeros(1,101)+1 1:-0.01:0];
- mycolor = mycolor';
- colormap(mycolor);
- grid on;
- xlabel('longitude(E)','FontSize',14);ylabel('latitude(S)','FontSize',14) ;
- set(gca,'fontsize',14);
- hold on;
- quiver(xx,yy,u_sf(:,:,4)'./size,v_sf(:,:,4)'./size,0,'color','k')';
- set(gca,'xtick',lonmin:5:lonmax,'fontsize',14);
- set(gca,'ytick',latmin:5:latmax,'fontsize',14);
- hold on;
- text(lonmin+2,latmax-2,'spring','fontsize',14);
- subplot(222)
- contourf(xx,yy,U_sm',66,'LineColor','none'); colorbar;caxis([0 1]);
- mycolor=zeros(3,202);
- mycolor(1,:)=[0:0.01:1 zeros(1,101)+1];
- mycolor(2,:)=[0:0.01:1 1:-0.01:0];
- mycolor(3,:)=[zeros(1,101)+1 1:-0.01:0];
- mycolor = mycolor';
- colormap(mycolor);
- grid on;
- xlabel('longitude(E)','FontSize',14);ylabel('latitude(S)','FontSize',14) ;
- set(gca,'fontsize',14);
- hold on;
- quiver(xx,yy,u_sf(:,:,7)'./size,v_sf(:,:,7)'./size,0,'color','k');
- set(gca,'xtick',lonmin:5:lonmax,'fontsize',14);
- set(gca,'ytick',latmin:5:latmax,'fontsize',14);
- hold on;
- text(lonmin+2,latmax-2,'summer','fontsize',14);
- subplot(223)
- contourf(xx,yy,U_fl',66,'LineColor','none'); colorbar;caxis([0 1]);
- mycolor=zeros(3,202);
- mycolor(1,:)=[0:0.01:1 zeros(1,101)+1];
- mycolor(2,:)=[0:0.01:1 1:-0.01:0];
- mycolor(3,:)=[zeros(1,101)+1 1:-0.01:0];
- mycolor = mycolor';
- colormap(mycolor);
- grid on;
- xlabel('longitude(E)','FontSize',14);ylabel('latitude(S)','FontSize',14) ;
- set(gca,'fontsize',14);
- hold on;
- quiver(xx,yy,u_sf(:,:,10)'./size,v_sf(:,:,10)'./size,0,'color','k');
- set(gca,'xtick',lonmin:5:lonmax,'fontsize',14);
- set(gca,'ytick',latmin:5:latmax,'fontsize',14);
- hold on;
- text(lonmin+2,latmax-2,'autumn','fontsize',14);
- subplot(224)
- contourf(xx,yy,U_wt',66,'LineColor','none'); colorbar;caxis([0 1]);
- mycolor=zeros(3,202);
- mycolor(1,:)=[0:0.01:1 zeros(1,101)+1];
- mycolor(2,:)=[0:0.01:1 1:-0.01:0];
- mycolor(3,:)=[zeros(1,101)+1 1:-0.01:0];
- mycolor = mycolor';
- colormap(mycolor);
- grid on;
- xlabel('longitude(E)','FontSize',14);ylabel('latitude(S)','FontSize',14) ;
- set(gca,'fontsize',14);
- hold on;
- quiver(xx,yy,u_sf(:,:,1)'./size,v_sf(:,:,1)'./size,0,'color','k');
- set(gca,'xtick',lonmin:5:lonmax,'fontsize',14);
- set(gca,'ytick',latmin:5:latmax,'fontsize',14);
- hold on;
- text(lonmin+2,latmax-2,'winter','fontsize',14);
复制代码
|
|
|
-
|