# References:
# ##########
#
# http://usjgofs.whoi.edu/mzweb/smppi/doney3.html
# http://usjgofs.whoi.edu/mzweb/smpdatadocs/moore_becmodel.html

# Refer to data as:
# Moore, J. K., S. C. Doney, and K. Lindsay.; U.S. JGOFS Synthesis & Modeling Project - Data.
# U.S. JGOFS. iPub: March 2005. 'date you accessed the data' http://usjgofs.whoi.edu/las/servlets/datasets [INVALID URL]

# Key reference:
# Moore, J. K., S. C. Doney, and K. Lindsay. 2004. Upper ocean ecosystem dynamics and iron cycling in a
# global three-dimensional model, Global Biogeochem. Cycles, 18, GB4028, doi:10.1029/2004GB002220.

# download the actual data files (beware: 96 MB each!)
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04450101.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440201.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440301.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440401.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440501.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440601.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440701.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440801.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04440901.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04441001.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04441101.nc
wget  http://usjgofs.whoi.edu/data/smp/moore/t.ECO3.33.04441201.nc

# extract actually required variables
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04450101.nc t.ECO3.33.01.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440201.nc t.ECO3.33.02.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440301.nc t.ECO3.33.03.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440401.nc t.ECO3.33.04.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440501.nc t.ECO3.33.05.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440601.nc t.ECO3.33.06.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440701.nc t.ECO3.33.07.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440801.nc t.ECO3.33.08.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04440901.nc t.ECO3.33.09.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04441001.nc t.ECO3.33.10.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04441101.nc t.ECO3.33.11.nc
ncks -v z_t,z_w,TLAT,TLONG,KMT,TAREA,HT,dz,REGION_MASK,TEMP,SALT,NO3,SiO3,O2,DIC,ALK,POC_FLUX_IN,CaCO3_FLUX_IN,SiO2_FLUX_IN,dust_FLUX_IN t.ECO3.33.04441201.nc t.ECO3.33.12.nc

# calculate annual averages
ncea t.ECO3.33.??.nc t.ECO3.33.nc

# discard the time dimension (which has only one entry with a meaningless value).
ncwa -O -a time t.ECO3.33.nc t.ECO3.33.nc


