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" 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 |
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