Set bottom + top ghost cells of ql to zero to prevent problems in
interpolation for masks
Create separate 2D mask for sampling surface quantities (e.g. projected
cloud cover). Not yet coupled to the statistics
Apply botom (surface) mask in time serie variables
Start at surface statistics. Todo: implement mask and add to NetCDF
Merge branch 'statmask' into surfacestat
Added surface time series (ustar, obl) including mask
Moved processing of time dependent boundary conditions to separate
function
Merge branch 'master' into surfacestat
Added generic script to read MicroHH statistics, and script to plot the
Gabls1 figures from Beare et al (2006)