Re: Netcdf output with non-numeric datatyep

Hi Wei,

This is very interesting idea, but unfortunately, there is something wrong with it.

First of all, the code you shown me can be ran errorlessly. However, the output "out.nc" cannot be recognised, i.g., when I try *ncdump out.nc*, the error message is:
ncdump: out.nc: NetCDF: Unknown file format

And, when I try to read *out.nc* in NCL using code like:
fi = addfile ("out.nc","r")
allmode = fi->allmode
 
The error is *Memory fault*.

This is bizarre as NCL can create the file, but the file cannot be recognised by itself.

Any idea about it?

Regards,
Liang

On 14 Apr 2014, at 21:00, Wei Huang <huangwei@ucar.edu> wrote:

> Liang,
>
> You may try something like:
>
> fout = "out.nc"
> system("/bin/rm -f "+fout)
> setfileoption( "nc", "format", "netcdf4" )
> ncdf = addfile(fout,"c")
>
> mode = (/"Original value","Component 1","Component 2","Component 3","Component 4"/)
> time = (/0, 1, 2/)
>
> mode@long_name = "mode"
> mode@units = "M"
>
> mode!0 = "mode"
> mode&mode = mode
>
> time!0 = "time"
> time&time = time
>
> time@long_name = "time"
> time@units = "day"
>
> print("mode=" + mode)
> print("time=" + time)
>
> allmode = new((/dimsizes(mode), dimsizes(time)/), float)
> allmode = 9.9
>
> allmode@long_name = "allmode"
> allmode@units = "mode/day"
>
> allmode!0 = "mode"
> allmode!1 = "time"
> allmode&mode = mode
> allmode&time = time
>
> ncdf->allmode = allmode
>
> delete(ncdf)
>
>
> Wei
>
>
> On Apr 7, 2014, at 8:18 AM, Liang Guo (Gmail) <liang.guo.81@gmail.com> wrote:
>
>> Dear NCL,
>>
>> I have a variable named *allmode(mode,time)*. For the dimension *mode*, I'd like to assign string value (/"Original value","Component 1","Component 2","Component 3","Component 4"/) instead of numeric value (/0,1,2,3,4/). Then, I try to write this variable into a Netcdf file. I wrote:
>> > fout = "out.nc"
>> > system("/bin/rm -f "+fout)
>> > setfileoption( "nc", "format", "netcdf4" )
>> > ncdf = addfile(fout,"c")
>> > ncdf->mode = allmode
>>
>> However, file can not be created successfully, with error information as:
>> warning:FileReadVarAtt: _FillValue attribute type differs from variable (mode) type in file (out), forcing type conversion; may result in overflow and/or loss of precision
>> ncclose: ncid 393216: NetCDF: Can't open HDF5 attribute
>>
>> Does any has ideas how to write this file out?
>>
>> Regards,
>> Liang
>>
>> _______________________________________________
>> 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 Fri May 02 05:29:13 2014