Making one more minor change in constants.

njoy · Sep 11, 2018 · 5c7c6f5 · 5c7c6f5
1 parent 549452a
commit 5c7c6f5
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 4 deletions.
diff --git a/src/gaminr.f90 b/src/gaminr.f90
@@ -125,8 +125,6 @@ subroutine gaminr
    real(kr),dimension(:,:,:),allocatable::ans
    real(kr),parameter::ekn=12.4e3_kr
    real(kr),parameter::emax=1.e12_kr
-!  real(kr),parameter::epair=.5110034e6_kr
-   real(kr),parameter::epair=amasse*amu*clight*clight/ev
    real(kr),parameter::zero=0
    integer,parameter::ncnt=9
    integer,parameter::nwpff=5000
@@ -1204,7 +1202,6 @@ subroutine gtff(e,enext,idisc,ff,nl,ng,iglo,nq,pff,nwpff)
    real(kr),parameter::c5=20.60744e0_kr
    real(kr),parameter::toler=1.e-6_kr
    real(kr),parameter::rndoff=1.0000001e0_kr
-   real(kr),parameter::epair=2._kr*amasse*amu*clight*clight/ev
    real(kr),parameter::emax=1.e12_kr
    real(kr),parameter::two=2.e0_kr
    real(kr),parameter::half=0.5e0_kr
@@ -1504,7 +1501,7 @@ subroutine gtff(e,enext,idisc,ff,nl,ng,iglo,nq,pff,nwpff)
             iglo=ig2pp
             nq=0
             ff(1,1)=2
-            ff(1,2)=e-epair
+            ff(1,2)=e-2.0*epair
          endif
          return
       endif

diff --git a/src/phys.f90 b/src/phys.f90
@@ -38,5 +38,6 @@ module physics
    real(kr),parameter,public::amassh=3.01493223469e0_kr !hellion (3)
    real(kr),parameter,public::amassa=4.0015061747e0_kr  !alpha
    real(kr),parameter,public::amasse=5.485799110e-4_kr  !electron
+   real(kr),parameter,public::epair=amasse*amu*clight*clight/ev
 end module physics