爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 24814|回复: 34

[脚本编辑] Grads计算温度平流

  [复制链接]

新浪微博达人勋

发表于 2015-6-2 11:38:16 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 shayingying 于 2015-6-2 11:39 编辑


      看了这个计算温度平流的帖子:
      http://bbs.06climate.com/forum.php?mod=viewthread&tid=2463

      以及官网上用cdiff计算温度平流的例子(例2):
     http://www.iges.org/grads/gadoc/gadocindex.html

      发现帖子里和官网介绍用cdiff计算温度平流的例子有一个地方不一样 贴主: pi=3.1416/180 这样下来cos函数里面就有问题 此外,dx也多除了一个180 只有将pi=3.1416的话才和官网是一样的,或者去掉式子中的180,即:


'define pi=3.1415926/180.0'

'define r=6.356e6'
'define dtx = cdiff(tpm,x)'
'define dty = cdiff(tpm,y)'
'define dx=cdiff(lon,x)*cos(lat*pi)*pi*r'
'define dy=cdiff(lat,y)*pi*r'
'define advT=-(upm*dtx/dx+vpm*dty/dy)'

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

新浪微博达人勋

发表于 2015-6-2 13:39:36 | 显示全部楼层
最近刚好在用,谢谢分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-7-30 15:34:13 | 显示全部楼层
请问计算出的温度平流单位是什么
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-9-22 17:44:12 | 显示全部楼层
你说的没错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-9-22 18:08:55 | 显示全部楼层

第六行为什么要设置
'set lev 1000 200'
屏蔽后图就变了
第十五行不是有
'set lev 850'
结果不一样,为什么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-30 10:45:57 | 显示全部楼层
楼主说的对,拿到程序应该先进行判断,确保正确后再使用。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-30 10:52:39 | 显示全部楼层
annalin0901 发表于 2015-9-22 18:08
第六行为什么要设置
'set lev 1000 200'
屏蔽后图就变了

屏蔽之后默认lev =1000,和你设置的850当然不一样
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-21 11:16:43 | 显示全部楼层
果然用这个程序得到了我想要的~谢谢分享~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-22 15:22:25 | 显示全部楼层
洋咩咩 发表于 2015-12-21 11:16
果然用这个程序得到了我想要的~谢谢分享~

好久没来气象家园了 不谢哈
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-22 15:22:40 | 显示全部楼层
右手边的小星星 发表于 2015-10-30 10:45
楼主说的对,拿到程序应该先进行判断,确保正确后再使用。

是啊 先看懂
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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