Re: fatal error: reading multi dimension data with one dimension fixed

Michael,

Guangshan has uploaded a file, I got the data and can repeat the problem

I have create a JIRA-1962 for this issue.
We will update you the progress on this issue.

Thanks,

Wei

On May 8, 2014, at 2:03 PM, Wei Huang <huangwei@ucar.edu> wrote:

> Michael, Guangshan,
>
> If you can provide us a file, we'd certainly want to take a look.
>
> Wei
>
> You can ftp larger datasets and files to:
>
>
> ftp ftp.cgd.ucar.edu
> anonymous
> <use your email address for the password>
> cd incoming
> put ...
> put ...
> .
> .
> .
> quit
> Note that for security reasons, you cannot list the contents of this directory, and neither can we. We need to know the exact names of the files in order to retrieve them.
>
>
>
>
> On May 8, 2014, at 1:38 PM, Michael Notaro <mnotaro@wisc.edu> wrote:
>
>> I looked at it with Guangshan and definitely something is wrong
>> in the newer versions of NCL. I included the ncdump from the file.
>> If you try to retrieve U(:,2,:,:) or U(:,:,2,:) or U(:,:,:,2), it will give an error,
>> but U(2,:,:,:) works ok.
>>
>>
>>
>>
>>
>>
>>
>>
>> [yslogin2]409%: ncdump -h isotag_t85_115.0k.cam2.h0.cel.plev.delta.611-650-ANN.nc
>> netcdf isotag_t85_115.0k.cam2.h0.cel.plev.delta.611-650-ANN {
>> dimensions:
>> time = UNLIMITED ; // (40 currently)
>> lat = 128 ;
>> lon = 256 ;
>> lev = 4 ;
>> variables:
>> float CLDHGH(time, lat, lon) ;
>> CLDHGH:time_rep = "Estimated from monthly modern calendar data" ;
>> CLDHGH:lon = 358.59375 ;
>> CLDHGH:lat = 88.9277353522959 ;
>> CLDHGH:units = "fraction" ;
>> CLDHGH:long_name = "Vertically-integrated high cloud" ;
>> CLDHGH:cell_method = "time: mean" ;
>> CLDHGH:_FillValue = 9.96921e+36f ;
>> float U(time, lev, lat, lon) ;
>> U:time_rep = "Estimated from monthly modern calendar data" ;
>> U:missing_value = 9.96921e+36f ;
>> U:lon = 49.21875 ;
>> U:lat = 48.3263610181882 ;
>> U:lev = 1000.f ;
>> U:cell_method = "time: mean" ;
>> U:long_name = "Zonal wind" ;
>> U:units = "m/s" ;
>> U:_FillValue = 9.96921e+36f ;
>>
>> On May 8, 2014, at 1:57 PM, Guangshan Chen wrote:
>>
>>> Dear NCL developers,
>>>
>>> I just notice that there is an error message popped out when I read multi dimension data with one dimension fixed.
>>>
>>> I could not remember it was shown in the old-version NCL.
>>>
>>> Here is what I did:
>>>
>>> Opening a file containing zonal wind and read zonal wind
>>> uwind = file_in->U(:,2,:,:)
>>>
>>> The following error will come out:
>>> Copyright (C) 1995-2014 - All Rights Reserved
>>> University Corporation for Atmospheric Research
>>> NCAR Command Language Version 6.2.0
>>> The use of this software is governed by a License Agreement.
>>> See http://www.ncl.ucar.edu/ for more details.
>>> fatal:Subscript out of range, error in subscript #0
>>>
>>> But if you read the wind as this: uwind = file_in->U(:,:,:,:), there will be no error message.
>>>
>>> However, even with the error message, it seems I still can get the wind at the level I want.
>>>
>>>
>>> Any idea why this happens?
>>>
>>> Thanks.
>>>
>>> - Guangshan
>>>
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> 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
Received on Thu May 8 16:11:45 2014