求助纠正错误！u，v风场叠加绘距平图！

月是故乡明

limoxi 发表于 2015-4-28 21:04
如图，明显太大了，距平都达到50！

可以尝试单独画出相减的两项看看问题出在哪

limoxi

如图，画出来的200hpa的气候平均态矢量风场不对，偏大很多。

1. 'define u20ave=ave(uwnd.1,t=1,t=360,12)'

复制代码

应该是ave命令不对，但用在高度场上没问题，难道是因为u，v分量有南北、东西方向，不能这样求平均？
单独画aveu20的等值线也是偏大，如下图

月是故乡明

limoxi 发表于 2015-4-29 12:18
如图，画出来的200hpa的气候平均态矢量风场不对，偏大很多。

应该是ave命令不对，但用在高度场上 ...

ave是可以用的，因为正负就已经代表了方向，你把ave里t+0,t=360统一成一种表示方法试一下，你单独画一个时次看看风场大小是否正确，确认一下你的数据无误

limoxi

如果不用ave就没有问题，画出来的风场正常 说明数据没问题
就是不能用ave，一用就大的夸张。

1. 'reinit''sdfopen f:/data/u20.nc''sdfopen f:/data/v20.nc''define u20ave=ave(uwnd.1,t=1,t=360,12)''define v20ave=ave(vwnd.2,t=1,t=360,12)''set lon 60 180''set lat 0 70''set t 414''set gxout vector''d u20ave;v20ave''gxprint d:/pic/w202.png';

'reinit''sdfopen f:/data/u20.nc''sdfopen f:/data/v20.nc''define u20ave=ave(uwnd.1,t=1,t=360,12)''define v20ave=ave(vwnd.2,t=1,t=360,12)''set lon 60 180''set lat 0 70''set t 414''set gxout vector''d u20ave;v20ave''gxprint d:/pic/w202.png';

就是不知道哪里不对劲

limoxi

'reinit'
'sdfopen f:/data/u20.nc'
'sdfopen f:/data/v20.nc'

'define u20ave=ave(uwnd.1,t=1,t=360,12)'
'define v20ave=ave(vwnd.2,t=1,t=360,12)'

'set lon 60 180'
'set lat 0 70'
'set t 414'
'set gxout vector'
'd u20ave;v20ave'
'gxprint d:/pic/w202.png'
;

limoxi

'reinit'
'sdfopen f:/data/u20.nc'
'sdfopen f:/data/v20.nc'
'set lon 60 180'
'set lat 0 70'
'set lev 200'
'set t 1 12'
'define u20ave=ave(uwnd.1,t+0,t=360,12)'
'set lon 60 180'
'set lat 0 70'
'set lev 200'
'set t 1 12'
'define v20ave=ave(vwnd.2,t+0,t=360,12)'
'set t 1 420'
'define u20=(uwnd.1-u20ave)'
'define v20=(vwnd.2-v20ave)'

'set lon 60 180'
'set lat 0 70'
'set t 414'
'set gxout vector'
*'set gxout contour'
'd u20;v20'
*'d v20'
*'gxprint d:/pic/awind201.png'
;



之前应该因为没有循环，maybe。不过现在好了！多谢！

1434272731

/猪头哼哼 发表于 2013-5-20 00:24
http://www.esrl.noaa.gov/psd/cgi-bin/data/composites/printpage.pl
楼主可以去这个网站在线合成试试~~~ ...

楼主，发现这个在线合成功能挺有用的，但是有缺陷的是好像没有合成分析的显著性检验

qq469015280

kongfeng0824 发表于 2013-5-20 05:50
对。可以吧网上的例子的数据和你画出来的对比一下既可以知道你的程序是否有问题。

请问风场距平场，是对u，v分别求距平得us，vs，再vector (us,vs)吗？？