problem with WRF output

Dear NCL Users, I have an output of WRF model and it 's needed to plot a parameter in name "EXTCOF55". when I got the plot, it's clear that there is a problem but I don't know where it is!!! my script is incorrect or my output is incorrect? you can see my script in follow and a sample image output in attach. can anyone help me please??!!! Dose anyone have any suggestion?! thanks in advance Marzie -------------------------------------------------- This is summary of "EXTCOF55" : ;float EXTCOF55(Time, bottom_top, south_north, west_east) ; ;EXTCOF55:FieldType = 104 ; ;EXTCOF55:MemoryOrder = "XYZ" ; ;EXTCOF55:description = "Extinction coefficients for .55um" ; ;EXTCOF55:units = "km^-1" ; ;EXTCOF55:stagger = "" ; ;EXTCOF55:coordinates = "XLONG XLAT" ; ----------------------------------------------------------- and this is my complete script: 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/WRF_contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" begin a = addfile("wrfout_d01_2009-07-05_10_00_00.nc","r") wks = gsn_open_wks("pdf","do3.11.00") gsn_define_colormap(wks,"sunshine_9lev") res= True res@lgTitleOn= True res@MainTitle= "REAL-TIME WRF" res@gsnMaximize= True res@gsnSpreadColors= True res@cnFillOn= True res@cnLinesOn= False res@cnLineLabelsOn= False pltres= True mpres= True mpres@gsnSpreadColors= True mpres@mpGeophysicalLineColor= "Black" mpres@mpNationalLineColor= "Black" mpres@mpUSStateLineColor= "Black" mpres@mpGridLineColor= "Black" mpres@mpLimbLineColor= "Black" mpres@mpPerimLineColor= "Black" mpres@mpGeophysicalLineThicknessF = 2 times= wrf_user_list_times(a); get times in the file   ntimes = dimsizes(times) do it = 0,ntimes-1,1; TIME LOOP print("Working on time: " + times(it) ) res = True res@TimeLabel = times(it); Set Valid time to use on plots ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; EXTCOF55= wrf_user_getvar(a,"EXTCOF55",it) res@lbTitleOn= True ;plotting for AOD opts=res pm= wrf_contour(a,wks,EXTCOF55(1,:,:),opts) delete(opts) plot = wrf_map_overlays(a,wks,(/pm/),pltres,mpres) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; end do; END OF TIME LOOP end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk