Hurricane Force Supercomputing: Petascale Simulations of Sandy

By Peter Johnsen, Mark Straka, Melvyn Shapiro, Alan Norton & Tom Galerneau

November 14, 2013

The devastation incurred by the landfall of Hurricane Sandy on the northeast coast of the United States just over one year ago exemplifies the need for further advances in accuracy and reliability in numerical weather prediction.  High resolution numerical weather simulations carried out on hundreds of thousands of processors on the largest supercomputers are providing these very insights.

The National Center for Atmospheric Research (NCAR) Weather Research and Forecasting (WRF) model has been employed on the largest yet storm prediction model using real data of over 4 billion points to simulate the landfall of Hurricane Sandy. Using an unprecedented 13,680 nodes (437,760 cores) of the Cray XE6 Blue Waters supercomputer at the National Center for Supercomputing Applications  at the University of Illinois, the team of Peter Johnsen from Cray, Inc., Mark Straka from NCSA, and Mel Shapiro, Alan Norton, and Tom Galarneau from NCAR achieved an unprecedented level of performance for any weather model.  The model used approximately 4 billion grid points at an extremely fine resolution of 500 meters.  Forecast data was written and analyzed by the NCAR team members using the NCAR VAPOR visualization suite.

The landfall of Hurricane Sandy along the New Jersey shoreline late on October 30th, 2012 produced a catastrophic storm surge extending from New Jersey to Rhode Island. The research highlighted here demonstrates the capability of the NCSA/Cray Blue Waters supercomputer to conduct a cloud-resolving WRF-ARW simulation of an intense cyclone over a relatively large domain at a very-fine spatial resolution.

The Blue Waters system is a Cray XE/XK hybrid machine composed of 362,240 AMD 6276 “Interlagos” processors and 4224 NVIDIA GK110 Kepler accelerators all connected by the Cray Gemini 3-D (24^3) torus interconnect. It provides sustained performance of 1 petaflop on a range of real-world science and engineering applications. Our motivation was to reduce time to solution as much as was under our control without major source code restructuring. The WRF version 3.3.1 source code was modified from the public distribution chiefly with concerns for reducing the I/O burden per MPI task and limiting the necessary information to a single MPI rank.

Topology Considerations Are Vital

Domain configuration and process layout using MPI rank ordering features of the Cray XE6 job scheduler (ALPS) form a cornerstone in efficiently using the XE6 3D torus interconnect and allowing WRF to scale this successfully.  We used the Cray grid_order utility to generate improved placement of the ranks for the primary communication pattern in the WRF solver, which is nearest neighbor halo exchanges. Reducing the number of neighbors communicating off-node is the primary goal. Using an alternate placement allows us to get 3 communication partners for most MPI ranks on the same node, instead of only 2, as would be with the default placement.  At very high scales, this strategy improves overall WRF performance by 18% or more.

We found the most effective way to run WRF on the AMD Bulldozer core-modules was to exploit WRF’s “hybrid” MPI/OpenMP structure, utilizing 2 OpenMP threads per MPI rank.  This puts 16 MPI ranks on each XE6 node.

The optimized placement we’ve employed also has the benefit of sending smaller east-west direction exchanges off-node and keeping as many larger north-south messages on-node as possible, resulting in 75% fewer bytes being sent over the network. We verified empirically the long-known tactic of decomposing WRF’s domain with many fewer MPI ranks in the X direction than the Y, as this leads to longer vectors on the inner compute loops.

Our simulations yielded an average Tflop count of 32.454 Tflops per second, per simulation time step. Parallel efficiency was still above 60% even on 13,680 XE6 nodes. Over 12 million off-node halo exchange messages totaling 280 GB were processed every WRF time step.

I/O Considerations at Scale

On the Blue Waters system, the Lustre file system was used for all file activity.

Two techniques were used to handle the large I/O requirements for the Sandy simulation –

  1. Parallel NetCDF (PnetCDF), jointly developed by Northwestern University and Argonne National Laboratory, was used where practical.  The MPICH library from Cray has a tuned MPI-IO implementation that aligns parallel I/O with the Lustre file system.  This format is required when post-processing tools are used.
  2. WRF has a multi-file option where each subdomain, or MPI rank, reads and writes unique files.  This was used for very large restart files and some of the pre-processing steps.  The Blue Waters Lustre file system was able to open and read 145,920 restart files in 18 seconds for a 4560 node case.

Additionally, use of WRF’s auxiliary history output options to select only the output fields of greatest interest, thus reducing the volume of output considerably, was of great utility in our work.

graph_1

Scalability of Hurricane Sandy run.  Sustained performance in Tflops/second (y-axis, left) and parallel efficiency over base run on 8,192 cores (y-axis, right) are shown.

Forecast Analysis and Validation

The following figures show a comparison of the maximum radar reflectivity (a surrogate for precipitation) from the simulations at 3-km and 500-m horizontal resolution. In both simulations, a broad region of heavy precipitation is located on the west and southwest side of Sandy, and is organized in a region where warm moist northeasterly flow intersects a northwesterly surge of cold continental air (not shown).

graph_2

Comparison of (a) 3-km and (b) 500-m horizontal resolution ARW simulations of maximum radar reflectivity (shaded according to the color bar in dBZ) verifying at 1500 UTC 29 October 2012.

The 500-m simulation is superior to that at 3-km because it shows the fine-scale linear structure of the convective precipitation bands, consistent with the available observations (not shown). The next images show a zoomed-in view of maximum radar reflectivity and 300-m wind speed within the inner-core of Sandy at 1800 UTC 29 October 2012. This zoomed perspective allows for examination of the full detail of the simulation, noting that the resolution of the simulation (7000×7000 grid points) exceeds the resolution of standard computer monitors by a factor of seven. Here we note the utility of ultra-advanced computational capability to represent the full range of scales spanning the storm-scale circulations down to fine-scale turbulent motions and individual cloud and precipitation systems.

graph_3

500-m ARW simulation of (a) maximum radar reflectivity (shaded according to the color bar in dBZ) and (b) 300-m wind speed (shaded according to the color bar in m s−1) verifying at 1800 UTC 29 October 2012.

The model accuracy for predicting such key output fields as rainfall, pressures, wind speeds, and storm track was graphically validated against actual atmospheric measurements from the storm using NCAR’s VAPOR software suite. Given recent advances in accessing and displaying large volume geophysical datasets as exemplified by the VAPOR software, it is now possible to view the full temporal evolution of numerical simulations and predictions of atmospheric and other geophysical systems. Examples of the advanced visualizations of Hurricane Sandy with VAPOR can be found at:

https://www.vapor.ucar.edu/sites/default/files/movies/sandy_SC13_web_0.mp4

The results of this research will be presented at the Supercomputing conference this month in November.  See the conference agenda here:

http://sc13.supercomputing.org/schedule/event_detail.php?evid=pap255

Cutting Edge Forecasting

NOAA has initiated the ten-year Hurricane Forecast Improvement Project (HFIP), which is evaluating a variety of modeling approaches, exploring the feasibility of real-time fine-scale hurricane projections. Its enhanced Hurricane WRF model (HWRF) is already being run in real time at a somewhat smaller scale. In a collaborative effort involving NOAA’s Hurricane Research Division and Environmental Modeling Center, Cray , NCSA, and NCAR, this code is already being run on Blue Waters to conduct performance studies at scale with grid nesting never before possible. Results are already promising that coming years’ hurricane seasons will be able to incorporate much finer detailed real-time forecasts generated by these simulations. The team is also exploring high resolution simulations with the Office of Naval Research ONR using the COAMPS model.

Research Team:

Peter Johnsen is a performance engineer and meteorologist with Cray, Inc.  Peter’s expertise is optimizing environmental applications on HPC systems.

Mark Straka specializes in performance analysis of scientific applications on the Blue Waters system at the National Center for Supercomputing Applications.

Melvyn Shapiro, Alan Norton, and Thomas Galarneau are research meteorologists with the National Center for Atmospheric Research and are studying many weather phenomena, including Hurricane Sandy’s unique nature.

Subscribe to HPCwire's Weekly Update!

Be the most informed person in the room! Stay ahead of the tech trends with industy updates delivered to you every week!

With New Owner and New Roadmap, an Independent Omni-Path Is Staging a Comeback

July 23, 2021

Put on a shelf by Intel in 2019, Omni-Path faced a uncertain future, but under new custodian Cornelis Networks, OmniPath is looking to make a comeback as an independent high-performance interconnect solution. A "significant refresh" – called Omni-Path Express – is coming later this year according to the company. Cornelis Networks formed last September as a spinout of Intel's Omni-Path division. Read more…

PEARC21 Panel Reviews Eight New NSF-Funded HPC Systems Debuting in 2021

July 23, 2021

Over the past few years, the NSF has funded a number of HPC systems to further supply the open research community with computational resources to meet that community’s changing and expanding needs. A review of these systems at the PEARC21 conference (July 19-22) highlighted... Read more…

Chameleon’s HPC Testbed Sharpens Its Edge, Presses ‘Replay’

July 22, 2021

“One way of saying what I do for a living is to say that I develop scientific instruments,” said Kate Keahey, a senior fellow at the University of Chicago and a computer scientist at Argonne National Laboratory, as s Read more…

PEARC21 Plenary Session: AI for Innovative Social Work

July 21, 2021

AI analysis of social media poses a double-edged sword for social work and addressing the needs of at-risk youths, said Desmond Upton Patton, senior associate dean, Innovation and Academic Affairs, Columbia University. S Read more…

Summer Reading: “High-Performance Computing Is at an Inflection Point”

July 21, 2021

At last month’s 11th International Symposium on Highly Efficient Accelerators and Reconfigurable Technologies (HEART), a group of researchers led by Martin Schulz of the Leibniz Supercomputing Center (Munich) presented a “position paper” in which they argue HPC architectural landscape... Read more…

AWS Solution Channel

Accelerate innovation in healthcare and life sciences with AWS HPC

With Amazon Web Services, researchers can access purpose-built HPC tools and services along with scientific and technical expertise to accelerate the pace of discovery. Whether you are sequencing the human genome, using AI/ML for disease detection or running molecular dynamics simulations to develop lifesaving drugs, AWS has the infrastructure you need to run your HPC workloads. Read more…

PEARC21 Panel: Wafer-Scale-Engine Technology Accelerates Machine Learning, HPC

July 21, 2021

Early use of Cerebras’ CS-1 server and wafer-scale engine (WSE) has demonstrated promising acceleration of machine-learning algorithms, according to participants in the Scientific Research Enabled by CS-1 Systems panel Read more…

With New Owner and New Roadmap, an Independent Omni-Path Is Staging a Comeback

July 23, 2021

Put on a shelf by Intel in 2019, Omni-Path faced a uncertain future, but under new custodian Cornelis Networks, OmniPath is looking to make a comeback as an independent high-performance interconnect solution. A "significant refresh" – called Omni-Path Express – is coming later this year according to the company. Cornelis Networks formed last September as a spinout of Intel's Omni-Path division. Read more…

Chameleon’s HPC Testbed Sharpens Its Edge, Presses ‘Replay’

July 22, 2021

“One way of saying what I do for a living is to say that I develop scientific instruments,” said Kate Keahey, a senior fellow at the University of Chicago a Read more…

Summer Reading: “High-Performance Computing Is at an Inflection Point”

July 21, 2021

At last month’s 11th International Symposium on Highly Efficient Accelerators and Reconfigurable Technologies (HEART), a group of researchers led by Martin Schulz of the Leibniz Supercomputing Center (Munich) presented a “position paper” in which they argue HPC architectural landscape... Read more…

PEARC21 Panel: Wafer-Scale-Engine Technology Accelerates Machine Learning, HPC

July 21, 2021

Early use of Cerebras’ CS-1 server and wafer-scale engine (WSE) has demonstrated promising acceleration of machine-learning algorithms, according to participa Read more…

15 Years Later, the Green500 Continues Its Push for Energy Efficiency as a First-Order Concern in HPC

July 15, 2021

The Green500 list, which ranks the most energy-efficient supercomputers in the world, has virtually always faced an uphill battle. As Wu Feng – custodian of the Green500 list and an associate professor at Virginia Tech – tells it, “noone" cared about energy efficiency in the early 2000s, when the seeds... Read more…

Frontier to Meet 20MW Exascale Power Target Set by DARPA in 2008

July 14, 2021

After more than a decade of planning, the United States’ first exascale computer, Frontier, is set to arrive at Oak Ridge National Laboratory (ORNL) later this year. Crossing this “1,000x” horizon required overcoming four major challenges: power demand, reliability, extreme parallelism and data movement. Read more…

Quantum Roundup: IBM, Rigetti, Phasecraft, Oxford QC, China, and More

July 13, 2021

IBM yesterday announced a proof for a quantum ML algorithm. A week ago, it unveiled a new topology for its quantum processors. Last Friday, the Technical Univer Read more…

ExaWind Prepares for New Architectures, Bigger Simulations

July 10, 2021

The ExaWind project describes itself in terms of terms like wake formation, turbine-turbine interaction and blade-boundary-layer dynamics, but the pitch to the Read more…

AMD Chipmaker TSMC to Use AMD Chips for Chipmaking

May 8, 2021

TSMC has tapped AMD to support its major manufacturing and R&D workloads. AMD will provide its Epyc Rome 7702P CPUs – with 64 cores operating at a base cl Read more…

Intel Launches 10nm ‘Ice Lake’ Datacenter CPU with Up to 40 Cores

April 6, 2021

The wait is over. Today Intel officially launched its 10nm datacenter CPU, the third-generation Intel Xeon Scalable processor, codenamed Ice Lake. With up to 40 Read more…

Berkeley Lab Debuts Perlmutter, World’s Fastest AI Supercomputer

May 27, 2021

A ribbon-cutting ceremony held virtually at Berkeley Lab's National Energy Research Scientific Computing Center (NERSC) today marked the official launch of Perlmutter – aka NERSC-9 – the GPU-accelerated supercomputer built by HPE in partnership with Nvidia and AMD. Read more…

Ahead of ‘Dojo,’ Tesla Reveals Its Massive Precursor Supercomputer

June 22, 2021

In spring 2019, Tesla made cryptic reference to a project called Dojo, a “super-powerful training computer” for video data processing. Then, in summer 2020, Tesla CEO Elon Musk tweeted: “Tesla is developing a [neural network] training computer called Dojo to process truly vast amounts of video data. It’s a beast! … A truly useful exaflop at de facto FP32.” Read more…

Google Launches TPU v4 AI Chips

May 20, 2021

Google CEO Sundar Pichai spoke for only one minute and 42 seconds about the company’s latest TPU v4 Tensor Processing Units during his keynote at the Google I Read more…

CentOS Replacement Rocky Linux Is Now in GA and Under Independent Control

June 21, 2021

The Rocky Enterprise Software Foundation (RESF) is announcing the general availability of Rocky Linux, release 8.4, designed as a drop-in replacement for the soon-to-be discontinued CentOS. The GA release is launching six-and-a-half months after Red Hat deprecated its support for the widely popular, free CentOS server operating system. The Rocky Linux development effort... Read more…

CERN Is Betting Big on Exascale

April 1, 2021

The European Organization for Nuclear Research (CERN) involves 23 countries, 15,000 researchers, billions of dollars a year, and the biggest machine in the worl Read more…

Iran Gains HPC Capabilities with Launch of ‘Simorgh’ Supercomputer

May 18, 2021

Iran is said to be developing domestic supercomputing technology to advance the processing of scientific, economic, political and military data, and to strengthen the nation’s position in the age of AI and big data. On Sunday, Iran unveiled the Simorgh supercomputer, which will deliver.... Read more…

Leading Solution Providers

Contributors

HPE Launches Storage Line Loaded with IBM’s Spectrum Scale File System

April 6, 2021

HPE today launched a new family of storage solutions bundled with IBM’s Spectrum Scale Erasure Code Edition parallel file system (description below) and featu Read more…

Julia Update: Adoption Keeps Climbing; Is It a Python Challenger?

January 13, 2021

The rapid adoption of Julia, the open source, high level programing language with roots at MIT, shows no sign of slowing according to data from Julialang.org. I Read more…

10nm, 7nm, 5nm…. Should the Chip Nanometer Metric Be Replaced?

June 1, 2020

The biggest cool factor in server chips is the nanometer. AMD beating Intel to a CPU built on a 7nm process node* – with 5nm and 3nm on the way – has been i Read more…

GTC21: Nvidia Launches cuQuantum; Dips a Toe in Quantum Computing

April 13, 2021

Yesterday Nvidia officially dipped a toe into quantum computing with the launch of cuQuantum SDK, a development platform for simulating quantum circuits on GPU-accelerated systems. As Nvidia CEO Jensen Huang emphasized in his keynote, Nvidia doesn’t plan to build... Read more…

Microsoft to Provide World’s Most Powerful Weather & Climate Supercomputer for UK’s Met Office

April 22, 2021

More than 14 months ago, the UK government announced plans to invest £1.2 billion ($1.56 billion) into weather and climate supercomputing, including procuremen Read more…

Q&A with Jim Keller, CTO of Tenstorrent, and an HPCwire Person to Watch in 2021

April 22, 2021

As part of our HPCwire Person to Watch series, we are happy to present our interview with Jim Keller, president and chief technology officer of Tenstorrent. One of the top chip architects of our time, Keller has had an impactful career. Read more…

Quantum Roundup: IBM, Rigetti, Phasecraft, Oxford QC, China, and More

July 13, 2021

IBM yesterday announced a proof for a quantum ML algorithm. A week ago, it unveiled a new topology for its quantum processors. Last Friday, the Technical Univer Read more…

Senate Debate on Bill to Remake NSF – the Endless Frontier Act – Begins

May 18, 2021

The U.S. Senate today opened floor debate on the Endless Frontier Act which seeks to remake and expand the National Science Foundation by creating a technology Read more…

  • arrow
  • Click Here for More Headlines
  • arrow
HPCwire