Integer ::j
       real :: T_ps(51), N_ps(51)
       open (41, file='inputs_T_n.dat')
       do j=1,51
       read (41,*) T_ps(j), N_ps(j)
       print*, 'Temp=',T_ps(j)
       print*, 'Density=',N_ps(j)
       end do
       close(41)
       end program