<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<font face="DejaVu Sans Mono" size="2">Thomas<br>
<br>
I have not tested, but have you looked into the "createopt"
parameter in r.out.gdal?<br>
<br>
From the manual:<br>
"createopt=string[,string,...]<br>
Creation option(s) to pass to the output format driver<br>
In the form of "NAME=VALUE", separate multiple entries with a
comma"<br>
<br>
So I believe that in principle you could pass the "BAND_NAMES"
creation option from GDAL.<br>
<br>
See:<br>
<br>
<a class="moz-txt-link-freetext" href="https://gdal.org/en/stable/drivers/raster/netcdf.html#creation-options">https://gdal.org/en/stable/drivers/raster/netcdf.html#creation-options</a><br>
<br>
<a class="moz-txt-link-freetext" href="https://grass.osgeo.org/grass85/manuals/r.out.gdal.html">https://grass.osgeo.org/grass85/manuals/r.out.gdal.html</a><br>
<br>
<br>
Hope this helps<br>
<br>
Hernán<br>
<br>
<br>
</font><br>
<div class="moz-cite-prefix">Den 2025-02-25 kl. 17:26, skrev Thomas
Adams via grass-user:<br>
</div>
<blockquote type="cite"
cite="mid:CAGxgkWj1aiLJB1yKoH2P4+WRb9B_rn1zq9-hcBevqOcV+1mhhA@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr">
<div>
<div>Hi all,</div>
<div><br>
</div>
<div>I really don't have a problem with exporting a netCDF
file from GRASS 8.4, but because of the intended use of the
output file, I need to have variable name defined
specifically, instead of the generic Band1:</div>
<div><br>
</div>
<div>ncdump -h sac_UZTWM.nc</div>
<div><br>
</div>
<div>netcdf sac_UZTWM {<br>
dimensions:<br>
x = 2698 ;<br>
y = 2427 ;<br>
variables:<br>
double Band1(y, x) ;<br>
Band1:long_name = "GDAL Band Number 1" ;<br>
Band1:_FillValue = -999. ;<br>
Band1:grid_mapping = "lambert_conformal_conic" ;<br>
Band1:COLOR_TABLE_RULES_COUNT = 255 ;<br>
Band1:COLOR_TABLE_RULE_RGB_0 = "0.000000e+00 5.496153e-01 68
1 84 68 2 86" ;<br>
Band1:COLOR_TABLE_RULE_RGB_1 = "5.496153e-01 1.099231e+00 68
2 86 69 4 87" ;<br>
Band1:COLOR_TABLE_RULE_RGB_2 = "1.099231e+00 1.648846e+00 69
4 87 69 5 89" ;<br>
Band1:COLOR_TABLE_RULE_RGB_3 = "1.648846e+00 2.199863e+00 69
5 89 70 7 90" ;<br>
Band1:COLOR_TABLE_RULE_RGB_4 = "2.199863e+00 2.749479e+00 70
7 90 70 8 92" ;<br>
Band1:COLOR_TABLE_RULE_RGB_5 = "2.749479e+00 3.299094e+00 70
8 92 70 10 93" ;<br>
Band1:COLOR_TABLE_RULE_RGB_6 = "3.299094e+00 3.848709e+00 70
10 93 70 11 94" ;<br>
Band1:COLOR_TABLE_RULE_RGB_7 = "3.848709e+00 4.398325e+00 70
11 94 71 13 96" ;<br>
<br>
</div>
<div>I need what I have below, using Climate Data Operators
(cdo)...</div>
<div><br>
</div>
<div>cdo -chname,Band1,sac_UZTWM sac_UZTWM.nc
sac_UZTWM.nc.test</div>
<div><br>
</div>
<div>I can change the variable name, which is exactly what I
need:</div>
<div><br>
</div>
<div>ncdump -h sac_UZTWM.nc<br>
netcdf sac_UZTWM {<br>
dimensions:<br>
x = 2698 ;<br>
y = 2427 ;<br>
variables:<br>
double x(x) ;<br>
x:standard_name = "projection_x_coordinate" ;<br>
x:long_name = "x coordinate of projection" ;<br>
x:units = "m" ;<br>
x:axis = "X" ;<br>
double y(y) ;<br>
y:standard_name = "projection_y_coordinate" ;<br>
y:long_name = "y coordinate of projection" ;<br>
y:units = "m" ;<br>
y:axis = "Y" ;<br>
char lambert_conformal_conic ;<br>
lambert_conformal_conic:grid_mapping_name =
"lambert_conformal_conic" ;<br>
lambert_conformal_conic:longitude_of_central_meridian = 74.
;<br>
lambert_conformal_conic:false_easting = 0. ;<br>
lambert_conformal_conic:false_northing = 0. ;<br>
lambert_conformal_conic:latitude_of_projection_origin = 23.
;<br>
lambert_conformal_conic:standard_parallel = 33., 45. ;<br>
lambert_conformal_conic:long_name = "CRS definition" ;<br>
lambert_conformal_conic:longitude_of_prime_meridian = 0. ;<br>
lambert_conformal_conic:semi_major_axis = 6378137. ;<br>
lambert_conformal_conic:inverse_flattening = 298.257223563 ;<br>
lambert_conformal_conic:crs_wkt =
"PROJCS[\"unknown\",GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"latitude_of_origin\",23],PARAMETER[\"central_meridian\",74],PARAMETER[\"standard_parallel_1\",33],PARAMETER[\"standard_parallel_2\",45],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]]"
;<br>
lambert_conformal_conic:spatial_ref =
"PROJCS[\"unknown\",GEOGCS[\"wgs84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"latitude_of_origin\",23],PARAMETER[\"central_meridian\",74],PARAMETER[\"standard_parallel_1\",33],PARAMETER[\"standard_parallel_2\",45],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]]"
;<br>
lambert_conformal_conic:GeoTransform = "-1567000 1000 0
2160000 0 -1000 " ;<br>
double sac_UZTWM(y, x) ;<br>
sac_UZTWM:long_name = "GDAL Band Number 1" ;<br>
sac_UZTWM:grid_mapping = "lambert_conformal_conic" ;<br>
sac_UZTWM:_FillValue = -999. ;<br>
sac_UZTWM:missing_value = -999. ;<br>
<br>
// global attributes:<br>
:CDI = "Climate Data Interface version 2.4.0 (<a
href="https://mpimet.mpg.de/cdi" moz-do-not-send="true"
class="moz-txt-link-freetext">https://mpimet.mpg.de/cdi</a>)"
;<br>
:Conventions = "CF-1.5" ;<br>
:GDAL = "GDAL 3.9.3, released 2024/10/07 (debug build)" ;<br>
:history = "Tue Feb 25 06:31:46 2025: cdo
-chname,Band1,sac_UZTWM sac_UZTWM.nc sac_UZTWM.nc.test\n",<br>
"Tue Feb 25 06:30:38 2025: GDAL Create(
/media/teaiii/research/pakistan/parameters/sac_UZTWM.nc, ...
)" ;<br>
:CDO = "Climate Data Operators version 2.4.0 (<a
href="https://mpimet.mpg.de/cdo" moz-do-not-send="true"
class="moz-txt-link-freetext">https://mpimet.mpg.de/cdo</a>)"
;<br>
}</div>
<div><br>
</div>
<div>But how using r.out.gdal can I do this; I tried something
that I thought should work, but I got an error</div>
<div><br>
</div>
r.out.gdal -c -m --overwrite input=uztwm@teaiii
output=/media/teaiii/research/pakistan/parameters/sac_UZTWM.nc
format=netCDF nodata=-999<br>
<br>
</div>
<div>Thank you,</div>
<div>Tom</div>
<div><br>
<div dir="ltr" class="gmail_signature"
data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr"><br>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre wrap="" class="moz-quote-pre">_______________________________________________
grass-user mailing list
<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a>
</pre>
</blockquote>
<br>
</body>
</html>