all_make_image.pro

pro all_make_image, phase=phase

  o=obj_new('cws_stats') ; ,scale=1.0
;  arr_pro = [] ; does not work 
;  arr_pro = ['cmb', 'cph', 'ctt', 'cth', 'ctp', 'cod', 'ref', 'lwp']   ; OK
  arr_pro = ['ctp']
  npro = n_elements(arr_pro)
; single=1
  o->set_group, 'GSF'

  h=12
  m=0
  day   = [13,17,18,22,03]
  month = [ 6, 6, 6, 6, 7]
; phase='water'

  for d = 2, 2 do begin 
;    for h=12,15 do begin
    for h=0,2 do begin
      for m=00,45,15 do begin 
;      for m=00,00,15 do begin 
        print, '=== TIME:', h, m
        o->set_date,2008,month[d],day[d],h,m
        for p = 0, npro-1 do begin 
          product = arr_pro[p]
          o->set_product, product  
          ;o->make_image, product=product, /single
          o->make_image, product=product
        endfor
      endfor
    endfor
 endfor

end