ugeom.F

Go to the documentation of this file.

      SUBROUTINE ugeom
      
* *** Define user geometry set up
*
#include "calor.inc"  
*
      dimension par(3)
*            
      dimension aair(2),zair(2),wair(2)
      dimension ah2o(2),zh2o(2),wh2o(2)
* 
* *** Air compound parameters
      DATA aair/14.01, 16.00/
      DATA zair/ 7.  ,  8.  /
      DATA wair/ 0.7 ,  0.3 /
*
* *** Air compound parameters
      DATA ah2o/ 1.01, 16.00/                                            
      DATA zh2o/ 1.  ,  8.  /
      DATA wh2o/ 2.  ,  1.  /
*
* *** Defines USER perticular materials
      CALL gsmixt( 1,'Air'      , aair ,zair, 1.29e-3, 2   , wair)
      CALL gsmate( 2,'H2 Liquid',  1.01,  1., 0.0708 , 865., 790., 0,0)
      CALL gsmixt( 3,'Water'    , ah2o ,zh2o, 1.0    ,-2   , wh2o)      
      CALL gsmate( 4,'Liquid Ar', 39.95, 18., 1.39   , 14.0, 84.0, 0,0)
      CALL gsmate( 5,'Aluminium', 26.98, 13., 2.7    ,  8.9, 37.2, 0,0)
      CALL gsmate( 6,'Iron'     , 55.85, 26., 7.87   , 1.76, 17.1, 0,0)
      CALL gsmate( 7,'Tungsten' ,183.85, 74., 19.30  , 0.35, 18.5, 0,0)
      CALL gsmate( 8,'Lead'     ,207.19, 82., 11.35  , 0.56, 18.5, 0,0)
      CALL gsmate( 9,'Uranium'  ,238.03, 92., 18.95  , 0.32, 12. , 0,0)
*      
      CALL gsmate(10,'Silicon'  , 28.09, 14.,  2.33  , 9.36, 45.5, 0,0)
      CALL gsmate(11,'Germanium', 72.61, 32., 5.323  , 2.30, 16.6, 0,0)      
*
* *** Defines USER tracking media parameters
      ifield = 0                                                             
      IF (field.GT.0.) ifield = 3
      fieldm = 10*field
      tmaxfd = 10.0
      stemax = 1000.
      if (stepmax.gt.0.) stemax = stepmax
      deemax = 0.20                                                            
      epsil  = 0.001                                                           
      stmin  = 0.010                                                           
*
      CALL gstmed( 1,'Container',imat, 0 ,ifield,fieldm,tmaxfd,
     *               stemax,deemax,epsil,stmin, 0 , 0 )
*
*
* *** Geometry
      par(1) = sizex /2. 
      par(2) = sizeyz/2.
      par(3) = sizeyz/2.
      CALL gsvolu('aBox','BOX ',1,par,3,ivol) 
*
      if (nlayers.gt.1) CALL gsdvn('layr','aBox',nlayers,1)       
* 
* *** Close geometry banks. (obligatory system routine) 
      CALL ggclos 
*
*
* *** dessin
      CALL gsatt('*','SEEN',1)
*
      DO ix = 1,3
        CALL gdopen(ix)
        scale =   18./sizex
        paxis =   0.
        saxis =   0.1*sizex
        CALL gdrawc('aBox',ix,0.,10.,9.3,scale,scale)
        CALL gdaxis(paxis,paxis,paxis,saxis)
        CALL gdscal(10. , 0.3)
        CALL gdclos
      END DO
*
                                                                               
      END