program bpconv

     ! (c) Dr. Johannes Orphal, 2000
     ! uses cross-section parameters provided by NASA 
     ! (file bp.par)
     
     real*8  lambda, sigma, T, c0, c1, c2
     integer i

     open(10,file='bp.par',status='old')
     open(11,file='bpconv.out')

     write(*,*) 'T in °C = '
     read(*,*) T

     do i=1,1956
        read(10,*) lambda,c0,c1,c2
        sigma = c0 + c1*T + c2 *T *T
        sigma = sigma * 1.0e-20
        write(11,111) lambda, sigma
     end do

     close(10)
     close(11)

 111 format(1x,1f7.3,1x,1e12.5)

     end