【提供代码获取方式】matlab使用风羽法画大气环流

以下部分代码：缺少关键性代码和函数！%% 本脚本利用wgrib2先把grib2文件转换为nc文件；在进行读取；画出位势高度场
% 步骤如下:
% 1 下载wgrib2并“安装”（为什么这里的“安装”打了引号呢？因为严格意义上来wgrib2并不需要安装，只需要在环境变量中进行设置即可）
% （1）下载链接：https://pan.baidu.com/s/1EFOBbGvHMPFg-iZS25zAwg 提取码：8bwn
% （2）wgrib2设置windows下wgrib/wgrib2和python读取grib2数据_meryl2001的博客-CSDN博客。
% 2 在cmd中用wgrib2软件将grib2文件转成nc文件
% （1）同时按下windows和R键，输入cmd，回车进入这个黑黑的界面。
% (2）转换命令：wgrib2 grib2file -netcdf ncfile。（grib2file为grib2文件的名字，ncfile为转换出的nc文件的名字）%%代码获取方式：见微信公众号
%% 公众号：海洋与大气科学
%% 作者mapm
clear;clc;close all;% ty_name 指的是台风号；year指的是台风发生的年份；
% month指的是台风发生的月份；%area台风所在范围；% inter1、2 是画图xlabel间隔；
ty_name='1515';year='2015';month='08';area=[90 180;-40 40];inter1=10;inter2=5;
days=[16;27];
%% 读取台风真实路径；
%% 本脚本读取台风路径：2018年24号”潭美“台风：
%% read txt
filename=['20151515.txt'];
fid=fopen(filename,'r');%open file
eng1='000000000';
row=0;%row at start is zero.
while ~feof(fid) %if eixt file ,then runrow=row+1; % read file according rowa=fgetl(fid);% data give a;num=length(a);eng1(row,1:num)=a;
end
fclose(fid);%close file ;
clear row a ans fid filename num;% 选取自己需要的；
%%
end1=eng1(2:end,:);
clear eng1;
time=str2num(end1(:,7:10));% read typhoon's time ;
tylon=str2num(end1(:,18:21))./10;% read typhoon's lon ;
tylat=str2num(end1(:,14:16))./10;% read typhoon's lat;
time1=str2num(end1(:,9:10));% read typhoon's time's hour
time2=str2num(end1(:,7:8));% read typhoon's time's day;
%% 对时间日期处理一下，补齐0；例如：1号补齐为01；小时6补齐为06；并且
%构造新的time_new,为完整日期，例如0706
for time1_i=1:length(time1)time111=time1(time1_i);time222=time2(time1_i);if length(num2str(time111))==1time1_new(time1_i,:)=['0',num2str(time111)];elsetime1_new(time1_i,:)=num2str(time111);endif length(num2str(time222))==1time2_new(time1_i,:)=['0',num2str(time222)];elsetime2_new(time1_i,:)=num2str(time222);end
end
time_new=[time2_new time1_new];
clear typhoon_name end1  time time1 time111 time1_i time1_new time2 time222 time2_new
%% 下面是fnl
filename='fnl_20150816_00_00.nc';
%% 读取风速和位势高度；500hpa;
lon = double(ncread(filename,'longitude'));
lat = double(ncread(filename,'latitude'));
GHT_500=double(ncread(filename,'HGT_500mb'));
U_500=double(ncread(filename,'UGRD_500mb'));
V_500=double(ncread(filename,'VGRD_500mb'));
%% 区域；
ln=find(lon>=area(1,1)&lon<=area(1,2));
la=find(lat>=area(2,1)&lat<=area(2,2));
lonx=lon(ln);laty=lat(la);
[x,y]=meshgrid(lonx,laty);
x=x';y=y';
%% 把区域变量找出来；
HGT=GHT_500(ln,la);
U=U_500(ln,la);
V=V_500(ln,la);
%%
% 读取颜色；
% colormore_1-49系列获取，回复”morecolor“
color=colormore_46(100);
num=46;
fig_h=figure;
set(gcf,'position',[50,50,1000,800],'color','w')
m_proj('miller','lon',[area(1,1) area(1,2)],'lat',[area(2,1) area(2,2)]);
hold on
[c,h]=m_contour(x,y,HGT,[5700:40:5920],'linecolor','k','ShowText','on','LineWidth',1.5);
clabel(c,h,'fontsize',14,'fontweight','bold','LabelSpacing',1000);
hold on
[c,h]=m_contour(x,y,HGT,[5880 5880],'linecolor','k','ShowText','on','LineWidth',1.5);
clabel(c,h,'fontsize',14,'fontweight','bold','LabelSpacing',1000);
hold on
ty_lon=tylon(9);ty_lat=tylat(9);
m_plot(ty_lon,ty_lat,'o','MarkerSize',10,'MarkerFaceColor','k','MarkerEdgeColor','k')%LineWidth——指定线宽
m_coast('linewidth',1.5,'color',[0.6,0.6,0.6]);
m_grid('ytick',[area(2,1):inter2:area(2,2)],'fontsize',16,'fontweight','bold','linewidth',3);
set(gca,'xtick',[area(1,1):inter1:area(1,2)],'xticklabel',[area(1,1):inter1:area(1,2)],'linewidth',2);
export_fig(fig_h,['colormore_',num2str(num),'大气环流彩图bar'],'-dpng','-r300');
close all

【提供代码获取方式】matlab使用风羽法画大气环流相关推荐

海神之光Matlab王者学习之路---代码获取方式
麦哥个人简介及代码获取方式
Bilibili视频主页:https://space.bilibili.com/1702212068
联邦学习算法介绍-FedAvg详细案例-Python代码获取
联邦学习算法介绍-FedAvg详细案例-Python代码获取一.联邦学习系统框架二.联邦平均算法(FedAvg) 三.联邦随梯度下降算法 (FedSGD) 四.差分隐私随联邦梯度下降算法 (DP- ...
java 读取流的字符编码格式_如何使用Java代码获取文件、文件流或字符串的编码方式...
标签: 今天通过网络资源研究了一下如何使用Java代码获取文件.文件流或字符串的编码方式,现将代码与大家分享: package com.ghj.packageoftool; import info.m ...
【Microsoft Azure 的1024种玩法】七十二.使用NodeJS手动编写调用Azure Cognitive Services Translator API代码获取其支持翻译的语言集
[简介] Azure Cognitive Services Translato 是一种基于云的神经网络机器翻译服务,同时也是 REST API 的 Azure 认知服务系列的一部分. 本篇文章将会主要 ...
Matlab：序列分析法MATLAB代码
Matlab:序列分析法MATLAB代码目录输出结果设计代码输出结果更新-- 设计代码 ###下面所有带代码中的n值需要以自己输入的数据为准###1.简单一次滑动平均法预测MATLAB程序代 ...
基于IMAGE法的房间回响模型创建、C++代码实现、matlab仿真
基于IMAGE法的房间回响模型创建.C++代码实现.matlab仿真 1.模型简介 \qquad在处理声音信号时,我们要对信号先进行采集.那么我们就必须要有,一个发出声音的声源,一个进行声音采集的传感 ...
MATLAB代码：基于雨流计数法的源-荷-储双层协同优化配置
MATLAB代码:基于雨流计数法的源-荷-储双层协同优化配置关键词:双层规划雨流计算法储能优化配置参考文档:<储能系统容量优化配置及全寿命周期经济性评估方法研究>第三章仿真平台: ...
CSDN代码和积分获取方式
主要有以下几种获取csdn代码方式,可私信博主了解更多代码或积分获取方式常规方式获取可用分 1.每天只要回复就可以获得10个可用分.注:回复后的第2天发放. 2.每周回复量大于10个帖子,将获得30 ...
黎明觉醒火种测试服务器维护,黎明觉醒火种测试什么时候上线黎明觉醒火种测试资格获取方式(图文)...
黎明觉醒是腾讯旗下的多人开放世界生存手游,对标的就是网易旗下的明日之后.在之前的曙光测试之后,这款游戏长时间来都没有传出过新消息,下面game234就来介绍一下黎明觉醒最新的火种测试什么时候上线,怎么 ...

【提供代码获取方式】matlab使用风羽法画大气环流

【提供代码获取方式】matlab使用风羽法画大气环流相关推荐

最新文章

热门文章