!$Id$
program test_avhrr
use avhrr_read_mdl

   type (avhrr_data_config) :: cnf
   type ( avhrr_data_out ) :: avh

   cnf % chan_on_rfl (1:6) = .true.
   cnf % chan_on_rad (4:6) = .true.
   cnf % chan_on_bt (4:6) = .true.
   
   cnf % file = 'NSS.GHRR.NN.D12122.S1049.E1223.B3579697.GC.gz'
   cnf % dir_1b = '/Users/awalther/Satellite_Input/n18/2012/122/'
   cnf % Ancil_data_dir ='/DATA/Ancil_Data/clavrx_ancil_data/'
   
   
    
  ! cnf % file = 'NSS.GHRR.NG.D88222.S0115.E0309.B0982931.GC.gz'
  ! cnf % dir_1b = '/Users/awalther/Satellite_Input/n10/'
   cnf % Ancil_data_dir ='/DATA/Ancil_Data/clavrx_ancil_data/'
   
   cnf % offset = [1,1]
   cnf % count = [100,400]
   cnf % use_intern_calib_refl = .true.
   cnf % use_intern_calib_thermal = .true.
   call get_avhrr_data ( cnf , avh )

   print*, avh % is_set
   
  ! print*, avh % band (5) % rad (1200:1210,200) 
   print*,'refl:', avh % band (1) % ref (1200:1410,200) 

end program test_avhrr