compass.landice.mesh.add_bedmachine_thk_to_ais_gridded_data

compass.landice.mesh.add_bedmachine_thk_to_ais_gridded_data(self, source_gridded_dataset, bedmachine_path)[source]

Copy BedMachine thickness to AIS reference gridded dataset. Replace thickness field in the compilation dataset with the one we will be using from BedMachine for actual thickness interpolation. There are significant inconsistencies between the masking of the two, particularly along the Antarctic Peninsula, that lead to funky mesh extent and culling if we use the thickness from 8km composite dataset to define the cullMask but then actually interpolate thickness from BedMachine. This function uses bilinear interpolation to interpolate from the 500 m resolution of BedMachine to the 8 km resolution of the reference dataset. It is not particularly accurate, but is fast and adequate for generating the flood filled mask for culling the mesh. Highly accurate conservative remapping is performed later for actually interpolating BedMachine thickness to the final MALI mesh.

Parameters:
  • source_gridded_dataset (str) – name of NetCDF file containing original AIS gridded datasets

  • bedmachine_path (str) – path to BedMachine dataset

Returns:

gridded_dataset_with_bm_thk (str) – name of NetCDF file with gridded dataset with BedMachine thk added