爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10664|回复: 14

发现wrf里奇怪的变量,是i,k维的,求解释。

[复制链接]

新浪微博达人勋

发表于 2013-11-27 23:45:42 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 green_tea789 于 2013-11-28 00:09 编辑

发现一个很奇观的现象,有一些变量是i,k维的。
只有module_cu_nsas.F, module_bl_gwdo.F, module_fddaobs_driver.F有(ims:ime,kms:kme)的变量。
我要用到module_bl_gwdo.F,输出地形重力波的一些参数,比如froude数。
为什么不是i,k,j维的呢?


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

新浪微博达人勋

发表于 2013-11-28 09:27:36 | 显示全部楼层
再往上看你就会发现J维的循环了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-28 13:24:38 | 显示全部楼层
的确没有。
在subroutine gwdo2d里的变量全都是(i,k)维的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-28 22:51:23 | 显示全部楼层
实在搞不明白。gwdo2d里的变量声明有的是(its:ite,kts:kte),而输出的变量dudt是(ims:ime,kms:kme),dudt居然在its,kts循环里被赋值。
然后通过call gwdo2d(dudt=rublten(ims,kms,j)  将2维的dudt变成3维的rublten输出.
its和ims的区别到底是什么啊,注释说its是start index for j in tile, ims是start index for j in memeory. 还是什么都不明白。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-12-26 16:31:01 | 显示全部楼层
楼主你好,我最近刚刚接触这个,请问如何求得你所说的froude数?WRF模式中可以直接计算出来吗?谢谢!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-1-10 13:48:36 | 显示全部楼层

这个可以在phys目录下查代码啊。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-12 11:22:13 | 显示全部楼层
green_tea789 发表于 2014-1-10 13:48
这个可以在phys目录下查代码啊。

谢谢楼主,您的意思就是查看Fortran程序吗?  还有个问题请教,Froude数应该是几维数据呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-1-14 22:26:48 | 显示全部楼层
牧场物语之主 发表于 2014-1-12 11:22
谢谢楼主,您的意思就是查看Fortran程序吗?  还有个问题请教,Froude数应该是几维数据呢?

是的。阿斯负担发送到发速度
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-1-15 11:16:17 | 显示全部楼层
green_tea789 发表于 2014-1-14 22:26
是的。阿斯负担发送到发速度

我没看懂您的意思,之前按照您说的那个方法一层层call出这个变量到 wrfout 结果变量是输出了,但是全是0....我觉得是调用的时候有问题呀。还望赐教。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-1-15 15:18:48 | 显示全部楼层
我也没成功呢,估计是参数化方案调用的问题。我还没时间搞这个,要弄清楚原因,在程序里加print就知道哪里出问题了。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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