# Data fetched via HTTP from ['https://dataverse.geus.dk/api/access/datafile/:persistentId?persistentId=doi:10.22008/FK2/F9P03L/FYNXIL']

# Data processed using the following steps:

gdalwarp -t_srs EPSG:3413 -r bilinear {input_dir}/geothermal_heat_flow_map_55km.nc {output_dir}/warped.tif

gdalwarp -cutline {assets_dir}/greenland_rectangle.geojson -crop_to_cutline -co COMPRESS=DEFLATE {input_dir}/warped.tif {output_dir}/geothermal_heat_flow_map_55km.tif

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

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