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

Erika,

Please split line:

  fout->time_bnds = f->time_bnds

to:

  time_bnds = f->time_bnds
  fout->time_bnds = time_bnds

and see what happened.

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 9:35 AM, Erika Folova <e.folova@gmail.com> wrote:

> Thank you Wei for your concern.. but I still got the same error!!
> "fatal:Dimension sizes of left hand side do not match right hand side
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 67 in file step2.ncl"
>
> Here I attached my scrip and printVarSummaries:
>
>
> Variable: data
> Type: float
> Total Size: 10212618240 bytes
> 2553154560 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 16060] x [level | 23] x [lat | 48] x [lon | 144]
> Coordinates:
> time: [1095..17154]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> lon: [ 0..357.5]
> Number Of Attributes: 3
> _FillValue : 9.96921e+36
> long_name : Deviation from mean
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
>
> Variable: climo
> Type: float
> Total Size: 232740864 bytes
> 58185216 values
> Number of Dimensions: 4
> Dimensions and sizes: [year_day | 366] x [level | 23] x [lat | 48] x [lon | 144]
> Coordinates:
> year_day: [1..366]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> lon: [ 0..357.5]
> Number Of Attributes: 5
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
> long_name : Deviation from mean
> _FillValue : 9.96921e+36
> information : Raw daily averages across all years
> smoothing : None
>
> Variable: climoHarm
> Type: float
> Total Size: 232740864 bytes
> 58185216 values
> Number of Dimensions: 4
> Dimensions and sizes: [year_day | 366] x [level | 23] x [lat | 48] x [lon | 144]
> Coordinates:
> year_day: [1..366]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> lon: [ 0..357.5]
> Number Of Attributes: 5
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
> long_name : Deviation from mean
> _FillValue : 9.96921e+36
> smoothing : FFT: 3 harmonics were retained.
> information : Smoothed daily climatological averages
>
> Variable: data
> Type: float
> Total Size: 10212618240 bytes
> 2553154560 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 16060] x [level | 23] x [lat | 48] x [lon | 144]
> Coordinates:
> time: [1095..17154]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> lon: [ 0..357.5]
> Number Of Attributes: 3
> long_name : Anomalies: Deviation from mean
> _FillValue : 9.96921e+36
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
> (0) 3556.132
> (0) -1506.295
>
> Variable: data
> Type: float
> Total Size: 10212618240 bytes
> 2553154560 values
> Number of Dimensions: 4
> Dimensions and sizes: [time | 16060] x [level | 23] x [lat | 48] x [lon | 144]
> Coordinates:
> time: [1095..17154]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> lon: [ 0..357.5]
> Number Of Attributes: 3
> long_name : Anomalies: Deviation from mean
> _FillValue : 9.96921e+36
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
> 0) 3541.621
> (0) -1521.45
>
> Variable: zm
> Type: float
> Total Size: 70920960 bytes
> 17730240 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 16060] x [level | 23] x [lat | 48]
> Coordinates:
> time: [1095..17154]
> level: [1000.. 1]
> lat: [0.9473684210526265..89.99999999999999]
> Number Of Attributes: 4
> long_name : Anomalies: Deviation from mean
> _FillValue : 9.96921e+36
> rmvmean_op_NCL : dim_rmvmean_n over dimension(s): lon
> average_op_ncl : dim_avg_n over dimension(s): lon
> fatal:Dimension sizes of left hand side do not match right hand side
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 67 in file step2.ncl
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> On Tue, Sep 10, 2013 at 4:41 PM, Wei Huang <huangwei@ucar.edu> wrote:
> 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
>
>
> <output.ncl>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 10 10:06:42 2013