爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 12766|回复: 3

[程序设计] 计算温度平流时出错

[复制链接]

新浪微博达人勋

发表于 2020-4-7 11:20:53 | 显示全部楼层 |阅读模式

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

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

x
lon=Data.lon;
lat=Data.lat;
t500=Data.Temperature(5,:,:);
u500=Data.U(4,:,:);
v500=Data.V(4,:,:);
PI=3.1415926;
R0=6371004;
dx=gradient(lon)*PI/180;
dy=gradient(lat)*PI/180;
T500=squeeze(t500);
U500=squeeze(u500);
V500=squeeze(v500);
[dtx500,dty500]=gradient(T500,1,1);
tp500=-((U500.*dtx500)./(dx.'*cos(lat.*PI/180))'+V500.*(dty500'*diag(1./dy))')/R0;
如上所示,在计算温度平流时最后这里(dx.'*cos(lat.*PI/180)维度不一致,查了很多公式都是这样的,lat是180*1,dx是360*1,所以运行一直报错,各位大神看看

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

新浪微博达人勋

发表于 2020-4-9 09:45:28 | 显示全部楼层
本帖最后由 wangzhaoying 于 2020-4-9 10:08 编辑

LAT和LON,meshgrid一下?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-19 19:23:07 | 显示全部楼层
楼主解决了吗?出现了同样的报错
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-19 21:06:18 | 显示全部楼层
开拓者 发表于 2021-4-19 19:23
楼主解决了吗?出现了同样的报错

http://bbs.06climate.com/forum.p ... &extra=page%3D6
用这个帖子的吧
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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