- 积分
- 2892
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-7-15
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-4-2 10:39:06
|
显示全部楼层
我还找到了一个官网的例子,Mark一下
; Example script to plot a field from a single metgrid file
; November 2008
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
;
a = addfile("./met_em.d01.2000-01-24_12:00:00.nc","r") ; Open a file
; We generate plots, but what kind do we prefer?
type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
wks = gsn_open_wks(type,"plt_metgrid_1")
res = True ; Set up some basic plot resources
res@MainTitle = "METGRID FILES"
res@Footer = False
pltres = True
mpres = True
lu = wrf_user_getvar(a,"LU_INDEX",0) ; Get land use from file
opts = res ; Set some plotting resources
opts@cnFillOn = True
opts@cnFillMode = "RasterFill"
opts@ContourParameters = (/ 1,16,1 /)
contour = wrf_contour(a,wks,lu,opts)
plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
end
|
|