PhysicalConstantsClass.f90 Source File

module PhysicalConstantsClass
    use iso_fortran_env
    implicit none

    
	type :: PhysicalConstants_
        real(real64) :: c =  299792458.0d0! speed of light in vacuum (m/s)
        real(real64) :: speed_of_light =  299792458.0d0! speed of light in vacuum (m/s)
        real(real64) :: Plank = 6.62607015*10.0d0**(-34.0) ! the Plank constant (J s)
        real(real64) :: g =9.80665d0  ! standard acceleration of gravity (m/s)
        real(real64) :: R =8.31446261815324d0 ! Molar gas constant (J/K/mol)
        real(real64) :: gas_constant = 8.31446261815324d0 ! Molar gas constant (J/K/mol)
        real(real64) :: Avogadro = 6.02214076*10.0d0**23 ! Avogadro constant (/mol)
        real(real64) :: k = 1.380649*10.0d0**(-23) ! Boltzmann constant (J/K)
        real(real64) :: Boltzmann = 1.380649*10.0d0**(-23) ! Boltzmann constant(J/K)
    contains
        procedure,public :: unit => unitPhysicalConstants
    end type
contains

function unitPhysicalConstants(obj,tgt) result(ret)
    class(PhysicalConstants_),intent(in) :: obj
    character(*),intent(in) :: tgt
    character(:),allocatable :: ret


    if(tgt=="c")then
        ret = "m/s"
        return
    endif
    if(tgt=="speed_of_light")then
        ret = "m/s"
        return
    endif
    if(tgt=="Plank")then
        ret = "J s"
        return
    endif
    if(tgt=="g")then
        ret = "m/s"
        return
    endif
    if(tgt=="R")then
        ret = "J/K/mol"
        return
    endif
    if(tgt=="gas_constant")then
        ret = "J/K/mol"
        return
    endif
    if(tgt=="Avogadro")then
        ret = "/mol"
        return
    endif
    if(tgt=="k")then
        ret = "J/K"
        return
    endif
    if(tgt=="Boltzmann")then
        ret = "J/K"
        return
    endif


end function

end module