# Data fetched via HTTP from ['https://arcticdata.io/metacat/d1/mn/v2/object/urn%3Auuid%3A61ff2294-4734-46ba-a0b0-845d69298131']

# Data processed using the following steps:

gdalwarp -cutline {assets_dir}/greenland_rectangle.geojson -crop_to_cutline -r bilinear -t_srs EPSG:3413 -srcnodata 0 -tr 1800 1800 -co COMPRESS=DEFLATE NETCDF:{input_dir}/percent_gris_g1800m_v3a_rcp_26_0_1000.nc:sftgif {output_dir}/extracted.tif

gdal_translate -co TILED=YES -co COMPRESS=DEFLATE -co PREDICTOR=3 {input_dir}/extracted.tif {output_dir}/compressed.tif

cp {input_dir}/compressed.tif {output_dir}/final.tif && gdaladdo -r average {output_dir}/final.tif 2 4 8 16