# 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} Icemask_Topo_Iceclasses_lon_lat_average_1km_GrIS.nc

gdal_calc.py --calc="numpy.where((B != 0), A, -9999)" --NoDataValue=-9999 --outfile={output_dir}/Topography.tif -A NETCDF:{input_dir}/Icemask_Topo_Iceclasses_lon_lat_average_1km_GrIS.nc:Topography -B NETCDF:{input_dir}/Icemask_Topo_Iceclasses_lon_lat_average_1km_GrIS.nc:Promicemask

cp {input_dir}/Topography.tif {output_dir}/edited.tif && gdal_edit.py -a_srs EPSG:3413 -a_ullr -639456.0 -655096.0 856544.0 -3355096.0 {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_Topography.tif && gdaladdo -r average {output_dir}/racmo_Topography.tif 2 4 8 16