This is the old phase!
Please see OCMIP Phase 2

1. Introduction

This file provides information to OCMIP participants describing how to
  1. retrieve data files as necessary to force carbon-cycle simulations that are to be run according to standard OCMIP protocol (as initially discussed at the first OCMIP meeting 20 January 1995, Hamburg, Germany), and
  2. send model results to LSCE where the standard analysis of OCMIP will take place.

1.1 File Transfer

All files mentioned below are available by ftp from LSCE. To transfer these files used for the standard OCMIP simulations, do the following:
  1. issue the ftp command,
  2. type the appropriate name and password (as prompted), and
  3. continue by typing only the underlined portions below, with a RETURN after each.
$ ftp
Name: ocmip
Password: xxxxxxxx
ftp> cd /home/geo/ocmip
ftp> prompt
ftp> binary
ftp> mget guide* README* *.f *.dat *.Z
ftp> quit

After transfer is complete, uncompress any files ending with ".Z" with

uncompress *.Z

The resulting invoice of files follows:

Filename Contents
README.satdat Explanation of satellite-derived data (from Boutin and Etcheto)
guidelines.ocmip General guidelines describing OCMIP and short-term plans
funcchem.f Fortran functions to calculate equilibrium constants and the solubility of carbon dioxide in sea water.
gridlab.f Includes documentation for transfer of model grids
trunlab.f Includes documentation for transfer of model results
c14equ.dat ,
Atmospheric Data used for OCMIP forcing
splco2.dat Atmospheric Data used for OCMIP forcing
stab.dat Atmospheric Data used for OCMIP forcing Fractional Sea Ice used to force model Gas transfer coefficient (i.e., solubility*piston velocity) Surface land mask used by Levitus

1.2 Timetable

Participants are asked to send model results for Phase A (alias Anything goes) as quickly as possible. We desperately need your results by the end of May, if the presentation at the First GAIM Science Conference in September is to be a success. In the spirit of Phase A, we ask that if you already have model results, please send them immediately. As for Phase B (where simulations must be forced with the standard model input specified below), that part of OCMIP does not officially start until just after the GAIM conference. Nonetheless, please begin Phase B as soon as possible and send your results when they become available. This will help avoid duplication of analytical efforts.

2. Model input

2.1 Subroutine (Fortran)

2.2 Atmospheric Records

2.3 Surface Fields

K = (b(t) + d * ( <v>+ <<u>2>)) * (Sc(t)/660)**(-0.5)

2.4 File Format

All files above were written as IEEE 32-bit binary; they are given on the standard "Levitus" grid (1 x 1 degree, having the first grid box (1,1) centered at -179.5 and -89.5 (i.e., 179.5W, 89.5S) and the last grid box (360,180) centered at +179.5, +89.5). Below is an example of the Fortran necessary to read OCMIP surface fields.
character*80 filename filename=''
open(unit= 60, file=filename, form='unformatted')
c For Cray only, uncomment following line to signal IEEE binary
c call asnfile(filename,'-N ieee -F f77',ier)
c Read data from January (mo=1) to December (mo=12)
do mo = 1,12
read(60) ((Kgw(i,j,mo), i=1,360), j=1,180)
end do

2.5 Plots of surface fields

PostScript plots of surface fields used for OCMIP can be found in the plots subdirectory. To visualize results, you'll need to first transfer them to your machine (by ftp as described in Section 1.1).
Each plot's filename has 5 fields: the first 3 fields describe the file contents as defined below; the last two fields do not change.
AAA = Data Type
BBB = Data Source
CCC = Grid
ps = PostScript
Z = Compressed file

Data Type:

Note: at the end of each "Data Type" specification, there are two additional pieces of information: Data Source: Grid:

3. Model output

This section is designed to answer the following question:

How do I organize my model results so that they can be analyzed at LSCE?

To facilitate the transfer of model results to LSCE for analysis, we avoid asking participants to produce data files according to a rigid structure. Instead, it appears much simpler if each participant just copies two subroutines (gridlab.f and trunlab.f), making only a minimum of changes (read statements, etc.) in order to pass his results as arguments to these two "standard" OCMIP routines.

A detailed explanation of the necessary changes are given as comments in each subroutine. Your modified versions of these subroutines will be employed at LSCE along with the files containing your model output, in order to transfer your results, iin a standard fashion, to the programs used for OCMIP analysis. After your model output is ready and you have modified the two subroutines mentioned above, please let me know ( ASAP. Then I will create a special directory into which you may transfer your results by ftp. Alternatively, we can transfer your results to LSCE (i.e., if you make them available by anonymous ftp).

As described in the OCMIP guidelines, all participants are requested to send their results on their original grid. For results from models that do NOT use a rectangular grid (i.e., only Hamburg and LODyC models, for the moment), we ask that you also duplicate all your results on a 2 x 2 degree rectangular grid (spanning from -180 to 180 Longitude, -90 to 90 Latitude); along with this additional output, please send modified versions of the data transfer subroutines We regret the extra effort this may require and hope that in the near future this will not be necessary.


Enting, I.G., T. M. L. Wigley, M. Heimann, 1994. Future Emissions and concentrations of carbon dioxide: key ocean / atmosphere / land analyses, CSIRO Aust. Div. Atmos. Res. Tech. Pap. No. 31, 118 pp.

Levitus, S., 1982. Climatological atlas of the World Ocean, NOAA Prof. Pap. 13, U.S. GPO., Washington, D.C., 173 pp. Walsh, J. 1978. A data set on northern hemisphere sea ice extent, 1953-1976. Glaciological Data, World Data Center for Glaciology (Snow and Ice), Report GD-2, 49-51.

Wanninkhof, R., 1992. Relationship between wind speed and gas exchange over the ocean, J. Geophys. Res., 97, 7373-7382 Weiss, R. F., 1974. Carbon dioxide in water and seawater: the solubility of a non-ideal gas, Marine Chem., 2, 203-215.

Zwally, H. J., J. Comiso, C. Parkinson, W. Campbell, F. Carsey, and P. Gloerson, 1983. Antarctic Sea Ice, 1973-1976: Satellite Passive Microwave Observations, NASA, 206 pp.