Re: [ferret_users] cancatenating netCDF files using descriptor file

[Andrew Wittenberg <Andrew.Wittenberg@xxxxxxxx>]

Use the NCO attribute editor to fix the time units in the original NetCDF files:

ncatted -h -a units,time,o,c,"days since 1970-01-01 00:00:00" old.nc new.nc

See the NCO documentation for details:

   http://nco.sourceforge.net/nco.html#ncatted-netCDF-Attribute-Editor

Andrew

jagadish karmacharya wrote:
> Dear all,
>
> I tried to cancatenate netCDF file using make_des. But the resultant 
> fill is not loaded properly.
>
> # ./make_des APHRO_MA_025deg_V0902.1971.nc APHRO_MA_025deg_V0902.1972.nc 
> APHRO_MA_025deg_V0902.1973.nc > new.des
>
> yes? use new.des
>  Warning: missing or improper date/time format for T0 in new
>
> I think the time axis is not proper in the source files:
>
> # ncdump -h APHRO_MA_025deg_V0902.1970.nc
> netcdf APHRO_MA_025deg_V0902.1970 {
> dimensions:
>         time = UNLIMITED ; // (365 currently)
>         longitude = 360 ;
>         latitude = 220 ;
> variables:
>         double time(time) ;
>                 time:units = "days since 1970-1-1 0" ;
>         double longitude(longitude) ;
>                 longitude:units = "degrees_east" ;
>                 longitude:long_name = "Longitude" ;
>         double latitude(latitude) ;
>                 latitude:units = "degrees_north" ;
>                 latitude:long_name = "Latitude" ;
>         float precip(time, latitude, longitude) ;
>                 precip:long_name = "daily precipitation analysis 
> interpolated into 0.25deg grids [mm/day]" ;
>                 precip:units = "mm/day" ;
>                 precip:level_description = "Earth surface" ;
>                 precip:missing_value = -99.9f ;
>         float rstn(time, latitude, longitude) ;
>                 rstn:long_name = "ratio of 0.05deg-grids with station [%]" ;
>                 rstn:units = "%" ;
>                 rstn:level_description = "Earth surface" ;
>                 rstn:missing_value = -99.9f ;
>
> // global attributes:
>                 :Conventions = "COARDS" ;
>                 :calendar = "standard" ;
>                 :comments = "APHRO_MA V0902 daily precipitation with 
> 0.25deg grids" ;
>                 :history = "Created by APHRODITE project team 08-Jun-2009" ;
> }
>
> # ncdump -h APHRO_MA_025deg_V0902.1971.nc
> netcdf APHRO_MA_025deg_V0902.1971 {
> dimensions:
>         time = UNLIMITED ; // (365 currently)
>         longitude = 360 ;
>         latitude = 220 ;
> variables:
>         double time(time) ;
>                 time:units = "days since 1971-1-1 0" ;
>                (I have deleted remaining lines)............................
>
>
> Please suggest me how to solve this problem.
>
> Best regards,
> Jagadish
>
>
>
>     <http://us.mc319.mail.yahoo.com/mc/showMessage?sMid=199&fid=Useful%2520ferret%2520mail&sort=date&order=down&startMid=175&filterBy=&.rand=801499036&midIndex=24&mid=1_117539_ACHOjkQAAXidQfA7kwyupT2I4BA>