请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 129398|回复: 146

[源程序] [原创]MATLAB流函数与势函数的数值计算程序

  [复制链接]

新浪微博达人勋

发表于 2014-3-2 20:12:40 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 沙颖凯 于 2014-11-4 22:27 编辑

MATLAB环境下流函数与势函数的数值计算程序,使用Lieberman迭代法,因为是迭代计算,为了节约计算资源,边界条件取0,所以边界值不可靠,被我用NaN覆盖了。

% [psi Upsi Vpsi]=psi_streamfunction(longitude, latitude, u, v);
% psi: Stream Fuction
% Upsi, Vpsi: Vorticity Wind
% [chi Uchi Vchi]=chi_potential(longitude, latitude, u, v);
% chi: Velosity Potential
% Uchi, Vchi: Divergence Wind

之前在[http://bbs.06climate.com/forum.php?mod=viewthread&tid=20665]贴过一些矢量场的数值计算的程序,这里的程序调用了里面计算涡度、散度、梯度和空间差分的子程序,所以要和上次的程序放在一起调用。

附一张流函数的测试图(500hPa)。

Test_Stream_Function

Test_Stream_Function



另外,我在[Endlich R. M., 1967: An iterative method for altering the kinematic properties of wind fields. J. Appl. Meteor., 6, 837--844]找到了一种同时计算流函数和势函数的逐步订正方法,没有按解Poission方程的思路求解。但我没有尝试成功。附件"Endlich"是失败的作品,目前不打算继续深究,不知道是否有人能了解或改进这种算法。
psi_streamfunction.m (1.7 KB, 下载次数: 635)

评分

参与人数 3金钱 +25 贡献 +1 收起 理由
冰雨星 + 8 很给力!
Minoz_me + 1 很给力!
wlzhongouc + 16 + 1

查看全部评分

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

新浪微博达人勋

发表于 2014-3-4 00:33:07 | 显示全部楼层
嗯,是这样。可以加一个判断既可以了:
lon1=lon(:,1);lat1=lat(:,1);
LLon=length(lon1);
【M,N】=size(lon);

if M==Llon,
lons=lons';lats=lats';
end
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2014-3-3 16:03:35 | 显示全部楼层
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2014-3-2 21:02:37 | 显示全部楼层
好东西抢沙发
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-2 22:05:06 | 显示全部楼层
感谢分享
有案例数据吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-2 22:08:45 | 显示全部楼层
kongfeng0824 发表于 2014-3-2 22:05
感谢分享
有案例数据吗?

NCEP/NCAR再分析数据里提取的

Mar_21_2009_00UTM.mat (71.12 KB, 下载次数: 141)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 09:14:11 | 显示全部楼层
牛 多谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 09:22:27 | 显示全部楼层
感谢分享!!!{:5_213:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-3 13:13:44 | 显示全部楼层
之前上传的势函数程序有错误,刚才更正了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 14:06:47 | 显示全部楼层
感谢分享 这两天正研究这个呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 14:07:22 | 显示全部楼层
是通过求解poisson方程算的 吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 14:07:34 | 显示全部楼层
是通过求解poisson方程算的 吗
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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