36a37
>       stop   ! HERE
57c58
<        Te    = Te_eV/8.617385d-5
---
>        Te    = Te_eV/8.617383d-5
3969d3969
<       integer(4),dimension(10000) ::istatwtup,istatwtlow
4000,4001d3999
<              istatwtlow(icnt)=ig(indx1glob)
<              istatwtup(icnt)=ig(indx2glob)
4019c4017
<        write(47,901)(radat(ii,j), j=1,8),istatwtlow(ii),istatwtup(ii)
---
>        write(47,901)(radat(ii,j), j=1,8)
4024c4022
< 901   FORMAT( 2x,1e11.5,2x,F11.5,2x,3(2x,e11.5),2(2x,f7.1),2x,1e11.5,i4,i4)
---
> 901   FORMAT( 2x,1e11.5,2x,F11.5,2x,3(2x,e11.5),2(2x,f7.1),2x,1e11.5)
4040,4041c4038,4039
<       h=6.6260755e-27 ! ergs*s
<       c=2.99792458d10  ! cm/s
---
>       h=6.626176e-27 ! ergs*s
>       c=2.997925e10  ! cm/s
4043c4041
<       const = c**2/(2*h) != 6.7818560099e+46 cm^2/(ergs*s^3) = 6.7819569724d46
---
>       const = c**2/(2*h) != 6.7818560099e+46 cm^2/(ergs*s^3) 
4054c4052
<           v=e*2.41799e14 ! eV -> Hz
---
>           v=e*2.41797e14 ! eV -> Hz