source: trunk/examples/extended/electromagnetic/TestEm4/geant3/src/uginit.F @ 807

Last change on this file since 807 was 807, checked in by garnier, 16 years ago

update

File size: 2.7 KB
Line 
1
2      SUBROUTINE UGINIT                                                         
3*                                                                               
4*     To initialise GEANT/USER  program and read data cards                     
5*
6      CHARACTER*20 filnam                                                                                                                                                                                                                             
7*                                                                               
8* *** Define the GEANT parameters                                                         
9      CALL GINIT
10     
11* *** read data cards
12      PRINT *, 'G3 > gives the filename of the data cards to be read:'
13      READ (*,'(A)') filnam
14      IF (filnam.EQ.' ') filnam = 'run01.dat'
15      OPEN (unit=5,file=filnam,status='unknown',form='formatted')                                                                         
16
17* *** read data cards                                                           
18      CALL GFFGO                                                                     
19*                                                                               
20      CALL GZINIT                                                               
21      CALL GPART
22*
23      CALL GDINIT                                                                     
24*                                                                               
25* *** Geometry and materials description                                       
26      CALL UGEOM                                                               
27*                                                                               
28* *** Energy loss and cross-sections initialisations                           
29      CALL GPHYSI                                                               
30*                                                                               
31      CALL GPRINT('MATE',0)                                                     
32      CALL GPRINT('TMED',0)                                                     
33      CALL GPRINT('VOLU',0)                                                     
34*                                                                                                                                                             
35* *** Define user histograms                                                   
36      CALL UHINIT                                                               
37*                                                                               
38      END                                                                       
Note: See TracBrowser for help on using the repository browser.