- 积分
- 1964
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-10-8
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2015-4-22 21:53:58
|
显示全部楼层
我的程序是这样的
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
begin
a = addfile("/public3/home/interior/careeri/huser037/WRF/WRFV3S/run/wrfout_d01_2010-06-29_00:00:00","r")
wks = gsn_open_wks("pdf","test")
; lake=wrf_user_getvar(a,"LAKEMASK",1)
lat = a->XLAT(1,:,:) ;
lon = a->XLONG(1,:,:)
; lon = wrf_user_getvar(a,"lon",1)
;time=ispan(16,264,1)
;print(time)
actt = wrf_user_getvar(a,"ctt",(/48,49/)) ; U component of the wind at 10m
printVarSummary(actt)
avg_ctt=dim_avg_n_Wrap(actt,0)
; Plotting options for Wind Vectors
res = True
res@cnFillOn = True
res@pmLabelBarOrthogonalPosF = -0.1
; res@ContourParameters = (/ 0.1 /)
; res@cnFillColors = (/"White","Chartreuse"/)
contour = wrf_contour(a,wks,avg_ctt,res)
delete(res)
pltres = True ; Set plot options
mpres = True ; Set map options
mpres@mpDataBaseVersion="MediumRes"
mpres@mpDataSetName="Earth..4"
mpres@mpOutlineSpecifiers=(/"China:states"/)
; MAKE PLOTS
plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
frame(wks)
end |
|