Re: trouble in plotting Pacific region

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Thu Feb 07 2013 - 10:37:28 MST

Hi Bo,
Assuming that your data's longitudes run from 0->360 (you can check that
by doing a printVarSummary(data)), you should specify lonL = 120 and
lonR = 320.. For example:

lonL = 120.
lonR = 320.
printVarSummary(data)
data_subset = data(:,:,{lonL:lonR}) ; data = time x lat x lon
printVarSummary(data_subset)

If your data's longitudes run from -180:180 and you wish to select
40:320, use the lonFlip function prior to subsetting..
Adam

On 02/07/2013 10:32 AM, Dong, Bo wrote:
>
> Hi Adam,
>
> I tried these setting, this time the base map is correctly plotted
> over the pacific region, but there isn't any color showing up on the map.
>
> It seems that the EOF still consider taking data from the region of
> 40W-120E, with the setting of lonL = -40. lonR = 120.
>
> just wondering how to let EOF take data from 120E-40W rather than
> 40W-120E? I also tried setting lonL = 120. lonR = -40. but failed...
>
> thanks much,
>
> Bo
>
> ------------------------------------------------------------------------
> *From:* ncl-talk-bounces@ucar.edu [ncl-talk-bounces@ucar.edu] on
> behalf of Adam Phillips [asphilli@ucar.edu]
> *Sent:* Thursday, February 07, 2013 12:16 PM
> *To:* ncl-talk@ucar.edu
> *Subject:* Re: [ncl-talk] trouble in plotting Pacific region
>
> Hi Bo,
> Try setting:
> res@mpMinLonF = 120.
> res@mpMaxLonF = 320.
> res@mpCenterLonF = (res@mpMinLonF+res@mpMaxLonF)/2.
> Adam
>
> On 02/07/2013 09:58 AM, Dong, Bo wrote:
>>
>> Dear All,
>>
>> I'm trying to plot EOF map for the 120E-40W pacific region. I set
>>
>> lonL = -40.
>>
>> lonR = 120.
>>
>> res@mpMinLonF <mailto:res@mpMinLonF> = lonL
>> res@mpMaxLonF <mailto:res@mpMaxLonF> = lonR
>>
>> However, the output plot shows the region of 40W-120E (i.e., left
>> from Atlantic across the meridian, and all the way west to Asia), not
>> 120E-40W Pacific region.
>>
>> I tried to set res@mpCenterLonF <mailto:res@mpCenterLonF> = 180, but
>> still doesn't work.
>>
>> Does anyone have idea how to plot this region correctly? Attached is
>> the eof script that I use.
>>
>> Thanks much!
>>
>> Bo
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> --
> ______________________________________________________________
> Adam Phillipsasphilli@ucar.edu
> NCAR/Climate and Global Dynamics Division (303) 497-1726
> P.O. Box 3000
> Boulder, CO 80307-3000http://www.cgd.ucar.edu/cas/asphilli

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Feb 7 10:37:37 2013

This archive was generated by hypermail 2.1.8 : Wed Feb 13 2013 - 09:25:58 MST