INTERNATIONAL Ocean Carbon Model Intercomparison Project
PHASE 2


MAIN CONTACT James Orr
ACCESS INFORMATION Message from Jim Orr (23 Sep 2002):

You may now access the OCMIP-2 model output archive via the DODS server at IPSL (http://dods.ipsl.jussieu.fr).

Recall that we agreed previously that everyone should have free access to our model output. I would also recommend that the data that we have used to evaluate models should also be released; however, I don't recall if we spoke about this specific issue. If there are objections, particularly from the data specialists, please let me know. Note we spent some time during OCMIP-2² to format the data so that it could be compared to the models.

To see the structure for the hierarchy of the OCMIP-2 model output archive start by pointing your browser at http://dods.ipsl.jussieu.fr/ocmip/phase2

Read further for an introduction to doing your own analysis on output in the OCMIP-2 DODS archive.

Some analysis packages, e.g. ferret and CDAT (http://esg.llnl.gov/cdat/), are DODS compatible. Thus in your analysis scripts, you can access the OCMIP model output archive remotely, thereby eliminating the need to transfer large files to your local machine. For instance, to make a plot of the bathymetry in the NCAR model, just execute the following 2 lines in Ferret

use "http://dods.ipsl.jussieu.fr/cgi-bin/nph-dods/ocmip/grids/NCAR_grid.nc"
shade bathy

For something a little more complicated, see the attached Ferret script (demo_cfc_2.jnl) written by Patrick Brockmann. On a single page, this script plots maps of the surface CFC-11 concentration in January for six models (MPIM, AWI, IPSL, SOC, UL, PRINCE); plotting takes about 40 seconds. For this script to work correctly, you must have already installed Ferret (v5.22) as well as FAST (Ferret Analysis Scripting Toolbox), which is also written by Pätrick. Amongst other things, the FAST package facilitates drawing pretty maps of all the OCMIP-2 models, even those with irregular grids.

To install FAST, go to http://www.ipsl.jussieu.fr/~brocksce/fast/FAST_install.html

Best regards, Jim

Contact James C. Orr
Laboratoire des Sciences du Climat et de l'Environnement
Unite Mixte de Recherche CEA-CNRS
LSCE, CEA Saclay
Bat. 709 - Orme
F-91191 Gif-sur-Yvette Cedex
FRANCE
mailto: orr@cea.fr
Phone: (33) (0)1 69 08 77 23
Fax: (33) (0)1 69 08 77 16
http://www.ipsl.jussieu.fr/~jomce

demo_cfc_2.jnl

go db_set http://dods.ipsl.jussieu.fr/cgi-bin/nph-dods/ocmip

go window_init
go page_size 640

set mode metafile demo_cfc_2.plt

go page_new -i 2 3 0 100 10 90

!----------------------
go logo_put ocmip

!----------------------
let new_a=CFC11[k=1,l=1]*1E+9
def sym qualifiers=/lev=(-1E-4)(1,9,0.5)(1E+4)

!----------------------
set v 1
def sym MODEL=MPIM
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
define view/xlim=0,1/ylim=0.,0.1 colorbar
set v colorbar
go colorbar_put -h 10 90 10 50 0.3 ($qualifiers)

!----------------------
set v 2
def sym MODEL=AWI
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
set v 3
def sym MODEL=IPSL
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
set v 4
def sym MODEL=SOC
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
set v 5
def sym MODEL=UL
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
set v 6
def sym MODEL=PRINCE
use "($DB)/phase2/($MODEL)/CFC/($MODEL)_CFC_1980.nc"
let a = if mask[k=1,d=($MODEL)_grid.nc] eq 1 then new_a
go map_cylin ($MODEL) 0 a 10 ($qualifiers)/title=($MODEL)/nokey earth_continents
cancel data/all

!----------------------
go textlayer_new
go textlayer_text_put 50 98 "CFC11\\nfrom OCMIP\\nDODS database"

!----------------------
cancel mode metafile