Create multivariable NetCDF from several Raster files

Is it possible to create a NetCDF with several variables from a collection of rasters (same spatial bounds)?

I've tried the Arcpy tool RasterToNetCDF, but it seems thats it`s only able to transform one raster each time.

Answer

I can find quite some information regarding Creating multiband netcdf with GDAL. Have a look here. If that doesn't help, have a look at gdal_merge.py, once the individual rasters are created something like this should do

gdal_merge.py image1.tif image2.tif image3.tif -o merged.nc.