# Data fetched via HTTP from ['ftp://ftp.ngdc.noaa.gov/geomag/wmm/wmm2020/shapefiles/2022/WMM_2022_all_shape_geographic.zip']

# Data processed using the following steps:

unzip -j -d {output_dir} {input_dir}/WMM_2022_all_shape_geographic.zip *X_SV_2022*

OGR_ENABLE_PARTIAL_REPROJECTION=TRUE ogr2ogr -lco ENCODING=UTF-8 -t_srs EPSG:3413 -clipdst {assets_dir}/latitude_shape_40_degrees.geojson -dialect sqlite -sql "Select Geometry, Contour, SIGN, \"INDEX\", CAST(Contour AS TEXT) || \" nT/y\" as label FROM X_SV_2022" {output_dir}/X_SV_2022.gpkg {input_dir}/X_SV_2022.shp