Re: dim_avg_Wrap(data)

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Mon Apr 11 2011 - 16:30:56 MDT

Hi Mark,
Try setting:

trXReverse = True

If that doesn't work let the group know..
Adam

On 04/11/2011 08:54 AM, Mark Chan wrote:
> Thanks very much, Dennis!
>
> I did like this (below), and got the correct result. Form here I know "bbb&lat"
> and "lat" are different. Thanks deeply again! By the way, how to modify the plot
> to be from 90N to 90S (from left to right)?
> plot(1) = gsn_csm_xy (wks, bbb&lat,bbb,zres)
>
> Mark
>
>
> ----- Original Message ----
> From: Dennis Shea<shea@ucar.edu>
> To: Mark Chan<cym263@yahoo.com>
> Cc: ncl-talk@ucar.edu
> Sent: Sat, April 9, 2011 9:40:26 PM
> Subject: Re: dim_avg_Wrap(data)
>
> What does the
>
> printVarSummary(data)
>
> indicate *before* you mangle the data?
>
> Why not use
>
> bbb=dim_avg_Wrap(data)
> printVarSummary(bbb)
>
> then
> iNH = ind(bbb&lat.ge.0)
> plot(1) = gsn_csm_xy (wks, bbb&lat(iNH:), bbb(iNH:), zres)
>
>
>
>
> On 4/8/11 5:03 PM, Mark Chan wrote:
>> NCL users,
>>
>> I have a global data file as "data(lat,lon)" and would make a zonal average as
>> below. I especially set southern hemisphere as missing data for distinguishing
>> between northern and southern hemisphere (see the attached plot).
>>
>>
>> plot=new(2,graphic)
>>
>> data!0 = "lat"
>> data!1 = "lon"
>> data&lat = lat(::-1)
>> data&lon = lon
>> data@_FillValue = msg_val
>>
>> clat=conform(data,lat,0)
>> data=where(clat.lt.0,data,-9999.0)
>> bbb=dim_avg_Wrap(data)
>>
>> plot(0) = gsn_csm_contour_map(wks,data,res)
>> plot(1) = gsn_csm_xy (wks, lat,bbb,zres)
>>
>>
>> (1) Why is the latitude not consistent for both spatial pattern (data) and
>> zonal
>> average (bbb)?
>>
>> (2) How to make the second plot starting from 90N to 90S instead of from 90S
> to
>> 90N?
>> (3) I could make the correct plot by "plot(1) = gsn_csm_xy (wks,
>> lat(::-1),bbb,zres)" but why didn't the "dim_avg_Wrap(data)" pass the
>> consistent
>> latitude?
>>
>> (4) By the way, can one explain a bit more about sfXArray and sfYArray? (extra
>> question, not related here :-) )
>>
>> Thanks a lot,
>> Mark
>>
>>
>>
>> _______________________________________________
>> 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

-- 
__________________________________________________
Adam Phillips 
asphilli@ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
Climate and Global Dynamics Division         fax: (303) 497-1333
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 Tue Apr 12 08:31:04 2011

This archive was generated by hypermail 2.1.8 : Tue Apr 19 2011 - 18:32:03 MDT