- 积分
- 16131
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-14
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-11-3 08:38:31
|
显示全部楼层
; Example script - plot terrain (using only basic NCL functions)
; Note no map background info is available
; November 2008
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
a = addfile("./geo_em.d01.nc","r") ; Open a file
; type = "x11"
; type = "pdf"
type = "ps"
; type = "ncgm"
wks = gsn_open_wks(type,"plt_geo_1") ; Create a plot workstation
opts = True
opts@MainTitle = "GEOGRID FIELDS"
opts@InitTime = False
opts@Footer = False
; ter = a->HGT_M(0,:,:) ; Read the variable to memory
ter = wrf_user_getvar(a,"HGT_M",0)
res = opts
res@cnFillOn = False
res@gsnSpreadColors = True
; res@cnLevelSelectionMode = "ManualLevels"
; res@cnMinLevelValF = 0.0
; res@cnMaxLevelValF = 6000.
; res@cnLevelSpacingF = 100.
;;
res@ContourParameters = (/0.,500.,100./)
contour = wrf_contour(a,wks,ter,res) ; Create plot
pltres = True
mpres = True
mpres@mpGeophysicalLineColor = "Black"
mpres@mpGeophysicalLineThicknessF = 1.0
mpres@mpGridLineColor = "Red"
mpres@mpLimbLineColor = "Orange"
mpres@mpNationalLineColor = "Blue"
mpres@mpNationalLineThicknessF = 1.0
mpres@mpPerimLineColor = "Purle"
mpres@mpDataSetName = "Earth..4"
mpres@mpDataBaseVersion = "Ncarg4_1"
mpres@mpOutlineSpecifiers = (/"China:states"/)
plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background
end
请指教 |
|