;
pro pcm_wv
;
ntime = 1188
y_1 = fltarr(ntime)
y_2 = fltarr(ntime)
y_3 = fltarr(ntime)
y_4 = fltarr(ntime)
y_5 = fltarr(ntime)
y_6 = fltarr(ntime)
;
data_1=dblarr(3,ntime)
openr,1,'pcm_B05_12_hpt_0001_monthly.dat'
readf,1,data_1
close,1
;
y_1=data_1[0,*]
y_1=transpose(y_1)
;
y_2=data_1[1,*]
y_2=transpose(y_2)
;
y_3=data_1[2,*]
y_3=transpose(y_3)
;
data_2=dblarr(3,ntime)
openr,1,'pcm_B05_12_hpt_0851_monthly.dat'
readf,1,data_2
close,1
;
y_4=data_2[0,*]
y_4=transpose(y_4)
;
y_5=data_2[1,*]
y_5=transpose(y_5)
;
y_6=data_2[2,*]
y_6=transpose(y_6)
;
dt = 1.0
time = findgen(ntime)*dt   ;*** create the time index
;
wave_1 = wavelet(y_1,dt,mother='paul',period=period_1,coi=coi_1,/verbose,signif=signif_1)
nscale_1 = n_elements(period_1)
;
wave_2 = wavelet(y_2,dt,mother='paul',period=period_2,coi=coi_2,/verbose,signif=signif_2)
nscale_2 = n_elements(period_2)
;
wave_3 = wavelet(y_3,dt,mother='paul',period=period_3,coi=coi_3,/verbose,signif=signif_3)
nscale_3 = n_elements(period_3)
;
wave_4 = wavelet(y_4,dt,mother='paul',period=period_4,coi=coi_4,/verbose,signif=signif_4)
nscale_4 = n_elements(period_4)
;
wave_5 = wavelet(y_5,dt,mother='paul',period=period_5,coi=coi_5,/verbose,signif=signif_5)
nscale_5 = n_elements(period_5)
;
wave_6 = wavelet(y_6,dt,mother='paul',period=period_6,coi=coi_6,/verbose,signif=signif_6)
nscale_6 = n_elements(period_6)
;
set_plot,'ps'
loadct,27
tvlct,redtab,greentab,bluetab,/get
redtab(1)=0
greentab(1)=0
bluetab(1)=0
tvlct,redtab,greentab,bluetab
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_h.eps'
contour,abs(wave_1)^2,time,period_1,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_1),min(period_1)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_1 = rebin(transpose(signif_1),ntime,nscale_1)
contour,abs(wave_1)^2/signif_1,time,period_1,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_1,color=180,thick=6,noclip=0 
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_p.eps'
contour,abs(wave_2)^2,time,period_2,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_2),min(period_2)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_2 = rebin(transpose(signif_2),ntime,nscale_2)
contour,abs(wave_2)^2/signif_2,time,period_2,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_2,color=180,thick=6,noclip=0 
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_t.eps'
contour,abs(wave_3)^2,time,period_3,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_3),min(period_3)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_3 = rebin(transpose(signif_3),ntime,nscale_3)
contour,abs(wave_3)^2/signif_3,time,period_3,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_3,color=180,thick=6,noclip=0 
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_me_h.eps'
contour,abs(wave_4)^2,time,period_4,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_4),min(period_4)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_4 = rebin(transpose(signif_4),ntime,nscale_4)
contour,abs(wave_4)^2/signif_4,time,period_4,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_4,color=180,thick=6,noclip=0 
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_me_p.eps'
contour,abs(wave_5)^2,time,period_5,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_5),min(period_5)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_5 = rebin(transpose(signif_5),ntime,nscale_5)
contour,abs(wave_5)^2/signif_5,time,period_5,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_5,color=180,thick=6,noclip=0 
;
device,/encapsulated,/color,bits=24,font_size=18,/helvetica,/bold,file='images/wv_me_t.eps'
contour,abs(wave_6)^2,time,period_6,color=1,charthick=4,ythick=6,xthick=6,xminor=1,yminor=0,xtickinterval=200,ytickinterval=12,xtitle='time (months)',ytitle='period (months)',yrange=[max(period_6),min(period_6)],nlevels=25,/fill,xstyle=1,ystyle=1,/ytype
signif_6 = rebin(transpose(signif_6),ntime,nscale_6)
contour,abs(wave_6)^2/signif_6,time,period_6,/overplot,color=1,charthick=6,thick=6,level=1.0,c_annot='95%'
plots,time,coi_6,color=180,thick=6,noclip=0 
;
return
end