!---------------------------------------------------------------------
!     Copyright (C) GFD Dennou Club, 2004. All rights reserved.
!---------------------------------------------------------------------
                                                                 !=begin
!= Module ChemPrm
!
!   * Developer: SUGIYAMA Ko-ichiro (sugiyama@gfd-dennou.org)
!   * Version: $Id: chemprm.f90,v 1.2 2005/01/28 12:41:21 sugiyama Exp $ 
!   * Tag Name: $Name:  $
!   * Change History: 
!
!== Overview 
!
!原子番号及び原子量を与えるモジュール 
!
!== Error Handling
!
!== Known Bugs
!
!== Note
!
!原子番号 20 までしか与えていないので注意
!
!== Future Plans
!
                                                                 !=end

module chemprm
  save
                                                                 !=begin
  !== Public Interface
  real(8)                    :: atmwt(20)       !原子量
  character(2)               :: atmsymbol(20)   !原子番号

  !== Public Data
  integer(8), parameter      :: amax = 20       !原子番号の最大値
                                                                 !=end  
contains
                                                                 !=begin
  !== Procedure Interface 
  !
  !=== Initialize module
  !
  !原子記号, 原子量を与える
  !
  subroutine chemperm_init

    !=== 原子記号
    atmsymbol = (/"H ", "He", &
         & "Li", "Be", "B ", "C ", "N ", "O ", "F ", "Ne", &
         & "Na", "Mg", "Al", "Si", "P ", "S ", "Cl", "Ar", &
         & "K ", "Ca"/)
    
    !=== 原子量 (kg/mol). 値は化学便覧改訂 4 版(丸善)より
    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/) 
                                                                 !=end    
  end subroutine chemperm_init
  
end module chemprm
