Re: draw xy plot from 3 dimensional variable

At the suggestion of an NCL user, I got the dataset
and ran the attached code. It generates graphics.

GPCC data are missing over the oceans.

Enough on this topic!

On 8/5/13 9:32 AM, Dennis Shea wrote:
> To minimize future questions to ncl-talk, I wrote the attached.
>
> *YOU* must show some effort.
>
> *Look* at the output from every printVarSummary.
>
> *Study* the documentation for "lonFlip"
>
>
>
> On 8/5/13 8:46 AM, Mary Haley wrote:
>> I'm not sure what kind of XY plot you want to create from a 3D
>> variable, but you can use the "ndtooned" to convert your 3D variable
>> to a 1D variable.
>>
>> To create an XY plot, see our examples page:
>>
>> http://www.ncl.ucar.edu/Applications/xy.shtml
>>
>> If you want a scatter (points) plot, see:
>>
>> http://www.ncl.ucar.edu/Applications/scatter.shtml
>>
>> --Mary
>>
>> On Aug 5, 2013, at 4:43 AM, Setareh Rahimi wrote:
>>
>>>
>>>
>>> ---------- Forwarded message ----------
>>> From: Setareh Rahimi <setareh.rahimi@gmail.com>
>>> Date: Sun, Aug 4, 2013 at 7:50 AM
>>> Subject: draw xy plot from 3 dimensional variable
>>> To: ncl-talk@ucar.edu
>>>
>>>
>>> Dear all,
>>>
>>> I know that my email have been seen by the member of the ncl groups,
>>> but since this is an ergent need I send this email again. Please
>>> kindlu advice me in this regard as soon as possible.
>>>
>>> After calculating SPI, I need to have a xy plot. My problem is that
>>> SPI has 3 dimensions, as you know to plot a xy plot 1 dimensional
>>> variables are needed. So my question is how can I convert 3
>>> dimensinal SPI to 1 dimensinal. I need to plot SPI for whole Africa
>>> cntinent, not specific lat/lon.
>>>
>>> Thank you in advance.
>>>
>>> ______________________________________________________________________________________________________________________________________
>>>
>>>
>>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>> load"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>
>>> begin
>>>
>>> f = addfile("precip.mon.total.v6.nc","r")
>>> prc2 = f->precip({62456:77035},{-35:37},{344:50})
>>> time = f->time({62456:77035} )
>>>
>>>
>>> yyyymm_i = cd_calendar(time, -2)
>>>
>>> ntim = dimsizes(time)
>>>
>>> dimp = dimsizes(prc2)
>>> ntim = dimp(0)
>>> nlat = dimp(1)
>>> mlon = dimp(2)
>>> len = (/3,6,12,24/)
>>> klen = dimsizes(len)
>>> spi = new((/klen,ntim,nlat,mlon/),float)
>>>
>>>
>>> do k=0,klen-1
>>> spi(k,:,:,:) = dim_spi_n(prc2, len(k), False, 0)
>>> end do
>>> end
>>>
>>>
>>>
>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>> _______________________________________________
>>> 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
>>
>
>
> _______________________________________________
> 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