- 积分
- 914
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我已经把网上说的调整经纬度间隔的方法通通试过一遍了,可是结局依然无用,跪求原因
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
;
a = addfile("1.nc","r")
wks = gsn_open_wks(“pdf”,"qctest")
; Set some basic resources
res = True
res@Footer = False
res@InitTime = False
res@vpWidthF = 0.8
res@vpHeightF =0.5
mpres = True ; Map resources
pltres = True ; Plot resources
times = wrf_user_getvar(a,"times",-1) ; get all times in the file
ntimes = dimsizes(times) ; number of times in the file
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
do it = 0,ntimes-1,24 ; TIME LOOP
print("Working on time: " + times(it) )
res@TimeLabel = times(it) ; Set Valid time to use on plots
qc = wrf_user_getvar(a,"QCLOUD",it)
qc = qc*1000.
qc@units = "g/kg"
level = 0
opts = res
opts@cnFillOn = True
opts@gsnSpreadColors = False
opts@cnMissingValPerimOn = True
opts@cnMissingValFillColor = 0
opts@cnMissingValFillPattern = 11
opts@tmXTOn = False
opts@tmYROn = False
opts@tmXBMode = "Explicit"
opts@tmXBValues = (/110,115,120,125,130/) ; Create tick marks
opts@tmXBLabels = (/"110","1150","120","125","130"/); Create labels
opts@tmXBLabelFontHeightF = 0.015
opts@tmYLMode = "Explicit"
opts@tmYLValues = (/30,40,50/) ; Create tick marks
opts@tmYLLabels = (/"30","40","50"/)
opts@cnLevelSelectionMode = "ExplicitLevels"
opts@cnLevels=(/0.0005,0.2,0.4,0.7/)
res@cnFillColors=(/"White","Blue","Green","Red","Grey"/)
contour = wrf_contour(a,wks,qc(level,:,:),opts)
plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
delete(contour)
end if
delete(opts)
end do ; END OF TIME LOOP
end
|
|