<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">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">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>