8 use,
intrinsic :: iso_fortran_env, only: dp => real64
10 real(dp),
parameter :: &
11 pi=3.1415926535897932384626433832795028841971693993751058209749445923&
12 &078164062862089986280348253421170679_dp,& !<pi
13 sigmab=5.67037321e-8_dp,&
14 kb=1.380648813e-23_dp,&
17 eulergamma=0.5772156649015328606&
18 &065120900824024310421_dp,&
20 hpc=6.6260695729e-34_dp,&
22 complex(dp),
parameter :: iu=(0.0e0_dp,1.0e0_dp)
23 character(len=80) :: &
25 real(dp),
dimension(4) :: &
26 tc=(/154.58_dp,126.19_dp,304.17_dp,511.7_dp/),&
27 pc=(/5.043e6_dp,3.396e6_dp,7.386e6_dp,45.1e5_dp/),&
28 mg=(/32e-3_dp,28e-3_dp,44e-3_dp,70e-3_dp/),&
29 tb=(/90.19_dp,77.36_dp,194.75_dp,322.4_dp/),&