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

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
Received on Tue Sep 10 06:48:18 2013