program read_CS_G15
! read the GEISA-2015 x-sections database
! compilation without specific options : ifort read_geisa_CS_G15.f90 -o read_geisa_CS_G15
! run : read_geisa_CS_G15 file_geisa_x-sections1, file_geisa_x-sections2, ...
!
integer nbr_pts
real wavenb, intensite
real T,P,resol,wavenb_min,wavenb_max
character*1 a1
character*15 reference
character*20 molecule
character*250 file_geisa_CS
nbfich=iargc()
do ifich=1,nbfich
   call getarg(ifich,file_geisa_CS)
   write(*,*) 'reading '//TRIM(file_geisa_CS)//' ...'
   open(10,file=file_geisa_CS,form='formatted',status='old')
   read(10,*)
   read(10,*)
   read(10,*)
   read(10,*)
   read(10,*)
   read(10,*) a1,molecule,T,P,resol,reference,nbr_pts,wavenb_min,wavenb_max
   write(*,'(A,1x,f6.2,1x,f9.2,1x,f8.4,A,i9,2f11.4)') molecule,T,P,resol,reference,nbr_pts,wavenb_min,wavenb_max
   read(10,*)
   read(10,*)
   read(10,*)
   do icpt=1,nbr_pts
      read (unit=10, fmt=*, IOSTAT=ios) wavenb, intensite
      if(icpt.eq.1) write(*,*) 'first point',wavenb, intensite
      if(icpt.eq.nbr_pts) write(*,*) 'last point',wavenb, intensite
   enddo
enddo

end