Re: Using different projections with wrf_contour/wrf_map_overlays

Mary,

I will take a look at your example.

Thanks,
Craig


On Mon, Jul 21, 2014 at 3:46 PM, Mary Haley <haley@ucar.edu> wrote:

> Hi Craig,
>
> The wrf_map_overlays function is not meant to be used for any projection
> other than the one provided on the WRF output file. That's why the first
> argument to this function is the file id, because it reads the map
> projection off the file, and sets this all internally.
>
> If you want to plot WRF output data on a different map projection you will
> need to:
>
> 1. Read the XLAT and XLONG arrays off the wrf output file and attach them
> to your data variable as special "lat2d" and "lon2d" attributes.
> 2. Set the appropriate map resources for zooming in on the map yourself.
> 3. Use the gsn_csm_contour_map function to create the plot.
>
> I've created a sample page showing how to do this using the 'HGT' variable
> on a WRF output file. I like this variable because it's easy to tell if
> the plot looks correct:
>
> http://www.ncl.ucar.edu/Applications/wrfgsn.shtml
>
> --Mary
>
>
>
> On Fri, Jul 18, 2014 at 3:38 PM, Craig Tierney - NOAA Affiliate <
> craig.tierney@noaa.gov> wrote:
>
>> I am plotting some wrf results from a 3km run on a Lambert grid. I would
>> like to get the grid to match some other results I have from a different
>> model, so I would like to plot on a Cylindrical Equidistant projection. I
>> tried setting the projection with:
>>
>> mpres@mpProjection = "CylindricalEquidistan"
>>
>> Then calling the appropriate NCL functions:
>>
>> contour = wrf_contour(cfn,wks,data,res)
>> ocontour=wrf_contour(cfn,wks,sp1,ores)
>> plot=wrf_map_overlays(cfn,wks,(/contour,ocontour/),pltres,mpres)
>>
>> This didn't seem to work. Is this possible?
>>
>> Thanks,
>> Craig
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>

Received on Mon Jul 21 09:48:32 2014