# Data fetched via HTTP from ['https://www.ncei.noaa.gov/thredds-ocean/fileServer/woa23/DATA/temperature/netcdf/decav91C0/0.25/woa23_decav91C0_t13_04.nc']

# Data processed using the following steps:

gdal_translate -b 25 NETCDF:{input_dir}/woa23_decav91C0_t13_04.nc:t_an {output_dir}/extracted.tif

gdalwarp -t_srs EPSG:3413 -r bilinear -tr 25000 25000 -te -5774572.727594967 -5774572.727594967 5774572.727594967 5774572.727594967 {input_dir}/extracted.tif {output_dir}/warped.tif

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

gdal_translate -co TILED=YES -co COMPRESS=DEFLATE -co PREDICTOR=3 {input_dir}/warped_and_cut.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