[ncl-talk] 3 dimensions to 2 dimensions
Md. Jalal Uddin
dmjalal90 at gmail.com
Tue Apr 30 21:26:16 MDT 2019
Is there any function that makes 3 dimensions to 2 dimensions array like
ndtooned
On Wed, May 1, 2019 at 11:23 AM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:
> What is it you really intend here -- a sum over time, or an average?
>
> I don't want to sum or average. I want to take the whole temporal value at
> each lat and lon and make a two dimensions Rain from three dimensions TRMM.
>
> On Wed, May 1, 2019 at 11:17 AM Rick Brownrigg <brownrig at ucar.edu> wrote:
>
>> Well, there are still two problems with that innermost nested loop:
>>
>> i) What is the value of the 1st (zero-th) index value into TRMM on the
>> first iteration of that loop?
>> do nt=0,ntim-1
>> Rain(nl,ml) = TRMM(nt-1,nl,ml) ; Subscript out of range, error in
>> subscript #0
>>
>> ii) again, this loop is overwriting the value of Rain(nl,ml) repeatedly,
>> for "ntim" times, only to end up with Rain(nl,ml) containing the value
>> TRMM(ntim-1,nl,ml) by the time the innermost loop exits. What is it you
>> really intend here -- a sum over time, or an average?
>>
>> I really don't know at all about the second set of errors from
>> copy_VarMeta -- it seems like some dimensionality requirement between
>> copy_from and copy_to is getting violated, but I don't have the full
>> picture. See the docs for discussion on requirements on the "leftmost
>> dimensions":
>>
>> http://ncl.ucar.edu/Document/Functions/Contributed/copy_VarMeta.shtml
>>
>> Perhaps get that inner loop correct and the rest may sort itself out....
>>
>> Rick
>>
>>
>> On Tue, Apr 30, 2019 at 8:42 PM Md. Jalal Uddin <dmjalal90 at gmail.com>
>> wrote:
>>
>>> Hi Rick,
>>> Now, I am getting the following errors (Highlight in red color).
>>>
>>> f = addfile("bob_vscs.nc","r")
>>> time =f->time
>>> lat =f->latitude
>>> lon =f->longitude
>>> ntim = dimsizes(time)
>>> nlat = dimsizes(lat)
>>> nlon = dimsizes(lon)
>>> TRMM = f->pcp(:,{0:30},{75:100}) ;TRMM(time, latitude,
>>> longitude)
>>>
>>> Rain =new((/nlat,nlon/),"float")
>>>
>>> do nl=0,nlat-1
>>> do ml=0,nlon-1
>>> do nt=0,ntim-1
>>> Rain(nl,ml) = TRMM(nt-1,nl,ml) ; Subscript out of range, error in
>>> subscript #0
>>> end do
>>> end do
>>> end do
>>> copy_VarMeta(TRMM(0,:,:),Rain) ; fatal:Coordinate variables must
>>> be the same dimension as their dimension fatal:No coordinate variable
>>> exists for dimension (longitude) in variable (var_to)
>>>
>>> Variable: TRMM
>>> Type: float
>>> Total Size: 18720000 bytes
>>> 4680000 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 390] x [latitude | 120] x [longitude |
>>> 100]
>>> Coordinates:
>>> time: [ 0..119208]
>>> latitude: [0.125..29.875]
>>> longitude: [75.125..99.875]
>>> Number Of Attributes: 7
>>> long_name : precipitation (mm/hr)
>>> _FillValue : -9999.9
>>> missing_value : -9999.9
>>>
>>>
>>>
>>>
>>>
>>> On Tue, Apr 30, 2019 at 8:48 PM Rick Brownrigg <brownrig at ucar.edu>
>>> wrote:
>>>
>>>> HI,
>>>>
>>>> It looks like there may be several problems going here:
>>>>
>>>> i) I speculate you mean something more like this:
>>>> time =f->time
>>>> lat =f->latitude
>>>> lon =f->longitude
>>>> ntim = dimsizes(time)
>>>> nlat = dimsizes(lat)
>>>> nlon = dimsizes(lon)
>>>>
>>>> ii) I'm not sure the intent here:
>>>>
>>>> Rain(nl,ml) = TRMM(nt,nl,ml)
>>>>
>>>> Each time through the inner loop, it will (re)assign the one by one
>>>> "nt" values of TRMM to the same element of Rain(nl, ml), exiting the loop
>>>> with it set to the "nt-1"th value. Its equivalent to foregoing the entire
>>>> inner loop and just saying Rain(nl,ml) = TRMM(nt-1,nl,ml)
>>>>
>>>> Hope that helps,
>>>> Rick
>>>>
>>>>
>>>>
>>>> On Tue, Apr 30, 2019 at 4:00 AM Md. Jalal Uddin <dmjalal90 at gmail.com>
>>>> wrote:
>>>>
>>>>> Hi all,
>>>>>
>>>>> I have tried to take all rainfall data for each latitude and
>>>>> longitude point and made 2 dimensions from 3 dimensions. I am getting the
>>>>> following errors.
>>>>> fatal:_NclBuildArray: each element of a literal array must have the
>>>>> same dimension sizes, at least one item doesn't
>>>>>
>>>>> The main code:
>>>>>
>>>>> f = addfile("bob_vscs.nc","r")
>>>>> ntim =f->time
>>>>> nlat =f->latitude
>>>>> nlon =f->longitude
>>>>> TRMM = f->pcp(:,{0:30},{75:100}) ;TRMM(time 390, latitude
>>>>> 120, longitude 100)
>>>>>
>>>>> Rain =new((/nlat,nlon/),"float")
>>>>>
>>>>> do nl=0,nlat-1
>>>>> do ml=0,nlon-1
>>>>> do nt=0,ntim-1
>>>>> Rain(nl,ml) = TRMM(nt,nl,ml)
>>>>> end do
>>>>> end do
>>>>> end do
>>>>>
>>>>> Could you help me to solve the problem, please?
>>>>>
>>>>> Best Regards,
>>>>> Jalal
>>>>> --
>>>>> *Md. Jalal Uddin*
>>>>> MSc in Applied Meteorology (English Language)
>>>>> Nanjing University of Information, Science and Technology, China
>>>>> B.Sc. in Disaster Management (Hons.)
>>>>> Patuakhali Science and Technology University, Bangladesh.
>>>>> Cell: +8613260859092, +8801792052662
>>>>> Web: www.dmjalal90.weebly.com
>>>>> Facebook: jalal.hossen.39
>>>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>>>> Twitter: dmjalal90
>>>>> Skype: dmjalal90
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>
>>>
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology (English Language)
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> Twitter: dmjalal90
>>> Skype: dmjalal90
>>>
>>
>
> --
> *Md. Jalal Uddin*
> MSc in Applied Meteorology (English Language)
> Nanjing University of Information, Science and Technology, China
> B.Sc. in Disaster Management (Hons.)
> Patuakhali Science and Technology University, Bangladesh.
> Cell: +8613260859092, +8801792052662
> Web: www.dmjalal90.weebly.com
> Facebook: jalal.hossen.39
> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
> Twitter: dmjalal90
> Skype: dmjalal90
>
--
*Md. Jalal Uddin*
MSc in Applied Meteorology (English Language)
Nanjing University of Information, Science and Technology, China
B.Sc. in Disaster Management (Hons.)
Patuakhali Science and Technology University, Bangladesh.
Cell: +8613260859092, +8801792052662
Web: www.dmjalal90.weebly.com
Facebook: jalal.hossen.39
LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
Twitter: dmjalal90
Skype: dmjalal90
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190501/b3ec7872/attachment.html>
More information about the ncl-talk
mailing list