# Data fetched via HTTP from ['https://www.ngdc.noaa.gov/geomag/data/poles/WMM2020_NP.xy']

# Data processed using the following steps:

sed "1i longitude latitude year" {input_dir}/WMM2020_NP.xy > {output_dir}/WMM2020_NP_with_header.xy

ogr2ogr -lco ENCODING=UTF-8 -t_srs EPSG:3413 -clipdst {assets_dir}/latitude_shape_40_degrees.geojson -makevalid -oo X_POSSIBLE_NAMES=longitude -oo Y_POSSIBLE_NAMES=latitude -s_srs EPSG:4326 {output_dir}/geomagnetic_north_pole.gpkg CSV:{input_dir}/WMM2020_NP_with_header.xy