Computer Programs
CCC-0650 DOORS3.2A.
last modified: 29-NOV-2022 | catalog | categories | new | search |

CCC-0650 DOORS3.2A.

DOORS3.2A, 1-,2-,3-dimensional discrete-ordinates system for deep-penetration neutron and photon transport

top ]


top ]

To submit a request, click below on the link of the version you wish to order. Rules for end-users are available here.

Program name Package id Status Status date
DOORS3.2A CCC-0650/04 Arrived 24-JUL-2007

Machines used:

Package ID Orig. computer Test computer
CCC-0650/04 Linux-based PC,PC Windows,DEC ALPHA W.S.,IBM RISC6000 WS,SUN W.S.
top ]

The DOORS3.2a discrete ordinates transport code system includes the most recent versions of CCC-0543/TORT-DORT, CCC-0254/ANISN-ORNL, CCC-0628/GBANISN and CCC-0351/FALSTF. It also includes the ISOPLOT code from the PSR-0155/DOGS package and various utility programs listed below which were previously included in the TORT-DORT package. In this release each module is a separate executable file. Several modules, as needed, can be run in a single job by using "jdos" to call the "drv" module which interprets the sequence specified in the input.


TORT calculates the flux of fluence of particles due to particles incident upon the system from extraneous sources or generated internally as a result of interaction with the system in two-or three-dimensional geometric systems, and DORT is used in one- or two-dimensional geometric systems. The principle application is to the deep-penetration transport of neutrons and photons. Reactor eigenvalue problems can also be solved. Numerous printed edits of the results are available, and results can be transferred to output files for subsequent analysis.


ANISN solves the one-dimensional Boltzmann transport equation for neutrons or gamma rays in slab, spherical, or cylindrical geometry. GBANISN is based on ANISN but was modified to allow randomizing of multibank fluxes within a group at all interfaces between dissimilar materials and a reduction in the number of outer iterations for problems involving neutron upscatter into higher energy groups. GBANISN requires fewer outer iterations than ANISN by performing "inner iterations" over energy groups within a "band" and converging those groups before moving to the next band. These "inner" iterations slightly resemble outer iterations in ANISN. Thus, a calculation with upscatter and no fission can be solved with one traditional outer iteration. GBANISN, like ANISN, includes a technique for handling general anisotropic scattering, pointwise convergence criteria, and alternate step function difference equations that effectively remove the oscillating flux distributions sometimes found in discrete ordinates solutions.


ISOPLOT uses the Sandia National Laboratory RSCORS graphical system. SNL's CMP system for code maintenance is used to build the Fortran source files for the target computer. The RSCORS library (librscors.a) can be used to build executables for BOT3P2.0 and BOT3P3.0 (Bologna Transport Analysis Pre-Post-Processors).


DOORS reads ANISN-format cross-section libraries, which are not included in the package. Users may choose from several available in RSICC's data library collection.


In October 2003, the package was updated to DOORS3.2a to meet increasing requests for support on Linux and Windows systems. The remainder of the package, other than relatively minor modifications to support PC architectures and current compilers, is unchanged.

top ]

The Boltzmann transport equation is solved using the method of discrete ordinates to treat the directional variable and weighted finite-difference methods, in addition to Linear Nodal and Linear Characteristic methods in TORT to treat spatial variables. Energy dependence is treated using a multigroup formulation. Time dependence is not treated. Starting in one corner of a mesh, at the highest energy, and with starting guesses for implicit sources, boundary conditions and recursion relationship are used to sweep into the mesh for each discrete direction independently. Integral quantities such as scalar flux are obtained  from weighted sums of the directional results. The calculation then proceeds to lower energy groups, one at a time.


Iterations are used to resolve implicitness caused by scattering between directions within a single energy group, by scattering from an energy group to another group previously calculated, by fission, and by certain types of boundary conditions. Methods are available to accelerate convergence of both inner and outer iterations. Anisotropic scattering is represented by a Legendre expansion of arbitrary order, and methods are available to mitigate the effect of negative scattering source estimates resulting from finite truncation of the expansion. Direction sets can be biased, concentrating work into directions of particular interest. Fixed sources can be specified at either external or internal mesh boundaries, or distributed within mesh cells.


Two multitasking options are available for solving the Linear Nodal method equations in TORT on UNICOS Cray platforms: Octant Parallel (OP) which concurrently performs the right-&left and left-&right sweeps for each in the mesh; Direction Parallel (DP) which concurrently sweeps each row in all directions within an octant. Wall-clock speedup generally scales up with problem size but is very sensitive to machine loading during execution. On a lightly loaded Y/MP speedup factor of 5 was achieved using 8 tasks with a problem as small as ~104,000 cells.

top ]

Penetration through large, non-scattering regions may become inaccurate due to ray effects. Problems with scattering ratios near unity or eigenvalue calculations with closely spaced eigenvalues may be quite time-consuming. Flexible dimensioning is used throughout so that no fixed limits on groups, problem size, etc., are applicable. External forces and nonlinear physical effects cannot be treated.

top ]

Central processor unit (CPU) time used by the flux sweep is roughly proportional to the number of flux calculations: spatial mesh cells * directions * energy groups * iterations/group. Times for running all test cases in the package range from about 1o minutes on an IBM RS/6000 Model 590 to several hours on slower machines.

top ]

The code system contains the following modules:

  • TORT     Three-dimensional neutron/photon transport

  • DORT     Two-dimensional neutron/photon transport

  • ANISN    One-dimensional neutron/photon transport

  • GBANISN  ANISN with "group band" option

  • JDOS     Script for performing calculations

  • DRV      DOORS driver module called by JDOS

  • ALC      Material cross-section library maintenance

  • GIP      Group cross-section preparation

  • TORSED   Couple R-Z DORT to X-Y-Z TORT calculation

  • TORSET   Couple primary TORT to secondary TORT calculation

  • VISA     Prepare a TORSED input file from a DORT output file

  • GRTUNCL  Uncolided flux and first collision source

  • FALSTF   Last collision source projection to point detector

  • BNDRYS   Convert internal boundary flux to internal boundary source

  • RTFLUM   Convert flux moment files between various formats

  • ISOPLOT  Display 2-D flux or response contours

  • XTORID   Prepare 2-D flux slices from TORT for display in ISPL3D

  • CMP      Sandia National Lab. code system maintenance and update

  • RSCORS   Sandia National Lab. subroutine library of graphical primitives

top ]
Package ID Status date Status
CCC-0650/04 24-JUL-2007 Masterfiled Arrived
top ]
  • S.L. Thompson: CMP Code Maintenance Package, SAND85-0825 Sandia National Laboratory Working Document (April 8, 1991)

  • S.L. Thompson: The RSCORS Graphics System, SAND99-XXXX Sandia National Laboratory Working Document (October 19, 1991)

  • S.L. Thompson: A Definition of the Basic Graphics Package (BGP) Intermediate File Format, Based on unpublished letter (March 1978)

CCC-0650/04, included references:
- W. A. Rhoades and D. B. Simpson:
The TORT Three-Dimensional Discrete Ordinates Neutron/Photon Transport Code,
ORNL/TM-13221 (October 1997)
- M. B. Emmett, W. A. Rhoades, R. L. Childs, and J. O. Johnson:
the DORT section of "A User's Manual for MASH 1.0 - A Monte Carlo Adjoint
Shielding Code System," ORNL/TM-11778 (March 1992)
- R. L. Childs
FALSTF Last-Flight Computer Program, ORNL/TM-12675 (Jan. 1996)
Input Modifications Required for ANISN and GBANISN (July 1998)
- W. W. Engle, Jr.:
ANISN, A One-Dimensional Discrete Ordinates Transport Code with Anisotropic
Scattering, K-1693 (March 1967)
- W. W. Engle, Jr., M. A. Boling, and B. W. Colston:
DTF-III, A One-Dimensional, Multigroup Neutron Transport Program, NAA-SR-10951
(March 1966)
- R. L. Childs and D. E. Cullen:
GROUPBAND-ANISN: A Code to Perform Multiband Calculations, Informal Notes (June
- W. A. Rhoades and M. B. Emmett:
DOS: The Discrete Ordinates System, ORNL/TM 8362 (September 1982)
- C. O. Slater:
The XTORID and ISPL3D Codes for Plotting TORT Activities, ORNL Internal Memo
(March 25, 1996)
- R. L. Childs:
GRTUNCL: First Collision Source Program, ORNL Informal Notes (1982)
- R. L. Childs, W. A. Rhoades:
Theoretical Basis of the Linear Nodal and Linear Characteristic Methods in the
TORT Computer Code, ORNL/TM-12246 (January 1993)
- W. A. Rhoades:
The TORSED Method for Construction of TORT Boundary Sources from External DORT
Flux Files, ORNL/TM-12359 (August 1993)
- W. A. Rhoades and D. B. Simpson:
Splicing and Bootstrapping Methods for DORT/TORT-to-TORT Coupling, (TORSET
Section), ORNL/TM-13350 (to be published)
- W. A. Rhoades:
The ALC1 Program for Cross-Section Library Management, ORNL/TM-4015 (December
top ]

DOORS3.2a modules run on IBM RS/6000, Sun, DEC Alpha and personal computers. The previous DOORS3.2 release ran on Cray, SGI, and Hewlett Packard. These computers are not available to test DOORS3.2a buy may still work. Reducing compiler optimization may be required.

top ]
Package ID Computer language
top ]

DOORS3.2a runs under Unix, Linux and Windows operating systems. Some numerical errors may occur when solving difficult problems on short-word configurations. Executables created with Portland Group, Inc. compilers are included both for Linux and Windows. All other systems require Fortran and C compilers. However, source code is only provided for rscors and sltutil programs.


DOORS 3.2a was tested at RSICC on the following systems:


  • PC Windows XP with Portland Group Fortran 4.0 2 and Portland Group C 4.0-2 - requires Cygwin environment to build (see

  • AMD Athlon running RedHat Linux 7.3 with Portland Group, Inc. 4.0-2 compiler

  • AMD Athlon running RedHat Linux 7.3 with GNU Fortran 0.5.26 and GNU gcc 2.96

  • PC Red Hat Linux 9.0 system with GNU Fortran 3.2.2 and GNU gcc 2.96/GNU gcc 3.2.2

  • DEC Alpha - Compaq Tru64 UNIX V5.0A - Compaq Fortran Compiler V5.5-1877-48BBF

  • DEC alpha - Digital UNIX V40.D (Rev. 878) - Compaq Fortran Compiler V5.5-1877-48BBF

  • IBM RS/6000 - AIX 5.1 system - IBM XL Fortran for AIX Version 08.01.0000.0000

  • IBM RS/6000 - AIX 4.3.3 XLF 6.1 > XLC/C++ 4.4

  • Sun SparcStation - Sun OS 5.6 - Sun WorkShop Compilers 5.0 98/12/15 FORTRAN 77 5.0

  • Sun SparcStation - Sun OS 5.7 - Sun WorkShop 6 2000/04/07 FORTRAN 77 5.1

  • UBUNTU Linux 7.04 with g77 [GNU Fortran (GCC) 3.4.6 (Ubuntu 3.4.6-5ubuntu1)]


The previous DOORS 3.2 release ran under the following systems that are no longer available for testing.


  • Cray Unicos 8.0.3, cft77 and Cray J90 Unicos, f90

  • SGI IP22, IRIX 5.3, F77 Ver 4.0.2

  • HP 9000/715, HP-UX 10.20, HP Fortran 77 Version B


Added May 2007: DOORS 3.2a will build successfully and run test cases under ubuntu linux 7.04 using g77. Detailed instructions can be found in entry 26 of the DOORS electronic notebook on RSICC's website. Only one small change to the mupopx11.lnx script is required to fix a problem with the X11 library.
Modify Install/mupopx11.lnx:
change line with X1


The distributed Windows executables can be run in a command prompt window (of WindowsXP or Windows2000) in a manner similar to UNIX executables (uses redirection for input and output). The included jdos utililty is a c shell script intended for UNIX and UNIX like systems. It can be used under the Cygwin environment on Windows. Without the availability of the jdos script, each of the codes can be run separately, or in a batch file. See Doc/doors-3.2a-notes.txt for PC information. Installation instructions are in file Doc/install.doc.


WinZIP 8.0 or newer is required to expand the distribution file under Windows.

top ]

Contributed by: Radiation Safety Information Computational Center
                Oak Ridge National Laboratory
                Oak Ridge, Tennessee, U. S. A.


Developed by:   Oak Ridge National Laboratory
                Oak Ridge, Tennessee, U. S. A.

top ]
Doc - install.doc and other information files and reports
Install - installation scripts
Outputs - developer's output for Cray and IBM
bin - executables for Windows and Linux
dort - dort, dos and falstf, etc. scripts
isoplot - isoplot source and test case
rscors -  source code
sltutil - source code, cmp and other Sandia National Lab. Utilities
tort - tort scripts
top ]
  • C. Static Design Studies
  • J. Gamma Heating and Shield Design
  • K. Reactor Systems Analysis

Keywords: complex geometry, cylindrical geometry, discrete ordinate method, gamma-ray, multigroup, neutron transport equation, one-dimensional, photon transport, shielding, slabs, spherical geometry, three-dimensional, two-dimensional.