# Data accessed manually by a human following instructions: RACMO data were obtained via a private data transer by Brice Noël. See the `scripts/preprocess-private-archive/racmo_qgreenla nd_jan2021/README.md` for additional information. # Data processed using the following steps: unzip {input_dir}/RACMO_QGreenland_Jan2021.zip -d {output_dir} subl.1958-2019.BN_RACMO2.3p2_FGRN055_1km.YY-mean.nc Icemask_Topo_Iceclasses_lon_lat_average_1km_GrIS.nc gdal_calc.py --calc="numpy.where((B != 0), A, -9999)" --NoDataValue=-9999 --outfile={output_dir}/subl.tif -A NETCDF:{input_dir}/subl.1958-2019.BN_RACMO2.3p2_FGRN055_1km.YY-mean.nc:subl -B NETCDF:{input_dir}/Icemask_Topo_Iceclasses_lon_lat_average_1km_GrIS.nc:Promicemask cp {input_dir}/subl.tif {output_dir}/edited.tif && gdal_edit.py -a_srs EPSG:3413 {output_dir}/edited.tif gdal_translate -co TILED=YES -co COMPRESS=DEFLATE -co PREDICTOR=3 {input_dir}/edited.tif {output_dir}/compressed.tif cp {input_dir}/compressed.tif {output_dir}/racmo_subl.tif && gdaladdo -r average {output_dir}/racmo_subl.tif 2 4 8 16