求助~wrfout_d02和wrfout_d03画图出错

油麦面疙瘩

我用ARWpost处理wrfout后，我的模拟用的是三层嵌套，运行出d01,d02,d03三个dat文件和ctl文件。
我要画的是一个区域内的累积降水量分布。下边是我的gs文件
'reinit'
'open f:\grib1\wrfout\wrfout-KF_d01_.ctl'
'enable print f:\grib1\wrfout\KF-rainc+rainnc01.gmf'
'set grid off'
'set grads off'
'set lon 92 109'
'set lat 32 44'
'set t 42'
'set z 1'
'set gxout shaded'
'set csmooth on'
'set cmin 5'
'cnbasemap RAINC+RAINNC'
'run cbarn.gs'
'set gxout contour'
'set csmooth on'
'set cmin 5'
'd RAINC+RAINNC'
'print'
'disable print'
;

d01画出来的图是这样的：

但是换成d02，图就变成这个样子了：

换成d03根本就没图了。。。


以下是我的namelist.ARWpost文件
&datetime
start_date = '2012-06-03_18:00:00',
end_date   = '2012-06-07_00:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/

&io
io_form_input  = 2,
input_root_name = '/public1/home/CAS/lzhpc36/WRFV3/run/wrfout_d03*'
output_root_name = '/public1/home/CAS/lzhpc36/dataout/wrfout-KF_d03*'
plot = 'all_list'
fields = 'height,pressure,tk,geopt,theta,tc,td,td2,rh,rh2,clfr,umet,vmet,wdir,wspd,wd10,ws10'
output_type = 'grads'
mercator_defs = .true.
/
split_output = .False.
frames_per_outfile = 1

output_type = 'grads'
output_type = 'v5d'

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'


&interp
interp_method = 1,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,


不知道哦有没有人和我遇到一样的问题？求大神指点~~~

油麦面疙瘩

{:eb302:}{:eb302:}求助啊 ~~~~~

guarenyouji

油麦面疙瘩 发表于 2014-2-20 14:42

                               
登录/注册后可看大图

求助啊 ~~~~~

我也小白，你是学校的学生吗？你的问题可能是用GS画第二重嵌套，以及第三重嵌套时，区域选大了，画第二重嵌套时，模拟区域肯定没有第一重嵌套范围大，你将脚本里'set lon 92 109''set lat 32 44'经纬度范围改小一些，同样画第三重嵌套时再小一些。能力有限，仅限参考，能否多交流，谢谢。

油麦面疙瘩

guarenyouji 发表于 2014-2-23 17:11

                               
登录/注册后可看大图

我也小白，你是学校的学生吗？你的问题可能是用GS画第二重嵌套，以及第三重嵌套时，区域选大了， ...

谢谢你帮我提出的意见，其实我模拟的时候第三重选的区域就是那么大的，只不过画d01的时候，因为我需要的就是那么大的区域，所以只截取了那一部分而已。
后来我又自己研究了研究，发现应该是我转模式的时候，namelist.wps中嵌套区域的位置选的有问题，现在改了，正在重新转，还没出结果，所以还不知道是不是就是这个原因。

guarenyouji

油麦面疙瘩 发表于 2014-2-24 10:57

                               
登录/注册后可看大图

谢谢你帮我提出的意见，其实我模拟的时候第三重选的区域就是那么大的，只不过画d01的时候，因为我需要的就 ...

请问下，你的脚本画图正常吗，我用你的脚本出图，gmf打开错误，不知道为啥

油麦面疙瘩

guarenyouji 发表于 2014-2-24 15:12

                               
登录/注册后可看大图

请问下，你的脚本画图正常吗，我用你的脚本出图，gmf打开错误，不知道为啥

我的脚本我画图的时候没问题啊，第一张图就是用这个脚本画的。

自导自演

不知道楼主的问题解决没有