PRINTSCRIPT; print $script_style; include "/var/www/html/core/partc"; $linkpage = <<< PRINTLINK gfdl homepage > people > cobweb homepage > people > v. balaji homepage > this page PRINTLINK; print $linkpage; // GFDL header include "/var/www/html/core/partd"; $titlepage = <<< TITLEPAGE Data analysis requiring detailed grid information TITLEPAGE; print $titlepage; // GFDL header include_once( '/var/lib/php/counter.inc' ); error_reporting(E_ERROR); require_once('../magpierss/rss_fetch.inc'); require_once('../magpierss/rss_utils.inc'); include "/var/www/html/core/parte"; $pagecontent = <<< ENDCONTENT

First of all, take a look at this 360x200 OM3 grid used by GFDL's IPCC ocean model:

While it is nominally described as a 1-degree ocean model, the resolution near the Equator is much finer. Similarly, knowing that an atmospheric model uses the sigma_coordinate and has 31 vertical levels does not answer the question, "Is the tropopause well-resolved in this model?" If you wanted to plot maps side by side from different models, you would want more precise information about the location of each model grid cell.

So the amount of information you can extract from IPCC documentation tables, as in our first use case, is quite minimal. It's quite remarkable, then, that scientists are carrying out analyses such as this example, from Held and Soden (2006), which shows a map of some quantity1 averaged across 20 different models from the IPCC archive.

How did Held and Soden pull this off? This is again through the magic of metadata standards. Besides the grid description that the producers of the model output have entered into the the IPCC Model Documentation Table, each dataset itself also contains some rather detailed information about each grid cell, recorded following the CF Conventions. For all the models in the IPCC archive,

  • The grid is required to be a lon_lat_grid.
  • Each variable location is associated with a specific geographic longitude and latitude.
  • Each variable location is associated with a cell, specified by a bounding box surrounding the nominal location.

Given that all datasets conform to these stipulations, it turns out to quite straightforward to perform analyses such as Held and Soden (2006). Metadata standards are the cornerstone of the scientific tour de force that is the IPCC archive.

But what if the grid were not a lon_lat_grid, but instead one of the more exotic ones in the grid menagerie?


1 It happens to be δ(P-E), though that's irrelevant to this example.
ENDCONTENT; print $pagecontent; $subj = 'grids'; $url = 'http://www.gfdl.noaa.gov/~vb/weblogs/journal.rdf'; $rss = fetch_rss($url); if( $rss ) { echo "\n"; foreach ($rss->items as $item) { $href = $item['link']; $title = $item['title']; if ( preg_match( "/\b$subj/i", $title ) ) { echo "\n"; } } } $url = 'http://www.gfdl.noaa.gov/~vb/weblogs/journal2007.rdf'; $rss = fetch_rss($url); if( $rss ) { foreach ($rss->items as $item) { $href = $item['link']; $title = $item['title']; if ( preg_match( "/\b$subj/i", $title ) ) { echo "\n"; } } } $pagecontent = <<< ENDCONTENT

emacs-muse-mode created by v. balaji (balaji@princeton.edu) in emacs using the emacs-muse mode.
ENDCONTENT; print $pagecontent; print "last modified: ". date( "d F Y", getlastmod() ); print "
this page visited: ".getCount(). " times "; include "/var/www/html/core/partf"; include "/var/www/html/core/partg";