PROGRAM read_fwd
      USE io_scene1
      IMPLICIT NONE

      INTEGER :: I, iu1, nMeasurData
      TYPE(MeasurData_type)  :: MeasurData
      CHARACTER*300     :: File1
      CHARACTER*300     :: DIR1
    
      iu1=1
!      DIR1='/data/data548/pub/yklee/mirs_20200122/data/TestbedData/DynamicData/fwd_analys/n20_atms/2019-02-15/'
!      File1=TRIM(DIR1)//'FWD_ECMWS_n20_d20190215_t0001520_e0002236_b06436_c20200410184711387322_nobc_ops.h5.HR'
!      File1=TRIM(DIR1)//'FWD_ECMWS_n20_d20190715_t1137173_e1137489_b08571_c20200410184934055027_noac_ops.bin'
      DIR1='/data/smcd8/yklee/mirs_nqs/data/TestbedData/DynamicData/fmsdr/nqs_atms/2023-08-30/'
      File1=TRIM(DIR1)//'FMSDR_TMS_nqs_d20230830_t0224026_e0224316_b29946_c20230830024754368000_oebc_ops.h5.HR'
!      DIR1='/data/smcd11/yklee/mirs_quicksounder/data/TestbedData/DynamicData/fmsdr/QuickSounder_atms/2023-08-30/'
!      File1=TRIM(DIR1)//'FMSDR_TMS_QuickSounder_d20230830_t0224026_e0224316_b29946_c20230830024754368000_oebc_ops.h5.HR'
!      File1=TRIM(DIR1)//'FWD_ECMWS_n20_d20190715_t1137173_e1137489_b08571_c20200410184934055027_noac_ops.h5.HR'

      OPEN(iu1,file=File1,form='unformatted')
      OPEN(11,FILE='quicksounder_clw_test_wvtimes10.dat')

      CALL ReadHdrMeasurmts(iu1,nMeasurData,MeasurData)
      DO I=1,nMeasurData
       CALL ReadMeasurmts(iu1,MeasurData)
       WRITE(11,*) nMeasurData
       WRITE(11,*) MeasurData%nchan
       WRITE(11,*) MeasurData%nPosScan,MeasurData%nScanLines
       WRITE(11,*) MeasurData%nqc
       WRITE(11,*) MeasurData%CentrFreq(1:MeasurData%nchan)
       WRITE(11,*) MeasurData%polar(1:MeasurData%nchan)
       WRITE(11,*) MeasurData%lat,MeasurData%lon,MeasurData%RelAziAngle,MeasurData%SolZenAngle  ! latitude and longitude
       WRITE(11,*) MeasurData%Node,MeasurData%iscanPos,MeasurData%iscanLine,MeasurData%Year,MeasurData%julDay,MeasurData%secs
       WRITE(11,*) MeasurData%angle(1:MeasurData%nchan)
       WRITE(11,*) MeasurData%Tb(1:MeasurData%nchan)   ! This is brightness temperature
       WRITE(11,*) MeasurData%qc(1:MeasurData%nqc)
       WRITE(11,*) MeasurData%secant_view(1:MeasurData%nchan)

      ENDDO

      STOP

END PROGRAM read_fwd