Dear NCL users,My name is Eun-A Kam.I'd like to make animation image and change 700x700 resolution.How can I fixing script?I would appreciate if anybody could help me how to fix it.
Here is my code :
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"
begin
; Open input file. in_file = addfile("../auxhist1_d01_2008-07-31_00:00:00.nc","r")
; Read time dimesion. times = wrf_user_list_times(in_file) ntimes=dimsizes(times)
type = "ncgm" wks = gsn_open_wks(type,"animation") gsn_define_colormap(wks,"wh-bl-gr-ye-re")
do it = 5, ntimes-1, 6 ; Time Loop start
if ( it .le. 9 ) then valtime = "0"+it else valtime = it end if
; type = "ps"; wks = gsn_open_wks(type,"d02_surface_20080731_"+valtime+"_precipitation"); gsn_define_colormap(wks,"wh-bl-gr-ye-re")
rainc = wrf_user_getvar(in_file,"RAINC",it) rainnc = wrf_user_getvar(in_file,"RAINNC", it) snownc = wrf_user_getvar(in_file,"SNOWNC",it) graupelnc = wrf_user_getvar(in_file,"GRAUPELNC",it) rainfall = rainc + rainnc + snownc + graupelnc rainfall_at_description = "Total Precipitation" rainfall_at_units = "mm"
res = True res_at_InitTime = False res_at_Footer = False res_at_gsnMaximize = False res_at_TimeLabel = times(it) opts= res opts_at_cnLinesOn = False opts_at_cnFillOn = True opts_at_cnLevelSelectionMode = "ExplicitLevels"
contour = wrf_contour(in_file,wks,rainfall,opts)
pltres = True mpres = True mpres_at_mpGeophysicalLineColor = "black" mpres_at_mpGridLineColor = "black"
plot = wrf_map_overlays(in_file,wks,(/contour/),pltres,mpres)
end do ; Time Loop endend
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 02 2009 - 06:12:22 MDT
This archive was generated by hypermail 2.2.0 : Thu Sep 03 2009 - 07:55:08 MDT