pro make_tropical_sf6_time_series dir = '/Volumes/Data_Gator/Data/' ;dir = '/Users/eray/Work/Data/' dyear = 1./12 ;0.0125 nt = 85.7/dyear years = findgen(nt)*dyear + 1939. + dyear/2. print,years[-1] ; Read in MOZART SF6. id = ncdf_open(dir+'Surface_Trace_Gas/SF6/RIGBY_2010_SF6_MOLE_FRACTION_1970_2008.nc') ncdf_varget,id,'longitude',lon ncdf_varget,id,'latitude',lat ncdf_varget,id,'level',level ncdf_varget,id,'time',time ncdf_varget,id,'SF6',sf6_mozart sf6_mozart_tropical_tpause = fltarr(n_elements(time)) ii = where(lat ge -10 and lat le 10) for t = 0, n_elements(time)-1 do begin sf6_mozart_tropical_tpause[t] = mean(sf6_mozart[*,ii,8,t]) endfor ; Extend the time series back 30 years to 1939. ;time = [indgen(31)+1939,time] ;sf6_mozart_tropical_tpause = [0.035 * exp((indgen(31))*(findgen(31)*.0015+0.02)),sf6_mozart_tropical_tpause] ;; Extend the time series forward 4 years to 2012. ;start_val = sf6_mozart_tropical_tpause[n_elements(time)-1] ;time = [time,indgen(4)+2009] ;time += 0.5 ;sf6_mozart_tropical_tpause = [sf6_mozart_tropical_tpause,start_val*exp((indgen(4)+1)*[0.04,0.042,0.044,0.046])] ; Read in the CMIP6 surface time series. dat = read_ascii(dir+'Surface_Trace_Gas/CMIP6_surface_tracer_time_series.txt',data_start=22) cmip_time = reform(dat.field01[0,*]) + 0.5 cmip_sf6 = reform(dat.field01[16,*]) ; Read in Garcia SF6. sf6_gar = read_ascii('~/Work/IDL/TLP_Model/SF6_growth_fit_Garcia_etal.txt',data_start=1) years_sf6_garcia = reform(sf6_gar.field1[0,*]) sf6_garcia = reform(sf6_gar.field1[1,*]) growth_garcia_sf6_ideal = [0.05,0.1,0.165,0.133,0.09,0.077,0.048,0.04] growth_years = [1939,1959,1963,1968,1984,1989,1997,2003,2013] sf6_tpause_garcia_ideal = fltarr(nt) for i = 0, n_elements(growth_garcia_sf6_ideal)-1 do begin ti = where(years ge growth_years[i] and years lt growth_years[i+1],nti) if ti[0] eq 0 then offset = 0.017 else offset = sf6_tpause_garcia_ideal[ti[0]-1] sf6_tpause_garcia_ideal[ti] = offset * exp((indgen(nti)+1)*dyear*growth_garcia_sf6_ideal[i]) endfor ; Read in Engel and Boenisch SF6. sf6_frankfurt = read_ascii(dir+'Surface_Trace_Gas/Unified_SF6_TTM.naf',data_start=53) years_sf6_frnk = reform(sf6_frankfurt.field1[0,*]) sf6_frnk = reform(sf6_frankfurt.field1[1,*]) years_sf6_frnk = [years_sf6_frnk,years[-1]] sf6_frnk = [sf6_frnk,sf6_tpause_garcia_ideal[-1]] ii = interpol(indgen(n_elements(years_sf6_frnk)),years_sf6_frnk,years) sf6_frnk_interp = interpolate(sf6_frnk,ii) ; Read in NOAA SF6. data = read_ascii(dir+'Surface_Trace_Gas/GMD/sf6/combined/GML_global_SF6.txt',data_start=68) sf6_year = reform(data.field01[0,*]) sf6_mon = reform(data.field01[1,*]) sf6_global = reform(data.field01[6,*]) sf6_sh = reform(data.field01[4,*]) sf6_nh = reform(data.field01[2,*]) sf6_cgo = reform(data.field01[26,*]) sf6_mlo = reform(data.field01[22,*]) sf6_smo = reform(data.field01[24,*]) sf6_hats_tropics = (sf6_mlo + sf6_smo) / 2. dt_sf6 = (sf6_hats_tropics - shift(sf6_hats_tropics,2)) / sf6_hats_tropics sf6_years = sf6_year + (sf6_mon-0.5)/12. ii = uniq(sf6_year) nyrs = n_elements(ii) sf6_global_yearly_avg = fltarr(nyrs) & sf6_sh_yearly_avg = sf6_global_yearly_avg & sf6_cgo_yearly_avg = sf6_global_yearly_avg for t = 0, nyrs-2 do begin sf6_global_yearly_avg[t] = mean(sf6_global[where(sf6_year eq sf6_year[ii[t]])]) sf6_sh_yearly_avg[t] = mean(sf6_sh[where(sf6_year eq sf6_year[ii[t]])]) sf6_cgo_yearly_avg[t] = mean(sf6_cgo[where(sf6_year eq sf6_year[ii[t]])]) endfor sf6_years_avg = sf6_year[ii] + 0.5 sf6_global_yearly_avg[-1] = sf6_global_yearly_avg[-2] + 0.3 ; Read in GMD marine boundary layer data for SF6 from https://gml.noaa.gov/ccgg/mbl/. dat = read_ascii(dir+'Surface_Trace_Gas/SF6/SF6_mbl_17S-17N.txt',data_start=72) sf6_mbl_time = reform(dat.field1[3,*]) sf6_mbl_tropics = reform(dat.field1[4,*]) ntm = n_elements(sf6_mbl_time) dat = read_ascii(dir+'Surface_Trace_Gas/SF6/SF6_mbl_30S-30N.txt',data_start=72) sf6_mbl_time2 = reform(dat.field1[3,*]) sf6_mbl_tropics2 = reform(dat.field1[4,*]) ; Read in latitude averaged mbl data from https://gml.noaa.gov/ccgg/mbl/. ny = 14 lat_grid_sf6 = indgen(ny)*10 - 60 sf6_mbl_tropics_lat = fltarr(ntm,ny) & ratio = sf6_mbl_tropics_lat & ratio_seas = fltarr(12,ny) & ratio_seas2 = ratio_seas & ratio_seas_interp = fltarr(1/dyear,ny) & diffs = sf6_mbl_tropics_lat ratio_seas_interp2 = ratio_seas_interp & diffs_avg = fltarr(ny) for y = 0, ny-1 do begin if lat_grid_sf6[y] lt 0 then hem = 'S' else hem = 'N' lat_string = strcompress(string(abs(lat_grid_sf6[y])),/r)+hem dat = read_ascii(dir+'Surface_Trace_Gas/SF6/SF6_mbl_'+lat_string+'.txt',data_start=72) sf6_mbl_time = reform(dat.field1[3,*]) sf6_mbl_tropics_lat[*,y] = reform(dat.field1[4,*]) ; ratio[*,y] = sf6_mbl_tropics_lat[*,y] / sf6_mbl_tropics ; diffs[*,y] = sf6_mbl_tropics_lat[*,y] - sf6_mbl_tropics ratio[*,y] = sf6_mbl_tropics_lat[*,y] / sf6_mbl_tropics2 diffs[*,y] = sf6_mbl_tropics_lat[*,y] - sf6_mbl_tropics2 mon_fracs = sf6_mbl_time - fix(sf6_mbl_time) for t = 0., 11. do begin ti = where(mon_fracs ge t/12. and mon_fracs lt (t+1)/12. and sf6_mbl_time lt 2003,nti) if nti gt 0 then ratio_seas[t,y] = mean(ratio[ti,y]) ti = where(mon_fracs ge t/12. and mon_fracs lt (t+1)/12. and sf6_mbl_time gt 2015,nti) if nti gt 0 then ratio_seas2[t,y] = mean(ratio[ti,y]) endfor diffs_avg[y] = mean(diffs[where(sf6_mbl_time lt 2003),y]) ii = interpol(indgen(12),findgen(12)/12.+1./24,findgen(1/dyear)*dyear+dyear/2.) ratio_seas_interp[*,y] = interpolate(ratio_seas[*,y],ii) ratio_seas_interp2[*,y] = interpolate(ratio_seas2[*,y],ii) endfor ; From Volk et al., 1997, Geller et al., 1997 global average from 1987-1996. time_volk = findgen(150)*1./12. + 1./24. + 1987. sf6_volk = 3.4361 + 0.2376*(time_volk-1996.) + 0.0049*(time_volk-1996.)^2 ; Cape Grim data sf6_cape_grim_years = indgen(31) + 1978.5 sf6_cape_grim = [0.59,0.68,0.78,0.89,0.99,1.11,1.23,1.37,1.53,1.69,1.87,2.05,2.22,2.41,2.6,2.81,3.03,3.26,3.51,3.75,3.98,4.2,4.4,4.6,4.8,5.02,5.25,5.47,5.7,5.95,6.23] sf6_global_cape_grim = 1.05 * sf6_cape_grim sf6_tropics_time_engel = [sf6_cape_grim_years[0:17],sf6_years_avg[1:nyrs-1]] sf6_tropics_engel = [1.03*sf6_global_cape_grim[0:1],1.02*sf6_global_cape_grim[2:15],1.01*sf6_global_cape_grim[16],1.015*sf6_global_cape_grim[17],0.99*sf6_global_yearly_avg[1:2], $ 0.995*sf6_global_yearly_avg[3:4],sf6_global_yearly_avg[5:nyrs-1]] years_sf6_garcia = [years_sf6_garcia,years[nt-1]] sf6_garcia = [sf6_garcia,7.465] ii = interpol(indgen(n_elements(years_sf6_garcia)),years_sf6_garcia,years) sf6_tpause_interp_garcia = interpolate(sf6_garcia,ii) sf6_tpause_interp_garcia = smooth(sf6_tpause_interp_garcia,4/dyear) sf6_tropics_time_engel = [years[0:38/dyear],sf6_tropics_time_engel] sf6_tropics_engel = [0.98*sf6_tpause_interp_garcia[0:30.5/dyear], $ indgen(38/dyear-30.5/dyear+2)*0.03/(38/dyear-30.5/dyear+1)+0.99*sf6_tpause_interp_garcia[30.5/dyear+1:38/dyear],sf6_tropics_engel] sf6_tropics_engel[0:15./dyear] = 0.00505 * exp(0.1*dyear*indgen(15./dyear+1)) ii = interpol(indgen(n_elements(sf6_tropics_time_engel)),sf6_tropics_time_engel,years) sf6_tpause_interp_engel = interpolate(sf6_tropics_engel,ii) ; Put together the best estimate for the tropical tpause SF6 time series. ; 1939 - 1978 : Use scaled WACCM SF6 at the equator and 149 hPa from Garcia et al., 2011. ; 1978 - 1996 : Use scaled Cape Grim SF6. ; 1996 - 2023 : Use average of NOAA MLO and SMO SF6. ;sf6_time_c = [years[0:39/dyear],sf6_cape_grim_years[0:1],years_sf6_frnk[where(years_sf6_frnk gt 1979.7 and years_sf6_frnk lt 1987)],time_volk[0:146],sf6_years[52:-1], $ ; (findgen(9)+1.)*1./12.+sf6_years[-1]] sf6_time_c = [years[0:39/dyear],sf6_cape_grim_years[0:1],years_sf6_frnk[where(years_sf6_frnk gt 1979.7 and years_sf6_frnk lt 1999.5)],sf6_years[58:-1],(findgen(1)+1.)*1./12.+sf6_years[-1]] ;print,sf6_time_c[-1] sf6_tropics = [0.00505 * exp(0.1*dyear*indgen(15./dyear+1)), $ ; From 1939-1954 do a 10% growth rate starting at .00505 ppt. 0.98*sf6_tpause_interp_garcia[15./dyear+1:30.5/dyear], $ ; From 1954-1969 scale the WACCM tpause SF6 by 0.98. indgen(39/dyear-30.5/dyear+1)*0.02/(39/dyear-30.5/dyear+1)+0.98*sf6_tpause_interp_garcia[30.5/dyear+1:39/dyear], $ ; From 1969-1978 linearly go from scaling WACCM by 0.98-1.0 1.05*sf6_cape_grim[0:1], $ ; From 1978-1980 use scaled annual Cape Grim data based on ratio of MLO-SMO avg/Cape Grim after 1996. sf6_frnk[where(years_sf6_frnk gt 1979.7 and years_sf6_frnk lt 1995.7)], $ ; From 1980-1995.6 use scaled monthly Cape Grim data from Engel et al. (2008). sf6_hats_tropics[4:-1]] ; From 1995.7-2012 use the monthly MLO-SMO average from NOAA. sf6_tropics_new = [0.5*sf6_tpause_interp_garcia[0:33/dyear] + 0.5*sf6_frnk_interp[0:33/dyear], $ ; From 1939-1972 take a mixture of the WACCM and Engel values. (indgen(39/dyear-33/dyear+1)*0.5/(39/dyear-33/dyear+1)+0.5)*sf6_tpause_interp_garcia[33/dyear+1:39/dyear] + $ ; From 1972-1978 change the ratio to the full WACCM value. (0.5-indgen(39/dyear-33/dyear+1)*0.5/(39/dyear-33/dyear+1))*sf6_frnk_interp[33/dyear+1:39/dyear], $ 1.05*sf6_cape_grim[0:1], $ ; From 1978-1980 use scaled annual Cape Grim data based on ratio of MLO-SMO avg/Cape Grim after 1996. ; sf6_frnk[where(years_sf6_frnk gt 1979.7 and years_sf6_frnk lt 1987)], $ ; From 1980-1987 use scaled monthly Cape Grim data from Engel et al. (2008). sf6_frnk[where(years_sf6_frnk gt 1979.7 and years_sf6_frnk lt 1999.5)], $ ; From 1980-1999 use scaled monthly Cape Grim data from Engel et al. (2008). ; sf6_volk[0:146], $ ; From 1987-1997 use Volk, Geller global average. ; sf6_hats_tropics[52:-1],1.035*sf6_hats_tropics[-12:-4]] ; From 1997-2023 use the monthly MLO-SMO average from NOAA. sf6_hats_tropics[58:-2],1.04*sf6_hats_tropics[-12:-11]] ; From 1997-2023 use the monthly MLO-SMO average from NOAA. sf6_time_1 = [sf6_time_c[0:39./dyear],sf6_cape_grim_years[0:16],sf6_years] sf6_tropics_1 = [sf6_tropics[0:39./dyear],1.05*sf6_cape_grim[0:16],sf6_hats_tropics] ii = interpol(indgen(n_elements(sf6_time_1)),sf6_time_1,years) sf6_tropics_interp = interpolate(sf6_tropics_1,ii) ii = interpol(indgen(n_elements(sf6_time_c)),sf6_time_c,years) sf6_tropics_interp = interpolate(sf6_tropics,ii) sf6_tropics_new_interp = interpolate(sf6_tropics_new,ii) sf6_tropics_surface = sf6_tropics_interp sf6_tropics_new_surface = sf6_tropics_new_interp sf6_time = years sf6_tropics = sf6_tropics_interp sf6_tropics_new = sf6_tropics_new_interp tii = interpol(indgen(ntm),sf6_mbl_time,sf6_time) ;tropical_sf6 = interpolate(sf6_mbl_tropics,tii) tropical_sf6 = interpolate(sf6_mbl_tropics2,tii) tropical_sf6_sm = smooth(sf6_tropics_new,48,/edge_truncate) deriv_yr = 12*sg_smooth(tropical_sf6_sm, NLEFT=12/2, NRIGHT=12/2-1, DERIV=1, DELTA=1.0, EDGE_PFCAST=1) sf6_gr = 100 * deriv_yr / tropical_sf6_sm sf6_gr_avg = mean(sf6_gr[-120:-12]) t1 = where(sf6_time lt sf6_mbl_time[0]+2) tropical_sf6[t1] = sf6_tropics_new[t1] t2 = where(sf6_time gt sf6_mbl_time[-1],nt2) ;tropical_sf6[t2] = sf6_tropics_new[t2] for t = 0, nt2-1 do tropical_sf6[t2[t]] = (1. + 1e-2*sf6_gr_avg) * tropical_sf6[t2[t]-12] print,sf6_mbl_time[-1],sf6_time[t2],sf6_tropics_new[t2] ;tropical_sf6 = sf6_tropics_new ;print,sf6_years[-12:-1],sf6_hats_tropics[-12:-1] tropical_sf6_lat = fltarr(nt,ny) ti = where(sf6_time lt sf6_mbl_time[0]+2,nti) ti2 = where(sf6_time ge sf6_mbl_time[0]+2,nti2) tii2 = interpol(indgen(ntm),sf6_mbl_time,sf6_time[ti2]) ti3 = where(sf6_time ge sf6_mbl_time[-1],nti3) ;print,sf6_time[ti3] ; Create additional time series that appends Engel time series before 1975. ii = interpol(indgen(n_elements(years_sf6_frnk)),years_sf6_frnk,sf6_time) sf6_frnk_interp = interpolate(sf6_frnk,ii) tif = where(sf6_time lt 1975,nti) sf6_tropics2 = [sf6_frnk_interp[tif],sf6_tropics[tif[-1]+1:-1]] for y = 0, ny-1 do begin tropical_sf6_lat[ti2,y] = interpolate(sf6_mbl_tropics_lat[*,y],tii2) ; for t = 0, nti-1 do tropical_sf6_lat[t,y] = ratio_seas_interp[t mod 1/dyear,y] * sf6_tropics_new[t] tropical_sf6_lat[ti,y] = diffs_avg[y] + tropical_sf6[ti] if nti3 gt 0 then for t = 0, nti3-1 do tropical_sf6_lat[ti3[t],y] = ratio_seas_interp2[t mod 1/dyear,y] * tropical_sf6[ti3[t]] endfor ;print,diffs_avg ; Shift the time series by 1 month to approximate the tropical tropopause value. ;sf6_tropics_noshift = sf6_tropics_new_interp ;sf6_tropics = shift(sf6_tropics_interp,7) ;sf6_tropics_new = shift(sf6_tropics_new_interp,7) ;sf6_tropics[0:6] = sf6_tropics[7] ;sf6_tropics_new[0:6] = sf6_tropics_new[7] ;tropical_sf6 = shift(tropical_sf6,7) ;tropical_sf6[0:6] = tropical_sf6[7] ;for y = 0, ny-1 do begin ; tropical_sf6_lat[*,y] = shift(tropical_sf6_lat[*,y],7) ; tropical_sf6_lat[0:6,y] = tropical_sf6_lat[7,y] ;endfor ; Adjust the 2015-16 data. ;ti = where(fix(years) eq 2015,nti) ;for t = 0, nti-1 do sf6_tropics[ti[t]] = 1.04*sf6_tropics[ti[t]-1./dyear] ;ti = where(fix(years) eq 2016,nti) ;for t = 0, nti-1 do sf6_tropics[ti[t]] = 1.04*sf6_tropics[ti[t]-1./dyear] ;ti = where(fix(years) eq 2015,nti) ;for t = 0, nti-1 do sf6_tropics_new[ti[t]] = 1.04*sf6_tropics_new[ti[t]-1./dyear] ;ti = where(fix(years) eq 2016,nti) ;for t = 0, nti-1 do sf6_tropics_new[ti[t]] = 1.04*sf6_tropics_new[ti[t]-1./dyear] save,sf6_time,tropical_sf6,tropical_sf6_lat,lat_grid_sf6,filename=dir+'Surface_Trace_Gas/SF6_tropical_time_series.sav' ; Write out SF6 in a text file. ;openw,1,dir+'Surface_Trace_Gas/SF6_tropical_time_series.txt' ;printf,1,'Tropical SF6 mixing ratios' ;printf,1,'1939-1972: Mixture of the WACCM (Garcia et al., 2011) and Engel et al. (2009) values.' ;printf,1,'1972-1978: Gradually change the ratio to the full WACCM value.' ;printf,1,'1978-1980: Scaled annual Cape Grim data based on ratio of MLO-SMO avg/Cape Grim after 1996.' ;printf,1,'1980-1995.6: Scaled monthly Cape Grim data from Engel et al. (2009). ;printf,1,'1995.7-2023: Monthly MLO-SMO average from NOAA. ;printf,1 ;printf,1,' Year SF6 (ppt)' ;for t = 2480, n_elements(sf6_time)-1, 4 do printf,1,sf6_time[t],sf6_tropics_noshift[t],format='(2F13.3)' ;close,1 gd = where(sf6_time ge 1945,ngd) ; Write out SF6 in a text file. ;openw,1,dir+'Surface_Trace_Gas/SF6_tropical_time_series_Harald.txt' ;openw,1,dir+'Surface_Trace_Gas/SF6_tropical_surface_time_series.txt' ;printf,1,'Tropical SF6 mixing ratios' ;printf,1,'1945-1972: Mixture of the WACCM (Garcia et al., 2011) and Engel et al. (2009) values.' ;printf,1,'1972-1978: Gradually change the ratio to the full WACCM value.' ;printf,1,'1978-1980: Scaled annual Cape Grim data based on ratio of MLO-SMO avg/Cape Grim after 1996.' ;printf,1,'1980-1995.6: Scaled monthly Cape Grim data from Engel et al. (2009). ;printf,1,'1995.7-2023: MBL data set from NOAA (https://gml.noaa.gov/ccgg/mbl/). ;printf,1 ;printf,1,' Year SF6 (ppt)' ;for t = 0, ngd-1 do printf,1,sf6_time[gd[t]],tropical_sf6[gd[t]],format='(2F13.3)' ;close,1 ; ;openw,1,dir+'Surface_Trace_Gas/SF6_tropical_surface_time_series_engel.txt' ;printf,1,'Tropical SF6 mixing ratios' ;printf,1,'1945-1975: Engel et al. (2009) values.' ;printf,1,'1972-1978: Gradually change the ratio to the full WACCM value.' ;printf,1,'1978-1980: Scaled annual Cape Grim data based on ratio of MLO-SMO avg/Cape Grim after 1996.' ;printf,1,'1980-1995.6: Scaled monthly Cape Grim data from Engel et al. (2009). ;printf,1,'1995.7-2023: MBL data set from NOAA (https://gml.noaa.gov/ccgg/mbl/). ;printf,1 ;printf,1,' Year SF6 (ppt)' ;for t = 0, ngd-1 do printf,1,sf6_time[gd[t]],sf6_tropics2[gd[t]],format='(2F13.3)' ;close,1 colors = ['sky blue','blue','purple','brown','pink','orange','red','lime green','green','tan','gold','medium orchid','lavendar','grey'] ;p = plot(indgen(2),/nodata,xrange=[1990,2025],yrange=[2,12],ytitle='ppt',title='Tropical Tropopause SF$_6$') ;;p = plot(sf6_time,sf6_tropics,thick=3,color='sky blue',/overplot) ;;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;;p = plot(sf6_mbl_time,sf6_mbl_tropics,thick=2,color='orange',/overplot) ;p = plot(sf6_time,tropical_sf6,thick=2,color='magenta',/overplot) ;p = plot(sf6_years,sf6_mlo,color='blue',thick=2,/overplot) ;p = plot(sf6_years,sf6_smo,color='green',thick=2,/overplot) p = plot(indgen(2),/nodata,xrange=[2020,2025],yrange=[10,12.5]) for y = 3, ny-4 do p = plot(sf6_time,tropical_sf6_lat[*,y],color=colors[y],thick=2,/overplot) ;p = plot(time_volk,sf6_volk,color='lime green',thick=2,/overplot) p = plot(sf6_time,tropical_sf6,thick=3,color='red',/overplot) ;p = plot(indgen(2),/nodata,xrange=[1945,1980],yrange=[0,0.9],ytitle='ppt',title='Tropical Average SF$_6$') ;p = plot(cmip_time,cmip_sf6,symbol='o',/sym_filled,color='sky blue',linestyle=6,/overplot) ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_frnk_interp,thick=2,color='green',/overplot) ;p = plot(sf6_time,sf6_tropics2,thick=2,color='lime green',/overplot) ;p.save,dir+'Plots/SF6_tropical_time_series_1945-80.png' ;p = plot(indgen(2),/nodata,xrange=[1990,2005],yrange=[0,6]) ;for y = 0, ny-1 do p = plot(sf6_time,tropical_sf6_lat[*,y],color=colors[y],/overplot) ;p = plot(sf6_time,tropical_sf6,thick=2,color='magenta',/overplot) ;p = plot(indgen(2),/nodata,xrange=[1980,2025]) ;for y = 0, ny-1 do p = plot(sf6_time,tropical_sf6_lat[*,y]-tropical_sf6,color=colors[y],/overplot) ;p = plot(indgen(2),/nodata,xrange=[1980,2025]) ;for y = 0, ny-1 do p = plot(sf6_time,1e2*(tropical_sf6_lat[*,y]-tropical_sf6)/tropical_sf6,color=colors[y],/overplot) ;p = plot(indgen(2),/nodata) ;for y = 0, ny-1 do p = plot(indgen(12)/12.+1./24,ratio_seas[*,y],thick=2,color=colors[y],/overplot) ;for y = 0, ny-1 do p = plot(indgen(1/dyear)*dyear,ratio_seas_interp[*,y],thick=2,color=colors[y],/overplot) ;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red',xrange=[1939,1970],yrange=[0,0.2]) ;p = plot(sf6_time,sf6_tropics,thick=2,color='blue',/overplot) ;p = plot(years,sf6_tpause_garcia_ideal,thick=2,color='orange',/overplot) ;p = plot(years,sf6_frnk_interp,thick=2,color='green',/overplot) ;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red') ;p = plot(sf6_time,sf6_tropics,thick=2,color='blue',/overplot) ;p = plot(years,sf6_tpause_garcia_ideal,thick=2,color='orange',/overplot) ;p = plot(years,sf6_frnk_interp,thick=2,color='green',/overplot) ;p = plot(sf6_years,sf6_mlo,color='blue',thick=2) ;p = plot(sf6_years,sf6_hats_tropics,color='sky blue',thick=2,/overplot) ;p = plot(sf6_years,sf6_smo,color='green',thick=2,/overplot) ;p = plot(sf6_years,sf6_cgo,color='orange',thick=2,/overplot) ;p = plot(sf6_time_1,sf6_tropics_1,color='blue',thick=2,xrange=[1970,1985],yrange=[0.1,2.0]) ;p = plot(sf6_cape_grim_years,sf6_global_cape_grim,color='orange',thick=2,/overplot) ;p = plot(sf6_time,sf6_tropics,color='green',thick=2,/overplot) ;p = plot(sf6_years,dt_sf6,thick=2,color='orange') ;p = plot(sf6_years,sf6_hats_tropics/sf6_cgo,thick=2,color='red',title='Ratio of Tropical (MLO and SMO avg) to Cape Grim NOAA HATS SF$_6$') ;p.save,'~/Work/Talks and Posters/SPARC BDC Workshop 2012/SF6_tropical_to_Cape_Grim_ratio.png' ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='plum') ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',/overplot) ;p = plot(sf6_time,sf6_tropics_surface,thick=3,color='blue',xrange=[1939,1956],yrange=[0,0.03],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',xrange=[1939,1972],yrange=[0,0.25],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=3,color='red',/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',xrange=[1963,1980],yrange=[0,0.8],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=3,color='red',/overplot) ;p = plot(time,sf6_mozart_tropical_tpause,thick=3,color='green',/overplot) ;;p = plot(years,sf6_tpause_interp_engel,thick=2,color='orange',/overplot) ;p = plot(sf6_cape_grim_years,1.05*sf6_cape_grim,thick=3,color='brown',linestyle=2,/overplot) ;;p = plot(sf6_cape_grim_years,1.03*sf6_cape_grim,thick=3,color='khaki',linestyle=2,/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_surface,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.66,'1.05 * Cape Grim',color='brown',/norm) ;t = text(0.17,0.7,'MOZART 10S-10N 132 hPa (Rigby et al. 2010)',color='green',/norm) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;;t = text(0.17,0.7,'Tropical Trop (Engel)',color='orange',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p.save,'~/Work/Talks and Posters/SPARC BDC Workshop 2012/SF6_tropical_tseries_1963-80.png' ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',xrange=[1963,1978],yrange=[0,0.62],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=3,color='red',/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',xrange=[1970,1985],yrange=[0.3,1.5],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=3,color='red',/overplot) ;p = plot(time,sf6_mozart_tropical_tpause,thick=3,color='green',/overplot) ;;p = plot(years,sf6_tpause_interp_engel,thick=2,color='orange',/overplot) ;p = plot(sf6_cape_grim_years,1.05*sf6_cape_grim,thick=3,color='brown',linestyle=2,/overplot) ;;p = plot(sf6_cape_grim_years,1.03*sf6_cape_grim,thick=3,color='khaki',linestyle=2,/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.66,'1.05 * Cape Grim',color='brown',/norm) ;t = text(0.17,0.7,'MOZART 10S-10N 132 hPa (Rigby et al. 2010)',color='green',/norm) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;;t = text(0.17,0.7,'Tropical Trop (Engel)',color='orange',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p.save,'~/Work/Talks and Posters/SPARC BDC Workshop 2012/SF6_tropical_tseries_1975-85.png' ;p = plot(sf6_time,sf6_tropics,thick=2,color='blue',xrange=[1983,2000],yrange=[1.0,4.5],ytitle='SF$_6$ (ppt)',title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red',/overplot) ;;p = plot(years,sf6_tpause_interp_engel,thick=2,color='orange',/overplot) ;p = plot(time,sf6_mozart_tropical_tpause,thick=2,color='green',/overplot) ;p = plot(sf6_cape_grim_years,sf6_global_cape_grim,thick=2,color='brown',linestyle=2,/overplot) ;;p = plot(sf6_cape_grim_years,1.03*sf6_cape_grim,thick=3,color='khaki',linestyle=2,/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.66,'1.05 * Cape Grim',color='brown',/norm) ;t = text(0.17,0.7,'MOZART 10S-10N 132 hPa (Rigby et al. 2010)',color='green',/norm) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;;t = text(0.17,0.7,'Engel Inferred',color='orange',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p.save,'~/Work/Talks and Posters/SPARC BDC Workshop 2012/SF6_tropical_tseries_1983-2000.png' ; ;p = plot(sf6_time,sf6_tropics,thick=2,color='blue',xrange=[1995,2013],yrange=[3.0,7.5],ytitle='SF$_6$ (ppt)', $ ; title='Tropical SF$_6$ Time Series') ;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red',/overplot) ;;p = plot(years,sf6_tpause_interp_engel,thick=2,color='orange',/overplot) ;p = plot(time,sf6_mozart_tropical_tpause,thick=2,color='green',/overplot) ;p = plot(sf6_cape_grim_years,sf6_global_cape_grim,thick=2,color='brown',linestyle=2,/overplot) ;p = plot(sf6_cape_grim_years,1.03*sf6_cape_grim,thick=3,color='khaki',linestyle=2,/overplot) ;p = plot(years_sf6_frnk,sf6_frnk,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics_surface,thick=2,color='sky blue',/overplot) ;t = text(0.17,0.66,'1.05 * Cape Grim',color='brown',/norm) ;t = text(0.17,0.7,'MOZART 10S-10N 132 hPa (Rigby et al. 2010)',color='green',/norm) ;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;;t = text(0.17,0.7,'Tropical Trop (Engel)',color='orange',/norm) ;t = text(0.17,0.82,'ERay (WACCM + Scaled Cape Grim + NOAA)',color='sky blue',/norm) ;t = text(0.7,0.82,'Tpause',color='blue',/norm) ;t = text(0.17,0.78,'Boenisch and Engel (2008)',color='orange',/norm) ;p.save,'~/Work/Talks and Posters/SPARC BDC Workshop 2012/SF6_tropical_tseries_1995-2013.png' ;p = plot(sf6_time,sf6_tropics,thick=2,color='blue',xrange=[1965,2025],yrange=[0,12],ytitle='ppt',title='Tropical Tropopause SF$_6$') ;;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red',/overplot) ;p = plot(years,sf6_tpause_interp_engel,thick=3,color='orange',/overplot) ;p = plot(sf6_time,sf6_tropics,thick=3,color='blue',/overplot) ;p = plot(sf6_time,sf6_tropics_new,thick=2,color='sky blue',/overplot) ;;p = plot(time,sf6_mozart_tropical_tpause,thick=2,color='green',/overplot) ;;p = plot(sf6_cape_grim_years,sf6_global_cape_grim,thick=2,color='brown',linestyle=2,/overplot) ;;t = text(0.17,0.82,'1.05 * Cape Grim',color='brown',/norm) ;;t = text(0.17,0.78,'MOZART 10S-10N 132 hPa (Rigby et al. 2010)',color='green',/norm) ;;t = text(0.17,0.74,'WACCM Eq 149 hPa (Garcia et al. 2011)',color='red',/norm) ;t = text(0.2,0.76,'Engel et al 2009',color='orange',/norm) ;t = text(0.2,0.8,'This study',color='blue',/norm) ;p.save,'~/Dropbox/Work/Journal Articles/Mine/TLP Age paper/SF6_time_series.png' ;p = plot(years,sf6_tpause_garcia_ideal,thick=2,color='indian red') ;p = plot(years,sf6_tpause_interp_garcia,thick=2,color='red',/overplot) end