program read_transitions_geisa2003
!PROGRAM FOR READ ASCII FORMAT GEISA 2003 DATABASE
implicit none
! FIELDS GEISA
     real(kind=8)                       :: waveNb                   ! A
     real(kind=8)                       :: intensity                ! B
     real                               :: collHW                   ! C
     real                               :: energyLowerLevel         ! D
     character(len=80)                  :: quant                    ! E
     real                               :: temperatureDependence    ! F
     integer                            :: isotopeId                ! G
     integer                            :: moleculeId               ! I
     character(len=3)                   :: internalGeisaCode        ! J

     integer(kind=2)                    :: hit_moleculeId           ! K
     integer(kind=2)                    :: hit_isotopeId            ! L
     real                               :: transition_probability   ! M
     real                               :: hit_selfBroadHW          ! N
     real                               :: hit_airBroadHW           ! O
     integer(kind=2)                    :: hit_HAI                  ! P
     integer                            :: hit_indLookUp            ! Q

     real                               :: temperatureDependenceAPS ! R
     real(kind=8)                       :: accWaveNb                ! A'
     real(kind=8)                       :: accIntensity             ! B'
     real                               :: accCollHW                ! C'
     real                               :: accTemperatureDependence ! F'
     real                               :: accAirBroadHW            ! O'
     real                               :: accTempDepAPS            ! R'
     real                               :: accSelfBroadHW           ! N'
     real                               :: temperatureDependenceSBH ! S
     real                               :: accTempDepSBH            ! S'
     real                               :: selfPressShift           ! T
     real                               :: accSelfPressShift        ! T'
     real                               :: temperatureDependenceSPS ! U
     real                               :: accTempDepSPS            ! U'
! OTHER DECLARATION
     integer                            :: n_arg,iargc,nbr_transitions
     integer                            :: unit_input_file
     character(len=250)                 :: prog,ascii_file_geisa

!
      unit_input_file=10
      nbr_transitions=0
!
      n_arg=iargc()
      call getarg(0,prog)
      if(n_arg.lt.1) then
         write(*,*)'Syntax error: '//trim(prog)// &
     & ' ascii format GEISA 2003 file '
         stop
      endif
      call getarg(1,ascii_file_geisa)
      open(unit_input_file, file=trim(ascii_file_geisa),form='formatted')
     
 10   read(unit_input_file,100,end=901) waveNb,intensity,collHW,energyLowerLevel,quant,temperatureDependence, &
                         isotopeId,moleculeId,internalGeisaCode,hit_moleculeId, &
                         hit_isotopeId,transition_probability,hit_selfBroadHW, &
                         hit_airBroadHW,hit_HAI,hit_indLookUp,temperatureDependenceAPS, &
                         accWaveNb,accIntensity,accCollHW,accTemperatureDependence, &
                         accAirBroadHW,accTempDepAPS,accSelfBroadHW,temperatureDependenceSBH, &
                         accTempDepSBH,selfPressShift,accSelfPressShift, &
                         temperatureDependenceSPS,accTempDepSPS

      nbr_transitions=nbr_transitions+1
      write(*,100) waveNb,intensity,collHW,energyLowerLevel,quant,temperatureDependence,   &
                         isotopeId,moleculeId,internalGeisaCode,hit_moleculeId, &
                         hit_isotopeId,transition_probability,hit_selfBroadHW, &
                         hit_airBroadHW,hit_HAI,hit_indLookUp,temperatureDependenceAPS, &
                         accWaveNb,accIntensity,accCollHW,accTemperatureDependence, &
                         accAirBroadHW,accTempDepAPS,accSelfBroadHW,temperatureDependenceSBH, &
                         accTempDepSBH,selfPressShift,accSelfPressShift, &
                         temperatureDependenceSPS,accTempDepSPS
! GOTO THE NEXT TRANSITION
      goto 10

  100 format(f12.6,1P,d11.4,0P,f6.4,f10.4,a36,f4.2,i3,i3,a3,i2,i1,1P,e10.3,0P,f5.4,f8.6,i3,i6, &
        & f6.4,f10.6,1P,d11.4,0P,f6.4,f4.2,f8.6,f6.4,f5.4,f4.2,f4.2,f8.6,f8.6,f4.2,f4.2)
  901 write(*,*) "End of ascii format GEISA 2003 file ", trim(ascii_file_geisa)
      write(*,*) "There are ",nbr_transitions," transitions in this file"

end