FUNCTION divgood,numer,denom
  result = 0.0*denom
  good   = where(denom NE 0.0)
  result(good) = numer(good)/denom(good)
return,result
END

PRO ratfilts

;  Read in all flat fields

   fits_read,'f140lp_a.fits',f140lp_a,hf140lp_a
   fits_read,'f140lp_b.fits',f140lp_b,hf140lp_b
   fits_read,'f140lp_b1.fits',f140lp_b1,hf140lp_b1
   fits_read,'f140lp_b2.fits',f140lp_b2,hf140lp_b2
   fits_read,'f140lp_b3.fits',f140lp_b3,hf140lp_b3
   fits_read,'f165lp.fits',f165lp,hf165lp
   fits_read,'f122m.fits',f122m,hf122m
   fits_read,'f150lp.fits',f150lp,hf150lp
   fits_read,'f150lp_coron.fits',f150lp_coron,hf150lp_coron
   fits_read,'f150lp_coron2.fits',f150lp_coron2,hf150lp_coron2
   fits_read,'f125lp.fits',f125lp,hf125lp
   fits_read,'f115lp.fits',f115lp,hf115lp
   fits_read,'pr130l.fits',pr130l,hpr130l
   fits_read,'pr110l.fits',pr110l,hpr110l
   fits_read,'krsplat.fits',krsplat,hkrsplat
   fits_read,'f115lp_s.fits',f115lps,hf115lps
   fits_read,'f115lp_s1.fits',f115lps1,hf115lps1
   fits_read,'f115lp_s2.fits',f115lps2,hf115lps2
   fits_read,'f115lp_s3.fits',f115lps3,hf115lps3
   fits_read,'f115lp_d2.fits',f115lpd2,hf115lpd2

;  Make some interesting ratios

   f165lp_f122m   = divgood(f165lp,f122m)
   f165lp_f115lp  = divgood(f165lp,f115lp)
   f122m_f115lp   = divgood(f122m,f115lp)
   rat_coron      = divgood(f150lp_coron2,f150lp)
   coront_first   = divgood(f150lp_coron2,f150lp_coron)
   pr130l_f125lp  = divgood(pr130l,f125lp)
   pr110l_f115lp  = divgood(pr110l,f115lp)
   f122m_f125lp   = divgood(f122m,f125lp)
   f140lpa_f150lp = divgood(f140lp_a,f150lp)
   f165lp_f150lp  = divgood(f165lp,f150lp)
   redblue        = 0.5*(f122m + f150lp)
   redblue_f115lp = divgood(redblue,f115lp)
   f140lpb_a      = divgood(f140lp_b,f140lp_a)
   f140lpb1_a     = divgood(f140lp_b1,f140lp_a)
   f140lpb2_a     = divgood(f140lp_b2,f140lp_a)
   f140lpb3_a     = divgood(f140lp_b3,f140lp_a)
   krsplat_f122m  = divgood(krsplat,f122m)
   f115lps_mma    = divgood(f115lps,f115lp)
   f115lps1_s2    = divgood(f115lps1,f115lps2)
   f115lps1_s3    = divgood(f115lps1,f115lps3)
   f115lps2_s3    = divgood(f115lps2,f115lps3)
   f115lpd2_s     = divgood(f115lpd2,f115lps)
   f115lpd2_mma   = divgood(f115lpd2,f115lp)

;  write results

   fits_write,'f165lp_f122m.fits',f165lp_f122m,hf165lp
   fits_write,'f165lp_f115lp.fits',f165lp_f115lp,hf165lp
   fits_write,'f122m_f115lp.fits',f122m_f115lp,hf122m
   fits_write,'rat_coron.fits',rat_coron,hf150lp_coron
   fits_write,'coront_first.fits',coront_first,hf150lp_coron2
   fits_write,'pr130l_f125lp.fits',pr130l_f125lp,hpr130l
   fits_write,'pr110l_f115lp.fits',pr110l_f115lp,hpr110l
   fits_write,'f122m_f125lp.fits',f122m_f125lp,hf122m
   fits_write,'f140lpa_f150lp.fits',f140lpa_f150lp,hf140lpa
   fits_write,'f165lp_f150lp.fits',f165lp_f150lp,hf165lp
   fits_write,'redblue_f115lp.fits',redblue_f115lp,hf115lp
   fits_write,'f140lpb_a.fits',f140lpb_a,hf140lp_b
   fits_write,'f140lpb1_a.fits',f140lpb1_a,hf140lp_b1
   fits_write,'f140lpb2_a.fits',f140lpb2_a,hf140lp_b2
   fits_write,'f140lpb3_a.fits',f140lpb3_a,hf140lp_b3
   fits_write,'krsplat_f22m.fits',krsplat_f122m,hkrsplat
   fits_write,'f115lps_mma.fits',f115lps_mma,hf115lps
   fits_write,'f115lps1_s2.fits',f115lps1_s2,hf115lps1
   fits_write,'f115lps1_s3.fits',f115lps1_s3,hf115lps1
   fits_write,'f115lps2_s3.fits',f115lps2_s3,hf115lps2
   fits_write,'f115lpd2_s.fits',f115lpd2_s,hf115lpd2
   fits_write,'f115lpd2_mma.fits',f115lpd2_mma,hf115lpd2

;  Look at results
   cube = [[[f165lp_f122m]], [[f165lp_f115lp]], [[f122m_f115lp]], $
           [[rat_coron]], [[coront_first]], [[pr130l_f125lp]], $
           [[pr110l_f115lp]], [[f122m_f125lp]], $
           [[f140lpa_f150lp]], [[f165lp_f150lp]], [[redblue_f115lp]], $
           [[f140lpb_a]], [[f140lpb1_a]], [[f140lpb2_a]], [[f140lpb3_a]], $
           [[krsplat_f122m]], [[f115lps_mma]], [[f115lps1_s2]], $
           [[f115lps1_s3]], [[f115lps2_s3]], [[f115lpd2_s]], $
           [[f115lpd2_mma]]]
   lab = ['f165lp_f122m', $
          'f165lp_f115lp', $ 
          'f122m_f115lp', $
          'rat_coron', $
          'coront_first', $
          'pr130l_f125lp', $
          'pr110l_f115lp', $
          'f122m_f125lp', $
          'f140lpa_f150lp', $
          'f165lp_f150lp', $
          'redblue_f115lp', $
          'f140lpb_a', $
          'f140lpb1_a', $
          'f140lpb2_a', $
          'f140lpb3_a', $
          'krsplat_f122m', $
          'f115lps_mma', $
          'f115lps1_s2', $
          'f115lps1_s3', $
          'f115lps2_s3', $
          'f115lpd2_s', $
          'f115lpd2_mma']
   mview,cube,LABELS=lab

END