TestEm2/geant3/src/uginit.F

Go to the documentation of this file.

      SUBROUTINE uginit                                                         
*                                                                               
*     To initialise GEANT/USER  program and read data cards                     
*                                                                                                                                                                                                                       
#include "pvolum.inc"
#include "celoss.inc"
*
      CHARACTER*20 filnam 
*                                                                                     
* *** Define the GEANT parameters                                                          
      CALL ginit 
                                                                          
* *** read data cards
      print *, 'G3 > gives the filename of the data cards to be read:'
      READ (*,'(A)') filnam
      IF (filnam.EQ.' ') filnam = 'run01.dat'
      OPEN (unit=5,file=filnam,status='unknown',form='formatted')
*                                              
* *** material definition                                                                                       
      CALL ffkey('MATE',imat,1,'INTEGER') 
*
* *** volumes and bins definition                                               
      CALL ffkey('BINS',nltot,4,'MIXED')
*      
* *** max allowed step size       
      CALL ffkey('STEPMX',stepmax,1,'REAL')      
*                                                         
* *** read data cards                                                           
      CALL gffgo 
*
* *** check size of arrays
      if (nltot.gt.nbin) then
        write (6,51) nltot, nbin
    nltot = nbin
      endif
      if (nrtot.gt.nbin) then
        write (6,52) nrtot, nbin
    nrtot = nbin
      endif         
51    FORMAT (/,5x,'warning (uginit): NLTOT= ',i3,' truncated to ',i3)
52    FORMAT (/,5x,'warning (uginit): NRTOT= ',i3,' truncated to ',i3)
*                                                                     
* *** achieve initialization                                                    
      CALL vzero(sel1,10*nbin+8)                                                  
*                                                                               
      CALL gzinit                                                               
      CALL gpart 
*
      CALL gdinit                                                                     
*                                                                               
* *** Geometry and materials description                                        
      CALL ugeom                                                                
*                                                                               
* *** Energy loss and cross-sections initialisations                            
      CALL gphysi                                                               
*                                                                               
      CALL gprint('MATE',0)                                                     
      CALL gprint('TMED',0)                                                     
      CALL gprint('VOLU',0)                                                     
*                                                                                                                                                              
* *** Define user histograms                                                    
      CALL uhinit                                                               
*                                                                               
      END