# Data fetched via HTTP from ['https://arcticdata.io/metacat/d1/mn/v2/object/urn%3Auuid%3Aed2d7235-2193-4ba3-a98f-f09d871199a1']

# 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_45_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