CFMIP Standard Output

Prepared by

 

Mark Webb and Karl E. Taylor

 

 

Based, originally, on AMIP Standard Output Tables developed by

PCMDI and the AMIP panel


A pdf version of this page can be downloaded from here.

  

26 September 2005


Change history: 

Added request for 3D fields to be on model levels rather than pressure levels and a more explicit distinction between stratiform, convective and total 3D cloud fields. Also added some extra diagnostics requested daily. KDW 26/09/05

Changed headings on monthly table to state preference for monthly means over climatological fields, and to allow for annual mean climatological CO2 forcing fields.  MJW 29/11/04

Contents:

·        Background and summary of format requirements

·        Model output summary

·        Tables of output fields

Monthly data and time-independent fields

Table CF1a: 2-d atmosphere or surface data (longitude, latitude, time)

Table CF1b: ISCCP simulator data (longitude, latitude, pressure2, tau, time)

Table CF1c: 3-d atmosphere data on pressure levels (longitude, latitude, pressure, time)

Table CF1d: 3-d data on model levels (longitude, latitude, model_level, time)

Table CF1e: time-independent 2-d surface data (longitude, latitude)

 

Annual mean data

Table CF2a: 2-d atmosphere or surface data (longitude, latitude, time)

Table CF2b: ISCCP simulator data (longitude, latitude, pressure2, tau, time)

 

Daily mean data

Table CF3a: 2-d atmosphere data (longitude, latitude, time)

Table CF3b: ISCCP simulator data (longitude, latitude, pressure2, tau, time)

Table CF3c: 3-d atmosphere data on pressure levels (longitude, latitude, pressure, time)

Table CF3d: 3-d atmosphere data on model levels (longitude, latitude, pressure, time)

 

Doubled CO2 radiative forcing

Table CF4: Monthly-mean 2-d radiative forcing data for doubled co2 (longitude, latitude, time).

 

·        Coordinate dimensions

 

Background and overview of format requirements

In order to facilitate participation in various community modeling experiments (e.g., IPCC, CMIP, AMIP, CFMIP), the Program for Climate Model Diagnosis and Intercomparison (with encouragement from the WGCM and WGNE and in collaboration with project leaders) is attempting to establish a more uniform set of model output requirements.  CFMIP was initiated before this effort had matured sufficiently, so initial CFMIP output has been collected in various structures and formats, but to facilitate analysis of future model output, any new CFMIP data should conform to the new requirements.

Modeling groups contributing output to the CFMIP database are henceforth asked to ensure that it meets rather strict format and metadata requirements.  These requirements yield files produced in network Common Data Form (netCDF; see http://www.unidata.ucar.edu/packages/netcdf), which has become the most popular form for exchanging ocean-atmosphere model output.  The files are "self-describing" and the metadata contained in the files conform with the NetCDF Climate and Forecast (CF) Metadata Conventions (see http://www.cgd.ucar.edu/cms/eaton/cf-metadata).  The CF conventions for netCDF data generalize and extend the Cooperative Ocean / Atmosphere Research Data Service (COARDS) conventions developed in the 1990s. Note that the CF convention establishes standard names for climate and weather variables, which identify the physical quantity. These standard names are given in the tables below. Note that more than one field can be associated with the same standard name because different fields sampled in different ways (e.g., surface air temperature vs. upper air temperature) refer to the same physical quantity.  Nevertheless, one can uniquely identify each stored field by considering additional metadata stored in the file (e.g., dimension information). Extended definitions of CF standard names, which basically answer the question "What do you mean precisely by this quantity?", may be found on the Web at http://www.cgd.ucar.edu/cms/eaton/cf-metadata/standard_name.html.

The standards for CFMIP contributions closely follow those established for IPCC/CMIP, which are given in the document http://www-pcmdi.llnl.gov/ipcc/IPCC_output_requirements.htm. This document should be read carefully before preparing contributions. Perhaps the easiest way to meet these requirements is to rewrite your model output through CMOR, a software library available from PCMDI and further described in the next paragraph. Put briefly, the requirements specify required metadata and file organization, along with acceptable coordinate systems.. Units and sign conventions of the data must conform to the tables below. Latitude-longitude grids must be rectilinear, i.e., have a unique set of longitudes that applies to all latitudes. Data on non-rectilinear grids must be interpolated to rectilinear grids before transmission to the PCMDI. With the exception of certain cloud-related fields (which as noted below should be provided on model levels), three-dimensional atmosphere variables must be interpolated to standard pressure levels given below.  Each output file should contain only a single output variable, though there may be many time points per file. This file organization contrasts with the typical model output history files, which contain all variables for a single time step.

To facilitate adherence to these standards, a PCMDI team has written (in FORTRAN 90) a standard output code called CMOR (pronounced "see more"), which is now available as open source. This code structures the data uniformly and writes netCDF files in full compliance with IPCC and CFMIP requirements. Use of CMOR is being encouraged (and in some cases required) by various other ongoing model intercomparison projects. A cmor_readme file (see http://www-pcmdi.llnl.gov/software/cmor/cmor_readme.html) and code documentation in pdf format (see http://www-pcmdi.llnl.gov/software/cmor/cmor_users_guide.pdf) are available. For further information, contact taylor13@llnl.gov.

The IPCC/CMIP requirements as defined by http://www-pcmdi.llnl.gov/ipcc/IPCC_output_requirements.htm should be satisfied by new submissions of CFMIP data except that in that document:

  1. All references to the IPCC Standard Output should refer instead to the CFMIP Standard Output (i.e., this document).
  2. The required global attribute "project_id" should be set to "CFMIP"
  3. As appropriate, the "experiment_id" should be set to one of the following:

'slab ocean control experiment'

'2xCO2 equilibrium experiment'

'+2K perpetual July experiment'

'0K perpetual July experiment'

'-2K perpetual July experiment'

4.                   The CFMIP files should be stored in the following directory structure:  A separate directory should be created containing all the variables in each model/experiment/standard_output_table category (e.g., CSIRO/Slabcntl/A1).  This will ensure that filenames within each directory are unique.  For uniformity, it is suggested that the trailing lower case letter, found on the tables below, should be dropped (e.g., 'CF1', not CF1a) and it is suggested that the experiments be abbreviated as follows:

Slabcntl (i.e., the slab ocean control experiment)

2xCO2 (i.e., the 2xCO2 equilibrium experiment)

Plus2KJul (i.e., the +2K perpetual July experiment)

CntlJul (i.e., the 0K perpetual July experiment)

Minus2KJul (i.e., the -2K perpetual July experiment)

 

5.                  For perpetual July experiments the calendar should be set to 'none', no matter what calendar is assumed by your model.  The basetime should be set to "days since 1-7-15".  This is the CF-convention method of indicating perpetual July.  The time coordinate values will normally record the number of days since the beginning of the experiment.  For the time mean, the bounds specify the beginning and end of the interval over which the time mean is formed.

6.                  For slab ocean monthly mean data, climatological values should be reported (averages of several January's, several February's, etc.).  The CF-convention specifies that time bounds information should be stored in an array identified by the "climatology" attribute (http://www.cgd.ucar.edu/cms/eaton/cf-metadata/CF-1.0.html#climatology), and the "cell_methods" attribute (http://www.cgd.ucar.edu/cms/eaton/cf-metadata/CF-1.0.html#cell-methods) provides further information concerning how the climatological values are computed.  For CFMIP, the CF conventions should be followed, but the complete set of climatological metadata accommodated by the conventions is not yet produced by CMOR.  It is therefore currently acceptable to omit the "climatology" attribute and simply record the cell bounds as if they were monthly means (and CMOR will indicate in the cell_methods attribute that a time-mean is stored, but not that it is a climatological mean).  In the future CMOR will be modified to accommodate climatological data.

The notes that appear in the tables below are meant to provide precise definitions of the requested fields.  Sometimes it may be impossible to satisfy the requests; in these cases, any deviations from the specifications below should be described in the "history" and/or "comment" attributes associated with the variable.


Model output summary

Various fields, sampled at various frequencies, are needed to carry out various proposed CFMIP analyses (see http://cfmip.metoffice.com/DR.html).  The tables below derive from tables (again see http://cfmip.metoffice.com/DR.html) that were originally modified from the AMIP standard output tables (see http://www-pcmdi.llnl.gov/projects/amip/OUTPUT/AMIP2/index.html) to meet the needs of CFMIP.  These "original" CFMIP tables have been restructured to be more consistent with IPCC standard output and to be consistent with the CMOR input tables that have been created for CFMIP.  The "analysis packages" indicated by letters in the fifth and sixth columns of each table refer to various proposed CFMIP analyses.  A few fields are not called for by any of the currently proposed projects, but have been requested as part of the IPCC standard model output, so they are included here because they may eventually be of some interest.

Tables CF1, CF2, and CF3 contain fields reported at monthly, annual, and daily frequencies, respectively.  Table CF4 contains doubled CO2 forcing fields. Tables CF1, CF2, and CF3 are further separated below into component tables, depending largely on whether the fields are spatially two or three dimensional and whether the three dimensional fields are stored on pressure levels or model levels.  Finally, a separate component table contains ISCCP simulator output.

CFMIP output is requested for different time periods, depending on the reporting frequency and the experiment.  In the case of "monthly mean" data:

  • For slab ocean experiments 12 x 20 months of monthly mean data should be reported after equilibrium is reached.
  • For perpetual-July experiments 24 months of monthly mean data (or longer if there is high variability) should be reported after equilibrium is reached.
  • For slab ocean runs, annual mean data should be reported for each year simulated (including spin-up).  No annual means are requested from prescribed SST experiments. 
  • For slab ocean experiments,  report daily values only  for the final 5 years at equilibrium.
  • For perpetual-July experiments, report daily values for the final 2 years at equilibrium.  

As noted below some of the tables parallel the IPCC standard output tables, but have been assigned a different table number.  Also a few minor differences between these tables and the original CFMIP tables exist. The current lists are preferable, however submission of data based on old CFMIP diagnostic lists are acceptable.


Tables of output fields

Monthly data and time-independent fields

For slab ocean experiments, 20*12 months of monthly mean data should be reported after equilibrium is reached.  If this is not possible, a monthly mean climatology is acceptable, although this will preclude much of the analysis in subproject 3. For perpetual-July experiments,  24 monthly means should be reported,  again after equilibrium is reached. Failing this, a climatological mean would be acceptable.

Table CF1a: monthly 2-d atmosphere or surface data (longitude, latitude, time).  The first 44 entries in this table reproduce IPCC Table A1a.  Sea ice concentration (sic, entry 48) appears in in IPCC Table O1c.

 

CF standard_name

output variable name

units

±2K  analysis pack-ages

Slab analysis pack-ages

notes

1

air_pressure_at_sea_level

psl

Pa

A, Z

E, Z

 

2

precipitation_flux

pr

kg m-2 s-1

A, Z

E, Z

includes both liquid and solid phases.

3

air_temperature

tas

K

A, D, Z

E, H, Z

near-surface (usually, 2 meter) air temperature; the CMOR singleton dimension default value of 2 m can be overridden, if absolutely necessary, by redefining axis "height1".

4

Moisture_content_of_soil_layer

mrsos

kg m-2

Z

Z

water in all phases in the upper 0.1 meters of soil, and averaged over the land portion of the grid cell (i.e., compute by dividing the total mass of water contained in the soil layer of the grid cell by the land area in the grid cell); report as "missing" or 0.0 where the land fraction is 0;  the CMOR singleton dimension default value of 0.1 m can be overridden, if absolutely necessary, by redefining axis "depth1".

5

soil_moisture_content

mrso

kg m-2

Z

E, Z

water in all phases summed over all soil layers, and averaged over the land portion of the grid cell (i.e., compute by dividing the total mass of water contained in the soil layer of the grid cell by the land area in the grid cell); report as "missing" or 0.0 where the land fraction is 0.

6

surface_downward_eastward_stress

tauu

Pa

Z

Z

 

7

surface_downward_northward_stress

tauv

Pa

Z

Z

 

8

surface_snow_thickness_where_snow

snd

m

 Z

this thickness when multiplied by the average area of the grid cell covered by snow yields the time-mean snow volume.  Thus, for time means, compute as the weighted sum of thickness (averaged over the snow-covered portion of the grid cell) divided by the sum of the weights, with the weights equal to the area covered by snow.  report as 0.0 in snow-free regions.

9

surface_upward_latent_heat_flux

hfls

W m-2

A,D,L,Z

E, L, Z

 

10

surface_upward_sensible_heat_flux

hfss

W m-2

A,D,L,Z

E,L,Z

 

11

surface_downwelling_longwave_flux_in_air

rlds

W m-2

A,D,L,Z

E,L,Z

 

12

surface_upwelling_longwave_flux_in_air

rlus

W m-2

A,D,L,Z

E,L,Z

 

13

surface_downwelling_shortwave_flux_in_air

rsds

W m-2

A,D,L,Z

E,L,Z

 

14

surface_upwelling_shortwave_flux_in_air

rsus

W m-2

A,D,L,Z

E,L,Z

 

15

surface_temperature

ts

K

A, D, Z

E, H, Z

"skin" temperature (i.e., SST for open ocean)

16

surface_air_pressure

ps

Pa

A, Z

E, Z

not mean sea-level pressure

17

snowfall_flux

prsn

kg m-2 s-1

Z

Z

 

18

convective_precipitation_flux

prc

kg m-2 s-1

Z

Z

 

19

atmosphere_water_vapor_content

prw

kg m-2

A, D, Z

E, Z

vertically integrated through the atmospheric column

20

soil_frozen_water_content

mrfso

kg m-2

Z

Z

summed over all soil layers, and averaged over the land portion of the grid cell (i.e., compute by dividing the total mass of frozen water contained in the soil layer of the grid cell by the land area in the grid cell); report as "missing" or 0.0 where the land fraction is 0.

21

surface_runoff_flux_where_land

mrros

kg m-2 s-1

Z

Z

compute as the total surface runoff leaving the land portion of the grid cell divided by the land area in the grid cell; report as "missing" or 0.0 where the land fraction is 0.

22

runoff_flux_where_land

mrro

kg m-2 s-1

Z

Z

compute as the total runoff (including "drainage" through the base of the soil model) leaving the land portion of the grid cell divided by the land area in the grid cell; report as "missing" or 0.0 where the land fraction is 0.

23

surface_snow_amount_where_land

snw

kg m-2

Z

E,Z

compute as the mass of surface snow on the land portion of the grid cell divided by the land area in the grid cell; report as "missing" or 0.0 where the land fraction is 0; exclude snow on vegetation canopy or on sea ice.

24

surface_snow_area_fraction_where_land

snc

%

Z

Z

fraction of grid cell covered by snow that lies on land; exclude snow that lies on sea ice.

25

surface_snow_melt_flux_where_land

snm

kg m-2 s-1

Z

Z

compute as the total surface melt water on the land portion of the grid cell divided by the land area in the grid cell; report as 0.0 for snow-free land regions; report as 0.0 or "missing" where the land fraction is 0.

26

eastward_wind

uas

m s-1

Z

Z

near-surface (usually, 10 meters) eastward component of wind; the CMOR singleton dimension default value of 10 m can be overridden, if absolutely necessary, by redefining axis "height2".

27

northward_wind

vas

m s-1

Z

Z

near-surface (usually, 10 meters) northward component of wind; the CMOR singleton dimension default value of 10 m can be overridden, if absolutely necessary, by redefining axis "height2".

28

specific_humidity

huss

1 (i.e., dimen-sionless fraction)

Z

Z

near-surface (usually, 2meters) specific humidity; the CMOR singleton dimension default value of 2 m can be overridden, if absolutely necessary, by redefining axis "height1".

29

toa_incoming_shortwave_flux

rsdt

W m-2

A, D, Z

E, H, Z

incident shortwave at the top of the atmosphere

30

toa_outgoing_shortwave_flux

rsut

W m-2

A, D, Z

E, H, Z

at the top of the atmosphere

31

toa_outgoing_longwave_flux

rlut

W m-2

A, D, Z

E, H, Z

at the top of the atmosphere (to be compared with satellite measurements)

32

net_downward_radiative_flux_at_ top_of_atmosphere_model

rtmt

W m-2

A, D, Z

E, H, Z

i.e., at the top of that portion of the atmosphere where dynamics are explicitly treated by the model.

33

net_downward_shortwave_flux_in_air

rsntp

W m-2