[ncl-talk] 3 dimensions to 2 dimensions
Md. Jalal Uddin
dmjalal90 at gmail.com
Tue Apr 30 20:41:57 MDT 2019
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190501/244e58ac/attachment.html>
More information about the ncl-talk
mailing list