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
Received on Sun Apr 10 12:40:35 2011