wrfout文件能否批量计算累计降水?

一顾留白

之前有通过WRF处理生成4天的wrfout文件，生成成功后用NCL作图出现问题，显示错误，见http://bbs.06climate.com/forum.php?mod=viewthread&tid=88871，
考虑是不是文件过大的问题，因此每六小时输出一次看能否处理，想问一下这样几个wrfout文件如何批量处理计算降水量？？

一顾留白

已解决，自问自答系列，修改namelist.input里frame相关的参数，使得文件分成几个小文件输出，运行最后一个时次的文件即可

yeraindrop

一顾留白 发表于 2019-2-28 15:37
已解决，自问自答系列，修改namelist.input里frame相关的参数，使得文件分成几个小文件输出，运行最后一个 ...

是修改后直接运行最后一个时次的wrfout文件？

一顾留白

yeraindrop 发表于 2020-2-11 14:33
是修改后直接运行最后一个时次的wrfout文件？

对，最后一个时次可直接计算总降水，不过我之前也有一个误区，cgwin的NCL对读取文件大小有限制，但是在linux系统下的并没有这个限制，所以建议还是在大型机上比较方便

yeraindrop

一顾留白 发表于 2020-2-11 19:40
对，最后一个时次可直接计算总降水，不过我之前也有一个误区，cgwin的NCL对读取文件大小有限制，但是在li ...

谢谢，楼主

1269737003

楼主可以分享下ncl程序吗，我也需要计算降水量，但我的wrfout文件不能读取

guoshengxin

dim_n_sum_warp,你试试这个函数呢

1269737003

一顾留白 发表于 2020-2-11 19:40
对，最后一个时次可直接计算总降水，不过我之前也有一个误区，cgwin的NCL对读取文件大小有限制，但是在li ...

请问为什么我的wrfout是六小时输出一个，但是用最后一个wrfout文件读出，降水量是0呀？

1269737003

楼主可以分享一下你的脚本文件吗？谢谢呀！

一顾留白

'open D:\CPAS_MERG_HLY_WORKING.lnx.ctl'
'define rain24=0.0'
i=1
while(i<=25)
'set t 'i''
'rain24=rain24+crain.1(t='i')'
i=i+1
endwhile

'set lon 106 122'
'set lat 20 28'
'set xlint 2'
'set ylint 2'

'set grads off'
'set grid off'
'set gxout shaded'


'set rgb 16 255 255 255'
'set rgb 17 173 255 47'
'set rgb 18 154 205 50'
'set rgb 19 0 255 0'
'set rgb 20 34 139 34'
'set rgb 21 255 255 0'
'set rgb 22 255 165 0'
'set rgb 23 255 0 0'
'set rgb 24 238 130 238'


'set ccols  0 16 17 18 19 20 21 22 23 24'
'set clevs  0 30 60 90 120 150 180 200'
'd rain24'
'cbar.gs'
'printim D:/2015hh.gif white'
'reinit'
我是做的这样的循环