Re: Question about overlaying WRF plot on non-WRF plot [Was: (no subject)]

Hi Merrill,

The wrf_vector plotting script simply creates a vector plot, and doesn't know anything about lat/lon coordinates. This is because usually you call wrf_vector as a precursor to calling wrf_map_overlays, which then does have the correct map information.

In order to get the overlay right, you need set some special vector resources for defining the lat/lon grid for the vector plot. Try this, and let me know if it doesn't work:

  wrflat = wrf_user_getvar(wf,"XLAT",0)
  wrflon = wrf_user_getvar(wf,"XLONG",0)
. . .
    opts@vfXArray = wrflon
    opts@vfYArray = wrflat
    vector = wrf_vector(wf,wks,u_plane(wtm,lv,:,:),v_plane(wtm,lv,:,:),opts)

--Mary

On Nov 14, 2012, at 12:06 PM, Merrill Bean wrote:

> I'm trying to overlay wrf wind barbs over chimere dust contours. I'm hacking together the raster3.ncl
> and the wrf_PressureLevel3.ncl exaples. I get warnings about resources that aren't valid for the wrf
> plot, but it does overlay. Just not correctly. Plotting the wrf data separately is different form the overlay.
> Also the titles disappeared when I overlay the plots.
> Questions:
> Should I do the overlays in ncl or the wrf options?
> What resources should I check to verify that I'm plotting on the same map (extents)? and getting labels?
> Any help is appreciated, I'll include the code.
> Thanks,
> Merrill
> (ncl noob)
> NCAR Command Language Version 6.1.0
> Ubuntu 10.4 2.6.32-45-server #99-Ubuntu SMP Tue Oct 16 16:41:38 UTC 2012 x86_64 GNU/Linux
> <combchimwrf.ncl>_______________________________________________
> 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 Wed Nov 14 14:54:02 2012