- 积分
- 1422
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
上图是我提取的降水序列用pdfx函数画 的,脚本如下:
ppt1_3 = asciiread("day_extreme_precipitation_amount_summer_wrf_rcp45_2006_2014.txt",-1,"float")
opt1=True
opt1@bin_min = 20
opt1@bin_max = 48
ppt1_4 = pdfx(ppt1_3, 0, opt1)
ppt2_3 = asciiread("day_extreme_precipitation_amount_summer_wrf_rcp45_2015_2030.txt",-1,"float")
opt2=True
opt2@bin_min = 20
opt2@bin_max = 48
ppt2_4 = pdfx(ppt2_3, 0, opt2)
nVar = 2
nBin = ppt1_4@nbins
xx = new ( (/nVar, nBin/), typeof(ppt1_4))
xx(0,:) = ppt1_4@bin_center
xx(1,:) = ppt2_4@bin_center
yy = new ( (/nVar, nBin/), typeof(ppt1_4))
yy(0,:) = (/ ppt1_4 /)
yy(1,:) = (/ ppt2_4 /)
pltDir = "./"
pltName = "pdf_extreme_pre_wrf_rcp45_2006_2014and2015_2030"
pltType = "png"
wks = gsn_open_wks(pltType, pltDir+pltName)
res = True
res@xyLineThicknesses = (/2.0,2.0/)
res@xyLineColors = (/"blue","red"/)
res@xyDashPatterns = "dash"
res@tiYAxisString = "PDF (%)"
plot = gsn_csm_xy (wks, xx, yy, res)
但是,我其实是想要这样的图:
还请指教
|
|