夏季500hPa平均位势高度场

冰凌

玄德来了 发表于 2016-4-14 08:51
用nc资料，http://www.esrl.noaa.gov/psd/data/gridded/
你进去这里看看，我用的是这里的资料

谢谢   我用逐月的数据做的平均

阿西达卡

多谢楼主分享

dgq957

冰凌 发表于 2016-4-14 10:41
谢谢   我用逐月的数据做的平均

能不能把你的gs 文件和CTL文件分享一下呀，我想做个冬季的平均，正愁得不会绘图呢

冰凌

dgq957 发表于 2016-4-27 00:20
能不能把你的gs 文件和CTL文件分享一下呀，我想做个冬季的平均，正愁得不会绘图呢

'sdfopen h:/ncep/hgt.mon.mean.nc'
'set lat 20 55'
'set lon 70 140'
'set xlint 10'
'set ylint 10'

'set mpdset cnworld'
'set grads off'
'set grid off'

*'set t 738'
'set lev 500'
**************1981-2010 summer  ave***************************
'define aa1=hgt(t=402)+hgt(t=414)+hgt(t=426)+hgt(t=438)+hgt(t=450)+hgt(t=462)+hgt(t=474)+hgt(t=486)+hgt(t=498)+hgt(t=510)+hgt(t=522)+hgt(t=534)+hgt(t=546)+hgt(t=558)+hgt(t=570)+hgt(t=582)+hgt(t=594)+hgt(t=606)+hgt(t=618)+hgt(t=630)+hgt(t=642)+hgt(t=654)+hgt(t=666)+hgt(t=678)+hgt(t=690)+hgt(t=702)+hgt(t=714)+hgt(t=726)+hgt(t=738)+hgt(t=750)'
'define aa2=hgt(t=403)+hgt(t=415)+hgt(t=427)+hgt(t=439)+hgt(t=451)+hgt(t=463)+hgt(t=475)+hgt(t=487)+hgt(t=499)+hgt(t=511)+hgt(t=523)+hgt(t=535)+hgt(t=547)+hgt(t=559)+hgt(t=571)+hgt(t=583)+hgt(t=595)+hgt(t=607)+hgt(t=619)+hgt(t=631)+hgt(t=643)+hgt(t=655)+hgt(t=667)+hgt(t=679)+hgt(t=691)+hgt(t=703)+hgt(t=715)+hgt(t=727)+hgt(t=739)+hgt(t=751)'
'define aa3=hgt(t=404)+hgt(t=416)+hgt(t=428)+hgt(t=440)+hgt(t=452)+hgt(t=464)+hgt(t=476)+hgt(t=488)+hgt(t=500)+hgt(t=512)+hgt(t=524)+hgt(t=536)+hgt(t=548)+hgt(t=560)+hgt(t=572)+hgt(t=584)+hgt(t=596)+hgt(t=608)+hgt(t=620)+hgt(t=632)+hgt(t=644)+hgt(t=656)+hgt(t=668)+hgt(t=680)+hgt(t=692)+hgt(t=704)+hgt(t=716)+hgt(t=728)+hgt(t=740)+hgt(t=752)'
'define ave1=(aa1+aa2+aa3)/90'
*'define ave1=aa1/30'
******end********1981-2010 summer  ave***************************

'define bb1=ave(hgt,t=738,t=740)'
'define aa4=bb1-ave1'

*'define aa4=hgt-ave1'

'set gxout contour'
'set ccolor 1'
*'set clevs 5600 5610 5620 5630 5640 5650 5660 5670 5680 5690 5700 5710 5720 5730 5740 5750 5760 5770 5780 5790 5800 5810 5820 5830 5840 5850 5860 5870 5880'
'd aa4'

'gxprint h:/ncep/2009_distance_hgt.png white'

冰凌

dgq957 发表于 2016-4-27 00:20
能不能把你的gs 文件和CTL文件分享一下呀，我想做个冬季的平均，正愁得不会绘图呢

月平均数据是nc格式的，不需要ctl文件

dgq957

冰凌 发表于 2016-4-28 10:40
'sdfopen h:/ncep/hgt.mon.mean.nc'
'set lat 20 55'
'set lon 70 140'

真是太谢谢你了，对我这个菜鸟受益匪浅，不过我运行了一下您的gs文件，为啥不保存图片呀，出来的怎么是黑色背景的图呀，不是设置的是白色背景吗？要是能用个循环就更好了。能不能把你的绘制风的gs也分享一下呀？

冰凌

dgq957 发表于 2016-4-28 22:47
真是太谢谢你了，对我这个菜鸟受益匪浅，不过我运行了一下您的gs文件，为啥不保存图片呀，出来的怎么是黑 ...

风也是一样一样的，哈哈哈，其实我也是菜鸟，那时候还没搞懂咋循环咧，所以办法比较笨

亚文

学习了！！！

Hi_Kid

正好用到了，用ncl可以直接读取nc资料，necp的再分析资料网站有很多，都可以免费下载

qinyuefan

https://www.earthsystemgrid.org/