Thanks Even!
The NETCDF file(s) I'm working with have hundreds of sub data sets.
Initially, I was breaking out all of them to .vrt's as you suggested,
but it got a bit messy. Counting the sub-datasets with gdalinfo
iteratively then using gdal_create -bands N was a lot cleaner.
Passing dozens
Scott,
# add another tif to band 2 and add meta data for it:
gdalwarp -q -srcband 1 -dstband 2 -t_srs EPSG:4326 src2.tif result.tif
BAND1 is lost at that stage. gdalwarp must override the metadata of the
output dataset with the one of src2.tif.
Presumably it should instead merge them, using
