From 3a9f2268def5f2354069f8d5ea9d78a4f76491ca Mon Sep 17 00:00:00 2001 From: walshjon Date: Wed, 10 Sep 2014 12:20:58 -0700 Subject: [PATCH] do not add 0 K energy points to union grid --- docs/source/usersguide/input.rst | 6 +++--- src/ace_header.F90 | 4 ++-- src/energy_grid.F90 | 3 --- src/physics.F90 | 5 ++--- 4 files changed, 7 insertions(+), 11 deletions(-) diff --git a/docs/source/usersguide/input.rst b/docs/source/usersguide/input.rst index f50c7ff217..c22008b793 100644 --- a/docs/source/usersguide/input.rst +++ b/docs/source/usersguide/input.rst @@ -274,7 +274,7 @@ attributes or sub-elements: of the nuclide to which a resonance scattering treatment is to be applied. The ``method`` attribute gives the type of resonance scattering treatment that is to be applied to the ``nuclide``. Acceptable inputs - none of - which are case-sensitive - for the ``method`` attribute are ``ARTS``, + which are case-sensitive - for the ``method`` attribute are ``ARES``, ``CXS``, ``WCM``, and ``DBRC``. Descriptions of each of these methods are documented here_. The ``xs_label`` attribute gives the label for the cross section data of the ``nuclide`` at a given temperature. The @@ -291,7 +291,7 @@ attributes or sub-elements: U-238 - ARTS + ARES 92238.72c 92238.00c 5.0e-6 @@ -317,7 +317,7 @@ attributes or sub-elements: be specified, each corresponding to a single nuclide at a single temperature. - *Defaults*: None (scatterer), ARTS (method), 0.01 eV (E_min), 1.0 keV (E_max) + *Defaults*: None (scatterer), ARES (method), 0.01 eV (E_min), 1.0 keV (E_max) ```` Element ---------------------- diff --git a/src/ace_header.F90 b/src/ace_header.F90 index b0fec99b6c..6d339891c9 100644 --- a/src/ace_header.F90 +++ b/src/ace_header.F90 @@ -159,14 +159,14 @@ module ace_header end type Nuclide !=============================================================================== -! NUCLIDE_0K temporarily contains all 0K cross section data and other parameters +! NUCLIDE0K temporarily contains all 0K cross section data and other parameters ! needed to treat resonance scattering before transferring them to NUCLIDE !=============================================================================== type Nuclide0K character(10) :: nuclide ! name of nuclide, e.g. U-238 - character(16) :: scheme = 'arts' ! target velocity sampling scheme + character(16) :: scheme = 'ares' ! target velocity sampling scheme character(10) :: name ! name of nuclide, e.g. 92235.03c character(10) :: name_0K ! name of 0K nuclide, e.g. 92235.00c real(8) :: E_min = 0.01e-6 ! lower cutoff energy for res scattering diff --git a/src/energy_grid.F90 b/src/energy_grid.F90 index a38a305d3e..66b95be0de 100644 --- a/src/energy_grid.F90 +++ b/src/energy_grid.F90 @@ -28,9 +28,6 @@ contains do i = 1, n_nuclides_total nuc => nuclides(i) call add_grid_points(list, nuc % energy) - - ! add 0K points to the grid - if (nuc % resonant) call add_grid_points(list, nuc % energy_0K) end do ! Set size of unionized energy grid diff --git a/src/physics.F90 b/src/physics.F90 index e0fbf80d25..04cc66afa2 100644 --- a/src/physics.F90 +++ b/src/physics.F90 @@ -2,8 +2,7 @@ module physics use ace_header, only: Nuclide, Reaction, DistEnergy use constants - use cross_section, only: elastic_xs_0K, find_energy_index, & - union_grid_index + use cross_section, only: elastic_xs_0K use endf, only: reaction_name use error, only: fatal_error, warning use fission, only: nu_total, nu_delayed @@ -885,7 +884,7 @@ contains if (.not. reject) exit end do - case ('arts') + case ('ares') E_red = sqrt((awr * E) / kT) E_low = (((E_red - 4.0_8)**2) * kT) / awr E_up = (((E_red + 4.0_8)**2) * kT) / awr