Re: always get ERROR to save "time_bnds" in a large netcdf file

From: Wei Huang <huangwei_at_nyahnyahspammersnyahnyah>
Date: Tue Sep 10 2013 - 08:41:53 MDT

Erika,

Can you try comment out line:

setfileoption("nc","Format","NetCDF4Classic")

and re-run your code?

Thanks,

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Sep 10, 2013, at 6:48 AM, Erika Folova <e.folova@gmail.com> wrote:

> Dear Community and NCL developer,
>
> Why do I always get error every time I need to save a "time_bnds" in a large netcdf file:
> the error is:
>
> fatal:Dimension sizes of left hand side do not match right hand side
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 66 in file step2.ncl
>
> the way I save the output in netcdf format is as follow:
>
> setfileoption("nc","Format","NetCDF4Classic")
> system ("/bin/rm -f heatflux.anom.nc") ; remove any pre-existing file
> fout = addfile("heatflux.anom.nc" ,"c") ; open output netCDF file
> fout->time_bnds = f->time_bnds ; double
> fout->anom = zm ; [time, level, lat], float
>
> end
>
>
> Information of input file is:
>
> netcdf hf_total_final {
> dimensions:
> time = 20440 ;
> level = 23 ;
> lat = 48 ;
> lon = 144 ;
> nb2 = 2 ;
> variables:
> float VT(time, level, lat, lon) ;
> VT:rmvmean_op_NCL = "dim_rmvmean_n over dimension(s): lon" ;
> VT:long_name = "Deviation from mean" ;
> VT:_FillValue = 9.96921e+36f ;
> double time(time) ;
> time:standard_name = "time" ;
> time:long_name = "time" ;
> time:bounds = "time_bnds" ;
> time:units = "days since 1955-01-01 00:00:00" ;
> time:calendar = "365_day" ;
> float level(level) ;
> level:standard_name = "air_pressure" ;
> level:long_name = "pressure level" ;
> level:units = "mb" ;
> level:positive = "down" ;
> level:axis = "Z" ;
> double lat(lat) ;
> lat:standard_name = "latitude" ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:axis = "Y" ;
> double lon(lon) ;
> lon:standard_name = "longitude" ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> lon:axis = "X" ;
> double time_bnds(time, nb2) ;
> time_bnds:units = "days since 1955-01-01 00:00:00" ;
> time_bnds:calendar = "365_day" ;
>
>
>
>
>
>
> Thank you
> -- Erika
> _______________________________________________
> 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 Tue Sep 10 08:42:07 2013

This archive was generated by hypermail 2.1.8 : Mon Sep 16 2013 - 13:43:52 MDT