Multi-dimensional NetCDF output using `OUTFILESFLAG`

[dprincz]

As requested by @mee067 in #115

it will be good to move from flat nc output which has only one variable (per layer, or gru), to a multi-dimensional file, where all layers, or grus are combined in a single file.