为什么我画不同时间的高度场距平是一样的（已解决）

瑾珂

'reinit'
'open L:\DJBY\20y_winter\2000_1\fnl_200001.ctl'
'set lev 500'
'define u1=ave(UGRdprs,t=1,t=4,1)'
'define v1=ave(VGRdprs,t=1,t=4,1)'
'define hgt1=HGTprs/9.8'
'define h1=ave(hgt1,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2001_01\fnl_200101.ctl'
'set lev 500'
'define u2=ave(UGRdprs,t=1,t=4,1)'
'define v2=ave(VGRdprs,t=1,t=4,1)'
'define hgt2=HGTprs/9.8'
'define h2=ave(hgt2,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2002_1\fnl_200201.ctl'
'set lev 500'
'define u3=ave(UGRdprs,t=1,t=4,1)'
'define v3=ave(VGRdprs,t=1,t=4,1)'
'define hgt3=HGTprs/9.8'
'define h3=ave(hgt3,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2003_1\fnl_200301.ctl'
'set lev 500'
'define u4=ave(UGRdprs,t=1,t=4,1)'
'define v4=ave(VGRdprs,t=1,t=4,1)'
'define hgt4=HGTprs/9.8'
'define h4=ave(hgt4,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2004_1\fnl_200401.ctl'
'set lev 500'
'define u5=ave(UGRdprs,t=1,t=4,1)'
'define v5=ave(VGRdprs,t=1,t=4,1)'
'define hgt5=HGTprs/9.8'
'define h5=ave(hgt5,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2005_1\fnl_200501.ctl'
'set lev 500'
'define u6=ave(UGRdprs,t=1,t=4,1)'
'define v6=ave(VGRdprs,t=1,t=420,1)'
'define hgt6=HGTprs/9.8'
'define h6=ave(hgt6,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2006_1\fnl_200601.ctl'
'set lev 500'
'define u7=ave(UGRdprs,t=1,t=4,1)'
'define v7=ave(VGRdprs,t=1,t=4,1)'
'define hgt7=HGTprs/9.8'
'define h7=ave(hgt7,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2007_1\fnl_200701.ctl'
'set lev 500'
'define u8=ave(UGRdprs,t=1,t=4,1)'
'define v8=ave(VGRdprs,t=1,t=4,1)'
'define hgt8=HGTprs/9.8'
'define h8=ave(hgt8,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2008_1\fnl_200801.ctl'
'set lev 500'
'define u9=ave(UGRdprs,t=1,t=4,1)'
'define v9=ave(VGRdprs,t=1,t=4,1)'
'define hgt9=HGTprs/9.8'
'define h9=ave(hgt9,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2009_1\fnl_200901.ctl'
'set lev 500'
'define u10=ave(UGRdprs,t=1,t=4,1)'
'define v10=ave(VGRdprs,t=1,t=4,1)'
'define hgt10=HGTprs/9.8'
'define h10=ave(hgt10,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2010_1\fnl_201001.ctl'
'set lev 500'
'define u11=ave(UGRdprs,t=1,t=4,1)'
'define v11=ave(VGRdprs,t=1,t=4,1)'
'define hgt11=HGTprs/9.8'
'define h11=ave(hgt11,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2011_1\fnl_201101.ctl'
'set lev 500'
'define u12=ave(UGRdprs,t=1,t=4,1)'
'define v12=ave(VGRdprs,t=1,t=4,1)'
'define hgt12=HGTprs/9.8'
'define h12=ave(hgt12,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2012_1\fnl_201201.ctl'
'set lev 500'
'define u13=ave(UGRdprs,t=1,t=4,1)'
'define v13=ave(VGRdprs,t=1,t=4,1)'
'define hgt13=HGTprs/9.8'
'define h13=ave(hgt13,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2013_1\fnl_201301.ctl'
'set lev 500'
'define u14=ave(UGRdprs,t=1,t=4,1)'
'define v14=ave(VGRdprs,t=1,t=4,1)'
'define hgt14=HGTprs/9.8'
'define h14=ave(hgt14,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2014_1\fnl_201401.ctl'
'set lev 500'
'define u15=ave(UGRdprs,t=1,t=4,1)'
'define v15=ave(VGRdprs,t=1,t=4,1)'
'define hgt15=HGTprs/9.8'
'define h15=ave(hgt15,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2015_1\fnl_201501.ctl'
'set lev 500'
'define u16=ave(UGRdprs,t=1,t=4,1)'
'define v16=ave(VGRdprs,t=1,t=4,1)'
'define hgt16=HGTprs/9.8'
'define h16=ave(hgt16,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2016_1\fnl_201601.ctl'
'set lev 500'
'define u17=ave(UGRdprs,t=1,t=4,1)'
'define v17=ave(VGRdprs,t=1,t=4,1)'
'define hgt17=HGTprs/9.8'
'define h17=ave(hgt17,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2017_1\fnl_201701.ctl'
'set lev 500'
'define u18=ave(UGRdprs,t=1,t=4,1)'
'define v18=ave(VGRdprs,t=1,t=4,1)'
'define hgt18=HGTprs/9.8'
'define h18=ave(hgt18,t=1,t=4,1)'


'open L:\DJBY\20y_winter\2018_1\fnl_201801.ctl'
'set lev 500'
'define u19=ave(UGRdprs,t=1,t=4,1)'
'define v19=ave(VGRdprs,t=1,t=4,1)'
'define hgt19=HGTprs/9.8'
'define h19=ave(hgt19,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2019_1\fnl_201901.ctl'
'set lev 500'
'define u20=ave(UGRdprs,t=1,t=4,1)'
'define v20=ave(VGRdprs,t=1,t=4,1)'
'define h20=ave(HGTprs/9.8,t=1,t=4,1)'

'open L:\DJBY\20y_winter\2020_1\fnl_202001.ctl'
'set lev 500'
'define u21=ave(UGRdprs,t=1,t=4,1)'
'define v21=ave(VGRdprs,t=1,t=4,1)'
'define hgt21=HGTprs/9.8'
'define h21=ave(hgt21,t=1,t=4,1)'


'define ah=(h1+h2+h3+h4+h5+h6+h7+h8+h9+h10+h11+h12+h13+h14+h15+h16+h17+h18+h19+h20+h21)/21'

'define au=(u1+u2+u3+u4+u5+u6+u7+u8+u9+u10+u11+u12+u13+u14+u15+u16+u17+u18+u19+u20+u21)/21'

'define av=(v1+v2+v3+v4+v5+v6+v7+v8+v9+v10+v11+v12+v13+v14+v15+v16+v17+v18+v19+v20+v21)/21'

'open L:\DJBY\ca_saf\201801\20180101.ctl'
'set lev 500'
'define u22=ave(UGRdprs,t=1,t=4,1)'
'define v22=ave(VGRdprs,t=1,t=4,1)'
'define hgt22=HGTprs/9.8'
'define h22=ave(hgt21,t=1,t=4,1)'


*'define dhh=h22-ah'

*'define du22=u22'
*'define dv22=v22'

'set lat 10 80'
'set lon 20 160'

'set grid off'
'set parea 1 9.5 0.5 8'
'set gxout contour'
'set mpdset cnworld'
'set map 9 1 1'
'draw map'
'set ccolor 1'
'set cint 4'
'set cthick 7'

'set gxout shaded'
'set grads off'
*'L:\DJBY\ca_saf\rb_matlab.gs'
'd ah'
'run cbarn_interp.gs '
*'L:\DJBY\ca_saf\cbar_matlab.gs'

'set gxout barb'
'set gxout vector'
'd skip(u22*2.5,3);skip(v22*2.5,3)'

'printim L:\DJBY\ca_saf\201801\jupin_500\fnl_500hpa_jan1_143.png white'

'disable gxprint'
；

占戈示申

理论上来说，你用的好像都是第一个file的数据。
比如：'define u2=ave(UGRdprs,t=1,t=4,1)'应该是'define u2=ave(UGRdprs.2,t=1,t=4,1)'……这个样子。

来自苹果APP客户端

瑾珂

占戈示申 发表于 2021-5-22 22:32
理论上来说，你用的好像都是第一个file的数据。
比如：'define u2=ave(UGRdprs,t=1,t=4,1)'应该是'define u ...

还是不行呀

瑾珂

问题已解决

屁屁呀

你好，请问可以加下我的扣扣1326818571，想请教下距平图，我实在弄不出来