- 积分
- 79
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-12-22
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2017-2-11 17:37:44
|
显示全部楼层
谢谢您的回复,我做的是一个差值图,下面是程序,烦请您看一下,谢谢
a = addfile("/home/zhy/program/test-wrf/WRFV3/test/em_real/output/wrfout_d01_2014-10-07-10.nc","r")
b = addfile("/home/zhy/program/test-wrf/WRFV3/test/em_real/output/wrfout_d01_2014-10-07-10.OW.nc","r")
type = "png"
wks = gsn_open_wks(type,"Meteorological")
res = True
res@MainTitle = "REAL-TIME WRF"
res@Footer = False
pltres = True
mpres = True
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
times = wrf_user_getvar(a,"times",-1) ; get all times in the file
ntimes = dimsizes(times) ; number of times in the file
f1 = wrf_user_getvar(a,"T2",-1)
f2 = wrf_user_getvar(b,"T2",-1)
t = f1-f2
do it = 0,ntimes-1,1 ; TIME LOOP
wks = gsn_open_wks("png","Meteorological"+times(it))
print("Working on time: " + times(it) )
res@TimeLabel = times(it) ; Set Valid time to use on plots
; Plotting options for T
opts = res
opts@cnFillOn = True
opts@gsnSpreadColors = True
opts@cnLevelSelectionMode ="Manuallevels"
opts@cnMinLevelValF=-0.6
opts@cnMaxLevelValF=0.2
opts@cnLevelSpacingF = 0.2
opts@cnLineLabelsOn = False ; do not use line labels
opts@cnLinesOn = False ; do not draw contour lines
contour_t = wrf_contour(a,wks,t(it,:,:),opts)
opts@cnInfoLabelOn = False
delete(opts)
plot1 = wrf_map_overlays(a,wks,(/contour_t/),pltres,mpres)
end do ; END OF TIME LOOP
end
|
|