This scenario is hypothetical, and by no means fully representative of all the processes involved in both the evolution of polar ice sheets, and the way sea level rises differently in various parts of the world. It does not for example account for local subsidence of the ground, flooding, erosion of the shorelines, glacial isostatic rebound, negative feedbacks from the atmosphere, etc. It is, however, meant to assess the sensitivity of local sea level rise to melting of polar regions.

• Digital elevation model for US coastal areas: Global Land One-km Base Elevation Project anisotropically interpolated onto unstructured mesh ranging from 500 m to 100 km resolution.
• US overlay: TrueMarble (4km resolution).

• Greenland mesh at 20 km resolution split into 3 areas (North, Center and South).
• Antarctica mesh at resolution ranging from 20 to 60 km, split into 4 areas (northwest, northeast, southwest and southeast).
• Computation of SLR using volume above hydrostatic equilibrium and conversion to global mean sea level rise from global ocean surface area. New US coastline computed using simple comparison between surface height and new sea level.
• Code to compute volume above floatation:

function VolumeAboveFloatation(md) { 
	let vol = 0;
	let di = md.materials.rho_ice / md.materials.rho_water;

	for (let i = 0; i < md.mesh.numberofvertices; ++i) {
		if (md.geometry.base[i] > 0){
			vol += (md.geometry.surface[i] - md.geometry.base[i]) * area[i];
		} else {
			let sh = (di - 1) / di * md.geometry.base[i];
			if (md.geometry.surface[i] > sh){
				vol += (md.geometry.surface[i] - sh) * area[i];
			}
		}
	}
	
	return vol;
}

• Code to convert to global mean sea level rise:

// Figure out the new volumes of the polar ice caps
let offset = 0.01;
let vg = VolumeAboveFloatation(mdgre);
let va = VolumeAboveFloatation(mdant);

// Overall volume change
let deltav = (vg - vgi) + (va - vai);

// Overall water equivalent volume change
deltav = deltav * di;

// equivalent slr
// 361 million square kilometers ocean surface 
// 10^6 * 10^6 (from etopo)
let newheight = offset - deltav / (361.934 * Math.pow(10,12));