| 1 |
|
|---|
| 2 | SUBROUTINE UGEOM
|
|---|
| 3 |
|
|---|
| 4 | * *** Define user geometry set up
|
|---|
| 5 | *
|
|---|
| 6 | #include "calor.inc"
|
|---|
| 7 | *
|
|---|
| 8 | DIMENSION PAR(3)
|
|---|
| 9 | *
|
|---|
| 10 | DIMENSION Aair(2),Zair(2),Wair(2)
|
|---|
| 11 | DIMENSION AH2O(2),ZH2O(2),WH2O(2)
|
|---|
| 12 | *
|
|---|
| 13 | * *** Air compound parameters
|
|---|
| 14 | DATA Aair/14.01, 16.00/
|
|---|
| 15 | DATA Zair/ 7. , 8. /
|
|---|
| 16 | DATA Wair/ 0.7 , 0.3 /
|
|---|
| 17 | *
|
|---|
| 18 | * *** Air compound parameters
|
|---|
| 19 | DATA AH2O/ 1.01, 16.00/
|
|---|
| 20 | DATA ZH2O/ 1. , 8. /
|
|---|
| 21 | DATA WH2O/ 2. , 1. /
|
|---|
| 22 | *
|
|---|
| 23 | * *** Defines USER perticular materials
|
|---|
| 24 | CALL GSMIXT( 1,'Air' , Aair ,Zair, 1.29E-3, 2 , Wair)
|
|---|
| 25 | CALL GSMATE( 2,'H2 Liquid', 1.01, 1., 0.0708 , 865., 790., 0,0)
|
|---|
| 26 | CALL GSMIXT( 3,'Water' , AH2O ,ZH2O, 1.0 ,-2 , WH2O)
|
|---|
| 27 | CALL GSMATE( 4,'Liquid Ar', 39.95, 18., 1.39 , 14.0, 84.0, 0,0)
|
|---|
| 28 | CALL GSMATE( 5,'Aluminium', 26.98, 13., 2.7 , 8.9, 37.2, 0,0)
|
|---|
| 29 | CALL GSMATE( 6,'Iron' , 55.85, 26., 7.87 , 1.76, 17.1, 0,0)
|
|---|
| 30 | CALL GSMATE( 7,'Tungsten' ,183.85, 74., 19.30 , 0.35, 18.5, 0,0)
|
|---|
| 31 | CALL GSMATE( 8,'Lead' ,207.19, 82., 11.35 , 0.56, 18.5, 0,0)
|
|---|
| 32 | CALL GSMATE( 9,'Uranium' ,238.03, 92., 18.95 , 0.32, 12. , 0,0)
|
|---|
| 33 | *
|
|---|
| 34 | CALL GSMATE(10,'Silicon' , 28.09, 14., 2.33 , 9.36, 45.5, 0,0)
|
|---|
| 35 | CALL GSMATE(11,'Germanium', 72.61, 32., 5.323 , 2.30, 16.6, 0,0)
|
|---|
| 36 | *
|
|---|
| 37 | * *** Defines USER tracking media parameters
|
|---|
| 38 | IFIELD = 0
|
|---|
| 39 | IF (Field.GT.0.) IFIELD = 3
|
|---|
| 40 | FIELDM = 10*Field
|
|---|
| 41 | TMAXFD = 10.0
|
|---|
| 42 | STEMAX = 1000.
|
|---|
| 43 | if (stepmax.gt.0.) STEMAX = stepmax
|
|---|
| 44 | DEEMAX = 0.20
|
|---|
| 45 | EPSIL = 0.001
|
|---|
| 46 | STMIN = 0.010
|
|---|
| 47 | *
|
|---|
| 48 | CALL GSTMED( 1,'Container',Imat, 0 ,IFIELD,FIELDM,TMAXFD,
|
|---|
| 49 | * STEMAX,DEEMAX,EPSIL,STMIN, 0 , 0 )
|
|---|
| 50 | *
|
|---|
| 51 | *
|
|---|
| 52 | * *** Geometry
|
|---|
| 53 | PAR(1) = SizeX /2.
|
|---|
| 54 | PAR(2) = SizeYZ/2.
|
|---|
| 55 | PAR(3) = SizeYZ/2.
|
|---|
| 56 | CALL GSVOLU('aBox','BOX ',1,PAR,3,IVOL)
|
|---|
| 57 | *
|
|---|
| 58 | if (nLayers.gt.1) CALL GSDVN ('layr','aBox',nLayers,1)
|
|---|
| 59 | *
|
|---|
| 60 | * *** Close geometry banks. (obligatory system routine)
|
|---|
| 61 | CALL GGCLOS
|
|---|
| 62 | *
|
|---|
| 63 | *
|
|---|
| 64 | * *** dessin
|
|---|
| 65 | CALL GSATT ('*','SEEN',1)
|
|---|
| 66 | *
|
|---|
| 67 | DO IX = 1,3
|
|---|
| 68 | CALL GDOPEN (IX)
|
|---|
| 69 | SCALE = 18./SizeX
|
|---|
| 70 | PAXIS = 0.
|
|---|
| 71 | SAXIS = 0.1*SizeX
|
|---|
| 72 | CALL GDRAWC ('aBox',IX,0.,10.,9.3,SCALE,SCALE)
|
|---|
| 73 | CALL GDAXIS (PAXIS,PAXIS,PAXIS,SAXIS)
|
|---|
| 74 | CALL GDSCAL (10. , 0.3)
|
|---|
| 75 | CALL GDCLOS
|
|---|
| 76 | END DO
|
|---|
| 77 | *
|
|---|
| 78 |
|
|---|
| 79 | END
|
|---|
