Difference between revisions of "Overview of the Generic PCM"

Revision as of 21:15, 31 August 2022

What is a GCM

A Global Climate Model (GCM) like the Generic PCM calculates the temporal evolution of the different variables (temperature, pressure, mass, mixing ratio of tracers/gas/clouds, winds, etc.) that control or describe the planetary meteorology and climate at different points of a 3D “grid” that covers the entire atmosphere.

To do that, the model solves 3D hydrodynamic equations that describe how atmospheric winds are created by various physical forcings (radiative heating/cooling, etc.) and how these winds transport and redistribute the various atmospheric fields (temperature, tracers, etc.). However, it can be shown that to a good level of approximation, the various forcings can be computed using only local conditions in a column of atmosphere, whereas long-range couplings between various points of the planet are mediated through atmospheric transport.

Physics / Dynamics interface

For these reasons, the 3D model is split in two components:

• a 3D dynamical core abble to integrate in time the general equations for atmospheric circulation when the various local forcings are known. This component (including source code) is common to all terrestrial-type atmospheres, and applicable in certain cases to the upper atmospheres of gas giant planets.
• a physical package that is specific to the planet in question (as it takes into account dedicated parametrizations representative of ongoing physical processes) and which computes atmospheric processes in each atmospheric column (which are sub-grid-scale from the point of view of the dynamical core).

In practice, we have several dynamical cores available with the model: LMDZ (global longitude/latitude grid), DYNAMICO (icosaedric grid), and WRF (rectangular box for local/regional studies)

From an initial state, the model calculates the temporal evolution of these variables, one timestep after another:

• At instant t, variables X_t (temperature for example) are known at every grid point of the atmosphere.
• The physical model calculates the local time derivatives of the various variables -- called the tendencies -- ((δX/δt)₁, (δX/δt)₂, etc.) due to each physical phenomenon: for example, heating due to absorption of solar radiation and infrared cooling, mixing due to turbulence in the planetary boundary layer, condensation/evaporation of water vapor, etc. These tendencies are calculated via parameterization of the processes involved that try to capture the physical behavior of the atmopshere at the subgrid scale.
• These tendencies are passed to the dynamical core that integrates the model forward in time to the next time step, t + δt, by solving the hydrodynamical equations.
• This new state can be used to recompute the tendencies for the next time step and so on.

Philosophy of the Generic PCM

The accuracy and predictive power of any GCM is often dictated by the choice made for the physical phenomena that are included and for the way they are parametrized. To be accurate, many Earth-centric climate models use very complex parametrization that are heavily tuned on Earth observations. Such parametrization can, however, behave very poorly far from the parameter space where they have been validated, thus reducing their predictive power for exotic atmospheres.

As the goal of the Generic PCM is to model the climate of very different planets with the same model, a particular emphasis has been put on using simple parametrizations based as much as possible on ab-initio physical equations and conservation principles to avoid unphysical behaviors in extreme regimes.

GCM inputs and outputs

input/output data

The input files are the files needed to initialize the model (state of the atmosphere at instant t0 as well as a dataset of boundary conditions). The output files are “historical files”, archives of the atmospheric flow history as simulated by the model: the “diagfi files”, the “stats files”, the daily averages, and so on. The "parameter files" are the files describing the physical parameters of the planet, as well as the physical parameterizations used in the model. The figure (on the right) shows the main architecture of the Generic PCM.

start and startfi files

GCM simulations needs a list of initial conditions/parameters to be run. This includes, for example, the initial temperature field, wind field, cloud field, position of water reservoirs on the surface, etc.

In the Generic PCM, initial conditions/parameters are stored in two files named "start.nc" and "startfi.nc". At each new run, the GCM reads these two files to set up the proper set of initial conditions.

start vs restart files

TO BE COMPLETED

parameter (.def) files

List of needed .def files:

• callphys.def: Physical parameters
• run.def: Numerical parameters
• gases.def: List of gases
• traceur.def: Numerical traceurs in the atmosphere
• z2sig.def: Pseudo-altitudes (in km) at which the user wants to set the vertical levels
• diagfi.def (optionnal): List of wanted outputs

input data

TO BE COMPLETED

output files

diagfi.nc diagspec_IR.nc / diagspec_VI.nc stats.nc etc.

TO BE COMPLETED

Our team of developers (by alphabetical order)

TO BE COMPLETED