isunlimited returns False on unlimited dimension in NetCDF4

Hi,

I have converted NetCDF3 model output to NetCDF4 with the command 'ncks
-4 ...'.

ncl's isunlimited function returns False on the the "time" dimension in
the resulting file, despite output from running ncdump -h on the file
indicating that the "time" dimension is unlimited.

The following ncl commands demonstrate this behavior, with ncl 6.1.2, as
well as with the nightly 6.2.0 build.

    case = "b40.20th.track1.1deg.005"
    wdir = "/glade/scratch/klindsay/"+case+"/ocn/hist/"

    fname = wdir+case+".pop.h2.1990-01-01.nc"
    f = addfile(fname, "r")
    print(isunlimited(f, "time")) ; prints True

    fname = wdir+case+".pop.h2.cdf4.1990-01-01.nc"
    f = addfile(fname, "r")
    print(isunlimited(f, "time")) ; prints False

Could this be looked into please?

Thanks,
Keith Lindsay
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 21 14:12:37 2013