source: trunk/examples/extended/electromagnetic/TestEm1/geant3/src/uginit.F @ 1358

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

update

File size: 3.1 KB
Line 
1
2      SUBROUTINE UGINIT                                                         
3*                                                                               
4*     To initialise GEANT/USER  program and read data cards                     
5*
6#include "calor.inc"
7#include "celoss.inc"                                                                                 
8*
9      CHARACTER*20 filnam                                                                                                                                                                                                     
10*                                                                               
11* *** Define the GEANT parameters                                                         
12      CALL GINIT
13                                                                         
14* *** read data cards
15      PRINT *, 'G3 > gives the filename of the data cards to be read:'
16      READ (*,'(A)') filnam
17      IF (filnam.EQ.' ') filnam = 'testem1.dat'
18      OPEN (unit=5,file=filnam,status='unknown',form='formatted')
19                                             
20*
21* *** Calor definition                                               
22      CALL FFKEY('CALOR',Imat,3,'MIXED')
23*                                                                   
24* *** read data cards                                                           
25      CALL GFFGO
26                                                                     
27*                                                                               
28      CALL GZINIT                                                               
29      CALL GPART
30*
31      CALL GDINIT                                                                     
32*                                                                               
33* *** Geometry and materials description                                       
34      CALL UGEOM                                                               
35*                                                                               
36* *** Energy loss and cross-sections initialisations                           
37      CALL GPHYSI                                                               
38*                                                                               
39      CALL GPRINT('MATE',0)                                                     
40      CALL GPRINT('TMED',0)                                                     
41      CALL GPRINT('VOLU',0)                                                     
42*                                                                                                                                                             
43* *** Define user histograms                                                   
44      CALL UHINIT
45*
46* *** some initialisation
47      NTRAK0 = 0
48      NTRAK1 = 0
49      NSTEP0 = 0
50      NSTEP1 = 0
51      CALL VZERO (NBCALL,12)                                                               
52*                                                                               
53      END                                                                       
Note: See TracBrowser for help on using the repository browser.