Previous Next Contents

2. Model runs

2.1 Gas exchange flux

For simulations of CFC-11 and CFC-12, we will directly model the finite air-sea flux F. In other words, surface CFC concentrations will NOT be set equal to temperature-derived equilibrium values determined from the solubility. Modelers must use the formulation for the standard OCMIP-2 air-to-sea flux,

(1a)   F = Kw (Csat - Csurf)


(1b)   Csat = alpha * pCFC * P/Po


All right hand terms, except P and Po, in equations (1a) and (1b) are different for CFC-11 and CFC-12, as well as other tracers.

2.2 The Piston Velocity Kw

For simulations of CFC-11 and CFC-12, modelers must use the standard OCMIP-2 formulation for the piston velocity Kw. The monthly climatology of Kw, to be interpolated linearly in time by each modeling group, is computed with the following equation adapted from Wanninkhof (1992, eq. 3):

(2)    Kw = (1 - Fice) [Xconv * a *(u2 + v)] (Sc/660)**-1/2


Practically speaking, to use equation (2) each group will interpolate the OCMIP-2 standard information to their own model grid. The standard information is provided by IPSL/LSCE as a monthly climatology on the 1 x 1 degree grid of Levitus (1982) in netCDF format (in file Gridded variables in that file include

For the variables Fice and xKw, continents on the 1 x 1 degree standard grid have been flooded with adjacent ocean values. Such an approach avoids discontinuities at land-sea boundaries during interpolation. See the Fortran program rgasx_ocmip2.f for an example of how to read the information in cfc1112.atm into your interpolation routines. After compilation, to link and use rgasx_ocmip2.f, one must have already installed netCDF.

The file may also be inspected with software that uses netCDF format, such as ncdump or Ferret. Ferret will be used for some of the analysis during OCMIP-2. We encourage participants to become familiar with Ferret now

After installation, one can visualize maps of the standard information in, by using the Ferret script vgasx_ocmip2.jnl.

After launching Ferret, simply issue the following command (at Ferret's "yes?" prompt)

yes? go vgasx_ocmip2.jnl

2.3 Oceanic and Atmospheric Components

Apart from Kw, there are two other terms in equation (1a). The ocean component Csurf [in mol/m^3] is computed by the model each timestep; the atmospheric component Csat is specified a priori via the three remaining terms:

  1. alpha: The CFC solubility alpha is to be computed using modeled SST and SSS, both of which vary in time at each grid point. For OCMIP-2 we use the solubility formulation determined by Warner and Weiss (1985, Table 5 for solubility in [mol/(l * atm)]). The function sol_cfc.f determines alpha accordingly, for both CFC-11 and CFC-12, but changes the units to [mol/(m^3 * picoatm)] so that model CFC concentrations can then be carried in SI units [mol/m^3].
  2. pCFC: Extrapolated records for observed CFC-11 and CFC-12 [in picoatm] constructed at 41S and 45N (Walker et al., pers. comm.). For OCMIP-2, each station will be treated as representative of its own hemisphere, except between 10S and 10N where those station values will be interpolated linearly as a function of latitude. Thus there are 3 zones:
  3. P: Is the total atmospheric pressure [atm] from the monthly mean climatology of Esbensen and Kushnir (1981). The latter, given originally on a 4 x 5 degree grid (latitude x longitude) in bars, is converted to atm by multiplying by (1/1.101325). Land and sea ice values in the original data set were filled with average values from adjacent ocean points. These monthly mean arrays were then linearly interpolated to the 1 x 1 degree grid of Levitus (see netCDF file

Technical notes:

  1. The Fortran subroutine cfc_interp.f interpolates pCFC spatially following the above algorithm. The code allows, in one pass, to spatially interpolate atmospheric pCFC-11 and pCFC-12, at a given timestep, to all grid points as a function of latitude.
  2. The ASCII file cfc1112.atm provides mid-year values of atmospheric pCFC-11 and pCFC-12 [in picoatm] at both stations for the period from 1931 to 1997. See the program read_cfcatm.f.
  3. Temporal interpolation of atmospheric pCFC-11 and pCFC-12 is to be made linearly for each time step, based on mid-year values (file cfc1112.atm).

Previous Next Contents