You can download the WOD files from the NODC website through an online query (http://www.nodc.noaa.gov/OC5/SELECT/dbsearch/dbsearch.html), you can also download the entire database if you look hard enough for the links (http://www.nodc.noaa.gov/OC5/WOD09/data09geo.html). I didn't feel like pointing and clicking myself to death so I wrote this script to download the entire thing. I also added a line in the script to compile a summary database of the CTD and XBT data for use in svp_tool.
#!/bin/bash
count=0
mkdir -p {OSD,MBT,CTD,XBT,PFL,MRB,DRB,APB,UOR,GLD}/OBS
mkdir -p {OSD,MBT,CTD,XBT,PFL,MRB,DRB,APB,UOR,GLD}/STD
echo "Downloading the WOD database"
for quadrant in 1 3 5 7
do
for latitude in 0 1 2 3 4 5 6 7 8
do
for longitude in 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17
do
let count++
for type in OSD MBT CTD XBT PFL MRB DRB APB UOR GLD
do
# Get the "observed" data
file=$type\O$quadrant$latitude$longitude.gz
wget http://data.nodc.noaa.gov/woa/WOD09/GEOGRAPHIC/$type/OBS/$file
mv $file $type/OBS
# Get the "standard level" data (observed is filtered and then
# interpolated to standard depth levels)
file=$type\S$quadrant$latitude$longitude.gz
wget http://data.nodc.noaa.gov/woa/WOD09/GEOGRAPHIC/$type/STD/$file
mv $file $type/STD
done
done
done
done
echo got total $count
# Now add a tmp directory for use by the WOD library for unzipping database files
# during extraction
mkdir tmp
# Now build database of CTD and XBT data
echo "Building a CTD/XBT summary database"
wod2db3 -out wod.db3 CTD/OBS/CTDO*gz XBT/OBS/XBTO*gz
echocount=0
mkdir -p {OSD,MBT,CTD,XBT,PFL,MRB,DRB,APB,UOR,GLD}/OBS
mkdir -p {OSD,MBT,CTD,XBT,PFL,MRB,DRB,APB,UOR,GLD}/STD
echo "Downloading the WOD database"
for quadrant in 1 3 5 7
do
for latitude in 0 1 2 3 4 5 6 7 8
do
for longitude in 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17
do
let count++
for type in OSD MBT CTD XBT PFL MRB DRB APB UOR GLD
do
# Get the "observed" data
file=$type\O$quadrant$latitude$longitude.gz
wget http://data.nodc.noaa.gov/woa/WOD09/GEOGRAPHIC/$type/OBS/$file
mv $file $type/OBS
# Get the "standard level" data (observed is filtered and then
# interpolated to standard depth levels)
file=$type\S$quadrant$latitude$longitude.gz
wget http://data.nodc.noaa.gov/woa/WOD09/GEOGRAPHIC/$type/STD/$file
mv $file $type/STD
done
done
done
done
echo got total $count
# Now add a tmp directory for use by the WOD library for unzipping database files
# during extraction
mkdir tmp
# Now build database of CTD and XBT data
echo "Building a CTD/XBT summary database"
wod2db3 -out wod.db3 CTD/OBS/CTDO*gz XBT/OBS/XBTO*gz
echo
echo "Add WOD_PATH="`pwd` "as an environment variable"
echo "Add WOD_DATABASE_NAME="`pwd`"/"wod.db3 "as an environment variable"
exit
No comments:
Post a Comment