We combine airborne gravity data collected in  2012 and 2016, with WISE and CECs radar sounder data from 2014 and 2015 to map bed elevation and ice thickness accross the Northern and Southern Patagonian Icefields. We perform a three dimensional inversion of the gravity field constrained by radar-derived thickness, ice boundary, and fjord bathymetry to infer bed elevation at 500-m spacing, with a precision of about 50 m. The enclosed netCDF contains the following fields:

- X: x UTM coordinate (meters)

-Y: y UTM coordinate (meters)

-bed: bedrock elevation (meters)

-misfit: Misfit after inversion (calculated minus observed gravity)

-misfitm: Converted misfit using a conversion factor of 5.8 mGal per 100m

-Source: Data source, 0=bathymetry , 1=gravity, 2=radar depth sounder

-Fulthick: Comprehensive Ice thickness calculated using gravity inversion inside the survey domain and Carrivick 2016 modified outside.

The projection is given as:

"PROJCS[\"WGS 84 / UTM zone 18S\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-75],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",10000000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"32718\"]]" ;