[ncl-talk] (no subject)

Modise Wiston modise.wiston at postgrad.manchester.ac.uk
Tue Feb 17 18:42:35 MST 2015


Hi Mary,

Thanks very much for your response please.. *sorry I took long to respond*
I managed to get the script working though (but in a different coding). However, I will go through the example (guide) you provided to see if I can do it the same way..

Thanks again please,

Wiston
________________________________
From: Mary Haley [haley at ucar.edu]
Sent: 13 February 2015 17:13
To: Modise Wiston
Cc: ncl-talk at ucar.edu
Subject: Re: [ncl-talk] (no subject)

Wiston,

If I understand your code correctly, then I think you need to use the special "ZoomIn" attribute that tells the NCL-WRF plotting scripts how to zoom in on a WRF plot.

It looks like you got the index locations of the lat/lon area of interest, but then you didn't do anything with this information:

       loc  = wrf_user_ll_to_ij(a,(/Min_Lon,Max_Lon/),(/Min_Lat,Max_Lat/),opts)

Once you get "loc", you can use this with the "ZoomIn" feature to tell the wrf_map_overlays procedure how to zoom in, via this code:

       ; loc(0,:) is west-east (x) ; loc(1,:) is south-north (y)
       ; subtract one since we want to use it as an index in NCL
         mpres at ZoomIn = True        ; set up map info for zoomed area
         mpres at Xstart = loc(0,0) - 1
         mpres at Xend   = loc(0,1) - 1
         mpres at Ystart = loc(1,0) - 1
         mpres at Yend   = loc(1,1) - 1


You will also want to comment or remove this code:

              ;;; Setting plotting options
;              mpres at mpMaxLatF  = Max_Lat
;              mpres at mpMinLatF  = Min_Lat
;           mpres at mpMinLonF  = Min_Lon
;              mpres at mpMaxLonF  = Max_Lon

To see a full example of using ZoomIn, go to:

http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/SPECIAL/olr.htm

and look at the "olr.ncl" script.

--Mary


On Wed, Feb 11, 2015 at 11:44 AM, Modise Wiston <modise.wiston at postgrad.manchester.ac.uk<mailto:modise.wiston at postgrad.manchester.ac.uk>> wrote:
Hi,

I want plot data from a larger domain to compare the output [low resolution] with results from a smaller grid [high resolution] for the same location..
The larger domain is the parent to the smaller domain (ratio 20:4). I want to extract data from the large grid just for the grid cells where the smaller domain is located and plot data to compare the results as if I am have results from the same grid at two different resolutions.

I defined the coordinates for the smaller domain (on the parent grid) and I am able to have the map background around the small grid only. However, my script plots data from across the whole parent grid *Not just for the grid cells defined* in the coordinates..

Any help on how I can do this please.. I want the plotted data for the small grid ONLY?
(See attached code in the script please)


Thank you,
Wiston

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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150218/7cf505e0/attachment.html 


More information about the ncl-talk mailing list