MM5 (weather model)

Last updated

The MM5 (short for Fifth-Generation Penn State/NCAR Mesoscale Model) is a regional mesoscale model used for creating weather forecasts and climate projections. It is a community model maintained by Penn State University and the National Center for Atmospheric Research. The MM5 is a limited-area, terrain-following sigma coordinate model that is used to replicate or forecast mesoscale and regional scale atmospheric circulation. [1] It has been updated many times since the 1970s to fix bugs, adapt to new technologies, and work on different types of computers and software. It is used in many different ways: for research and for weather prediction. In research, it is used to compare it to other models, to see what works and what does not work. It is also used for air quality models. [2]

Contents

Active development on the model ended with version 3.7.2 in 2005, and it has been largely superseded by the Weather Research and Forecasting (WRF) model. [3]

Features

MM5 is globally relocatable, which helps support different latitudes, terrain types, elevations, soil types, etc. . The model can be either hydrostatic or non-hydrostatic, depending on the desired outcome. The fact that the model is regional implies that it requires initial conditions and lateral boundary conditions. This means that each boundary (there are four) has initialized wind speeds, temperatures, pressure and moisture fields. [4] Thus, gridded data is needed for this program. This model takes and then analyzes its data based on pressure surfaces. However, these surfaces must first be interpolated by a specific vertical coordinate before it can be analyzed. [4] This vertical coordinate, sigma, is computed and then used throughout the program. Σ is defined as: Σ = (p-pt)/p*, p* = ps-pt, Where p is pressure, ps is surface pressure, and pt is the pressure at the top of the model. [5] When Σ is close to the ground, the program follows the actual terrain, but when Σ is higher up, the program looks at isobaric surfaces. Σ ranges from 0 to 1. [4] [6] It has adaptable and multiple nesting capabilities, which allows multiple programs to run at once, while utilizing 2-way nesting. MM5 features inputs from actual data, which is helpful because routine observations can be used. Then, data can be compared and used in context with other models. [6] MM5 also features terrain-following vertical coordinates and four-dimensional data assimilation (FDDA). [6] FDDA is used when there is a lot of data that was taken over a longer period of time. Then this data that needed to be taken over a longer period of time gets placed into FDDA. It is also utilized for dynamical initialization and four-dimensional data sets. [4] Most importantly, MM5 is well documented and has many places for user support.

Functions

There are many pieces of the MM5 modeling system, each with different functions. These pieces include: TERRAIN, REGRID, LITTLE_R, NESTDOWN, INTERPF, INTERPB, and GRAPH/RIP. In the basic sense, the model begins by taking the information and then makes a model of the land. This is done by using the piece TERRAIN. [7] Then, the model guesses the different pressures that exist in the atmosphere throughout this plot of land, which is done through REGRID. [7] Then, the model takes these guesses, and with the help of some surface observations, can create an objective analysis, using LITTLE_R. [8] RAWINS is an older version of LITTLE_R and thus LITTLE_R is more commonly used. Then INTERPF takes that data from RAWINS/LITTLE_R and REGRID to interpolate that data to the sigma coordinate explained above. [7] MM5 itself computes the time combination. The NESTDOWN function has the ability to change the vertical sigma levels. INTERPB generates the first guess for RAWINS, as well as gathers the files used to make guesses for REGRID. Finally, GRAPH/RIP generates the plots that users can view from all of the data given from the other programs. [7] RIP stands for Read/Interpolate/Plot, and it uses NCAR Graphics to help envisage the output from the MM5. [9]

Flow chart of the functions of the MM5 Mm5.gif
Flow chart of the functions of the MM5

Memory and code basics

MM5 is written in FORTRAN. These FORTRAN programs must be compiled on a local computer, and some need to be recompiled each time a model configuration is changed. [10] The program uses pointers to assign variables to values. These pointers go into parts of the memory to assign specific values to the desired variables. [11] MM5 can also perform multiple tasks at once. Specifically, two different tasks can be executed on different processors simultaneously and MM5 uses this as much as possible. [11] This multi-tasking also uses nesting, and MM5 allows up to nine domains (processes) running at the same time and they interact throughout the whole process. [4] The model utilizes two-way nesting, which occurs when the input from one nest's coarse mesh, which is a high density of cells in one area, comes from one of the four boundaries, but the feedback to the coarser mesh happens through the nest interior. [4] Each domain gathers information from its parent domain each timestep, then it runs three timesteps, and then delivers the information back to its parent domain. [4] There are three different ways to do two-way nesting: Nest interpolation, nest analysis input, and nest terrain input. [4] Interpolation occurs when the terrain is smooth, such as water. There is no input that is required for this type of two-way nesting. Nest input requires a file called MMINPUT, and that file contains the meteorological and the terrain information so that initially, a better analysis can be done. Finally terrain input requires a TERRAIN file. Then the meteorological fields are interpolated. [4] When multi-tasking occurs, the variables must be marked as either shared or private. Shared implies that the processors all have access to the same part of memory, while private implies that each processor must have its own private copy of an array with its personal memory location. [11] The multi-tasking occurs specifically in the subroutines Solve1, Solve3, and Sound.

Requirements

MM5 can be run in many different ways, depending on the computer. The model can be run on a single processor computer, shared memory architecture or distributed memory architecture. It can be run on many different platforms as well such as IBM, SGI origin 200/2000, CRAYs (J90, C90, T3E), DEC Alphas, Alphas running Linux, Sun, and more. The computer it is run on must have a Fortran 90 and 77 compilers and a C compiler. Additionally, it can have NCAR graphics and MPI tools for running the MPI MM5 model. However, this is not necessary. The MM5 model takes up at least half a gigabyte of memory and a few gigabytes of disk space. [6]

Source code

The MM5 code itself has over 220 subroutines, with more than 55,000 lines of code. [12] It uses Standard Fortran 77 with “Cray” pointers. A list of the subroutines of all of the functions within the MM5, with a short description and the routines that are called within it can be found here. [13]

Use

The MM5 had been adapted for use in many different types of atmospheric simulations:

Developments and improvements

TERRAIN's latest update uses USGS global 30 sec terrain elevation data, which looks at the entire world, where before was only used in the continental United States. [21] This new update also allows for a better quality terrain image that can be used in MM5. REGRID's improvement makes it simpler for its users to input the data and it also is more portable. [21] LIITLE_R was developed in 2001 to replace RAWINS. LITTLE_R's improvement again makes it easier for users to input data. [21]

See also

Related Research Articles

<span class="mw-page-title-main">Fortran</span> General-purpose programming language

Fortran is a third generation, compiled, imperative programming language that is especially suited to numeric computation and scientific computing.

In computer programming, a parameter or a formal argument is a special kind of variable used in a subroutine to refer to one of the pieces of data provided as input to the subroutine. These pieces of data are the values of the arguments with which the subroutine is going to be called/invoked. An ordered list of parameters is usually included in the definition of a subroutine, so that, each time the subroutine is called, its arguments for that call are evaluated, and the resulting values can be assigned to the corresponding parameters.

IDL, short for Interactive Data Language, is a programming language used for data analysis. It is popular in particular areas of science, such as astronomy, atmospheric physics and medical imaging. IDL shares a common syntax with PV-Wave and originated from the same codebase, though the languages have subsequently diverged in detail. There are also free or costless implementations, such as GNU Data Language (GDL) and Fawlty Language (FL).

<span class="mw-page-title-main">Numerical weather prediction</span> Weather prediction using mathematical models of the atmosphere and oceans

Numerical weather prediction (NWP) uses mathematical models of the atmosphere and oceans to predict the weather based on current weather conditions. Though first attempted in the 1920s, it was not until the advent of computer simulation in the 1950s that numerical weather predictions produced realistic results. A number of global and regional forecast models are run in different countries worldwide, using current weather observations relayed from radiosondes, weather satellites and other observing systems as inputs.

GRIB is a concise data format commonly used in meteorology to store historical and forecast weather data. It is standardized by the World Meteorological Organization's Commission for Basic Systems, known under number GRIB FM 92-IX, described in WMO Manual on Codes No.306. Currently there are three versions of GRIB. Version 0 was used to a limited extent by projects such as TOGA, and is no longer in operational use. The first edition is used operationally worldwide by most meteorological centers, for Numerical Weather Prediction output (NWP). A newer generation has been introduced, known as GRIB second edition, and data is slowly changing over to this format. Some of the second-generation GRIB is used for derived products distributed in the Eumetcast of Meteosat Second Generation. Another example is the NAM model.

Data assimilation is a mathematical discipline that seeks to optimally combine theory with observations. There may be a number of different goals sought – for example, to determine the optimal state estimate of a system, to determine initial conditions for a numerical forecast model, to interpolate sparse observation data using knowledge of the system being observed, to set numerical parameters based on training a model from observed data. Depending on the goal, different solution methods may be used. Data assimilation is distinguished from other forms of machine learning, image analysis, and statistical methods in that it utilizes a dynamical model of the system being analyzed.

Automatic parallelization, also auto parallelization, or autoparallelization refers to converting sequential code into multi-threaded and/or vectorized code in order to use multiple processors simultaneously in a shared-memory multiprocessor (SMP) machine. Fully automatic parallelization of sequential programs is a challenge because it requires complex program analysis and the best approach may depend upon parameter values that are not known at compilation time.

<span class="mw-page-title-main">National Center for Atmospheric Research</span> US federally funded research and development center

The US National Center for Atmospheric Research is a US federally funded research and development center (FFRDC) managed by the nonprofit University Corporation for Atmospheric Research (UCAR) and funded by the National Science Foundation (NSF). NCAR has multiple facilities, including the I. M. Pei-designed Mesa Laboratory headquarters in Boulder, Colorado. Studies include meteorology, climate science, atmospheric chemistry, solar-terrestrial interactions, environmental and societal impacts.

<span class="mw-page-title-main">Tropical cyclone forecast model</span> Computer program that uses meteorological data to forecast tropical cyclones

A tropical cyclone forecast model is a computer program that uses meteorological data to forecast aspects of the future state of tropical cyclones. There are three types of models: statistical, dynamical, or combined statistical-dynamic. Dynamical models utilize powerful supercomputers with sophisticated mathematical modeling software and meteorological data to calculate future weather conditions. Statistical models forecast the evolution of a tropical cyclone in a simpler manner, by extrapolating from historical datasets, and thus can be run quickly on platforms such as personal computers. Statistical-dynamical models use aspects of both types of forecasting. Four primary types of forecasts exist for tropical cyclones: track, intensity, storm surge, and rainfall. Dynamical models were not developed until the 1970s and the 1980s, with earlier efforts focused on the storm surge problem.

<span class="mw-page-title-main">Atmospheric model</span> Mathematical model of atmospheric motions

In atmospheric science, an atmospheric model is a mathematical model constructed around the full set of primitive, dynamical equations which govern atmospheric motions. It can supplement these equations with parameterizations for turbulent diffusion, radiation, moist processes, heat exchange, soil, vegetation, surface water, the kinematic effects of terrain, and convection. Most atmospheric models are numerical, i.e. they discretize equations of motion. They can predict microscale phenomena such as tornadoes and boundary layer eddies, sub-microscale turbulent flow over buildings, as well as synoptic and global flows. The horizontal domain of a model is either global, covering the entire Earth, or regional (limited-area), covering only part of the Earth. The different types of models run are thermotropic, barotropic, hydrostatic, and nonhydrostatic. Some of the model types make assumptions about the atmosphere which lengthens the time steps used and increases computational speed.

CALPUFF is an advanced, integrated Lagrangian puff modeling system for the simulation of atmospheric pollution dispersion distributed by the Atmospheric Studies Group at TRC Solutions.

<span class="mw-page-title-main">Weather Research and Forecasting Model</span> Numerical weather prediction system

The Weather Research and Forecasting (WRF) Model is a numerical weather prediction (NWP) system designed to serve both atmospheric research and operational forecasting needs. NWP refers to the simulation and prediction of the atmosphere with a computer model, and WRF is a set of software for this. WRF features two dynamical (computational) cores, a data assimilation system, and a software architecture allowing for parallel computation and system extensibility. The model serves a wide range of meteorological applications across scales ranging from meters to thousands of kilometers.

The MEMO model is a Eulerian non-hydrostatic prognostic mesoscale model for wind-flow simulation. It was developed by the Aristotle University of Thessaloniki in collaboration with the Universität Karlsruhe. The MEMO Model together with the photochemical dispersion model MARS are the two core models of the European zooming model (EZM). This model belongs to the family of models designed for describing atmospheric transport phenomena in the local-to-regional scale, frequently referred to as mesoscale air pollution models.

NetCDF is a set of software libraries and self-describing, machine-independent data formats that support the creation, access, and sharing of array-oriented scientific data. The project homepage is hosted by the Unidata program at the University Corporation for Atmospheric Research (UCAR). They are also the chief source of netCDF software, standards development, updates, etc. The format is an open standard. NetCDF Classic and 64-bit Offset Format are an international standard of the Open Geospatial Consortium.

In computer programming, a function, procedure, method, subroutine, routine, or subprogram is a callable unit of software logic that has a well-defined interface and behavior and can be invoked multiple times.

<span class="mw-page-title-main">Sigma coordinate system</span> Coordinate system used in fluid dynamics

The sigma coordinate system is a common coordinate system used in computational models for oceanography, meteorology and other fields where fluid dynamics are relevant. This coordinate system receives its name from the independent variable used to represent a scaled pressure level.

<span class="mw-page-title-main">NCAR-Wyoming Supercomputing Center</span> High performance computing center in Wyoming, US

The NCAR-Wyoming Supercomputing Center (NWSC) is a high-performance computing (HPC) and data archival facility located in Cheyenne, Wyoming, that provides advanced computing services to researchers in the Earth system sciences.

PROSE was the mathematical 4GL virtual machine that established the holistic modeling paradigm known as Synthetic Calculus. A successor to the SLANG/CUE simulation and optimization language developed at TRW Systems, it was introduced in 1974 on Control Data supercomputers. It was the first commercial language to employ automatic differentiation (AD), which was optimized to loop in the instruction-stack of the CDC 6600 CPU.

<span class="mw-page-title-main">Regional Ocean Modeling System</span> Free-surface, terrain-following, primitive equations ocean model

Regional Ocean Modeling System (ROMS) is a free-surface, terrain-following, primitive equations ocean model widely used by the scientific community for a diverse range of applications. The model is developed and supported by researchers at the Rutgers University, University of California Los Angeles and contributors worldwide.

The Model for Prediction Across Scales (MPAS) is an Earth system modeling software that integrates atmospheric, oceanographic, and cryospheric modeling across scales from regional to planetary. It includes climate and weather modeling and simulations that were used initially by researchers in 2013. The atmospheric models were created by the Earth System Laboratory at the National Center for Atmospheric Research and the oceanographic models were created by the Climate, Ocean, and Sea Ice Modeling Group at Los Alamos National Laboratory. The software has been used to model real-time weather as well as seasonal forecasting of convection, tornadoes and tropical cyclones. The atmospheric modeling component of the software can be used with other atmospheric modeling software including the Weather Research and Forecasting Model, the Global Forecast System, and the Community Earth System Model.

References

  1. "MM5 Community Model Homepage". .mmm.ucar.edu. Retrieved 2022-04-20.
  2. http://www2.mmm.ucar.edu/mm5/workshop/ws03/sessionJ1/Klausmann.pdf [ bare URL PDF ]
  3. Dudhia, Jimy (June 2005). MM5 Version 3.7 (The Final Version). WRF/MM5 Users' Workshop – June 2005. Boulder, Colorado: National Center for Atmospheric Research.
  4. 1 2 3 4 5 6 7 8 9 "Chapter 1: Introduction". .mmm.ucar.edu. Retrieved 2022-04-20.
  5. http://www2.mmm.ucar.edu/mm5/documents/mm5-code-pdf/sec1.pdf [ bare URL PDF ]
  6. 1 2 3 4 "MM5 Community Model Homepage".
  7. 1 2 3 4 "MM5 Community Model Homepage". .mmm.ucar.edu. Retrieved 2022-04-20.
  8. "MM5 OnLine Tutorial : LITTLE_R".
  9. "MM5 OnLine Tutorial : RIP".
  10. "Chapter 2: Getting Started".
  11. 1 2 3 http://www2.mmm.ucar.edu/mm5/documents/mm5-code-pdf/sec2.pdf [ bare URL PDF ]
  12. ww2.mmm.ucar.edu
  13. http://www2.mmm.ucar.edu/mm5/documents/mm5-code-pdf/sec5.pdf [ bare URL PDF ]
  14. "SUNY Stony Brook MM5 Mesoscale Forecasts for 2000030612". cheget.msrc.sunysb.edu. Archived from the original on 6 March 2000. Retrieved 6 June 2022.
  15. "Real-Time Mm5 Forecasts". www.atmos.umd.edu. Archived from the original on 5 October 2002. Retrieved 6 June 2022.
  16. "University of Wisconsin MM5 Real Time Forecast and Forecast Sensitivity Page". helios.aos.wisc.edu. Archived from the original on 9 March 2001. Retrieved 6 June 2022.
  17. "Archived copy". Archived from the original on 2008-06-13. Retrieved 2008-06-25.{{cite web}}: CS1 maint: archived copy as title (link)
  18. "Archived copy". Archived from the original on 2008-06-13. Retrieved 2008-06-29.{{cite web}}: CS1 maint: archived copy as title (link)
  19. "Echam5-Mm5".
  20. "Polar Meteorology Group Home Page". Archived from the original on 2008-06-07. Retrieved 2008-06-29.
  21. 1 2 3 "MM5 Community Model Homepage".