[ncl-talk] write out in HDF5 format by NCL

[Dennis Shea]

Hello,NCL 6.1.2 was released February 7, 2013 ... 4.5 years ago

**Significant* *upgrades have been made to NCL's  internal data structures.
These are needed to handle features of HDF5.

Please upgrade to NCL 6.4.0 or, possibly, (beta) NCL 6.5.0

http://www.ncl.ucar.edu/Download/

ncl-install at ucar.edu

Cheers



On Tue, Aug 15, 2017 at 5:32 PM, Yi-Chih Huang <dscpln at gmail.com> wrote:

> David,
>
>      I am working on Version 6.1.2.  I created one variable to output in
> my NCL script instead of the fn file.  Thus, I commented all the fn and fi
> below for consistency.  I got the error message "NclFileAddGrp: Old File
> Structure DO NOT Support Group" over and over again and wonder how to fix
> this error?
>
>     Thanks,
>
>                     Yi-Chih
>
> ; fn = "/Users/huangwei/myncl/nc4_work/wrtgroup/nc4uvt.nc"
>
> ; fi = addfile(fn, "r")
>
> ;  fAtt at source_file   = fn
>
> ; time = fi->time
> ; lev = fi->lev
> ; lat = fi->lat
> ; lon = fi->lon
>
> ; t = fi->T
> ; u = fi->U
> ; v = fi->V
>
>
> On Wed, Aug 16, 2017 at 3:09 AM, David Brown <dbrown at ucar.edu> wrote:
>
>> Hi Yi-Chih,
>>
>> I tried running the script you pointed to using NCL 6.4.0 and there
>> were no errors. I am not sure what you mean
>> by "commented the three lines below". Which lines specifically are you
>> referring to? Also which version of NCL are you running?
>>  -dave
>>
>>
>> On Sat, Aug 12, 2017 at 3:10 AM, Yi-Chih Huang <dscpln at gmail.com> wrote:
>> > Hello,
>> >
>> >     I am trying to write out one variable in HDF5 format in an NCL
>> script
>> > following the commands on the web page
>> > https://www.ncl.ucar.edu/Applications/hdf5.shtml.  Because the output
>> data
>> > are calculated in the NCL script, I commented the three lines below.
>> > However, I got the error message "NclFileAddGrp: Old File Structure DO
>> NOT
>> > Support Group." as follows.  Could anyone indicate how to fix this
>> error?
>> >
>> >     Thanks much,
>> >
>> >                    Yi-Chih
>> >
>> >
>> > fatal:["FileSupport.c":1638]:_NclFileAddGrp: Old File Structure DO NOT
>> > Support Group.
>> > ^Mfatal:["FileSupport.c":1638]:_NclFileAddGrp: Old File Structure DO
>> NOT
>> > Support Group.
>> > ^Mfatal:["FileSupport.c":1638]:_NclFileAddGrp: Old File Structure DO
>> NOT
>> > Support Group.
>> > ^Mfatal:["Execute.c":8128]:Execute: Error occurred at or near line 340
>> in
>> > file q_ln-ln_pre_q01.ncl
>> >
>> > ##########################################
>> > ; fn = "/Users/huangwei/myncl/nc4_work/wrtgroup/nc4uvt.nc"
>> >
>> > ; fi = addfile(fn, "r")
>> >
>> > ;  fAtt at source_file   = fn
>> >
>> > ##########################################     as on the web page
>> > https://www.ncl.ucar.edu/Applications/hdf5.shtml
>> >
>> >       Write HDF5 file with group:
>> >
>> >  setfileoption("nc", "FileStructure", "Advanced")
>> >
>> >  fn = "/Users/huangwei/myncl/nc4_work/wrtgroup/nc4uvt.nc"
>> >
>> >  fi = addfile(fn, "r")
>> >
>> > ;printVarSummary(fi)
>> > ;print(fi)
>> >
>> >  time = fi->time
>> >  lev = fi->lev
>> >  lat = fi->lat
>> >  lon = fi->lon
>> >
>> >  t = fi->T
>> >  u = fi->U
>> >  v = fi->V
>> >
>> > ;printVarSummary(t)
>> > ;printVarSummary(u)
>> > ;printVarSummary(v)
>> >
>> > ;print("t(0,0,0,0) = " + t(0,0,0,0))
>> > ;print("u(0,1,1,1) = " + u(0,1,1,1))
>> > ;print("v(0,2,2,2) = " + v(0,2,2,2))
>> >
>> > ;print("t&lat(0) = " + t&lat(0))
>> > ;print("t at units = " + t at units)
>> > ;print("u at units = " + u at units)
>> > ;print("v at units = " + v at units)
>> >
>> > ;u1 = fi->U(::2)
>> > ;v1 = fi->V(3:121:3)
>> >
>> >  ntim  = dimsizes(time)       ; get dimension sizes
>> >  nlev  = dimsizes(lev)
>> >  nlat  = dimsizes(lat)
>> >  nlon  = dimsizes(lon)
>> >
>> > ;------------------------------------------------------------
>> >  setfileoption("h5", "FileStructure", "Advanced")
>> >
>> >  fon = "ncl_wrt_uvt.h5"
>> >  system("/bin/rm -f " + fon) ; remove if exists
>> >  fo = addfile(fon, "c")
>> >
>> > ;===================================================================
>> > ; explicitly declare file definition mode. Improve efficiency.
>> > ;===================================================================
>> > ;setfileoption(fo,"CompressionLevel", 6)
>> > ;setfileoption(fo,"CacheSize", 3200000)
>> > ;setfileoption(fo,"CacheNelems", 1027)
>> > ;setfileoption(fo,"CachePreemption", 0.25)
>> >
>> > ..............
>> >
>> >
>> > _______________________________________________
>> > ncl-talk mailing list
>> > ncl-talk at ucar.edu
>> > List instructions, subscriber options, unsubscribe:
>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>> >
>>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170815/5b1d52e7/attachment.html