求助！！wrf输出变量和Arwpost输出变量对比

schLltt

如图 为namelist.Arwpost截图

请问 wspd,wir与U,V有什么区别？画风场的时候应该用哪组数据？
       u10m,v10m与U10，V10有什么区别？
多谢！！

lhaikun@163.com

U;V分别为风的两个方向的分量，为矢量，画矢量风的时候用；
wspd=wind speed;wdir= wind direction;是标量，画风速大小分布之类的可以用。
wspd=sqrt(U**2+V**2);wdir也是根据UV算的；
你说的u10m；v10m和U10;V10应该是一回事吧，你给的图片中也没有那个东西啊

schLltt

lhaikun@163.com 发表于 2015-8-23 12:37
U;V分别为风的两个方向的分量，为矢量，画矢量风的时候用；
wspd=wind speed;wdir= wind direction;是标量 ...

首先多谢你的解答~~~~~~~
U10;V10是wrf的输出变量，u10m;v10m是经过Arwpost处理后的变量，这两组有没有区别啊~~

lhaikun@163.com

schLltt 发表于 2015-8-23 15:37
首先多谢你的解答~~~~~~~
U10;V10是wrf的输出变量，u10m;v10m是经过Arwpost处理后的变量，这两组有没有 ...

没有..

schLltt

lhaikun@163.com 发表于 2015-8-23 16:40
没有..

大神不要走啊

schLltt

schLltt 发表于 2015-8-23 16:56
大神不要走啊

'reinit'
'open I:/d02.ctl'
'set lon 108 120'
'set lat 25 30'
'set grads off'
'set grid off'
i=1
while(i<2)
'set t 'i''
'set map 1 1 10'
'set gxout shaded'
'd T2'
'cbarn'
'set gxout contour'
'd T2'
'printim I:/2014.11.'math_int((i-1)/24)+6'.'math_mod(i-1,24)*1'.png x1000 y800 white'
'c'
i=i+1
endwhile
;

这个是画2m温度场的gs脚本，可是为什么画出来的的图是这个样子的？换成海平面气压(d slp)，效果也是一样的，请问哪里出问题了，该怎么解决？

schLltt

lhaikun@163.com 发表于 2015-8-23 12:37
U;V分别为风的两个方向的分量，为矢量，画矢量风的时候用；
wspd=wind speed;wdir= wind direction;是标量 ...

'reinit'
'open I:/d02.ctl'
'set lon 108 120'
'set lat 25 30'
'set grads off'
'set grid off'
i=1
while(i<2)
'set t 'i''
'set map 1 1 10'
'set gxout shaded'
'd T2'
'cbarn'
'set gxout contour'
'd T2'
'printim I:/2014.11.'math_int((i-1)/24)+6'.'math_mod(i-1,24)*1'.png x1000 y800 white'
'c'
i=i+1
endwhile
;

这个是画2m温度场的gs脚本，可是为什么画出来的的图是这个样子的？换成海平面气压(d slp)，效果也是一样的，请问哪里出问题了，该怎么解决？

lhaikun@163.com

你arwpost输出结果有问题，或者你的arw输出的ctl有问题，你再检查下

schLltt

lhaikun@163.com 发表于 2015-8-23 21:00
你arwpost输出结果有问题，或者你的arw输出的ctl有问题，你再检查下

图片为Arwpost运行提示
！！！！！！！！以下是ctl文件！！！！！！！！！！！！
dset ^d02.dat
undef 1.e30
title  OUTPUT FROM WRF V3.3.1 MODEL
xdef  309 linear  107.31900   0.04054054
ydef  260 linear   23.65130   0.04054054
zdef   19 levels  
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef  720 linear 00Z01NOV2014      60MN      
VARS    6
U             19  0  x-wind component (m s-1)
V             19  0  y-wind component (m s-1)
T2             1  0  TEMP at 2 M (K)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
slp            1  0  Sea Levelp Pressure (hPa)
ENDVARS
！！！！！！！！！！以下为namelist.Arwpost!!!!!!!!!!
&datetime
start_date = '2014-11-01_00:00:00',
end_date   = '2014-12-01_00:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/home/lttt/Program/WRFV3/run/wrfout_d02_2014-11-01_00:00:00'
output_root_name = '/home/lttt/Program/ARWpost/d02'
plot = 'list'
fields = 'slp,U,V,T2,U10,V10'
mercator_defs = .false.
/
split_output = .true.
output_type = 'grads'
frames_per_outfile = 1000

&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,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
能帮看看哪儿有问题吗~~~~~~~~~~~~

lhaikun@163.com

你的namelist.ARWpost看起来没什么问题，只能提供一些可能的问题排查方法：
1、你plot一下t=n(n>1)时候随便一张，看能不能出图
2、你生成的dat文件多大，很大的话考虑一下用上split_output = .true.
frames_per_outfile = 1000 这两句
3、将output_type = 'grads' 放到运行的语法中去（虽然default为它），尝试一下是否能解决
4、你生成的ctl中我发现没有pdef的设置和options的设置，不知道是不是这样的问题，我在服务器上做的arw生成的ctl是自带如上两个设置的