im=readfits('m911226',h) 
 plot_nar,'26-dec-1991'