PROGRAM read_scene
      USE io_dep
      IMPLICIT NONE

      INTEGER :: nprf1, nprf2, I, iu1, iu2, K
      TYPE(dep_type)  :: dep1,dep2
      CHARACTER*300     :: File1,FILE2,fileout_lat,fileout_lon,fileout_tpw,fileout_sfctyp
    
      iu1=1
      iu2=2
!ykl      File1='NWP_ECMWS_n21_d20200201_t0001253_e0001569_b11416_c20200201000808262491_noac_ops.h5.HR'
!ykl      File2='../../n20_atms/2020-02-01/NWP_ECMWS_n20_d20200201_t0001253_e0001569_b11416_c20200201000808262491_noac_ops.h5.HR'
      File1='./2021-08-17/DEP_TMS_n20_d20210817_t0927306_e0928023_b19409_c20210817095735741992_oeac_ops.h5.HR.ORB'
      fileout_lat='lat.bin'
      fileout_lon='lon.bin'
      fileout_tpw='tpw.bin'
      fileout_sfctyp='sfctyp.bin'
      OPEN(21,FILE=fileout_lat,access='direct',recl=4)
      OPEN(22,FILE=fileout_lon,access='direct',recl=4)
      OPEN(23,FILE=fileout_tpw,access='direct',recl=4)
      OPEN(24,FILE=fileout_sfctyp,access='direct',recl=4)

      OPEN(iu1,file=File1,form='unformatted')
!ykl      OPEN(iu2,file=File2,form='unformatted')
      OPEN(11,FILE='n21.dat')
!ykl      OPEN(12,FILE='n20.dat')

      CALL ReadHdrDep(iu1,dep1,nprf1)
!ykl      CALL ReadHdrdep(iu2,dep2,nprf2)
      K=0
      DO I=1,nprf1
       CALL ReadDep(iu1,dep1)
!ykl       CALL Readdep(iu2,dep2)

       !---Surface-level paremeters
       WRITE(11,*) dep1%iTypAtm,dep1%TPW,dep1%Clw,dep1%RWP,DEP1%SWP,      &
         DEP1%IWP,DEP1%GWP,DEP1%RR,DEP1%SFR,DEP1%CldTop,DEP1%CldBase, &
         DEP1%CldThick,DEP1%PrecipType,DEP1%RainFlag,DEP1%LWP
       !---QC variables
       WRITE(11,*) dep1%qc(1:4)
       !---Positioning variables
       WRITE(11,*) dep1%lat,dep1%lon,dep1%node,dep1%scanUTC,dep1%scanYear,dep1%scanDay,dep1%iscanPos,dep1%iscanLine, &
                   DEP1%angle,DEP1%RelAziAngle,DEP1%SolZenAngle

       K=K+1
       WRITE(21,rec=K) dep1%lat
       WRITE(22,rec=K) dep1%lon
       WRITE(23,rec=K) dep1%TPW
       WRITE(24,rec=K) REAL(dep1%iTypAtm)

       GOTO   2001

       !---Surface-level paremeters
       WRITE(12,*) dep2%iTypAtm,dep2%TPW,dep2%Clw,dep2%RWP,DEP2%SWP,      &
         DEP2%IWP,DEP2%GWP,DEP2%RR,DEP2%SFR,DEP2%CldTop,DEP2%CldBase, &
         DEP2%CldThick,DEP2%PrecipType,DEP2%RainFlag,DEP2%LWP
       !---QC variables
       WRITE(12,*) dep2%qc(1:4)
       !---Positioning variables
       WRITE(12,*) dep2%lat,dep2%lon,dep2%node,dep2%scanUTC,dep2%scanYear,dep2%scanDay,dep2%iscanPos,dep2%iscanLine, &
                   DEP2%angle,DEP2%RelAziAngle,DEP2%SolZenAngle


2001  CONTINUE
      ENDDO

      STOP

END PROGRAM read_scene