Re: problem with WRF output

Dear Marzie,

It's hard to tell if the problem is with the data or the plot, since it's being overlaid on a map.

First, you should look at the data itself. What does:

print("min/max data = " + min(EXTCOF55(1,:,:)) + "/" + max(EXTCOF55(1,:,:)))

report?

What about just:

print("min/max data = " + min(EXTCOF55) + "/" + max(EXTCOF55))

Usually, when I'm trying to plot something that doesn't look correct, I like to look at just the contour plot so I can make sure it isn't the map overlay that is the problem.

To look at the contour plot: after you call wrf_contour, add the two lines:

   draw(pm)
   frame(wks)

Finally, maybe you want to look at "raster" contours rather than smooth contours. Try adding:

opts@cnFillMode = "RasterFill"

--Mary

On Oct 21, 2012, at 8:23 AM, marziyeh dadizadeh wrote:

> 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
> <aod.pdf>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 23 09:56:21 2012