[R-sig-Geo] Rasterbrick to netCDF

[ppanday]

Hello

I have an output from R-script in the form of a rasterbrick with 3
dimensions (time, latitude, longitude). I want to be able to write it to a
netCDF format, add another dimension (time, level, latitude, longitude), and
be able to edit dimension attributes.
Below you will see the properties when I write my rasterbrick to a netCDF
format using writeRaster.
I would like to change 'value' dimension to time, and edit time:units as
well.
Alternatively, is there a way I can copy dimensions and attributes from
another netCDF template and put my rasterbrick in that empty netCDF file?
Any help is greatly appreciated.

Sincerely,
Prajjwal


dimensions:
	longitude = 720 ;
	latitude = 360 ;
	value = UNLIMITED ; // (12 currently)
variables:
	double longitude(longitude) ;
		longitude:units = "degrees_east" ;
		longitude:long_name = "longitude" ;
	double latitude(latitude) ;
		latitude:units = "degrees_north" ;
		latitude:long_name = "latitude" ;
	int value(value) ;
		value:units = "unknown" ;
		value:long_name = "value" ;
	float prec(value, latitude, longitude) ;
		prec:_FillValue = -3.4e+38 ;
		prec:missing_value = -3.4e+38 ;
		prec:long_name = "prec" ;
		prec:projection = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
;
		prec:projection_format = "PROJ.4" ;
		prec:min = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0. ;
		prec:max = 663.033325195312, 553.867553710938, 595.225219726562,
577.460388183594, 688.026123046875, 1178.05859375, 1098.39819335938,
1269.41442871094, 1184.70092773438, 1013.62524414062, 755.313537597656,
633.667541503906 ;

// global attributes:
		:Conventions = "CF-1.4" ;
		:created_by = "R, packages ncdf and raster (version 2.0-31)" ;
		:date = "2013-01-24 15:58:09" ;







--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Rasterbrick-to-netCDF-tp7582376.html
Sent from the R-sig-geo mailing list archive at Nabble.com.