31:T629,# EMC Horizontal Slicer 1.3
## Beta #### Quick links [→ EMC—Earth Model Collaboration](http://ds.iris.edu/ds/products/emc/)
[→ Repository Explorer—search and discover EMC models](https://observablehq.com/@earthscope/emc-repository-explorer?collection=@earthscope/models)
[→ netCDF Reader—learn how to read netCDF files in Observable](https://observablehq.com/@earthscope/emc-netcdf-reader) #### Overview The Horizontal Slicer is a flexible tool for quickly visualizing horizontal slices from the EMC-compatible models. The model could be either inputted from the user's local disk or forwarded by the Repository Explorer notebook. Optional auxiliary data such as global active faults, plate boundaries, earthquake distribution and volcano locations could also be superimposed on the model plots. This notebook currently supports 2D and 3D models under the netCDF 3 format. To use this notebook: * *Load a model from the EMC repository using the [**Explorer notebook**](https://observablehq.com/@earthscope/earth-model-collection-explorer?collection=@earthscope/models)*. **OR** * *Load an EMC-compatible model from your computer using the **Local Model** menu below.* ***Note:*** Initially, the **${defaultModel.name}** file is loaded as the **Default Model**. * *If necessary, change the X, Y, and Z parameter settings for the model.* * *Adjust the **Auxiliary Data** display parameters.* * *Tune opacity and grid sampling parameters as needed.* ***Note:*** Observable notebooks update reactively as you modify the parameters.32:T496,title = { if (zValues.length > 0) { return html`

${modelDescription}

${selectedVar} (${ selections.selectedVarInfo[1] }, ${varUnits[selectedVar]}) at the depth of ${d3.format(",")( selections.selectedZ )} km model depth range (${zValues[0].toFixed(1)} to ${zValues[ zValues.length - 1 ].toFixed(1)} km)

`; } else { return html`

${modelId}

${modelDescription}

${selectedVar} (${selections.selectedVarInfo[1]}, ${varUnits[selectedVar]}) from a 2D model

`; } }33:T609,{ switch (colormap.selectedCmap) { case "Turbo": return legend({ color: d3.scaleSequential( [dataExtent[0], dataExtent[1]], d3.interpolateTurbo ), title: selections.selectedVarInfo + " (" + varUnits[selectedVar] + ")", invert: colormap.invertCmap }); break; case "Plasma": return legend({ color: d3.scaleSequential( [dataExtent[0], dataExtent[1]], d3.interpolatePlasma ), title: selections.selectedVarInfo + " (" + varUnits[selectedVar] + ")", invert: colormap.invertCmap }); break; case "rainbow": return legend({ color: d3.scaleSequential( [dataExtent[0], dataExtent[1]], d3.interpolateRainbow ), title: selections.selectedVarInfo + " (" + varUnits[selectedVar] + ")", invert: colormap.invertCmap }); break; case "Viridis": return legend({ color: d3.scaleSequential( [dataExtent[0], dataExtent[1]], d3.interpolateViridis ), title: selections.selectedVarInfo + " (" + varUnits[selectedVar] + ")", invert: colormap.invertCmap }); break; case "Cividis": return legend({ color: d3.scaleSequential( [dataExtent[0], dataExtent[1]], d3.interpolateCividis ), title: selections.selectedVarInfo + " (" + varUnits[selectedVar] + ")", invert: colormap.invertCmap }); break; } }34:T1908,chart = { if (chartType == "Interactive Globe") { const context = DOM.context2d(width, height); context.font = "14px roboto"; const path = d3.geoPath(globeProjection, context).pointRadius(3); context.cursorStyle = "crosshair"; function render({ lofi = false }) { context.clearRect(0, 0, width, height); context.globalAlpha = 1; context.beginPath(); path(sphere); context.lineWidth = 1.5; context.fillStyle = "#cfe0ff"; context.strokeStyle = "#000"; context.fill(); context.stroke(); context.beginPath(); path(d3.geoGraticule10()); context.lineWidth = 0.5; context.strokeStyle = "#aaa"; context.stroke(); // land context.beginPath(); path(land); context.fillStyle = "#faf8f5"; context.fill(); context.stroke(); // countries context.beginPath(); path(countriesGeometry); context.fillStyle = "#faf8f5"; context.fill(); context.stroke(); if (!lofi) { context.globalAlpha = sliceOpacity; path.pointRadius(3); // Sample sample.forEach((d) => { let feature = { type: "Point", coordinates: [+d.longitude, +d.latitude] }; context.beginPath(); path(feature); context.fillStyle = d[selections.selectedVarInfo[0]] ? colorScale(d[selections.selectedVarInfo[0]]) : "rgba(0, 0, 200, 0)"; // don't draw NaN values context.fill(); }); context.globalAlpha = misc.auxDataOpacity; path.pointRadius(2); if (auxData.includes("earthquakes")) { // Earthquakes context.beginPath(); path({ type: "MultiPoint", coordinates: earthquakeSamples.map((d) => [ +d.Longitude, +d.Latitude ]) }); context.fillStyle = auxScale("earthquakes"); context.fill(); } if (auxData.includes("volcanoes")) { // Volcanoes context.beginPath(); path({ type: "MultiPoint", coordinates: volcanoSamples.map((d) => [+d.long, +d.lat]) }); context.fillStyle = auxScale("volcanoes"); context.fill(); context.strokeStyle = "black"; context.stroke(); } if (auxData.includes("active fault")) { // Active Faults context.beginPath(); path({ type: "MultiLineString", coordinates: activeFaults.features.map( (d) => d.geometry.coordinates ) }); context.strokeStyle = auxScale("active fault"); context.lineWidth = 1; context.stroke(); } if (auxData.includes("divergent")) { // Divergent context.beginPath(); path({ type: "MultiPoint", coordinates: divergentSamples.map((d) => [ +d.longitude, +d.latitude ]) }); context.strokeStyle = auxScale("divergent"); context.lineWidth = 0.5; context.stroke(); } if (auxData.includes("convergent")) { // Convergent context.beginPath(); path({ type: "MultiPoint", coordinates: convergentSamples.map((d) => [ +d.longitude, +d.latitude ]) }); context.strokeStyle = auxScale("convergent"); context.lineWidth = 0.5; context.stroke(); } if (auxData.includes("transform")) { // Transform context.beginPath(); path({ type: "MultiPoint", coordinates: transformSamples.map((d) => [ +d.longitude, +d.latitude ]) }); context.strokeStyle = auxScale("transform"); context.lineWidth = 0.5; context.stroke(); } context.globalAlpha = 1; } } let T; function zoomed(e) { T = e.transform; mutable currentZoom = T; render({ lofi: true }); // Render fewer details while user moves map. } context.canvas.onmousemove = (event) => { if (T) { let coordinate = globeProjection.invert(d3.pointer(event)); mutable mouse = d3.pointer(event); } render({ lofi: false }); }; return d3 .select(context.canvas) .call( zoom(globeProjection) .on("zoom.render", zoomed) .on("end.render", () => render({ lofi: false })) ) .call(render) .node(); } else { return Plot.plot({ projection: { type: "mercator", domain: { type: "MultiPoint", coordinates: sample.map((d) => [d.longitude, d.latitude]) } }, color: { type: "linear", scheme: colormap.selectedCmap.toLowerCase(), domain: [dataExtent[0], dataExtent[1]], range: [0, 1], reverse: colormap.invertCmap }, width: height, marks: [ Plot.geo(land, { clip: true }), Plot.geo(countriesGeometry, { clip: true }), Plot.geo(states, { clip: true }), Plot.raster( sample.filter((d) => !isNaN(d[selections.selectedVarInfo[0]])), { clip: true, x: (d) => d.longitude, y: (d) => d.latitude, fill: (d) => d[selections.selectedVarInfo[0]], opacity: sliceOpacity, interpolate: "barycentric", blur: 2 } ), Plot.geo(activeFaults, { stroke: auxScale("active fault") }), Plot.dot( volcanoSamples.map((d) => [+d.long, +d.lat], { fill: auxScale("volcanoes") }) ), Plot.dot( divergentSamples.map((d) => [+d.longitude, +d.latitude]), { fill: auxScale("divergent") } ), Plot.dot( earthquakeSamples.map((d) => [+d.Longitude, +d.Latitude]), { fill: auxScale("earthquakes"), r: 2 } ), Plot.dot( convergentSamples.map((d) => [+d.longitude, +d.latitude]), { fill: auxScale("convergent") } ) ] }); } }35:T51e,earthquakePlot = { let circleOpacity = 0; if (eqDisplay.includes("markers")) { circleOpacity = 0.3; } let contourOpacity = 0; if (eqDisplay.includes("contours")) { contourOpacity = 0.9; } if (auxData.includes("earthquakes")) { return Plot.plot({ x: { tickFormat: "0.1f", tickRotate: 45, label: "magnitude", domain: [magValues[0], magValues[magValues.length - 1]] }, inset: 50, width, y: { label: "Depth (km)", domain: [depthQuantile, 0], tickFormat: "0.0f" }, color: { scheme: "ylgnbu" }, grid: true, marks: [ Plot.dot( earthquakeSamplesQuantile, Plot.binY( { r: "count" }, { x: (d) => +d["Magnitude"], y: (d) => +d["Depth/km"], opacity: circleOpacity, fill: "maroon" } ) ), Plot.density( earthquakeSamplesQuantile, Plot.binY( { r: "count" }, { x: (d) => +d["Magnitude"], y: (d) => +d["Depth/km"], opacity: contourOpacity, bandwidth } ) ) ] }); } else { return html``; } }36:T787,PARAMS = ({ // find them when we want to get them directly from the model file. // These are the parameters that mostly come from the model metadata. I am defining them here so we can eqNumberMax: 2000, eqNumberDefault: 1000, eqMagMax: 12, eqMagMinDefault: 5, projections: [ ["equal-earth", "equal-earth"], ["albers", "albers"], ["albers-usa", "albers-usa"], ["mercator", "mercator"], ["equirectangular", "equirectangular"], ["azimuthal-equidistant", "azimuthal-equidistant"], ["transverse-mercator", "transverse-mercator"], ["gnomonic", "gnomonic"], ["azimuthal-equal-area", "azimuthal-equal-area"], ["stereographic", "stereographic"], ["orthographic", "orthographic"], ["equirectangular", "equirectangular"] ], projectionDefault: { value: "equal-earth", label: "Projection" }, invertCmapDefault: { value: true, label: "Invert Colormap" }, gridResampleDefault: { value: 1, label: "Resample grids by", step: 1 }, gridCountDefault: 200000, colorMaps: [ ["Turbo, Turbo, An Improved Rainbow", "Turbo"], ["rainbow", "rainbow"], ["Viridis", "Viridis"], ["Cividis", "Cividis"], ["Plasma", "Plasma"] ], emcRepositoryURL: "https://ds.iris.edu/files/products/emc/emc-files/", colorMapDefault: { value: 0.2, label: "Colormap" }, sliceOpacityDefault: { value: 0.15, label: "Slice Opacity", step: 0.05 }, mapOpacityDefault: { value: 0.6, label: "Basemap Opacity", step: 0.05 }, mapMarkerRadiusDefault: { value: 5, label: "Map Marker Radius", step: 1 }, auxData: [ ["Earthquakes¹", "earthquakes"], ["Volcanoes²", "volcanoes"], ["Divergent P.B.³", "divergent"], ["Convergent P.B.³", "convergent"], ["Transform P.B.³", "transform"], ["Active Fault⁴", "active fault"] ], auxDataValue: [ "volcanoes", "earthquakes", "divergent", "convergent", "transform", "active fault" ] })37:[{"id":"9c175be6367ed138205fe43bc00f32dc610b9e48dc75e8e5993876855978970e4e831968df1f28174f3969bbd55ba01e3b98f7a613b80b066389445f8146a473","href":"https://static.observableusercontent.com/files/9c175be6367ed138205fe43bc00f32dc610b9e48dc75e8e5993876855978970e4e831968df1f28174f3969bbd55ba01e3b98f7a613b80b066389445f8146a473","name":"transform.csv","type":"text/csv","lastModified":1672161109043.097,"size":132706},{"id":"d39298ea9b7ff20be4c0cd05c4f3f9c8ce0e8065613a38fb69436a43cb265d26c83d38b2bcf239f969445acb55a84dc0a2501eb8544d4a5e18f56eebaaf5c20e","href":"https://static.observableusercontent.com/files/d39298ea9b7ff20be4c0cd05c4f3f9c8ce0e8065613a38fb69436a43cb265d26c83d38b2bcf239f969445acb55a84dc0a2501eb8544d4a5e18f56eebaaf5c20e","name":"convergent.csv","type":"text/csv","lastModified":1672161109043.097,"size":98750},{"id":"f022ee10b969e07156b5fc54f2941a409e62030ce814b4a0a9ff4c6b14d6abffd9ddfcc6cf7bb09fd2f9db3d09ccf1015c872c784d7813cd677c3909a08bc912","href":"https://static.observableusercontent.com/files/f022ee10b969e07156b5fc54f2941a409e62030ce814b4a0a9ff4c6b14d6abffd9ddfcc6cf7bb09fd2f9db3d09ccf1015c872c784d7813cd677c3909a08bc912","name":"divergent.csv","type":"text/csv","lastModified":1672161109043.097,"size":102155},{"id":"5620cd7159029727afd64f5b6739bc4d5de7b8c2292883d475640eed4efe7ca2c1077eb5ca598c92751ee35d7e427e76ec706f50339facd70601bc69e06ca79b","href":"https://static.observableusercontent.com/files/5620cd7159029727afd64f5b6739bc4d5de7b8c2292883d475640eed4efe7ca2c1077eb5ca598c92751ee35d7e427e76ec706f50339facd70601bc69e06ca79b","name":"WOVOdat_volcano_locations.csv","type":"text/csv","lastModified":1672161109043.097,"size":99904},{"id":"f4433a4508d35df7867681877371a172dcba6ea982f56f8309822fe4bed13a50907447ba71c596e041c76269e4934937cd878236325589f96e2cff3245f76e8d","href":"https://static.observableusercontent.com/files/f4433a4508d35df7867681877371a172dcba6ea982f56f8309822fe4bed13a50907447ba71c596e041c76269e4934937cd878236325589f96e2cff3245f76e8d","name":"plot.umd.min (2).js","type":"application/javascript","lastModified":1672161109043.097,"size":137943},{"id":"5da3f659e07ce6c96d368d02f1eb49ee05f10740530e1bdffe56cd4d462260377608071694661a83b8937e6a74c37f2d804ac2e90cceb209ef8a83ce899538d5","href":"https://static.observableusercontent.com/files/5da3f659e07ce6c96d368d02f1eb49ee05f10740530e1bdffe56cd4d462260377608071694661a83b8937e6a74c37f2d804ac2e90cceb209ef8a83ce899538d5","name":"CUSRA2021.r0.0.nc","type":"application/octet-stream","lastModified":1673476172786.462,"size":2077612},{"id":"f75ca3dc7c0b65cf225cea300e01e5e3cb5abf4ad75592936a2b6c79b797e933a208355d31d5b160f5b1db2a7de61fa402fe279d036a052211cd09462f524cad","href":"https://static.observableusercontent.com/files/f75ca3dc7c0b65cf225cea300e01e5e3cb5abf4ad75592936a2b6c79b797e933a208355d31d5b160f5b1db2a7de61fa402fe279d036a052211cd09462f524cad","name":"land-110m.json","type":"application/json","lastModified":1674240166862.074,"size":55207},{"id":"5b7bb44872469fc6f9f4256d902fba909728666ad0d5263654ce95f4470ec90527fdc39a1d1b70490f457d65c31795d651a4afc0745f72351916e715fddc7352","href":"https://static.observableusercontent.com/files/5b7bb44872469fc6f9f4256d902fba909728666ad0d5263654ce95f4470ec90527fdc39a1d1b70490f457d65c31795d651a4afc0745f72351916e715fddc7352","name":"countries-110m@1.json","type":"application/json","lastModified":1674240407036.148,"size":107761},{"id":"4f4339611a8deff730136068de8f6cb87f342b1183691c8618ac209804e5848159d5d48d523bed28146a075667326e125608638bec2b099b6693610f91c0a9f6","href":"https://static.observableusercontent.com/files/4f4339611a8deff730136068de8f6cb87f342b1183691c8618ac209804e5848159d5d48d523bed28146a075667326e125608638bec2b099b6693610f91c0a9f6","name":"gem_active_faults.json","type":"application/json","lastModified":1674505969141.114,"size":12302802},{"id":"4918abcd2b8f4440c28b9e3f88330557d07ccfda7a1177633c154fcc137f7303b181b7b32f2e6711b03337ffacda82df428be3bf7401e43729fa94c90228ad33","href":"https://static.observableusercontent.com/files/4918abcd2b8f4440c28b9e3f88330557d07ccfda7a1177633c154fcc137f7303b181b7b32f2e6711b03337ffacda82df428be3bf7401e43729fa94c90228ad33","name":"states-10m.json","type":"application/json","lastModified":1674662103299.225,"size":114554},{"id":"5874aaadeac5469822c1486e416d406e06cb993f08b62e999dc7362ec30ae2ea08b39bf8b7b64a873d7e56200b796b255cc2fe3d303c0ebe3f3e63d984be305a","href":"https://static.observableusercontent.com/files/5874aaadeac5469822c1486e416d406e06cb993f08b62e999dc7362ec30ae2ea08b39bf8b7b64a873d7e56200b796b255cc2fe3d303c0ebe3f3e63d984be305a","name":"wUS-SH-2010_percent.nc","type":"application/octet-stream","lastModified":1678685114158.673,"size":1728968}] 2f:{"title":"EMC Horizontal Slicer 1.3","cells":[{"id":0,"value":"$31","pinned":false,"hidden":false,"mode":"md"},{"id":4880,"value":"#### Settings","pinned":false,"hidden":false,"mode":"md"},{"id":2710,"value":"viewof model = Inputs.form({\n defaultModel_file: Inputs.text({\n label: \"Default Model\",\n value: defaultModel.name,\n disabled: true\n }),\n localModel: Inputs.file({\n label: \"Local Model\",\n required: false,\n value: \"\",\n reset: true\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":371,"value":"**${modelId}** (${Object.keys(dimensions).length}D): ${modelDescription}\n(**read more about this model**)
","pinned":false,"hidden":false,"mode":"md"},{"id":1658,"value":"viewof coordinates = Inputs.form({\n xVariable: Inputs.select(modelDimensionVars, {\n label: \"X-Variable\",\n value: \"longitude\"\n }),\n yVariable: Inputs.select(modelDimensionVars, {\n label: \"Y-Variable\",\n value: \"latitude\"\n }),\n zVariable: Inputs.select(modelDimensionVars, depthOption)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":705,"value":"viewof colormap = Inputs.form({\n selectedCmap: Inputs.select(\n new Map(PARAMS.colorMaps),\n PARAMS.colorMapDefault\n ),\n invertCmap: Inputs.toggle(PARAMS.invertCmapDefault)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":971,"value":"viewof misc = Inputs.form({\n minMagnitude: Inputs.range([0, PARAMS.eqMagMax], {\n label: \"Min. Magnitude\",\n step: 0.1,\n value: PARAMS.eqMagMinDefault\n }),\n maxEq: Inputs.range([0, PARAMS.eqNumberMax], {\n value: PARAMS.eqNumberDefault,\n step: 10,\n label: \"Max. Earthquakes\"\n }),\n mapOpacity: Inputs.range([0, 1], {\n value: (PARAMS.mapOpacityDefault.value * widthFactor).toFixed(2),\n label: PARAMS.mapOpacityDefault.label,\n step: PARAMS.mapOpacityDefault.step\n }),\n auxDataOpacity: Inputs.range([0, 1], {\n label: \"Aux. Data Opacity\",\n step: 0.01,\n value: 1\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":3844,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":817,"value":"viewof auxData = Inputs.checkbox(new Map(PARAMS.auxData), {\n value: PARAMS.auxDataValue.filter(\n (d) => d !== \"volcanoes\" && d !== \"earthquakes\"\n ),\n label: \"Auxiliary Data\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":861,"value":"

${auxData.map((d) => htl.html`

${d.toUpperCase()}

`)}

","pinned":false,"hidden":false,"mode":"html"},{"id":3345,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":3100,"value":"*Click, drag, and scroll to interact with the globe.*","pinned":false,"hidden":false,"mode":"md"},{"id":4622,"value":"*Open the image cell menu on the left to download image as a PNG.*","pinned":false,"hidden":false,"mode":"md"},{"id":1145,"value":"viewof selections = Inputs.form({\n selectedVarInfo: Inputs.select(varInfoItems, {\n label: \"Variable\",\n step: 1\n }),\n // selectedZ: Inputs.select(zValues, {\n // label: coordinates.zVariable + \" (\" + varUnits[coordinates.zVariable] + \")\",\n // step: 1,\n // value: zValues[1]\n // })\n selectedZ: Inputs.range([0, zValues.length - 1], {\n label: coordinates.zVariable + \" (\" + varUnits[coordinates.zVariable] + \")\",\n step: 1,\n format: (d) => zValues[d],\n value: 0\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4835,"value":"## View Options:","pinned":false,"hidden":false,"mode":"md"},{"id":4820,"value":"viewof chartType = Inputs.radio([\"Interactive Globe\", \"Rasterized Model\"], {\n label: \"Select viewer\",\n value: \"Interactive Globe\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2979,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":3153,"value":"viewof sliceOpacity = Inputs.range([0, 1.0], {\n value: (PARAMS.sliceOpacityDefault.value * widthFactor).toFixed(2),\n label: PARAMS.sliceOpacityDefault.label,\n step: PARAMS.sliceOpacityDefault.step\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":3879,"value":"viewof grids = Inputs.form({\n stats: Inputs.text({\n label: \"\",\n value:\n \"✅ \" +\n d3.format(\",\")(data.length) +\n \" grid points loaded across \" +\n d3.format(\",\")(zValues.length) +\n \" depths\",\n disabled: true\n }),\n decimation: Inputs.range([1, 15], {\n value: 1, // Math.ceil(data.length / PARAMS.gridCountDefault),\n label: PARAMS.gridResampleDefault.label,\n step: PARAMS.gridResampleDefault.step\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4814,"value":"The plot below is a rasterized interpolation of the selected model value across **${d3.format(\",\")(data.length)}** grid points.","pinned":false,"hidden":false,"mode":"md"},{"id":4618,"value":"*Open the image cell menu on the left to download plots as SVG or PNG.*","pinned":false,"hidden":false,"mode":"md"},{"id":4083,"value":"distributionOfEarthquakes = {\n if (auxData.includes(\"earthquakes\")) {\n return md`**Distribution of earthquakes** in the area as a function of magnitude and depth. The marker size is proportional to the number of event at that depth with the corresponding magnitude. Adjust the sliders to include and exclude *outliers* and scale your view.`;\n } else {\n return html`No earthquake data requested!`;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":4287,"value":"viewof quantile = Inputs.range([0.01, 1], {\n label: \"Depth Quantile\",\n step: 0.01,\n value: 0.99\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4275,"value":"viewof bandwidth = Inputs.range([0, 40], {\n label: \"Bandwidth\",\n value: 20,\n step: 0.2\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4385,"value":"viewof eqDisplay = {\n let value = [];\n if (earthquakes.length >= 100) {\n value = [\"markers\", \"contours\"];\n } else if (earthquakes.length >= 1) {\n value = [\"markers\"];\n }\n return Inputs.checkbox([\"markers\", \"contours\"], {\n value: value,\n label: \"Display type\"\n });\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3987,"value":"$35","pinned":false,"hidden":false,"mode":"ojs"},{"id":640,"value":"{\n if (Object.keys(dimensions).length == 3) {\n return md`**Distribution of the ${selectedVar} values** across the ${d3.format(\n \",\"\n )(selections.selectedZ)} ${varUnits[coordinates.zVariable]} ${\n coordinates.zVariable\n } slice:`;\n } else {\n return md`**Distribution of the ${selectedVar}** values:`;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":631,"value":"Plot.plot({\n width: width,\n height: 200,\n color: {\n domain: dataExtent,\n scheme: colormap.selectedCmap,\n legend: false,\n reverse: colormap.invertCmap\n },\n marks: [\n Plot.rectY(\n sample,\n Plot.binX(\n { y: \"count\", label: \"Count\" },\n {\n x: selectedVar,\n stroke: plotData,\n fill: plotData,\n label: `${selectedVar} (${varUnits[selectedVar]})`\n }\n )\n ),\n Plot.ruleY([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4583,"value":"## References","pinned":false,"hidden":false,"mode":"md"},{"id":4588,"value":"* ¹ [USGS Earthquake Hazards Program - Earthquake Catalog](https://earthquake.usgs.gov/fdsnws/event/1/)\n* ² [The World Organization of Volcano Observatories, WOVOdat](https://www.wovodat.org/populate/convertie/Volcano_zone/main.php?data_type=zone_index)\n* ³ [The University of Texas Institute for Geophysics (UTIG) PLATES program](https://ig.utexas.edu/marine-and-tectonics/plates-project/)\n* ⁴ [GEM Global Active Faults Database](https://github.com/GEMScienceTools/gem-global-active-faults)","pinned":false,"hidden":false,"mode":"md"},{"id":62,"value":"---\n## Appendix","pinned":false,"hidden":false,"mode":"md"},{"id":4614,"value":"## Checks","pinned":false,"hidden":false,"mode":"md"},{"id":4613,"value":"yExtentCheck = {\n let range = d3.extent(yValues);\n let err = false;\n let message = \" \";\n if (range[0] < -90) {\n message = \"Invalid minimum: \" + range[0];\n err = true;\n }\n if (range[1] > 90) {\n message = message + \" Invalid maximum: \" + range[1];\n err = true;\n }\n if (err) {\n return htl.html`

${coordinates.yVariable} extent: ${message}

`;\n } else {\n return htl.html`

${coordinates.yVariable} extent valid

`;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":4612,"value":"xExtentCheck = {\n let range = d3.extent(xValues);\n let err = false;\n let message = \" \";\n if (range[0] < -180) {\n message = \" Invalid minimum: \" + range[0];\n err = true;\n }\n if (range[1] > 180) {\n message = message + \" Invalid maximum: \" + range[1];\n err = true;\n }\n if (err) {\n return htl.html`

${coordinates.xVariable} extent: ${message}

`;\n } else {\n return htl.html`

${coordinates.xVariable} extent valid

`;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":4611,"value":"## Definitions","pinned":false,"hidden":false,"mode":"md"},{"id":4684,"value":"sphere = ({\n type: \"Sphere\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":4670,"value":"mutable currentZoom = null","pinned":false,"hidden":false,"mode":"ojs"},{"id":2989,"value":"height = Math.max(640, width / 2)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3601,"value":"widthFactor = {\n if (width >= 1200) {\n return 5;\n } else if (width >= 900) {\n return 4;\n } else if (width >= 700) {\n return 3;\n } else if (width >= 500) {\n return 2;\n } else {\n return 1;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3843,"value":"depthOption = {\n if (Object.keys(dimensions).length == 3) {\n return { label: \"Z-Variable\", value: \"depth\" };\n } else {\n return { label: \"Z-Variable\", value: \"N/A\" };\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3011,"value":"mutable mouse = []","pinned":false,"hidden":false,"mode":"ojs"},{"id":2324,"value":"defaultModel = FileAttachment(\"wUS-SH-2010_percent.nc\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":1489,"value":"modelId = {\n if (\"model\" in metadata) {\n return metadata.model;\n } else if (\"id\" in metadata) {\n return metadata.id;\n }\n return inputModel.name;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3799,"value":"url = modelURL(reader)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1498,"value":"modelDescription = {\n if (metadata.title) {\n return metadata.title;\n }\n return modelId;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2307,"value":"sample = {\n var arr = [];\n\n for (var i = 0; i < plotData.length; i++) {\n const dict = {};\n dict[coordinates.yVariable] = plotY[i];\n dict[coordinates.xVariable] = plotX[i];\n dict[selectedVar] = plotData[i];\n arr[i] = dict;\n }\n return arr;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1258,"value":"dimensions = getDimensions(reader)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1270,"value":"modelDimensionVars = {\n let dim = Object.keys(dimensions);\n if (dim.length == 3) {\n return dim;\n } else if (dim.length == 2) {\n dim.push(\"N/A\");\n return dim;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1274,"value":"modelVariables = {\n return getVariables(reader).model_variables;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1276,"value":"variables = getVariables(reader)","pinned":true,"hidden":true,"mode":"ojs"},{"id":1460,"value":"varInfoItems = Object.entries(varLongName)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3266,"value":"modelCenter = {\n let dict = {};\n dict[coordinates.xVariable] = Math.round(\n xValues.reduce((total, current) => total + current) / xValues.length\n );\n dict[coordinates.yVariable] = Math.round(\n yValues.reduce((total, current) => total + current) / yValues.length\n );\n // Center the glibal models on the US.\n if (dict[coordinates.xVariable] == 0 && dict[coordinates.yVariable] == 0) {\n dict[coordinates.xVariable] = -110;\n dict[coordinates.yVariable] = 40;\n }\n return dict;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1630,"value":"selectedVar = selections.selectedVarInfo[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":1303,"value":"varAttributes = {\n return variables.model_variables;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2388,"value":"varMissingValue = {\n if (modelVariables[selectedVar].missing_value) {\n return modelVariables[selectedVar].missing_value;\n }\n return NaN;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2149,"value":"selectedVarDimensions = [\n modelDimensionVars[selectedVarAttr.dimensions[0]],\n modelDimensionVars[selectedVarAttr.dimensions[1]],\n modelDimensionVars[selectedVarAttr.dimensions[2]]\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":2164,"value":"{\n return d3.zip(data).map((d) => ({\n longitude: xValues[0] > 180 ? -360 + xValues[0] : xValues[0],\n latitude: yValues[0],\n var: d[0]\n }));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1619,"value":"selectedVarAttr = selectDictElement(reader.variables, selectedVar)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1410,"value":"varLongName = {\n let dict = {};\n for (const [key, value] of Object.entries(variables.model_variables)) {\n dict[key] = value.long_name;\n }\n return dict;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1339,"value":"varUnits = {\n let dict = {};\n for (const [key, value] of Object.entries(variables.model_variables)) {\n dict[key] = value.units;\n }\n\n // TODO (brooks)\n dict[\"depth\"] = \"km\";\n\n return dict;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3251,"value":"rawData = readRawData(reader)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1510,"value":"data = {\n // var data = extractParameterData(reader, selectedVar, zVariable, selectedZ);\n //.filter(function (el) {\n //return el[zVariable] == selectedZ;\n //});\n var data = rawData[selectedVar].values;\n //reader.getDataVariable(selectedVar));\n if (varMissingValue != NaN) {\n for (let i = 0; i < data.length; i++) {\n if (data[i] == \"\" || data[i] == varMissingValue) {\n data[i] = NaN;\n }\n }\n }\n return data;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2193,"value":"plotData = {\n let dataSegment = data.slice(dataRange[\"start\"], dataRange[\"end\"]);\n dataSegment = getEveryNth(dataSegment, grids.decimation);\n return dataSegment;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1761,"value":"dataExtent = d3.extent(plotData)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2173,"value":"dataRange = {\n const dict = {};\n const startIndex = zIndex * xValues.length * yValues.length;\n const endIndex = startIndex + xValues.length * yValues.length - 1;\n dict[\"start\"] = startIndex;\n dict[\"end\"] = endIndex;\n dict[\"length\"] = xValues.length * yValues.length;\n return dict;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2279,"value":"d3.extent(plotData)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2202,"value":"xIndex = {\n const key = Object.keys(selectedVarDimensions).find(\n (key) => selectedVarDimensions[key] === coordinates.xVariable\n );\n return eval(key);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1747,"value":"xAttr = selectDictElement(reader.variables, coordinates.xVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1543,"value":"xValues = {\n return rawData[coordinates.xVariable].values;\n}\n//reader.getDataVariable(xVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2232,"value":"plotX = {\n let x;\n let out = [];\n if (xIndex > yIndex) {\n x = repeat(xValues, yValues.length);\n\n x = getEveryNth(x, grids.decimation);\n return x;\n } else {\n for (const x of xValues) {\n const t = repeat([x], yValues.length);\n out.push(t);\n }\n let merge = out.flat(1);\n merge = getEveryNth(merge, grids.decimation);\n return merge;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2281,"value":"d3.extent(plotX)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1764,"value":"xExtent = {\n let range = d3.extent(xValues);\n if (range[0] < -180) {\n range[0] = -180;\n }\n if (range[1] > 180) {\n range[1] = 180;\n }\n return range;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2225,"value":"yIndex = {\n const key = Object.keys(selectedVarDimensions).find(\n (key) => selectedVarDimensions[key] === coordinates.yVariable\n );\n return eval(key);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1751,"value":"yAttr = selectDictElement(reader.variables, coordinates.yVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1541,"value":"yValues = {\n return rawData[coordinates.yVariable].values;\n}\n//reader.getDataVariable(yVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1551,"value":"yExtent = {\n let range = d3.extent(yValues);\n if (range[0] < -90) {\n range[0] = -90;\n }\n if (range[1] > 90) {\n range[1] = 90;\n }\n return range;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2269,"value":"plotY = {\n let y;\n let out = [];\n if (xIndex > yIndex) {\n for (const y of yValues) {\n const t = repeat([y], xValues.length);\n out.push(t);\n }\n let merge = out.flat(1);\n merge = getEveryNth(merge, grids.decimation);\n return merge;\n } else {\n y = repeat(yValues, xValues.length);\n\n y = getEveryNth(y, grids.decimation);\n return y;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2283,"value":"d3.extent(plotY)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1755,"value":"zAttr = selectDictElement(reader.variables, coordinates.zVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1546,"value":"zValues = {\n if (coordinates.zVariable == \"N/A\") {\n return [];\n } else {\n return rawData[coordinates.zVariable].values;\n }\n}\n//reader.getDataVariable(zVariable)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1769,"value":"zExtent = d3.extent(zValues)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2169,"value":"zIndex = selections.selectedZ","pinned":false,"hidden":false,"mode":"ojs"},{"id":4031,"value":"eqDepthValues = {\n let dict = {};\n earthquakeSamples.map((d) => (dict[+d[\"Depth/km\"]] = true));\n return Object.keys(dict)\n .map((d) => +d)\n .sort(d3.ascending);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1052,"value":"magValues = {\n let dict = {};\n earthquakeSamples.map((d) => (dict[+d.Magnitude] = true));\n return Object.keys(dict)\n .map((d) => +d)\n .sort(d3.ascending);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1104,"value":"// https://observablehq.com/@fil/netcdf\nnetcdf = import(\"https://cdn.skypack.dev/netcdfjs@2.0.1?min\").then(\n (d) => d.NetCDFReader\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3690,"value":"shareableSuffix = {\n let url = \"?\"; // Everything before the query\n\n let strippedParams = {};\n\n for (const [key, value] of Object.entries(exportConfiguration)) {\n if (value) {\n strippedParams[key] = value;\n }\n }\n\n return url + new URLSearchParams(strippedParams).toString();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3691,"value":"baseUrl = document.baseURI.split(\"?\")[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":3689,"value":"shareableUrl = baseUrl.split(\"?\")[0] + shareableSuffix","pinned":false,"hidden":false,"mode":"ojs"},{"id":3830,"value":"location","pinned":false,"hidden":false,"mode":"ojs"},{"id":3688,"value":"// Every property has a string as its value.\nexportConfiguration = {\n let obj = {\n model: metadata.id\n };\n\n for (const property in obj) {\n obj[property] = JSON.stringify(obj[property]);\n }\n\n return obj;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3687,"value":"mutable importConfiguration = {\n let dict = {};\n let i = 0;\n\n // Hack: Fixes apparent bug with Observable passing location.search down from parent frame.\n while (!location.search && i < 2) {\n i++;\n yield Promises.tick(1000, {});\n }\n\n if (location.search) {\n let queries = location.search.length\n ? location.search\n .substring(1)\n .split(\"&\")\n .map((d) => d.split(\"=\"))\n : [];\n\n queries.forEach((d) => {\n if (d[1]) {\n try {\n dict[d[0]] = JSON.parse(decodeURIComponent(d[1]));\n } catch (error) {\n console.log(\"Unable to process query parameter: \", d[0], d[1]);\n }\n }\n });\n }\n\n yield dict;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":4394,"value":"## Functions","pinned":false,"hidden":false,"mode":"md"},{"id":4676,"value":"// default project that should remain when user changes values that re-evaluate the chart cell\n// that renders the globe.\nmutable globeProjection = {\n let projection = d3.geoOrthographic().fitExtent(\n [\n [1, 1],\n [width - 1, height - 1]\n ],\n sphere\n );\n\n // Center the projection on the model bounding box.\n projection.rotate([\n -modelCenter[coordinates.xVariable],\n -modelCenter[coordinates.yVariable]\n ]);\n\n return projection;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3069,"value":"function colorScale(value) {\n let domainStart = dataExtent[0];\n let domainEnd = dataExtent[1];\n if (colormap.invertCmap) {\n domainStart = dataExtent[1];\n domainEnd = dataExtent[0];\n }\n switch (colormap.selectedCmap) {\n case \"Turbo\":\n return d3.interpolateTurbo(\n d3.scaleLinear().domain([domainStart, domainEnd]).range([0, 1])(value)\n );\n break;\n case \"Plasma\":\n return d3.interpolatePlasma(\n d3.scaleLinear().domain([domainStart, domainEnd]).range([0, 1])(value)\n );\n break;\n case \"rainbow\":\n return d3.interpolateRainbow(\n d3.scaleLinear().domain([domainStart, domainEnd]).range([0, 1])(value)\n );\n break;\n case \"Viridis\":\n return d3.interpolateViridis(\n d3.scaleLinear().domain([domainStart, domainEnd]).range([0, 1])(value)\n );\n break;\n case \"Cividis\":\n return d3.interpolateCividis(\n d3.scaleLinear().domain([domainStart, domainEnd]).range([0, 1])(value)\n );\n break;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":988,"value":"quakeRadius = {\n const scale = d3\n .scaleSqrt()\n .domain([0, 100])\n .range([magValues[0], magValues[magValues.length - 1]]);\n return function (quake) {\n return scale(Math.exp(quake[\"Magnitude\"]));\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1608,"value":"// Return a dictionary from a list of dictionaries for a given name.\nfunction selectDictElement(dictList, name) {\n const selected = dictList.find((element) => element.name === name);\n return selected;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2114,"value":"auxScale = d3.scaleOrdinal(\n [\n \"earthquakes\",\n \"volcanoes\",\n \"convergent\",\n \"divergent\",\n \"transform\",\n \"active fault\"\n ],\n [\"red\", \"orange\", \"magenta\", \"green\", \"blue\", \"black\"]\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2540,"value":"// https://bobbyhadz.com/blog/javascript-get-every-nth-element-of-array\nfunction getEveryNth(arr, nth) {\n const result = [];\n\n for (let i = 0; i < arr.length; i += nth) {\n result.push(arr[i]);\n }\n\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2244,"value":"function repeat(arr, n) {\n var repeated = new Array(n).fill(arr).flat();\n return repeated;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2922,"value":"## Data and Parameters","pinned":false,"hidden":false,"mode":"md"},{"id":3656,"value":"### The notebook parameters and defaults","pinned":false,"hidden":false,"mode":"md"},{"id":905,"value":"$36","pinned":false,"hidden":false,"mode":"ojs"},{"id":3676,"value":"### Model, data, and metadata","pinned":false,"hidden":false,"mode":"md"},{"id":1296,"value":"reader = {\n return new netcdf(await inputModel.arrayBuffer());\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3722,"value":"importConfiguration.model","pinned":true,"hidden":false,"mode":"ojs"},{"id":3685,"value":"PARAMS.emcRepositoryURL.concat(importConfiguration.model)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3183,"value":"inputModel = {\n if (model.localModel != null) {\n return model.localModel;\n } else if (\"model\" in importConfiguration) {\n if (importConfiguration.model.length > 3) {\n return fetch(PARAMS.emcRepositoryURL.concat(importConfiguration.model));\n }\n } else {\n return defaultModel;\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3177,"value":"metadata = getMetadata(reader)","pinned":false,"hidden":false,"mode":"ojs"},{"id":3648,"value":"### Other Data:","pinned":false,"hidden":false,"mode":"md"},{"id":3652,"value":"#### - Point data","pinned":false,"hidden":false,"mode":"md"},{"id":4309,"value":"depthQuantile = d3.quantile(\n earthquakeSamples.map((d) => d[\"Depth/km\"]),\n quantile\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":4295,"value":"earthquakeSamplesQuantile = {\n return earthquakeSamples.filter((d) => d[\"Depth/km\"] < depthQuantile);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3929,"value":"earthquakeSamples = {\n let samp = earthquakes.filter(\n (d) =>\n +d.Longitude >= metadata.geospatial_lon_min &&\n +d.Longitude <= metadata.geospatial_lon_max &&\n +d.Latitude >= metadata.geospatial_lat_min &&\n +d.Latitude <= metadata.geospatial_lat_max\n );\n return samp;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":956,"value":"earthquakes = {\n if (auxData.includes(\"earthquakes\")) {\n return d3\n .text(\n \"https://earthquake.usgs.gov/fdsnws/event/1/query?starttime=1900-01-01T00:00:00&minmag=\" +\n misc.minMagnitude +\n \"&orderby=magnitude&format=text&limit=\" +\n misc.maxEq +\n \"&maxlat=\" +\n yExtent[1] +\n \"&minlon=\" +\n xExtent[0] +\n \"&maxlon=\" +\n xExtent[1] +\n \"&minlat=\" +\n yExtent[0] +\n \"&nodata=404\"\n )\n .then((text) => d3.dsvFormat(\"|\").parse(text));\n } else {\n return [];\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3944,"value":"volcanoSamples = {\n let samp = volcanoes.filter(\n (d) =>\n +d.long >= metadata.geospatial_lon_min &&\n +d.long <= metadata.geospatial_lon_max &&\n +d.lat >= metadata.geospatial_lat_min &&\n +d.lat <= metadata.geospatial_lat_max\n );\n return samp;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":773,"value":" volcanoes = d3\n .dsvFormat(\"|\")\n .parse(removeComments(await FileAttachment(\"WOVOdat_volcano_locations.csv\").text()))","pinned":false,"hidden":false,"mode":"ojs"},{"id":4632,"value":"activeFaults = FileAttachment(\"gem_active_faults.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":3635,"value":"#### - Line data (land)","pinned":false,"hidden":false,"mode":"md"},{"id":4806,"value":"states = topojson.feature(states10m, states10m.objects.states)","pinned":false,"hidden":false,"mode":"ojs"},{"id":4805,"value":"states10m = FileAttachment(\"states-10m.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":4549,"value":"land = topojson.feature(world, world.objects.land)","pinned":false,"hidden":false,"mode":"ojs"},{"id":107,"value":"world = FileAttachment(\"land-110m.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":4558,"value":"countriesGeometry = topojson.feature(countries, countries.objects.countries)","pinned":false,"hidden":false,"mode":"ojs"},{"id":4553,"value":"countries = FileAttachment(\"countries-110m@1.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":3670,"value":"#### - Line data (plate boundaries)","pinned":false,"hidden":false,"mode":"md"},{"id":3954,"value":"function pbFilter(pb) {\n return pb;\n let samp = pb.filter(\n (d) =>\n +d.longitude >= metadata.geospatial_lon_min &&\n +d.longitude <= metadata.geospatial_lon_max &&\n +d.latitude >= metadata.geospatial_lat_min &&\n +d.latitude <= metadata.geospatial_lat_max\n );\n return samp;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":3952,"value":"divergentSamples = pbFilter(divergent)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2820,"value":"divergent = FileAttachment(\"divergent.csv\").csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":3958,"value":"convergentSamples = pbFilter(convergent)","pinned":true,"hidden":false,"mode":"ojs"},{"id":2857,"value":"convergent = FileAttachment(\"convergent.csv\").csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":3960,"value":"transformSamples = pbFilter(transform)","pinned":true,"hidden":false,"mode":"ojs"},{"id":2859,"value":"transform = FileAttachment(\"transform.csv\").csv({ typed: true })","pinned":false,"hidden":false,"mode":"ojs"},{"id":2897,"value":"### Dependencies","pinned":false,"hidden":false,"mode":"md"},{"id":1548,"value":"d3 = require(\"d3@6\", \"d3-geo-projection@3\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":16,"value":"import { removeComments } from \"@unavco/earthscope-notebook-utilities\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":3175,"value":"import {\n readRawData,\n getMetadata,\n getDimensions,\n getVariables,\n modelURL,\n extractParameterData\n} from \"@unavco/emc-notebook-netcdf-utilities\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2981,"value":"import { zoom } from \"@d3/versor-zooming\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":3347,"value":"import { Legend, Swatches } from \"@d3/color-legend\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":3459,"value":"import { legend } from \"@johnhaldeman/color-legend-invertable\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":4920,"value":"plotUmdMin2 = FileAttachment(\"plot.umd.min (2).js\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W37","stdlib":"1","titling":"auto"}