;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 prof

poc=19  ;zadat pocet profilu
ran1=08*3600+20*60+20 ;082020-zacatek casoveho intervalu 
ran2=08*3600+21*60+20 ;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
window,0,xsize=600,ysize=600
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

y=make_array(s,value=ips)



;-----srovnani plot_3Dbox a surface---------
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_3Dbox,x,y,z,/noerase, $
           yrange=[ran1,ran2],xrange=[-3,2],zrange=[0,1.2], $
           ytickformat='yticks', $
           xtitle='wavelength',ytitle='time',ztitle='intensity', $
           title='PROFILES', charsize=2.5, $
           xmargin=[4,4], ymargin=[4,4], zmargin=[4,4]
           az=15
endfor



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

p=100
pole=congrid(pole,p,p)

surface,pole,/noerase, $
        zrange=[0,1.2],yrange=[0,p], xrange=[0,p], $
        xcharsize=0, ycharsize=0,zcharsize=2.5, $
        xmargin=[10,10],ymargin=[10,10], $
         xmargin=[15,15],ymargin=[15,15], $
         ;xticks=1,yticks=1
         az=15       
       


;a=tvrd()
;write_gif,'int.gif',a     ;doplnit jmeno.gif



end