plot_choice='x'

; read data

if (n_elements(rdflg) eq 0) then rdflg=0
if (rdflg eq 0) then begin

dir='/gen0s/hida/fmt/1999/990216/ha/'

file0=file_list(dir,'*.024*ha.*')
file=file0
file0=file_list(dir,'*.025*ha.*')
file=[file,file0]
;file0=file_list(dir,'*.030*ha.*')
;file=[file,file0]
rd_fmt,file,indexfma,datafma
index2map,indexfma,datafma,map
mapfma=rot_map(map,roll=0)

dir='/gen0s/hida/fmt/1999/990216/hm/'

file0=file_list(dir,'*.024*hm.*')
file=file0
file0=file_list(dir,'*.025*hm.*')
file=[file,file0]
;file0=file_list(dir,'*.030*hm.*')
;file=[file,file0]
rd_fmt,file,indexfmm,datafmm
index2map,indexfmm,datafmm,map
mapfmm=rot_map(map,roll=0)

dir='/gen0s/hida/fmt/1999/990216/hp/'

file0=file_list(dir,'*.024*hp.*')
file=file0
file0=file_list(dir,'*.025*hp.*')
file=[file,file0]
;file0=file_list(dir,'*.030*hp.*')
;file=[file,file0]
rd_fmt,file,indexfmp,datafmp
index2map,indexfmp,datafmp,map
mapfmp=rot_map(map,roll=0)


rdflg=1
endif

; process data
if (n_elements(prflg) eq 0) then prflg=0
if (prflg eq 0) then begin
prflg=1
endif

; choice of device
!x.margin=[5.,2.]
!y.margin=[2.,2.]
charsize=2
if (    (plot_choice eq 'ps')  $
     or (plot_choice eq 'cps') $
     or (plot_choice eq 'eps') $
     or (plot_choice eq 'ceps')) then begin
  set_plot,'ps'
  !p.font=0
  encapsulated=0
  if ((plot_choice eq 'cps') or (plot_choice eq 'ceps')) then $
    device,/color,bits_per_pix=5
  if ((plot_choice eq 'eps') or (plot_choice eq 'ceps')) then encapsulated=1
  device,encapsulated=encapsulated
  device,/portrait
  full_size_a4,encapsulated=encapsulated,/portrait ;,yfactor=0.6
  !p.charsize=0.8
  !x.margin=[5.,2.]
  !y.margin=[2.,2.]
   if ((plot_choice eq 'ps') or (plot_choice eq 'eps')) then begin
   endif
endif else begin
  set_plot,'x'
endelse

; plot

mx=n_elements(mapfma)

movie=bytarr(480,480,mx)

if 1 then begin
xrange=[100,400] & yrange=[-480,-180]

x0=xrange(0) & y0=yrange(0)
dx=xrange(1)-xrange(0) & dy=yrange(1)-yrange(0)

color_max=!d.n_colors-4
yellow=10 & green=8 & blue=5 & red=14


for m=0,mx-1 do begin
mfm=m

plot_map,mapfma(mfm),/log $
  ,xrange=xrange,yrange=yrange $
  ,bottom=17,/notitle,/nolabels


plot_map,mapfmm(mfm),/over,lcolor=blue,levels=[80,100,120],lthick=3
plot_map,mapfmp(mfm),/over,levels=[140,160],lcolor=red,lthick=3

   infor='FMT '+strmid(indexfma(mfm).date_obs,11,8)
   xyouts,x0+dx*0.05,y0+dy*0.03,infor,charsize=charsize,color=color_max,thick=2

img=tvrd(0,0)
movie(*,*,m)=img

endfor


endif


; back to default

resetty

if (!d.name ne 'X') then begin
   device,/close
   set_plot,'x'
endif

if (plot_choice eq 'png') then begin
   img=tvrd(0,0)
   mkpng,img,'.'
endif

end