;plot_CO2_LBC_WACCM.pro
;
; This routine will read in the lower boundary montly global mean valuse for CO2 used in the 1955-2010 run of WACCM-X and plot for 1955-2010 and
; 1970-2000
;
;  Open and read the CO2 from the LBC file
;
fid = ncdf_open('/glade/p/cesm/cseg/inputdata/atm/waccm/lb/LBC_1765-2500_1.9x2.5_CMIP5_RCP45_za_c120204.nc')

ncdf_varget,fid,'CO2_LBC_mean',CO2_LBC_Monthly

ncdf_varget,fid,'date',date
;
;  Plot the monthly CO2 for 1955-2010 and 1970-2000.  Field date goes from 1765 to 2500 monthly.  
;  Jan 1, 1955 is index 2280, Jan 1, 2011 is index 2952 and  Jan 1, 1970 is index 2460
;

!P.THICK = 3
!x.thick = 3
!y.thick = 3

SET_PLOT, 'PS'

psfile = 'CO2_LBC_Monthly_56yr_GlobalMean_WACCMX.ps'

DEVICE, /COLOR, FILENAME = psfile
loadct, 39, ncolors=10

CYears = ['1955','1960','1965','1970','1975','1980','1985','1990','1995','2000','2005','2010']

yTitle = 'CO2 (ppm)'
xTitle = 'Year'
Title = 'CO2 1955-2010 Surface Monthly Global Mean'

plot, CO2_LBC_Monthly[2280:2952]*1.E06,title=Title,xtitle=xTitle,ytitle=yTitle,xstyle=1,ystyle=1,charsize=1.3 ,xticks=11,xtickname=CYears,yrange=[300.,400.],/nodata
oplot, CO2_LBC_Monthly[2280:2952]*1.E06,color=40
;plot, date[2460:29

perCO2Increase = (CO2_LBC_Monthly[2952] - CO2_LBC_Monthly[2460]) / MEAN(CO2_LBC_Monthly[2280:2952]) * 100.
cperCO2Increase = STRTRIM(perCO2Increase,2)
perCO2IncreaseDec = perCO2Increase / 5.6
cperCO2IncreaseDec = STRTRIM(perCO2IncreaseDec,2)


XYOUTS,20,390,'CO2 % Increase 1955-2010 = '+cperCO2Increase+'% ('+cperCO2IncreaseDec+' %/decade)',/data,color=230,charsize=1.1, charthick=3

print, 'CO2 % Increase 1955-2010 ', (CO2_LBC_Monthly[2952] - CO2_LBC_Monthly[2460]) / MEAN(CO2_LBC_Monthly[2280:2952]) * 100.

print, 'CO2 % Increase Per Decade 1955-2010 ', ( (CO2_LBC_Monthly[2952] - CO2_LBC_Monthly[2460]) / MEAN(CO2_LBC_Monthly[2280:2952]) * 100. ) /5.6

device,/close

END