matlab ncep,Matlab全球风场图—以ERA5、NCEP和ECMWF再分析数据为例-程序员宅基地

技术标签: matlab ncep  

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x

2 M6 B* g4 Y4 D* \! X8 r& L一、matlab读取NCEP再分析数据并绘制风场

/ U- h7 c) Q% g. [4 Z9 g

) k  P( w- N) ?, O" l%该程序用于求水汽通量散度

# z: T5 d6 Y) x3 T9 Y%注意!我们这里读到的u是四维矩阵,分别是lon*lat*level*time,1 e, |( l# I* B5 `: A- p" N6 g

clc;clear;close all

( Q; w5 a5 Y. {# P& S& k7 xf_hgt = 'ps_level_20170121_0130.nc';! k& g! b8 M: ~2 w4 b9 c2 I  K, @

% ncdisp(f_hgt);* X8 }% C4 z( U, p& L

time=ncread(f_hgt,'time');/ m" p4 L9 m1 h2 u5 r8 h

level=ncread(f_hgt,'level');+ j0 E$ W  D* ]& ~: X- z

lon=ncread(f_hgt,'longitude');8 @0 u2 e( _' d0 |

lat=ncread(f_hgt,'latitude');

1 ]" ]2 d/ z7 u5 s7 e%%%%%%时间转换

4 |' [9 T. j+ T3 Ltime  = double(time);

' M4 e1 \  C% z4 }  N- H! eformat = 'mm dd, yyyy HH:MM:SS.FFF AM';%转换格式% y# x6 g# Q$ i! h7 ]

dstr = datestr((datenum('1900-01-01') + time./24),format);%转换后时间字符串存储

5 E/ A6 ]5 z8 B+ ATM = datevec(dstr);%将时间字符数组转化为数值数组

6 E- }9 ^+ Z% ?5 Xtidx=find(TM(:,2)==1 & TM(:,3)==28 & TM(:,4)==00);%筛选7月25日08时(世界时加8)/ ~% m. U( r* K( O- V5 G! Y! w

ps_lev=find(level ==850);%%删选出850hPa高度

6 P. H! B: R( U- O! Y( V; wstart=[1,1,ps_lev,tidx];%所指定变量的每一维的开始读取的位置<

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_39761573/article/details/115826803

智能推荐

【CesiumforUnreal插件】UE5 快速构建Cesium场景 快速入门!!!_cesium 在unreal中建筑物是白色-程序员宅基地

文章浏览阅读972次,点赞2次,收藏4次。文章作者技术和水平有限,如果文中出现错误,希望大家能指正,同时有问题的话,欢迎大家留言讨论。_cesium 在unreal中建筑物是白色

2023软件设计师中级备考经验分享(文中有资料链接分享)_2023软件设计师考试学习资料-程序员宅基地

文章浏览阅读2.6k次,点赞10次,收藏24次。上午题:推荐把软考通APP里的历年真题刷3-4遍,直接刷真题,然后背解析就行,没必要分章节练习。每次考试的真题会有几个往年的题,其他题目的话都是万变不离其宗,把历年真题刷几遍,背好解析。上午题几乎不是问题!下午题:推荐把近5年真题做一遍,把一个up的下午题真题解析课程听1遍就差不多了,2011年-2018年的真题解析课也得听一遍,考试的时候会有往年题出现。B站up主资料分享:【2023下半年中级软件设计师备考资料分享-哔哩哔哩】 https://b23.tv/C2kJjDm。我自己都没想到会这么好的成绩。_2023软件设计师考试学习资料

PS封装码流的解析_码流解析-程序员宅基地

文章浏览阅读4.2k次,点赞4次,收藏35次。1.解析PS码流的流程图2.获取数据2.1 读取源码文件2.2 冗余数据处理3.解析数据3.1 解析PS Header3.2 解析 PS system header(可有可无)3.3 解析PS system map(I帧才会有)案例:3.4 解析PES包3.5 保存裸数据3.6 保存时间戳3.7 保存帧类型..._码流解析

投影仪工作原理细说-程序员宅基地

文章浏览阅读2.7w次。项目中用到了投影系统,本科上课、开组会的时候对投影仪相当的不爽,感觉用着很麻烦。因为要调很多参数,色彩、位置balabala的。然并卵,一旦项目中采用了投影仪方案,呵呵,冤家路窄,好吧,解说一下自己对投影仪原理的认识~投影机的成像原理 投影机是一种用来放大显示图像的投影装置。 说到投影机显示图像的原理,基本上所有类型的投影

ADX345获取加速度在LCD显示-程序员宅基地

文章浏览阅读2k次。ADX345使用IIC读取数据IIC配置方法一:使用IO模拟IIC操作彩屏显示使用之前修改好的ili9325的程序(lcd.c,lcd.h,gui.c,gui.h),在main.c中直接添加lcd_init()初始化函数,就可以使用lcd和gui中的绘图或者显示字符的函数了,不过需要注意要重定向printf输出到串口。IIC使用IO模拟,创建IIC初始化和读写函数ADX34_adx345

KGCN---pytorch代码(2)---aggregator_kgcn代码实现-程序员宅基地

文章浏览阅读334次,点赞4次,收藏7次。将当前节点的向量(self_vectors)与邻居节点的向量(neighbor_vectors)+邻居关系(neighbor_relations)+以及用户嵌入(user_embeddings)+act(激活函数)结合。在得分张量的最后添加一个维度,将其形状从 [batch_size, -1, n_neighbor] 调整为 [batch_size, -1, n_neighbor, 1]结果是一个形状为 [batch_size, -1, n_neighbor] 的张量,表示每个邻居对当前节点的重要性得分。_kgcn代码实现

随便推点

DashScope灵积模型服务 java testcase - 特色功能 模型监督学习_灵积模型声音训练java-程序员宅基地

文章浏览阅读923次,点赞16次,收藏15次。DashScope 测试 和模型比较代码 DashScope灵积模型服务建立在“模型即服务”(Model-as-a-Service,MaaS)的理念基础之上,围绕AI各领域模型,通过标准化的API提供包括模型推理、模型微调训练在内的多种模型服务。_灵积模型声音训练java

ffmpeg调整音频倍速_ffmpeg 音频变速-程序员宅基地

文章浏览阅读581次。ffmpeg命令行调整音频倍速,解释各个参数的含义。_ffmpeg 音频变速

MyBatis拦截器四种类型和自定义拦截器的使用流程_mybatis 拦截器-程序员宅基地

文章浏览阅读1.3k次,点赞35次,收藏17次。MyBatis拦截器四种类型和自定义拦截器的使用流程_mybatis 拦截器

Flutter基础_flutter appbar actions 布局中-程序员宅基地

文章浏览阅读2.4k次。Fluter基础Flutter环境搭建Android Studio开发环境搭建Flutter目录结构,入口函数以及基础组件目录结构文件夹 作用android android 平台相关代码ios ios 平台相关代码lib flutter 相关代码,我们主要编写的代 码就在这个文件夹test 用于存放测试代码pubspec.yaml 配置文件,一般存放一些第三方库的依 赖。入口函数void main_flutter appbar actions 布局中

HDU-1789 Doing Homework again 动态规划 Or 贪心-程序员宅基地

文章浏览阅读58次。http://acm.hdu.edu.cn/showproblem.php?pid=1789题义是给定一个作业序列,求如何分配使得得到的分数最多。设 dp[i][j] 代表截止到第i个作业,第j天所能够完成的最多分数。dp方程为:if (1<= j <= e[i].time) dp[i][j] = max(dp[i-1][j], d[i-1][j-1]+e[i].s..._leetcode doing homework 贪心算法

ThreadX 内部系统时钟服务_threadx如何确定系统时钟-程序员宅基地

文章浏览阅读2.5k次。ThreadX中,有两个函数可以获取和设置内部系统时钟服务,tx_time_get 获取当前时间。tx_time_set 设置当前时间。tx_time_get获取当前时间。原型ULONG tx_time_get(VOID)描述这项服务用来返回当前的内部系统时钟。每个定时节拍使内部系统时钟加1。系统时钟在初始化时被设置为0,同样可以被tx_time_set服务_threadx如何确定系统时钟

推荐文章

热门文章

相关标签