- 积分
- 209
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-6
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近用gsn_csm_pres_hgt画几张ERA5数据的气压-时间剖面图,全程没有报错,但是发现使用trYReverse=True压根没用,不加trYReverse=True纵坐标Y轴顶端也是200hPa,但数据中(0,0)对应的200hPa数据却没有画到顶端,还是画在了底端,不论怎么调都换不过去,最后还是用了gsn_contour,可算能自由颠倒Y轴值了,求助各位大神,这是怎么回事?
下边是我的作图部分代码:
a=addfile(filename1,"r")
time=a->time
ddhh_zuobiao=ispan(0, 72, 1) ;画图横坐标设置
z=a->level(0:10)
d1=short2flt(a->d(48:120,0:10,20,36))
d1=d1*(10^5)
d12=transpose(d1)
res1=True
res1@cnLinesOn=False
res1@cnFillOn = True
res1@cnLineLabelsOn =False
res1@gsnPresHgtHeightLabelOn=False
res1@gsnAddCyclic = False
res1@sfXArray=ddhh_zuobiao
res1@sfYArray=z
res1@tiYAxisString="pressure(hPa)"
res1@tiXAxisString="time"
res1@tmYLMode = "Explicit"
res1@tmYLValues = z
res1@tmYLLabels=(/"200","","300"," ","400"," ","500"," ","600"," ","700"/)
res1@cnLevelSelectionMode="ExplicitLevels"
res1@cnLevels =(/-20,-15,-10,-5,0,5,10,15,20/)
res1@cnFillColors=(/5,6,7,8,9,12,13,14,15,16/)
res1@lbLabelBarOn =False
plot=gsn_csm_pres_hgt(wks, d12, res1)
用来画图的数据格式如下:Variable: d12
Dimensions and sizes: [level | 11] x [time | 73]
Coordinates:
level: [200..700]
time: [1080048..1080120]
Number Of Attributes: 9
longitude : 99
latitude : 40
units : s**-1
long_name : Divergence
standard_name : divergence_of_wind
|
|