tstEnergyResource.f90

./Core/re/f/tests/tstEnergyResource.f90

program tstenergyresource

#include "ScaleSTL/FortranTestMacros.h"

 use origen_energyresource_m

 implicit none

 type(origen_energyresource) :: fiscap
 real(C_DOUBLE) :: kcap, kfiss


 call fiscap%initialize

 !test behavior for a bad zzzaaai, i.e. -5000.
 insist(.not.fiscap%getKappaValues(-5000,kcap,kfiss),'found ZZZAAAI -5000 when did not expect to')
 insist(abs(kcap/5.0d0-1.0d0) < 1.0d-10,'By default, kappa-capture should be 5 MeV')
 insist(abs(kfiss/200.0d0-1.0d0) < 1.0d-10,'By default, kappa-fission should be 200 MeV')

 !test behavior for a good zzzaaai, i.e. 92235.
 insist(fiscap%getKappaValues(92235,kcap,kfiss),'should have found ZZZAAAI 92235')
 insist(abs(kcap/6.5451d0-1.0d0) < 1.0d-10,'92235 kappa-capture should be 6.5451 MeV')
 insist(abs(kfiss/194.02d0-1.0d0) < 1.0d-10,'92235 kappa-fission should be 194.02 MeV')


 call fiscap%destroy

end program