@ANNOTC.PRO
close,15
swit=0
READ, swit, PROMPT='Enter Station Switch: 1 for Pechora, 2 for Tura, 3 for Resolute, 4 for Eureka'
MONTH=['          JAN','          FEB','          MAR','         APR','         MAY','         JUN',' ']
if(swit eq 1)then tit_stat='pechora'
if(swit eq 2)then tit_stat='tura'
if(swit eq 3)then tit_stat='resolute'
if(swit eq 4)then tit_stat='eureka'
np=0 & AA=' '
openr,15,'col_o3.'+tit_stat
      readf,15,AA ;+', Latitude '+alat+', Longitude '+alon
      for j=1,3 do readf,15,format='(10x)'
      while ~ EOF(15) do begin
         readf,15,AA
         np=np+1
      endwhile
close,15
timer=0.d0 & colo3_1r=0.d0 & colo3_2r=0.d0 & colo3_3r=0.d0 & col_OMIr=0.d0 & col_sbuvr=0.d0
openr,15,'col_o3.'+tit_stat
      readf,15,tit_stat ;+', Latitude '+alat+', Longitude '+alon
      for j=1,3 do readf,15,format='(10x)'
      time=dblarr(np) & colo3_1=time & colo3_2=time & colo3_3=time & col_omi=time & col_sbuv=time
      for j=0,np-1 do begin
          readf,15,timer,colo3_1r,colo3_2r,colo3_3r,col_OMIr,col_sbuvr
          ;print,timer,colo3_1r,colo3_2r,colo3_3r,col_OMIr,col_sbuvr
          time(j)=timer & colo3_1(j)=colo3_1r & colo3_2(j)=colo3_2r & colo3_3(j)=colo3_3r & col_OMI(j)=col_OMIr & col_sbuv(j)=col_sbuvr
      endfor
close,15
nlen=strlen(tit_stat)
ja=0 & while(strmid(tit_stat,ja,1) eq ' ')do ja=ja+1
ja1=ja & while(strmid(tit_stat,ja1,1) ne ',')do ja1=ja1+1
;print,nlen,ja,ja1,strmid(tit_stat,ja,ja1-ja)
;stop
AAM=strmid(tit_stat,ja,ja1-ja)
nyy=1 & nzz=2
!P.MULTI=[0,nyy,nzz]
!X.TICKS=6
!P.BACKGROUND=-1
!P.COLOR=1
dev=2
Invert=1
Encaps=1
IF (DEV EQ 0) THEN BEGIN
   DEVICE, RETAIN=2
   WINDOW, 0, XSIZE=850, YSIZE=800, RETAIN=2
ENDIF
IF (DEV EQ 1) THEN BEGIN
   SET_PLOT,'PS'
   !P.FONT = -1
   if(Encaps eq 1)then begin
      DEVICE, /ENCAPSULATED
      EPSfile=AAM+'.eps'
      DEVICE, FILE=EPSFile
      DEVICE, /PORTRAIT, XSIZE=20, YSIZE=22, YOFFSET=8, XOFFSET=1
   endif else begin
      PSfile=AAM+'.ps'
      DEVICE, FILE=PSFile
      DEVICE, /PORTRAIT, XSIZE=20, YSIZE=22, YOFFSET=2, XOFFSET=1
   endelse
   DEVICE, BITS_PER_PIXEL=8
ENDIF
IF (DEV EQ 2) THEN BEGIN
   SET_PLOT,'PS'
   phaser, /portrait
   !P.FONT = -1
   if(Encaps eq 1)then begin
      DEVICE, /ENCAPSULATED
      EPSfile=AAM+'.eps'
      DEVICE, FILE=EPSFile
      DEVICE, /PORTRAIT, XSIZE=20, YSIZE=22, YOFFSET=8, XOFFSET=3
   endif else begin
      PSfile=AAM+'.ps'
      DEVICE, FILE=PSFile
      DEVICE, /PORTRAIT, XSIZE=20, YSIZE=22, YOFFSET=2, XOFFSET=2
   endelse
   DEVICE, BITS_PER_PIXEL=8
ENDIF
LOADCT,26
    check,colo3_1 ;& check,col_sbuv(0:ich-1)
    check,colo3_2 ;& check,col_sbuv(0:ich-1)
    check,colo3_3 ;& check,col_sbuv(0:ich-1)
    yminv=fltarr(5) & ymaxv=yminv
    yminv(0)=min(colo3_1) & ymaxv(0)=max(colo3_1)
    yminv(1)=min(colo3_2) & ymaxv(1)=max(colo3_2)
    yminv(2)=min(colo3_3) & ymaxv(2)=max(colo3_3)
jsbuv=0
dcol_sbuv=fltarr(365)
col_sbuv_av=col_sbuv & time_sbuv_av=time
for jm=0,180-1 do if(col_sbuv(jm) gt 0.)then begin
    col_sbuv_av(jsbuv)=col_sbuv(jm)
    time_sbuv_av(jsbuv)=time(jm)
    dcol_sbuv(jsbuv)=colo3_3(jm)-col_sbuv(jm)
    jsbuv=jsbuv+1
endif
yminv(4)=min(col_sbuv_av(0:jsbuv-1)) & ymaxv(4)=max(col_sbuv_av(0:jsbuv-1))
check,col_sbuv_av(0:jsbuv-1)
jomi=0
col_omi_av=col_omi & time_omi_av=time
for jm=0,180-1 do if(col_omi(jm) gt 0.)then begin
    col_omi_av(jomi)=col_omi(jm)
    time_omi_av(jomi)=time(jm)
    jomi=jomi+1
endif
yminv(3)=min(col_omi_av(0:jomi-1)) & ymaxv(4)=max(col_omi_av(0:jomi-1))
check,col_omi_av(0:jomi-1)
    ymin=fix(min(yminv)*0.1)*10 & ymax=fix(max(ymaxv)*0.1+1)*10
    ymax=500
    xtck=year_fin-year_beg & xb=year_fin-0.38*xtck & xf=xb+0.1*xtck
    plot,time_omi_av(0:jomi-1),col_OMI_av(0:jomi-1),xrange=[year_beg,year_fin],xstyle=1,XMARGIN=[8,1],  $
    title=strmid(tit_stat,ja-1,nlen-ja+1),charthick=THK,  $
    CHARSIZE=1.5,xticks=6,yrange=[ymin,ymax],ytitle='!3Dobson Units',     $
    ymargin=[2,2],ystyle=1,XTICKNAME = MONTH
    oplot,time_omi_av(0:jomi-1),col_OMI_av(0:jomi-1),linestyle=0,col=2,thick=THK
    ;oplot,time_ftir,oco_ftir,linestyle=0,col=100,thick=THK
    ;oplot,time_sbuv_av(0:jsbuv-1),col_sbuv(0:jsbuv-1),linestyle=0,col=100,thick=THK
    oplot,time(0:np-1),colo3_1(0:np-1),linestyle=0,col=36,thick=THK
    oplot,time(0:np-1),colo3_2(0:np-1),linestyle=0,col=72,thick=THK
    oplot,time(0:np-1),colo3_3(0:np-1),linestyle=0,col=150,thick=THK
    ;oplot,time(0:np-1),col_sm(0:np-1),linestyle=2,col=150,thick=THK
    ;oplot,time(0:np-1),colo3_4(0:np-1),linestyle=0,col=250,thick=THK
    oplot,time_sbuv_av(0:jsbuv-1),col_sbuv_av(0:jsbuv-1),linestyle=0,col=100,thick=THK
    dy=(ymax-ymin)*0.045 & TCOR=2
    ANNOTC,xb,ymax-2*dy,xf,fn1,0,36,THK-TCOR
    ANNOTC,xb,ymax-3*dy,xf,fn2,0,72,THK-TCOR
    ANNOTC,xb,ymax-4*dy,xf,fn3,0,150,THK-TCOR
    ;ANNOTC,xb,ymax-4*dy,xf,fn4,0,250,THK-TCOR
    ;ANNOTC,xb,ymax-4*dy,xf,fn3,0,150,THK
    ;ANNOTC,xb,ymax-5*dy,xf,fn4,0,250,THK
    ANNOTC,xb,ymax-5*dy,xf,'!3OMI',0,2,THK-TCOR
    ANNOTC,xb,ymax-1*dy,xf,'SBUV 1979-2019',0,100,THK-TCOR
    dy=(ymaxd-ymind)*0.05
    xb=year_fin-0.4*xtck & xf=xb+0.1*xtck
    plot,time(0:np-1),colo3_1(0:np-1)-colo3_2(0:np-1),xrange=[year_beg,year_fin],xstyle=1, $
    XMARGIN=[8,1],title='!3Model difference',charthick=THK,  $
    CHARSIZE=1.5,xticks=6,yrange=[ymind,ymaxd],ytitle='!3Dobson Units',     $
    ymargin=[2,2],ystyle=1,XTICKNAME = MONTH
    ;oplot,time(0:np-1),colo3_1(0:np-1)-colo3_4(0:np-1),linestyle=0,col=50,thick=THK
    oplot,time(0:np-1),colo3_1(0:np-1)-colo3_3(0:np-1),linestyle=0,col=150,thick=THK
    oplot,time(0:np-1),colo3_1(0:np-1)-colo3_2(0:np-1),linestyle=0,col=16,thick=THK
    oplot,time_sbuv_av(0:jsbuv-1),dcol_sbuv(0:jsbuv-1),linestyle=0,col=100,thick=THK
    oplot,time_sbuv_av(0:jsbuv-1),fltarr(jsbuv),linestyle=2,thick=THK-2
    ANNOTC,xb,ymind+1.5*dy,xf,fn3+' vs ',0,100,THK-TCOR
    ANNOTC,xf,ymind+0.5*dy,xf,'SBUV 1979-2019',0,300,THK-TCOR
    ANNOTC,xb,ymind+2.55*dy,xf,fn1+' vs '+fn2,0,16,THK-TCOR
;    ANNOTC,xb,ymind+3*dy,xf,fn1+' vs '+fn4,0,50,THK-TCOR
    ANNOTC,xb,ymind+3.5*dy,xf,fn1+' vs '+fn3,0,150,THK-TCOR
;endfor
if(dev ne 0)then DEVICE, /CLOSE
end