Zubal Head Directory README.  September 30, 2005
c Modified Zubal Phantom Input Deck by Jeff Evans. Revisions by Chenguan Li.
c   Department of Mechanical Engineering, The Ohio State University, Columbus 43210, USA.
c
c Based on the Zubal Phantom. ref - 
c   I. G. Zubal, C. R. Harrell, E. O. Smith, and A. L. Smith, 
c   "Two dedicated software, voxel-based, anthropomorphic (torso and head) phantoms," 
c   in Proceedings of the International Workshop, National Radiological Protection Board, 
c   Chilton, UK, on 6 and 7 July 1995, edited by P. J. Dimbylow (NRPB, Chilton, 1996), pp. 105-111. 
c
c This geometric model has been used in:  J.F. Evans, T. E. Blue, N. Gupta,
c   "Absorbed dose estimates to structures of the brain and head using a high-resolution 
c   voxel-based head phantom." Med Phys. 2001 May;28(5):780-6. 
c  
c Most of the head is 85 x 109 x 120 lattice of voxels, where each voxel is 2.2 x 2.2 x 1.4 mm3. 
c   29 critical structures of the head are identified by their individual universes.  It is
c   attached to a simplistic model of the neck and torso.
c
c Other citations in this input deck:
c   J. E. Woollard, T. E. Blue, N. Gupta, and R. A. Gahbauer, "Development and application of 
c   neutron field optimization parameters for an accelerator-based neutron source for boron 
c   neutron capture therapy," Nucl. Technol. 115, 100-113 (1996). 
c
c   ICRU 46, "Photon, electron, proton, and neutron interaction data for body tissues," 
c   International Commission on Radiation Units and Measurements, Bethesda, MD, 1992.

This input deck is quite descriptive and uses a lattice-based geometry.  While MCNP5 will run this
input deck, the following changes will reduce set-up times from hours or days to minutes.  
The following changes to the source code have not been checked for accuracy and may change tally results.
These changes must be used with caution, and they require that the user re-compile the MCNP source
to create their own version of MCNP.

Comment out all lines (except continue lines) between 590 continue and 670 continue in checkcs.F90


