原子番号及び原子量を与えるモジュール
原子番号 20 までしか与えていないので注意
real(8) :: atmwt(20) !原子量 character(2) :: atmsymbol(20) !原子番号
integer(8), parameter :: amax = 20 !原子番号の最大値
原子記号, 原子量を与える
subroutine chemperm_init
atmsymbol = (/"H ", "He", & & "Li", "Be", "B ", "C ", "N ", "O ", "F ", "Ne", & & "Na", "Mg", "Al", "Si", "P ", "S ", "Cl", "Ar", & & "K ", "Ca"/)
atmwt = (/1.00794d-3, 4.002602d-3, & & 6.941d-3, 9.012128d-3, 10.811d-3, 12.0107d-3, & & 14.00674d-3, 15.9994d-3, 18.9984032d-3, 20.1797d-3, & & 22.989770d-3, 24.3050d-3, 26.981538d-3, 28.0855d-3, & & 30.973761d-3, 32.066d-3, 35.4527d-3, 39.948d-3, & & 39.0983d-3, 40.078d-3/)