[GRASS-SVN] r52065 - in grass/trunk/raster3d/r3.out.netcdf: . test_suite
svn_grass at osgeo.org
svn_grass at osgeo.org
Thu Jun 14 02:20:11 PDT 2012
Author: huhabla
Date: 2012-06-14 02:20:10 -0700 (Thu, 14 Jun 2012)
New Revision: 52065
Modified:
grass/trunk/raster3d/r3.out.netcdf/main.c
grass/trunk/raster3d/r3.out.netcdf/r3.out.netcdf.html
grass/trunk/raster3d/r3.out.netcdf/test_suite/test.r3.out.netcdf.sh
Log:
Code indention and documentation update
Modified: grass/trunk/raster3d/r3.out.netcdf/main.c
===================================================================
--- grass/trunk/raster3d/r3.out.netcdf/main.c 2012-06-14 06:35:44 UTC (rev 52064)
+++ grass/trunk/raster3d/r3.out.netcdf/main.c 2012-06-14 09:20:10 UTC (rev 52065)
@@ -1,11 +1,11 @@
/****************************************************************************
*
- * MODULE: r3.out.netcdf
+ * MODULE: r3.out.netCDF
*
* AUTHOR(S): Soeren Gebbert
*
- * PURPOSE: Export a 3D raster map as netcdf file
+ * PURPOSE: Export a 3D raster map as netCDF file
*
* COPYRIGHT: (C) 2012 by the GRASS Development Team
*
@@ -51,7 +51,7 @@
#define UNITS "units"
#define DEGREES_EAST "degrees_east"
#define DEGREES_NORTH "degrees_north"
-#define HISTORY_TEXT "GRASS GIS 7 NetCDF export of r3.out.netcdf"
+#define HISTORY_TEXT "GRASS GIS 7 netCDF export of r3.out.netcdf"
#define CF_SUPPORT "CF-1.5"
#define ERR(e) {fatalError(nc_strerror(e));}
@@ -94,7 +94,7 @@
param.output = G_define_standard_option(G_OPT_F_OUTPUT);
param.output->key = "output";
- param.output->description = _("Name for netcdf output file");
+ param.output->description = _("Name for netCDF output file");
param.null = G_define_option();
param.null->key = "null";
@@ -103,7 +103,7 @@
param.null->multiple = NO;
param.null->description =
_
- ("The value to be used for null values, default is the NetCDF standard");
+ ("The value to be used for null values, default is the netCDF standard");
param.proj = G_define_flag();
param.proj->key = 'p';
@@ -584,9 +584,9 @@
G_gisinit(argv[0]);
module = G_define_module();
G_add_keyword(_("raster3d"));
- G_add_keyword(_("netcdf"));
+ G_add_keyword(_("netCDF"));
G_add_keyword(_("export"));
- module->description = _("Export a 3D raster map as netcdf file.");
+ module->description = _("Export a 3D raster map as netCDF file.");
/* Get parameters from user */
setParams();
@@ -616,7 +616,7 @@
G_fatal_error(_("Error opening 3d raster map <%s>"),
param.input->answer);
- /* Create netcdf file */
+ /* Create netCDF file */
if ((retval = nc_create(param.output->answer, NC_CLOBBER, &ncid)))
ERR(retval);
@@ -647,7 +647,7 @@
if (!Rast3d_close(map))
fatalError(_("Unable to close 3D raster map"));
- /* Close the netcdf file */
+ /* Close the netCDF file */
if ((retval = nc_close(ncid)))
ERR(retval);
Modified: grass/trunk/raster3d/r3.out.netcdf/r3.out.netcdf.html
===================================================================
--- grass/trunk/raster3d/r3.out.netcdf/r3.out.netcdf.html 2012-06-14 06:35:44 UTC (rev 52064)
+++ grass/trunk/raster3d/r3.out.netcdf/r3.out.netcdf.html 2012-06-14 09:20:10 UTC (rev 52065)
@@ -1,76 +1,80 @@
<h2>DESCRIPTION</h2>
-The module <em>r3.out.netcdf</em> exports a 3D raster map as NetCDF file.
+The module <em>r3.out.netcdf</em> exports a 3D raster map as netCDF file.
Maps are valid 3D raster maps in the current mapset search path. The
-<em>output</em> parameter is the name of a NetCDF file that will be
+<em>output</em> parameter is the name of a netCDF file that will be
written in the current working directory. The module is sensitive
to mapset region settings (set with <em>g.region</em>) and the vertical unit settings
of the 3D raster map. The vertical unit can be set with <a href="r3.support.html">r3.support</a>.
<h2>Note</h2>
-The resulting NetCDF file will have time as third dimension in case a temporal
+The resulting netCDF file will have time as third dimension in case a temporal
vertical unit is specified for the input map, otherwise the third dimension is assumed as spatial.
If the 3D raster map has a time stamp with absolute time, the start time will be
-used as starting point for the NetCDF time series. If the absolute time stamp is missing
+used as starting point for the netCDF time series. If the absolute time stamp is missing
the date "Jan. 01. 1900 00:00:00" will be used as default. In case of relative time stamp no start time
-is set in the NetCDF file. The number of depths of the 3D raster map are the number of time slices
-in the NetCDF file. You can use <em>t.rast.to.rast3</em> to convert a
+is set in the netCDF file. The number of depths of the 3D raster map are the number of time slices
+in the netCDF file. You can use <em>t.rast.to.rast3</em> to convert a
space time raster dataset into a 3D raster map representing a space time voxel
-cube, then export the raster time series as NetCDF file.
+cube, then export the raster time series as netCDF file.
<p>
-Spatial coordinates are exported as cell centered coordinates. The projection is stored
-in the metadata as <em>crs</em> attributes using a WKT parameter string.
-The export of projection parameters is suppressed when the XY-projection is set.
-The NetCDF projection metadata storage follows the spatial_ref GDAL/NetCDF suggestion
+Spatial coordinates are exported as cell centered coordinates. The projection can be optionally stored
+in the metadata as <em>crs</em> attributes .
+The netCDF projection metadata storage follows the spatial_ref GDAL/netCDF suggestion
<a href="http://www.gdal.org/frmt_netcdf.html">here</a>
-and the NetCDF CF 1.6 convention
+and the netCDF CF 1.6 convention
<a href="http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.6/cf-conventions.html#coordinate-system">here</a>
using WKT projection information. Additional a PROJ.4 string is
stored in the <em>crs</em> attribute section.
+The export of projection parameters is suppressed when the XY-projection is set.
<p>
-The NetCDF library must be installed on the system and activated at configuration time to compile this module.
+The range of the 3D raster map is set in the netCDF output file. Optionally a the netCDF missing value and
+_FillValue can be set using the <em>null</em> option. As default the netCDF floating point _FillValues are used
+to fill empty cell.
+<p>
+The netCDF library must be installed on the system and activated at configuration time to compile this module.
<h2>Examples</h2>
-<h3>Using lat/lon projection without vertical unit specification</h3>
+<div class="code"><pre>
+g.region s=-90 n=90 w=-180 e=180 b=0 t=5 res=10 res3=10 tbres=1 -p3
-<p>Lets create a 3D raster map without the specification of a vertical unit</p>
+#####################################################################
+# We create a simple volume map with floating point values
+#####################################################################
-<div class="code"><pre>
-g.region s=-90 n=90 w=-180 e=180 b=0 t=50 res=10 res3=2 -p3
-r3.mapcalc expr="volume_double = double(col() + row() + depth())"
-r3.info
-
+r3.mapcalc --o expr="volume_float = float(col() + row() + depth())"
+r3.info volume_float
+----------------------------------------------------------------------------+
- | Layer: volume_double Date: Wed May 16 20:08:22 2012 |
+ | Layer: volume_float Date: Thu Jun 14 08:40:56 2012 |
| Mapset: PERMANENT Login of Creator: soeren |
| Location: TestLL |
| DataBase: /1/soeren/grassdata |
- | Title: volume_double |
+ | Title: volume_float |
| Units: none |
| Vertical unit: units |
| Timestamp: none |
|----------------------------------------------------------------------------|
| |
| Type of Map: 3d cell Number of Categories: 0 |
- | Data Type: DCELL |
- | Rows: 90 |
- | Columns: 180 |
- | Depths: 25 |
- | Total Cells: 405000 |
- | Total size: 2835864 Bytes |
- | Number of tiles: 108 |
- | Mean tile size: 26258 Bytes |
- | Tile size in memory: 32400 Bytes |
- | Number of tiles in x, y and z: 6, 6, 3 |
- | Dimension of a tile in x, y, z: 30, 15, 9 |
+ | Data Type: FCELL |
+ | Rows: 18 |
+ | Columns: 36 |
+ | Depths: 5 |
+ | Total Cells: 3240 |
+ | Total size: 10114 Bytes |
+ | Number of tiles: 1 |
+ | Mean tile size: 10114 Bytes |
+ | Tile size in memory: 12960 Bytes |
+ | Number of tiles in x, y and z: 1, 1, 1 |
+ | Dimension of a tile in x, y, z: 36, 18, 5 |
| |
| Projection: Latitude-Longitude (zone 0) |
- | N: 90N S: 90S Res: 2 |
- | E: 180E W: 180W Res: 2 |
- | T: 50 B: 0 Res: 2 |
- | Range of data: min = 3 max = 295 |
+ | N: 90N S: 90S Res: 10 |
+ | E: 180E W: 180W Res: 10 |
+ | T: 5 B: 0 Res: 1 |
+ | Range of data: min = 3 max = 59 |
| |
| Data Source: |
| |
@@ -80,147 +84,171 @@
| generated by r3.mapcalc |
| |
| Comments: |
- | r3.mapcalc expression="volume_double = double(col() + row() + depth(\ |
- | ))" |
+ | r3.mapcalc expression="volume_float = float(col() + row() + depth())" |
| |
+----------------------------------------------------------------------------+
-</pre></div>
-<p>Export the map and show the NetCDF header</P>
+#####################################################################
+# We use the netCDF ncdump tool to have a look at the header
+# and coordinates of the exported netCDF file
+#####################################################################
-<div class="code"><pre>
-r3.out.netcdf input=volume_double output=test_double.nc
-ncdump -h test_double.nc
-
-netcdf test_double {
+r3.out.netcdf --o input=volume_float output=test_float.nc
+ncdump -c test_float.nc
+
+netcdf test_float {
dimensions:
- longitude = 180 ;
- latitude = 90 ;
- z = 25 ;
+ longitude = 36 ;
+ latitude = 18 ;
+ z = 5 ;
variables:
- char crs ;
- crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
- crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
- crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
- float longitude(longitude) ;
- longitude:units = "degrees_east" ;
- longitude:long_name = "Longitude values" ;
- longitude:standard_name = "longitude" ;
- longitude:axis = "X" ;
- float latitude(latitude) ;
- latitude:units = "degrees_north" ;
- latitude:long_name = "Latitude values" ;
- latitude:standard_name = "latitude" ;
- latitude:axis = "Y" ;
- float z(z) ;
- z:units = "meter" ;
- z:long_name = "z coordinate of projection" ;
- z:standard_name = "projection_z_coordinate" ;
- z:positive = "up" ;
- z:axis = "Z" ;
- double volume_double(z, latitude, longitude) ;
- volume_double:grid_mapping = "crs" ;
+ float longitude(longitude) ;
+ longitude:units = "degrees_east" ;
+ longitude:long_name = "Longitude values" ;
+ longitude:standard_name = "longitude" ;
+ longitude:axis = "X" ;
+ float latitude(latitude) ;
+ latitude:units = "degrees_north" ;
+ latitude:long_name = "Latitude values" ;
+ latitude:standard_name = "latitude" ;
+ latitude:axis = "Y" ;
+ float z(z) ;
+ z:units = "meter" ;
+ z:long_name = "z coordinate of projection" ;
+ z:standard_name = "projection_z_coordinate" ;
+ z:positive = "up" ;
+ z:axis = "Z" ;
+ float volume_float(z, latitude, longitude) ;
+ volume_float:valid_min = 3.f ;
+ volume_float:valid_max = 59.f ;
// global attributes:
- :Conventions = "CF-1.5" ;
- :history = "GRASS GIS 7 NetCDF export of r3.out.netcdf" ;
+ :Conventions = "CF-1.5" ;
+ :history = "GRASS GIS 7 netCDF export of r3.out.netcdf" ;
+data:
+
+ longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
+ -65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
+ 105, 115, 125, 135, 145, 155, 165, 175 ;
+
+ latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
+ -65, -75, -85 ;
+
+ z = 0.5, 1.5, 2.5, 3.5, 4.5 ;
}
-</pre></div>
-<p>Now we create a 3D raster map with a temporal vertical unit</p>
+#####################################################################
+# Exporting the projection specific settings and a null value
+#####################################################################
-<div class="code"><pre>
-g.region s=-90 n=90 w=-180 e=180 b=0 t=50 res=10 res3=2 -p3
-r3.mapcalc expr="volume_double = double(col() + row() + depth())"
-r3.timestamp map=volume_time_double date='1 Jan 2001'
-r3.support map=volume_time_double vunit="days"
-r3.info
+r3.out.netcdf --o -p null=-1 input=volume_float output=test_float.nc
+ncdump -c test_float.nc
- +----------------------------------------------------------------------------+
- | Layer: volume_time_double Date: Wed May 16 20:26:28 2012 |
- | Mapset: PERMANENT Login of Creator: soeren |
- | Location: TestLL |
- | DataBase: /1/soeren/grassdata |
- | Title: volume_time_double |
- | Units: none |
- | Vertical unit: days |
- | Timestamp: 1 Jan 2001 |
- |----------------------------------------------------------------------------|
- | |
- | Type of Map: 3d cell Number of Categories: 0 |
- | Data Type: DCELL |
- | Rows: 90 |
- | Columns: 180 |
- | Depths: 25 |
- | Total Cells: 405000 |
- | Total size: 2835864 Bytes |
- | Number of tiles: 108 |
- | Mean tile size: 26258 Bytes |
- | Tile size in memory: 32400 Bytes |
- | Number of tiles in x, y and z: 6, 6, 3 |
- | Dimension of a tile in x, y, z: 30, 15, 9 |
- | |
- | Projection: Latitude-Longitude (zone 0) |
- | N: 90N S: 90S Res: 2 |
- | E: 180E W: 180W Res: 2 |
- | T: 50 B: 0 Res: 2 |
- | Range of data: min = 3 max = 295 |
- | |
- | Data Source: |
- | |
- | |
- | |
- | Data Description: |
- | generated by r3.mapcalc |
- | |
- | Comments: |
- | r3.mapcalc expression="volume_time_double = double(col() + row() + d\ |
- | epth())" |
- | |
- +----------------------------------------------------------------------------+
+netcdf test_float {
+dimensions:
+ longitude = 36 ;
+ latitude = 18 ;
+ z = 5 ;
+variables:
+ char crs ;
+ crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
+ crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
+ crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
+ float longitude(longitude) ;
+ longitude:units = "degrees_east" ;
+ longitude:long_name = "Longitude values" ;
+ longitude:standard_name = "longitude" ;
+ longitude:axis = "X" ;
+ float latitude(latitude) ;
+ latitude:units = "degrees_north" ;
+ latitude:long_name = "Latitude values" ;
+ latitude:standard_name = "latitude" ;
+ latitude:axis = "Y" ;
+ float z(z) ;
+ z:units = "meter" ;
+ z:long_name = "z coordinate of projection" ;
+ z:standard_name = "projection_z_coordinate" ;
+ z:positive = "up" ;
+ z:axis = "Z" ;
+ float volume_float(z, latitude, longitude) ;
+ volume_float:valid_min = 3.f ;
+ volume_float:valid_max = 59.f ;
+ volume_float:missing_value = -1.f ;
+ volume_float:_FillValue = -1.f ;
+ volume_float:grid_mapping = "crs" ;
-</pre></div>
-<p>Export the map and show the NetCDF header</P>
+// global attributes:
+ :Conventions = "CF-1.5" ;
+ :history = "GRASS GIS 7 netCDF export of r3.out.netcdf" ;
+data:
-<p>NetCDF header of an exported 3D raster map with vertical temporal unit in a Lat/Lon projection:</P>
+ longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
+ -65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
+ 105, 115, 125, 135, 145, 155, 165, 175 ;
-<div class="code"><pre>
-r3.out.netcdf --o input=volume_time_double output=test_time_double.nc
-ncdump -h test_time_double.nc
+ latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
+ -65, -75, -85 ;
-netcdf test_time_double {
+ z = 0.5, 1.5, 2.5, 3.5, 4.5 ;
+}
+
+#####################################################################
+# Assigning time as vertical unit and setting an absolute time stamp
+#####################################################################
+
+r3.timestamp map=volume_float date='1 Jan 2001/5 Jan 2001'
+r3.support map=volume_float vunit="days"
+r3.out.netcdf --o -p null=-1 input=volume_float output=test_float.nc
+ncdump -c test_float.nc
+
+netcdf test_float {
dimensions:
- longitude = 180 ;
- latitude = 90 ;
- time = 25 ;
+ longitude = 36 ;
+ latitude = 18 ;
+ time = 5 ;
variables:
- char crs ;
- crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
- crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
- crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
- float longitude(longitude) ;
- longitude:units = "degrees_east" ;
- longitude:long_name = "Longitude values" ;
- longitude:standard_name = "longitude" ;
- longitude:axis = "X" ;
- float latitude(latitude) ;
- latitude:units = "degrees_north" ;
- latitude:long_name = "Latitude values" ;
- latitude:standard_name = "latitude" ;
- latitude:axis = "Y" ;
- int time(time) ;
- time:units = "days since 2001-01-01 00:00:00" ;
- time:long_name = "Time in days" ;
- time:calendar = "gregorian" ;
- time:positive = "up" ;
- time:axis = "T" ;
- double volume_time_double(time, latitude, longitude) ;
- volume_time_double:grid_mapping = "crs" ;
+ char crs ;
+ crs:crs_wkt = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
+ crs:spatial_ref = "GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]]" ;
+ crs:crs_proj4 = " +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000" ;
+ float longitude(longitude) ;
+ longitude:units = "degrees_east" ;
+ longitude:long_name = "Longitude values" ;
+ longitude:standard_name = "longitude" ;
+ longitude:axis = "X" ;
+ float latitude(latitude) ;
+ latitude:units = "degrees_north" ;
+ latitude:long_name = "Latitude values" ;
+ latitude:standard_name = "latitude" ;
+ latitude:axis = "Y" ;
+ int time(time) ;
+ time:units = "days since 2001-01-01 00:00:00" ;
+ time:long_name = "Time in days" ;
+ time:calendar = "gregorian" ;
+ time:positive = "up" ;
+ time:axis = "T" ;
+ float volume_float(time, latitude, longitude) ;
+ volume_float:valid_min = 3.f ;
+ volume_float:valid_max = 59.f ;
+ volume_float:missing_value = -1.f ;
+ volume_float:_FillValue = -1.f ;
+ volume_float:grid_mapping = "crs" ;
// global attributes:
- :Conventions = "CF-1.5" ;
- :history = "GRASS GIS 7 NetCDF export of r3.out.netcdf" ;
+ :Conventions = "CF-1.5" ;
+ :history = "GRASS GIS 7 netCDF export of r3.out.netcdf" ;
+data:
+
+ longitude = -175, -165, -155, -145, -135, -125, -115, -105, -95, -85, -75,
+ -65, -55, -45, -35, -25, -15, -5, 5, 15, 25, 35, 45, 55, 65, 75, 85, 95,
+ 105, 115, 125, 135, 145, 155, 165, 175 ;
+
+ latitude = 85, 75, 65, 55, 45, 35, 25, 15, 5, -5, -15, -25, -35, -45, -55,
+ -65, -75, -85 ;
+
+ time = 0, 1, 2, 3, 4 ;
}
+
</pre></div>
<h2>SEE ALSO</h2>
Modified: grass/trunk/raster3d/r3.out.netcdf/test_suite/test.r3.out.netcdf.sh
===================================================================
--- grass/trunk/raster3d/r3.out.netcdf/test_suite/test.r3.out.netcdf.sh 2012-06-14 06:35:44 UTC (rev 52064)
+++ grass/trunk/raster3d/r3.out.netcdf/test_suite/test.r3.out.netcdf.sh 2012-06-14 09:20:10 UTC (rev 52065)
@@ -13,21 +13,22 @@
r3.mapcalc --o expr="volume_double = double(col() + row() + depth())"
r3.mapcalc --o expr="volume_time_double = double(col() + row() + depth())"
r3.mapcalc --o expr="volume_time_float = float(col() + row() + depth())"
-r3.timestamp map=volume_time_double date='1 Jan 1900/5 Jan 1900'
+r3.timestamp map=volume_time_double date='1 Jan 2001/5 Jan 2001'
r3.support map=volume_time_double vunit="days"
r3.timestamp map=volume_time_float date='5 seconds/10 seconds'
r3.support map=volume_time_float vunit="seconds"
# @test
r3.out.netcdf --o input=volume_float output=test_float.nc
-#r3.info volume_float
-#ncdump -h test_float.nc
+r3.info volume_float
+ncdump -h test_float.nc
r3.out.netcdf --o null=-100 input=volume_double output=test_double.nc
-#r3.info volume_double
-#ncdump -h test_double.nc
+r3.info volume_double
+ncdump -h test_double.nc
r3.out.netcdf --o -p input=volume_time_double output=test_time_double.nc
-#r3.info volume_time_double
-#ncdump -h test_time_double.nc
+r3.info volume_time_double
+ncdump -h test_time_double.nc
r3.out.netcdf --o -p null=-1000 input=volume_time_float output=test_time_float.nc
-#r3.info volume_time_float
-#ncdump -h test_time_float.nc
+r3.info volume_time_float
+ncdump -h test_time_float.nc
+g.remove rast3=volume_float,volume_double,volume_time_double,volume_time_float
\ No newline at end of file
More information about the grass-commit
mailing list