;nastaveni formatu casove osy

function yticks,axis,index,value

   hour=long(value)/3600
   minute=long(value-3600 * hour) / 60
   sec=value mod 60
   
   return,string(hour, minute, sec, $
          format="(i2.2,':',i2.2,':',i2.2)")
end

;-----------------------------------------------------


pro prof5

poc=19  ;zadat pocet profilu

ran11=08 & ran12=20 & ran13=20  ;zadat hodiny, minuty, vteriny
ran21=08 & ran22=21 & ran23=20

ran1=ran11*3600+ran12*60+ran13 ;082020-zacatek casoveho intervalu 
ran2=ran21*3600+ran22*60+ran23 ;082120-konec casoveho intervalu


for i=1,poc do begin

case i of

1: begin 
   a='/home1/pkotrc/jura/pol/1.fts'
   ips='082023'
   end
2: begin 
   a='/home1/pkotrc/jura/pol/2.fts'
   ips='082037'
   end
3:begin 
   a='/home1/pkotrc/jura/pol/3.fts'
   ips='082038'
   end
4:begin 
   a='/home1/pkotrc/jura/pol/4.fts'
   ips='082046'
   end
5: begin 
   a='/home1/pkotrc/jura/pol/5.fts'
   ips='082050'
   end
6: begin 
   a='/home1/pkotrc/jura/pol/51.fts'
   ips='082051'
   end
7: begin 
   a='/home1/pkotrc/jura/pol/6.fts'
   ips='082053'
   end
8: begin 
   a='/home1/pkotrc/jura/pol/54.fts'
   ips='082054'
   end
9: begin 
   a='/home1/pkotrc/jura/pol/55.fts'
   ips='082055'
   end
10: begin 
   a='/home1/pkotrc/jura/pol/7.fts'
   ips='082057'
   end
11: begin 
   a='/home1/pkotrc/jura/pol/58.fts'
   ips='082058'
   end
12: begin 
   a='/home1/pkotrc/jura/pol/8.fts'
   ips='082100'
   end
13: begin 
   a='/home1/pkotrc/jura/pol/9.fts'
   ips='082103'
   end
14: begin 
   a='/home1/pkotrc/jura/pol/10.fts'
   ips='082104'
   end
15: begin 
   a='/home1/pkotrc/jura/pol/11.fts' 
   ips='082105'
   end
16: begin 
   a='/home1/pkotrc/jura/pol/12.fts'
   ips='082107'
   end
17: begin 
   a='/home1/pkotrc/jura/pol/13.fts'
   ips='082110'
   end
18: begin 
   a='/home1/pkotrc/jura/pol/14.fts'
   ips='082113'
   end
19: begin 
   a='/home1/pkotrc/jura/pol/15.fts'
   ips='082114'
   end

endcase

;------------------------------------------

b=readfits(a)


x=b(*,0)-6562.808
z=b(*,1)


if i eq 1 then begin
   okno=700
   window,0,xsize=okno,ysize=okno

   s=size(z)
   s=s(1)        ;velikost z
   pole=make_array(s,ran2-ran1)  ;pole pouzite v surface
endif


ips1=strmid(ips,0,2) & ips2=strmid(ips,2,2) & ips3=strmid(ips,4,2)
ips=ips1*3600+ips2*60+ips3   ;cas pro dany profil
y=make_array(s,value=ips)

proc1=100*(x(0)-floor(x(0)))/5
proc2=100*(ceil(x(s-1))-x(s-1))/5
zac=proc1*s/100 & kon=proc2*s/100
zpol=congrid(z,s-kon-zac)
pole(0:zac-1,ips-ran1)=0 & pole(s-kon:s-1,ips-ran1)=0
prvni=0
for j=zac,s-kon-1 do begin
pole(j,ips-ran1)=zpol(prvni)
prvni=prvni+1
endfor
;-----------------------
plot,z,y,/noerase, $
     ytickformat='yticks', yrange=[ran1,ran2],xrange=[0,1.2], $
     position=[0.75,0.2,0.95,0.8],ytitle='time',xtitle='intensity', $
     charsize=1.5,xticks=3,xstyle=1,yticks=3,ystyle=1
     

endfor

;-------------------------------
q=(ran2-ran1)/3           ;3 - hodnota u yticks 

cas=[ran1,ran1+q,ran1+2*q,ran2] ;y hodnoty
lambda=findgen(6)-3          ;x hodnoty


p1=0.17 & p2=0.2 & p3=0.55 & p4=0.8  ;pozice grafu v okne

sizex=(p3-p1)*okno 
sizey=(p4-p2)*okno
posx=p1*okno
posy=p2*okno


pole=congrid(pole,sizex,sizey)
pole=bytscl(pole)
tv,pole,posx,posy


x=make_array(5,value=0) & y=make_array(5,value=0)

plot,x,y,ytickformat='yticks', background=0,color=255, /noerase,$
      yrange=[ran1,ran2],xrange=[-3,2], $
      yticks=3,ytickv=cas,xticks=5,xtickv=lambda, $
      position=[p1,p2,p3,p4], $
      charsize=1.5,xtitle='wavelength',ytitle='time'

xyouts,300,600,'PROFILES',color=255,charsize=3,/device

a=tvrd()
write_gif,'5.gif',a


end