爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 11212|回复: 10

[分享资料] 涡度计算法

[复制链接]

新浪微博达人勋

发表于 2012-8-13 14:34:40 | 显示全部楼层 |阅读模式

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

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

x
各位大侠 下面的涡度计算 怎么用grads写哦 非常感谢! 1.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-13 15:15:43 | 显示全部楼层
你这样的公式让我想了网格插值积分的程序,大批量的数据计算需要循环,这样算出来的数存到一个变量不好实现。感觉用fortran数组循环处理起来更容易实现。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-13 21:19:26 | 显示全部楼层
有点象数值分析里的插值方法..
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-17 11:07:21 | 显示全部楼层
本帖最后由 waylee2008 于 2012-8-17 11:16 编辑

木头版主 这个olr.mon.mean.nc涡度怎么用grads写
我的想法是 固定中心点 用4个插值循环(x+1,y),(x-1,y) ,(x,y-1),(x,y+1)   再赋变量计算
这样能行吗
问题是插值有重复 不知道对不对
纠结中 不知道怎么写版主 帮帮忙!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-17 14:06:11 | 显示全部楼层
你的公式求SE,SE后面的坐标是怎么样的我看不清楚 ,公式能手写不要截图吗?循环的话我感觉两个循环就够了,x、y格点变化,不过每次输出一个点(SE(X,Y))的值,这样没有办法放数组里,只能是循环一次出一个值输出出来。不知道你准备用这数据画图还是干嘛
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-17 16:10:46 | 显示全部楼层
本帖最后由 waylee2008 于 2012-8-17 16:12 编辑

3.JPG

麻烦木头版主 帮帮忙看看 非常感谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-17 17:48:18 | 显示全部楼层

你的公式我晓得了,如果s格点为360*180的话,得出的涡度值格点只有359*179,因为你得保证每个涡度周围都有四个点,起点是(2,2),结束是(359.179)。计算的数据输出到二进制文件中,用两层循环,每次循环输出一个值,循环控制经纬度,用格点形式,周围四个点表示比如:s(x+1),s(y+1),s(x-1),s(y-1)这样的形式,然后定义公式就好了,具体的你思考下,自己试试~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-26 22:13:40 | 显示全部楼层
麻烦各位用fortran编写下 谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-24 19:36:08 | 显示全部楼层
应该是要用fortran了,grads里没有那么多计算
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-24 22:57:00 | 显示全部楼层
grads里有函数数hcurl计算涡度
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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