8 use,
intrinsic :: iso_fortran_env, only: dp => real64
10 real(dp),
parameter :: pi=3.141592653589793238462_dp
11 real(dp),
parameter :: sigmab=5.67037321e-8_dp
12 real(dp),
parameter :: kb=1.380648813e-23_dp
13 real(dp),
parameter :: rg=8.31446218_dp
14 real(dp),
parameter :: na=rg/kb