1 September 2010, eaton

Processing to create 1995-2105 dataset for RCP6.0 transient runs from
decadal avg data provided by Jean-Francois Lamarque.

For the RCP6.0 case I only produced one version of the aerosol and
ozone datasets corresponding to what we called "version 1" in the other RCP
datasets.  This version begins with the years 1995 and 2000 from the 20th
century dataset and appends the new datasets using mid-decade timestamps,
i.e., 2005 through 2105.

1. Copy input datasets to bluefire disk:

% cd /ptmp/eaton/inputdata/rcp60
% msrcp 'mss:/LAMAR/IPCC_AR5/Climatologies_20100806/RCP60/ozone*' .

2. Copy existing rcp26(v1) datasets for ozone.  These will be used to extract
   the 1994, 1995 and 2000 data that will be prepended to the RCP60 data.
   The 1994 data is a copy of 1995 to supply wrap points so the run can start
   on Jan 01 of 1995.

% cp /fis/cgd/cseg/csm/inputdata/atm/cam/ozone/ozone_rcp26_v1_1.9x2.5_L26_1995-2105_c100316.nc .


3. Convert decadal avg files from r8 to r4 data, and remove the variable T.
   This step is done using the cam interpolation routine from
   tools/interpic_new/.  Then concatenate these files...

% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2000-2009.nc /ptmp/eaton/inputdata/rcp60/ozone_2000-2009.nc ozone_2000-2009.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2010-2019.nc /ptmp/eaton/inputdata/rcp60/ozone_2010-2019.nc ozone_2010-2019.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2020-2029.nc /ptmp/eaton/inputdata/rcp60/ozone_2020-2029.nc ozone_2020-2029.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2030-2039.nc /ptmp/eaton/inputdata/rcp60/ozone_2030-2039.nc ozone_2030-2039.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2040-2049.nc /ptmp/eaton/inputdata/rcp60/ozone_2040-2049.nc ozone_2040-2049.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2050-2059.nc /ptmp/eaton/inputdata/rcp60/ozone_2050-2059.nc ozone_2050-2059.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2060-2069.nc /ptmp/eaton/inputdata/rcp60/ozone_2060-2069.nc ozone_2060-2069.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2070-2079.nc /ptmp/eaton/inputdata/rcp60/ozone_2070-2079.nc ozone_2070-2079.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2080-2089.nc /ptmp/eaton/inputdata/rcp60/ozone_2080-2089.nc ozone_2080-2089.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2090-2099.nc /ptmp/eaton/inputdata/rcp60/ozone_2090-2099.nc ozone_2090-2099.r4.nc
% /blhome/eaton/bin/interpic -e T -p 4 -t /ptmp/eaton/inputdata/rcp60/ozone_2100-2109.nc /ptmp/eaton/inputdata/rcp60/ozone_2100-2109.nc ozone_2100-2109.r4.nc

% ncrcat ozone_2000-2009.r4.nc ozone_2010-2019.r4.nc ozone_2020-2029.r4.nc \
ozone_2030-2039.r4.nc ozone_2040-2049.r4.nc ozone_2050-2059.r4.nc \
ozone_2060-2069.r4.nc ozone_2070-2079.r4.nc ozone_2080-2089.r4.nc \
ozone_2090-2099.r4.nc ozone_2100-2109.r4.nc ozone_2000-2109.nc.516970 

4. Strip off 1994, 1995, and 2000 data from the existing rcp.

% ncks -dtime,0,35 ozone_rcp26_v1_1.9x2.5_L26_1995-2105_c100316.nc ozone_1994-2000.nc

5. Prepend 1994-2000 to 2000-2109:

% ncrcat ozone_1994-2000.nc ozone_2000-2109.nc.516970 ozone_1994-2109.nc

6. overwrite the time, date, datesec vars with the time coordinate
   extracted from a previous rcp file.

% ncks --append -v time,date,datesec time-date_1995-2105.nc ozone_1994-2109.nc

7. zonally average the data and use ncview to make sure it looks reasonable...

% ncwa -a lon ozone_1994-2109.nc ozone_1994-2109_za.nc

8. Move to final location:

% mv ozone_1994-2109.nc ozone_rcp60_1.9x2.5_L26_1995-2105_c100901.nc
% cp ozone_rcp60_1.9x2.5_L26_1995-2105_c100901.nc /fis/cgd/cseg/csm/inputdata/atm/cam/ozone