根据风场数据,画出剖面和断面海洋流场
剖面
clc;
clear;
dep=load('C:\Users\风轻叶落\Desktop\xu\xiashi\bs.depth');
u_file=dir(['C:\Users\风轻叶落\Desktop\xu\xiashi\u*']);
v_file=dir(['C:\Users\风轻叶落\Desktop\xu\xiashi\v*']);
w_file=dir(['C:\Users\风轻叶落\Desktop\xu\xiashi\w*']);
for i=1:9
u(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\xiashi\',u_file(i).name]);
v(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\xiashi\',v_file(i).name]);
w(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\xiashi\',w_file(i).name]);
end
% dep=load('C:\Users\风轻叶落\Desktop\xibeifeng\bs.depth');
% u_file=dir(['C:\Users\风轻叶落\Desktop\xibeifeng\u*']);
% v_file=dir(['C:\Users\风轻叶落\Desktop\xibeifeng\v*']);
% w_file=dir(['C:\Users\风轻叶落\Desktop\xibeifeng\w*']);
% for i=1:9
% u(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xibeifeng\',u_file(i).name]);
% v(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xibeifeng\',v_file(i).name]);
% w(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xibeifeng\',w_file(i).name]);
% end
row=24; %行 28 28 24
col=10; %列 18 38 10
lat=37.167 +(row-1)/12; %纬度
lon=117.583 +(col-1)/12; %经度
mydep=dep(row,col); %对应经纬度的水深
%sigma坐标下的参数*水深
h1=0*mydep;
h2=-0.021*mydep;
h3=-0.042*mydep;
h4=-0.083*mydep;
h5=-0.167*mydep;
h6=-0.333*mydep;
h7=-0.500*mydep;
h8=-0.667*mydep;
h9=-0.833*mydep;
h10=-1*mydep;
h=[h1,h2,h3,h4,h5,h6,h7,h8,h9];
for i=1:9
myu(i)=u(row,col,i);
myv(i)=v(row,col,i);
myw(i)=w(row,col,i);
plo=quiver3(0,0,h(i),myu(i),myv(i),myw(i),'r');
hold on
plo2=quiver3(0,0,h(1),0,2,0,'b','LineWidth',2);
set(plo2,'maxheadsize',1);%风向不符合实际,前面风场好像也不对?
% set(gca,'xtick',[-6:1:6]);
% set(gca,'ytick',[-6:1:6]);
hold on
plo=quiver3(0,0,h(9),myu(i),myv(i),0,'b');
end
axis([-5 10 -3 4]);
%修饰图
xx=[0,0,0,0,0,0,0,0,0];
yy=xx;
zz=h;
plot3(xx,yy,zz,'-oblack');
plot3(myu,myv,h+myw,'green')
h0(1:9)=h9;
plot3(myu,myv,h0,'green')
for i=1:9
plot3([myu(i),myu(i)],[myv(i),myv(i)],[h9,zz(i)],'--oblack');
end
text(0,2,0,'风向')
text(-3,3,0,num2str(lat,'纬度=%.4f°'))
% text(-3,3,0,'纬度:39.5837°')
text(-3,3,-3,num2str(lon,'经度=%.4f°'))
title('渤海夏季南风流速剖面','FontSize',20);
xlabel('u(cm/s)');
ylabel('v(cm/s)');
zlabel('水深(m)');
断面
clc;
clear;
dep=load('C:\Users\风轻叶落\Desktop\xu\dongshi\bs.depth');
u_file=dir(['C:\Users\风轻叶落\Desktop\xu\dongshi\u*']);
v_file=dir(['C:\Users\风轻叶落\Desktop\xu\dongshi\v*']);
w_file=dir(['C:\Users\风轻叶落\Desktop\xu\dongshi\w*']);
for i=1:9
u(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\dongshi\',u_file(i).name]);
v(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\dongshi\',v_file(i).name]);
w(:,:,i)=load(['C:\Users\风轻叶落\Desktop\xu\dongshi\',w_file(i).name]);
end
row0=1;
row1=46;
for row =row0:row1
% row=43; %行
col=38; %列
lat(row)=37.167 +(row-1)/12; %纬度
lon=117.583 +(col-1)/12; %经度
mydep=dep(row,col); %对应经纬度的水深
%sigma坐标下的参数*水深
h1=0*mydep;
h2=-0.021*mydep;
h3=-0.042*mydep;
h4=-0.083*mydep;
h5=-0.167*mydep;
h6=-0.333*mydep;
h7=-0.500*mydep;
h8=-0.667*mydep;
h9=-0.833*mydep;
h10=-1*mydep;
h=[h1,h2,h3,h4,h5,h6,h7,h8,h9]
for i=1:9
wH=1;
myu(row,i)=u(row,col,i);
myv(row,i)=v(row,col,i);
myw(row,i)=w(row,col,i)*wH;
% myv(row,i)=min(myv);
% myw(row,i)=min(myw);
% plo=quiver(lat(row),h(i),myv(row,i),myw(row,i),0.1,'r');
% plo=quiver(lat(row),-i,myv(row,i),myw(row,i),0.08,'r');
% set(plo,'maxheadsize',2)
% hold on
% plo2=quiver3(0,0,h(1),0,2,0,'b','LineWidth',2);
% set(plo2,'maxheadsize',1);%风向不符合实际,前面风场好像也不对?
% set(gca,'xtick',[-6:1:6]);
% set(gca,'ytick',[-6:1:6]);
% hold on
% plo=quiver3(0,0,h(9),myu(i),myv(i),0,'b');
end
% axis([-3 3 -6 4]);
% set(gca,'ytick',[5:-1:-60]);
end
for i=1:46
for j=1:9
if (myw(i,j)>0&myw(i,j)<1)
myw(i,j)=1;
elseif (myw(i,j)<0&myw(i,j)>-1)
myw(i,j)=-1;
elseif (myw(i,j)>10)
myw(i,j)=10;
elseif(myw(i,j)<-10)
myw(i,j)=-10;
end
end
end
% for i=1:46
% for j=1:9
% if (myv(i,j)>0&myv(i,j)<1)
% myv(i,j)=1;
% elseif (myv(i,j)<0&myv(i,j)>-1)
% myv(i,j)=-1;
% elseif (myv(i,j)>10)
% myv(i,j)=10;
% elseif(myv(i,j)<-10)
% myv(i,j)=-10;
% end
% end
% end
hold on
zS=0.05;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%即使w v在同一量级,x y坐标尺度不对称也会影响效果,调整z轴跨度
for row =row0:row1
col=38; %列
lat(row)=37.167 +(row-1)/12; %纬度
lon=117.583 +(col-1)/12; %经度
mydep=dep(row,col); %对应经纬度的水深
%sigma坐标下的参数*水深
h1=0*mydep;
h2=-0.021*mydep;
h3=-0.042*mydep;
h4=-0.083*mydep;
h5=-0.167*mydep;
h6=-0.333*mydep;
h7=-0.500*mydep;
h8=-0.667*mydep;
h9=-0.833*mydep;
h10=-1*mydep;
h=[h1,h2,h3,h4,h5,h6,h7,h8,h9]
for i=1:9
plo=quiver(lat(row),h(i)*zS,myv(row,i),myw(row,i),0.03,'r');
set(plo,'maxheadsize',1)
end
set(gca,'xtick',[min(lat):1:max(lat)]);
end
hold on;
plot(lat(row0:row1),-dep(row0:row1,col)*zS)
%构造闭合边界
plot([min(lat(row0:row1)),min(lat(row0:row1))],[-max(dep(row0:row1,col))*zS,0])
plot([max(lat(row0:row1)),max(lat(row0:row1))],[-max(dep(row0:row1,col))*zS,0])
plot([min(lat(row0:row1)),max(lat(row0:row1))],[-max(dep(row0:row1,col))*zS,-max(dep(row0:row1,col))*zS])
laa=lat;
laa(1,47:48)=[lat(row1),lat(row0)]; %补充闭合边界的点到矩阵中
dee=dep(row0:row1,col);
dee(47:48)=[max(dep(row0:row1,col)),max(dep(row0:row1,col))];
dee=-dee*zS;
c=[1:48];
patch(laa,dee,c)
% set(gca,'XTick',-3.5:0.5:0,'XTickLabel',{'-60''-50''-40''-30''-20''-10''0'});
set(gca,'YTick',-3.5:0.5:0.5,'YTickLabel',{'-70','-60','-50','-40','-30','-20','-10','0','10'});
% set(gca,'XTickLabel',0:-10:-60);
title('渤海冬季断面流速矢量图 (经度:120.6663°)','FontSize',16);
xlabel('纬度(°N)');
ylabel('深度(m)');
根据风场数据,画出剖面和断面海洋流场相关推荐
- Excel如何直接根据数据画出统计图
给定一组数据,如何来用excel画出其柱状分布图呢? 如下数据: 120 121 122 125 130 132 135 137 138 139 140 150 140 141 125 如上数据根据数 ...
- 爬取NBA球员数据画出球员的雷达图
目录 爬取NBA球员数据并进行清洗 画出两个球员的雷达图 创建时间轴组件自动以时间轴播放 代码整理: 总结 爬取NBA球员数据并进行清洗 这里我是从https://nba.hupu.com/playe ...
- matlab导入数据画二维云图,matlab中用xyz三组数据画出2维云图
各位Matlab的高手,你们好. 我最近遇到一个技术问题,久未能解决,还望知道的能指教下,谢谢.我的问题是, 我有三组数据xyz.其中,xy分别为二维面的坐标,z为对于坐标的温度值.由于试件发生了变形 ...
- 通过导入txt数据画出python turtle图形_【Python】txt文件读取绘画
import turtle as t #绘画库 t.title('自动轨迹绘制') #标题 t.setup(800,600,0,0) t.pencolor('red') t.pensize(5) #数 ...
- python画出K线图及技术指标
目录 安装mplfinance及Ta-lib 画图 安装mplfinance及Ta-lib mplfinance是基于matplotlib的金融数据可视化分析模块,前身是mpl_finance,比起m ...
- Arcgis画出断面线并批量导出
ARCGIS生成河道断面文件 (1)准备好河道中心线,高程点,面文件 (2)利用高程点插值为TIN或者DEM文件: (3)按照以下步骤点击工具条(3D Analyst) 第一个按钮画断面线,第二个按钮 ...
- 【解决办法】pandas画出时序数据(股票数据)横轴不是时间
简述 遇到了这个问题,被坑了很久. 首先我们要假设我们一直认为index是时间数据.然后我们发现没有看到横轴为时间 (如果不是的这么认为的话,就记得先把index设置为时间数据) 可能性 遇到这个问题 ...
- matlab导入txt生成曲面,求助:怎样将txt的数据导入到Matlab中并根据参数画出波形...
导师让我把波形画出来,这下完全不会啊!希望各位高手帮帮在下-- txt中部分数据如下: ID DDD HH:MM:SS.mmmuuun PARA1 PARA2 CH RISE ...
- df python 增加数据_美国确诊超100万!教你用Python画出全球疫情动态图(附数据下载)...
CDA数据分析师 出品 [导语]:今天我们教你用Python绘制全球疫情动态图,技术部分请看第二部分. 获取数据:扫描下方公众号回复关键字"全球疫情" 目前,全球新冠疫情还十 ...
最新文章
- JS设计模式——3.封装与信息隐藏
- oracle23290,oracle标题
- Cortex-M3中的寄存器组
- “敏捷版”全链路压测
- boost实用工具:assign库了解学习
- 在RelativeLayout中对控件的显示隐藏操作问题
- 凭什么,Netty能成为最流行的NIO框架?
- 无心剑中译马塞尔·普鲁斯特《追忆似水年华》
- 2021牛客寒假算法基础集训营3,签到题DGHIJ
- Repository HDU - 2846 (trie)
- 关于scrollTop
- mysql用身份证号判断男女_如何根据身份证号码辨别性别呢
- jtag接口定义 jtag接口的主要作用是什么
- Python入门学习笔记(三)
- 网页图片透明html代码,HTML_让css使网页图片半透明,你用过css么?当然,我是指你 - phpStudy...
- 关于嵌入式的bin、hex、axf、map
- java使用jacob将word,excel,ppt转成html
- 想调节—下父元素的透明度,但是又不影响子元素的透明度,怎么破?
- SpringBoot JPA多对一 持久化是报错object references an unsaved transient instance - save the transient instanc
- 阿里云实人认证 人脸识别