爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5773|回复: 6

[程序设计] 读取nc文件风场数据以及画风场和风应力旋度

[复制链接]

新浪微博达人勋

发表于 2015-8-21 14:23:07 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
首先读取nc文件里的风场的x、y,再画风应力矢量图,最后生成风应力旋度等值面图为底。
addpath D:\datatrain
ncid=netcdf.open('AVHRR_SST.nc','NOWRITE');%打开nc文件
ncdisp('AVHRR_SST.nc');
TimeData=ncread('AVHRR_SST.nc','TIME');
LonData=ncread('AVHRR_SST.nc','LON3186_3755');
LatData=ncread('AVHRR_SST.nc','LAT1081_1309');
SSTData=ncread('AVHRR_SST.nc','SST');
[X,Y]=meshgrid(LatData,LonData);
contourf(Y,X,SSTData(:,:,9));%画9月等值面图
m_proj('miller','lat',[0,25],'lon',[95,125]);
m_contourf(x,y,w,[0:0.1:11]);
shading flat;
hold on;
h=m_quiver(x,y,u1,v1,0);
colorbar
m_coast('patch',[.6 .6 .6],'edgecolor','none');
m_grid('box','fancy','linestyle','none');

评分

参与人数 1金钱 +6 贡献 +4 收起 理由
Aires + 6 + 4

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-12-27 16:50:51 | 显示全部楼层
m_contourf(x,y,w,[0:0.1:11]); 我知道这一行你要画风应力旋度,但是程序压根就没有代码计算出风应力旋度!
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-8-21 18:00:13 | 显示全部楼层
赞一个,好东西
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-9-29 18:19:44 | 显示全部楼层
没看明白,不知道风应力旋度是哪个变量……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-10 10:31:45 | 显示全部楼层
学习了,谢谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-3 17:01:15 | 显示全部楼层
就想问问您w怎么算的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-3-31 15:40:26 | 显示全部楼层
谢谢分享!!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表