* * HDDSgeant3 - fortran geometry definition package * for the Hall D experiment. * * WARNING: DO NOT EDIT THIS FILE * * This file was generated automatically from the * HDDS xml geometry definition by the hdds-geant * translator. Any changes made to this file will * disappear as soon as it is regenerated from the * xml source. To introduce Geant3 optimizations, * see the subroutine Goptimize() in goptimize.F. * subroutine tgr_HDDSgeant3 implicit none integer imate character*20 chnama,namate real a,z,dens,radl,absl,ubuf(99) integer nwbuf real amat(99),zmat(99),wmat(99) integer nlmat integer itmed character*20 natmed integer nmat,isvol,ifield real fieldm,tmaxfd,stemax,deemax,epsil,stmin character*4 chname,chshap,chmoth integer nmed,npar,ivolu real par(99) integer irot real theta1,phi1,theta2,phi2,theta3,phi3 integer nr,ndiv,iaxis,numed,ndvmax real step,c0 real x,y character*4 chonly imate = 2 chnama = 'Hydrogen' a = 1.00797 z = 1 dens = 8.38e-05 radl = 731265 absl = 606205 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 3 chnama = 'Oxygen' a = 15.9995 z = 8 dens = 0.00133 radl = 25744.4 absl = 68421.1 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 4 chnama = 'Aluminum' a = 26.9815 z = 13 dens = 2.7 radl = 8.89259 absl = 39.4074 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 5 namate = 'Calcium' a = 40.08 z = 20 dens = 1.55 nlmat = 1 wmat(1) = 1 call gsmixt(imate,namate,a,z,dens,nlmat,wmat) imate = 6 chnama = 'Silicon' a = 28.086 z = 14 dens = 2.33 radl = 9.36481 absl = 45.4936 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 7 chnama = 'Iron' a = 55.847 z = 26 dens = 7.87 radl = 1.75858 absl = 16.7598 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 1 namate = 'Concrete' wmat(1) = 0.004 call gfmate(2,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf) wmat(2) = 0.509 call gfmate(3,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf) wmat(3) = 0.034 call gfmate(4,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf) wmat(4) = 0.07 call gfmate(5,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf) wmat(5) = 0.345 call gfmate(6,chnama,amat(5),zmat(5),dens,radl,absl,ubuf,nwbuf) wmat(6) = 0.038 call gfmate(7,chnama,amat(6),zmat(6),dens,radl,absl,ubuf,nwbuf) dens = 2.35 nlmat = 6 call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat) itmed = 1 natmed = 'HILL Concrete' nmat = 1 isvol = 0 ifield = 0 fieldm = 0 tmaxfd = 0 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'HILL' chshap = 'BOX ' nmed = 1 npar = 3 par(1) = 2500 par(2) = 2500 par(3) = 5000 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.1) stop 'consistency check #1 failed' imate = 9 chnama = 'Nitrogen' a = 14.0067 z = 7 dens = 0.00116 radl = 32750 absl = 75689.7 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 10 chnama = 'Argon' a = 39.948 z = 18 dens = 0.00166 radl = 11777.1 absl = 70602.4 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) imate = 8 namate = 'Air' wmat(1) = 0.7494 call gfmate(9,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf) wmat(2) = 0.2369 call gfmate(3,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf) wmat(3) = 0.0129 call gfmate(10,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf) wmat(4) = 0.0008 call gfmate(2,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf) dens = 0.001214 nlmat = 4 call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat) itmed = 2 natmed = 'AREA Air' nmat = 8 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'AREA' chshap = 'BOX ' nmed = 2 npar = 3 par(1) = 350 par(2) = 225 par(3) = 975 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.2) stop 'consistency check #1 failed' chname = 'AREA' nr = 1 chmoth = 'HILL' x = -130 y = 50 z = 900 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) imate = 11 chnama = 'Vacuum' a = 0 z = 0 dens = 0 radl = 1e+30 absl = 1e+30 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) itmed = 3 natmed = 'GONI Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'GONI' chshap = 'BOX ' nmed = 3 npar = 3 par(1) = 30 par(2) = 30 par(3) = 30 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.3) stop 'consistency check #1 failed' chname = 'GONI' nr = 1 chmoth = 'AREA' x = 130 y = -50 z = -900 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 4 natmed = 'GOBX Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'GOBX' chshap = 'PGON' nmed = 4 npar = 10 par(1) = -45 par(2) = 360 par(3) = 4 par(4) = 2 par(5) = -30 par(6) = 29.4 par(7) = 30 par(8) = 30 par(9) = 29.4 par(10) = 30 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.4) stop 'consistency check #1 failed' chname = 'GOBX' nr = 1 chmoth = 'GONI' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 5 natmed = 'GOBE Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'GOBE' chshap = 'PGON' nmed = 5 npar = 10 par(1) = -45 par(2) = 360 par(3) = 4 par(4) = 2 par(5) = -0.3 par(6) = 2 par(7) = 30 par(8) = 0.3 par(9) = 2 par(10) = 30 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.5) stop 'consistency check #1 failed' chname = 'GOBE' nr = 1 chmoth = 'GONI' x = 0 y = 0 z = -29.7 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'GOBE' nr = 2 chmoth = 'GONI' x = 0 y = 0 z = 29.7 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) imate = 12 chnama = 'Carbon' a = 12.011 z = 6 dens = 2.265 radl = 18.8521 absl = 38.1015 nwbuf = 0 call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf) itmed = 6 natmed = 'XTAL Carbon' nmat = 12 isvol = 1 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'XTAL' chshap = 'BOX ' nmed = 6 npar = 3 par(1) = 0.25 par(2) = 0.25 par(3) = 0.0001 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.6) stop 'consistency check #1 failed' chname = 'XTAL' nr = 1 chmoth = 'GONI' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 7 natmed = 'QUAV Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 0.1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'QUAV' chshap = 'TUBE' nmed = 7 npar = 3 par(1) = 0 par(2) = 2 par(3) = 137.1 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.7) stop 'consistency check #1 failed' chname = 'QUAV' nr = 1 chmoth = 'AREA' x = 130 y = -50 z = -733.985 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 8 natmed = 'QUAP Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 0.1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'QUAP' chshap = 'TUBE' nmed = 8 npar = 3 par(1) = 1.8 par(2) = 2 par(3) = 137 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.8) stop 'consistency check #1 failed' chname = 'QUAP' nr = 1 chmoth = 'QUAV' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 9 natmed = 'QUAD Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'QUAD' chshap = 'TUBE' nmed = 9 npar = 3 par(1) = 2 par(2) = 12 par(3) = 30 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.9) stop 'consistency check #1 failed' chname = 'QUAD' nr = 1 chmoth = 'AREA' x = 130 y = -50 z = -748.985 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 10 natmed = 'SHOE Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'SHOE' chshap = 'BOX ' nmed = 10 npar = 3 par(1) = 22.5 par(2) = 10.5 par(3) = 153.231 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.10) stop 'consistency check #1 failed' irot = 1 theta1 = 84.1 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 5.9 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'SHOE' nr = 1 chmoth = 'AREA' x = 119.815 y = -38 z = -444.49 irot = 1 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 2 theta1 = 84.1 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 5.9 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'SHOE' nr = 2 chmoth = 'AREA' x = 119.815 y = -62 z = -444.49 irot = 2 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 11 natmed = 'YOK1 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'YOK1' chshap = 'BOX ' nmed = 11 npar = 3 par(1) = 24 par(2) = 22.5 par(3) = 153.231 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.11) stop 'consistency check #1 failed' irot = 3 theta1 = 84.1 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 5.9 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOK1' nr = 1 chmoth = 'AREA' x = 181.984 y = -50 z = -438.065 irot = 3 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 12 natmed = 'YOKC Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'YOKC' chshap = 'BOX ' nmed = 12 npar = 3 par(1) = 54.5 par(2) = 24 par(3) = 153.231 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.12) stop 'consistency check #1 failed' irot = 4 theta1 = 84.1 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 5.9 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOKC' nr = 1 chmoth = 'AREA' x = 151.646 y = -3.5 z = -441.201 irot = 4 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 5 theta1 = 84.1 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 5.9 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOKC' nr = 2 chmoth = 'AREA' x = 151.646 y = -96.5 z = -441.201 irot = 5 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 6 theta1 = 83.208 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.792 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'SHOE' nr = 3 chmoth = 'AREA' x = 72.443 y = -38 z = -96.697 irot = 6 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 7 theta1 = 83.208 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.792 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'SHOE' nr = 4 chmoth = 'AREA' x = 72.443 y = -62 z = -96.697 irot = 7 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 13 natmed = 'YOK2 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'YOK2' chshap = 'BOX ' nmed = 13 npar = 3 par(1) = 24 par(2) = 22.5 par(3) = 153.231 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.13) stop 'consistency check #1 failed' irot = 8 theta1 = 83.208 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.792 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOK2' nr = 1 chmoth = 'AREA' x = 134.504 y = -50 z = -89.3054 irot = 8 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 14 natmed = 'PXHO Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'PXHO' chshap = 'TUBE' nmed = 14 npar = 3 par(1) = 0 par(2) = 1.5 par(3) = 154.5 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.14) stop 'consistency check #1 failed' irot = 9 theta1 = 96.79 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.79 phi3 = 0 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'PXHO' nr = 1 chmoth = 'YOK2' x = -4.536 y = 0 z = 0 irot = 9 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 10 theta1 = 83.208 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.792 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOKC' nr = 3 chmoth = 'AREA' x = 104.218 y = -3.5 z = -92.9125 irot = 10 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 11 theta1 = 83.208 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.792 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'YOKC' nr = 4 chmoth = 'AREA' x = 104.218 y = -96.5 z = -92.9125 irot = 11 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 15 natmed = 'VBS1 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBS1' chshap = 'TRAP' nmed = 15 npar = 11 par(1) = 148.743 par(2) = 7.35 par(3) = 180 par(4) = 3.8 par(5) = 38.7425 par(6) = 38.7425 par(7) = 0 par(8) = 3.8 par(9) = 52.5245 par(10) = 52.5245 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.15) stop 'consistency check #1 failed' chname = 'VBS1' nr = 1 chmoth = 'AREA' x = 101.663 y = -50 z = -448.743 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 16 natmed = 'VBF1 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBF1' chshap = 'TRAP' nmed = 16 npar = 11 par(1) = 148.743 par(2) = 7.35 par(3) = 180 par(4) = 0.5 par(5) = 38.7425 par(6) = 38.7425 par(7) = 0 par(8) = 0.5 par(9) = 52.5245 par(10) = 52.5245 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.16) stop 'consistency check #1 failed' chname = 'VBF1' nr = 1 chmoth = 'VBS1' x = 0 y = -3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'VBF1' nr = 2 chmoth = 'VBS1' x = 0 y = 3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 17 natmed = 'VBB1 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBB1' chshap = 'BOX ' nmed = 17 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 149.25 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.17) stop 'consistency check #1 failed' irot = 12 theta1 = 85.27 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 4.73 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBB1' nr = 1 chmoth = 'VBS1' x = 45.133 y = 0 z = 0 irot = 12 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 18 natmed = 'VBE1 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBE1' chshap = 'BOX ' nmed = 18 npar = 3 par(1) = 38.7425 par(2) = 2.8 par(3) = 0.5 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.18) stop 'consistency check #1 failed' chname = 'VBE1' nr = 1 chmoth = 'VBS1' x = 19.187 y = 0 z = -148.242 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 19 natmed = 'VBX1 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBX1' chshap = 'TUBE' nmed = 19 npar = 3 par(1) = 0 par(2) = 1.5 par(3) = 0.525 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.19) stop 'consistency check #1 failed' chname = 'VBX1' nr = 1 chmoth = 'VBE1' x = 9.149 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 20 natmed = 'VBW1 Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBW1' chshap = 'BOX ' nmed = 20 npar = 3 par(1) = 0.05 par(2) = 2.8 par(3) = 151.011 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.20) stop 'consistency check #1 failed' irot = 13 theta1 = 80.06 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.94 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBW1' nr = 1 chmoth = 'VBS1' x = -45.334 y = 0 z = 0 irot = 13 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 21 natmed = 'VBS2 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBS2' chshap = 'TRAP' nmed = 21 npar = 11 par(1) = 180 par(2) = 8.15 par(3) = 180 par(4) = 3.8 par(5) = 44.0245 par(6) = 44.0245 par(7) = 0 par(8) = 3.8 par(9) = 55.583 par(10) = 55.583 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.21) stop 'consistency check #1 failed' chname = 'VBS2' nr = 1 chmoth = 'AREA' x = 48.196 y = -50 z = -120 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 22 natmed = 'VBF2 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBF2' chshap = 'TRAP' nmed = 22 npar = 11 par(1) = 180 par(2) = 8.15 par(3) = 180 par(4) = 0.5 par(5) = 44.0245 par(6) = 44.0245 par(7) = 0 par(8) = 0.5 par(9) = 55.583 par(10) = 55.583 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.22) stop 'consistency check #1 failed' chname = 'VBF2' nr = 1 chmoth = 'VBS2' x = 0 y = -3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'VBF2' nr = 2 chmoth = 'VBS2' x = 0 y = 3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 23 natmed = 'VBB2 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBB2' chshap = 'BOX ' nmed = 23 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 181.339 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.23) stop 'consistency check #1 failed' irot = 14 theta1 = 83.03 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.97 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBB2' nr = 1 chmoth = 'VBS2' x = 51.304 y = 0 z = 0 irot = 14 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 24 natmed = 'VBW2 Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBW2' chshap = 'BOX ' nmed = 24 npar = 3 par(1) = 0.05 par(2) = 2.8 par(3) = 182.745 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.24) stop 'consistency check #1 failed' irot = 15 theta1 = 80.06 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.94 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBW2' nr = 1 chmoth = 'VBS2' x = -49.504 y = 0 z = 0 irot = 15 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 25 natmed = 'VBS3 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBS3' chshap = 'TRAP' nmed = 25 npar = 11 par(1) = 230 par(2) = 11.76 par(3) = 180 par(4) = 3.8 par(5) = 45.083 par(6) = 45.083 par(7) = 0 par(8) = 3.8 par(9) = 29.993 par(10) = 29.993 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.25) stop 'consistency check #1 failed' chname = 'VBS3' nr = 1 chmoth = 'AREA' x = -35.953 y = -50 z = 290 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 26 natmed = 'VBF3 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBF3' chshap = 'TRAP' nmed = 26 npar = 11 par(1) = 230 par(2) = 11.76 par(3) = 180 par(4) = 0.5 par(5) = 45.083 par(6) = 45.083 par(7) = 0 par(8) = 0.5 par(9) = 29.993 par(10) = 29.993 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.26) stop 'consistency check #1 failed' chname = 'VBF3' nr = 1 chmoth = 'VBS3' x = 0 y = -3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'VBF3' nr = 2 chmoth = 'VBS3' x = 0 y = 3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 27 natmed = 'VBB3 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBB3' chshap = 'BOX ' nmed = 27 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 236.582 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.27) stop 'consistency check #1 failed' irot = 16 theta1 = 76.45 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 13.55 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBB3' nr = 1 chmoth = 'VBS3' x = 37.038 y = 0 z = 0 irot = 16 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 28 natmed = 'VBW3 Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBW3' chshap = 'BOX ' nmed = 28 npar = 3 par(1) = 0.05 par(2) = 2.8 par(3) = 233.508 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.28) stop 'consistency check #1 failed' irot = 17 theta1 = 80.06 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.94 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBW3' nr = 1 chmoth = 'VBS3' x = -37.238 y = 0 z = 0 irot = 17 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 29 natmed = 'VBS4 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBS4' chshap = 'TRAP' nmed = 29 npar = 11 par(1) = 60.2715 par(2) = 3.55 par(3) = 180 par(4) = 3.8 par(5) = 29.993 par(6) = 29.993 par(7) = 0 par(8) = 3.8 par(9) = 8.761 par(10) = 8.761 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.29) stop 'consistency check #1 failed' chname = 'VBS4' nr = 1 chmoth = 'AREA' x = -87.565 y = -50 z = 580.271 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 30 natmed = 'VBF4 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBF4' chshap = 'TRAP' nmed = 30 npar = 11 par(1) = 60.2715 par(2) = 3.55 par(3) = 180 par(4) = 0.5 par(5) = 29.993 par(6) = 29.993 par(7) = 0 par(8) = 0.5 par(9) = 8.761 par(10) = 8.761 par(11) = 0 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.30) stop 'consistency check #1 failed' chname = 'VBF4' nr = 1 chmoth = 'VBS4' x = 0 y = -3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'VBF4' nr = 2 chmoth = 'VBS4' x = 0 y = 3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 31 natmed = 'VBB4 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBB4' chshap = 'BOX ' nmed = 31 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 61.9585 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.31) stop 'consistency check #1 failed' irot = 18 theta1 = 76.6 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 13.4 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBB4' nr = 1 chmoth = 'VBS4' x = 18.877 y = 0 z = 0 irot = 18 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 32 natmed = 'VBW4 Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBW4' chshap = 'BOX ' nmed = 32 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 60.662 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.32) stop 'consistency check #1 failed' irot = 19 theta1 = 96.53 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 6.53 phi3 = 0 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBW4' nr = 1 chmoth = 'VBS4' x = -19.077 y = 0 z = 0 irot = 19 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 33 natmed = 'VBS5 Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBS5' chshap = 'BOX ' nmed = 33 npar = 3 par(1) = 8.5225 par(2) = 3.8 par(3) = 2.03 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.33) stop 'consistency check #1 failed' irot = 20 theta1 = 76.6 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 13.4 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBS5' nr = 1 chmoth = 'AREA' x = -91.308 y = -50 z = 640.543 irot = 20 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 34 natmed = 'VBF5 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBF5' chshap = 'BOX ' nmed = 34 npar = 3 par(1) = 8.5225 par(2) = 0.5 par(3) = 2.03 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.34) stop 'consistency check #1 failed' chname = 'VBF5' nr = 1 chmoth = 'VBS5' x = 0 y = -3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'VBF5' nr = 2 chmoth = 'VBS5' x = 0 y = 3.3 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 35 natmed = 'VBB5 Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBB5' chshap = 'BOX ' nmed = 35 npar = 3 par(1) = 0.5 par(2) = 2.8 par(3) = 2.03 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.35) stop 'consistency check #1 failed' chname = 'VBB5' nr = 1 chmoth = 'VBS5' x = 8.023 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 36 natmed = 'VBW5 Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBW5' chshap = 'BOX ' nmed = 36 npar = 3 par(1) = 3.5225 par(2) = 2.8 par(3) = 0.5 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.36) stop 'consistency check #1 failed' chname = 'VBW5' nr = 1 chmoth = 'VBS5' x = -5 y = 0 z = 1.53 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 37 natmed = 'VBXF Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBXF' chshap = 'BOX ' nmed = 37 npar = 3 par(1) = 0.5 par(2) = 2.5 par(3) = 500 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.37) stop 'consistency check #1 failed' irot = 21 theta1 = 80.06 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.94 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'VBXF' nr = 1 chmoth = 'AREA' x = -24.238 y = -50 z = 6.0974 irot = 21 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 38 natmed = 'VBXS Air' nmat = 8 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'VBXS' chshap = 'BOX ' nmed = 38 npar = 3 par(1) = 0.5 par(2) = 1.5 par(3) = 490 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.38) stop 'consistency check #1 failed' chname = 'VBXS' nr = 1 chmoth = 'VBXF' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 39 natmed = 'EDUT Iron' nmat = 7 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'EDUT' chshap = 'BOX ' nmed = 39 npar = 3 par(1) = 5 par(2) = 2.5 par(3) = 600 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.39) stop 'consistency check #1 failed' irot = 22 theta1 = 76.6 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 13.4 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) x = -130 y = 50 z = 900 chname = 'EDUT' nr = 1 chmoth = 'HILL' x = x-227.401 y = y-50 z = z+1227 irot = 22 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 40 natmed = 'EDIT Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'EDIT' chshap = 'BOX ' nmed = 40 npar = 3 par(1) = 4 par(2) = 2 par(3) = 600 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.40) stop 'consistency check #1 failed' chname = 'EDIT' nr = 1 chmoth = 'EDUT' x = -0.5 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 41 natmed = 'MSCP Carbon' nmat = 12 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'MSCP' chshap = 'BOX ' nmed = 41 npar = 3 par(1) = 10 par(2) = 5 par(3) = 60 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.41) stop 'consistency check #1 failed' irot = 23 theta1 = 80.056 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.944 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCP' nr = 1 chmoth = 'AREA' x = -15.2456 y = -50 z = -108.783 irot = 23 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 42 natmed = 'MSCI Air' nmat = 8 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'MSCI' chshap = 'BOX ' nmed = 42 npar = 3 par(1) = 9.5 par(2) = 4 par(3) = 59 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.42) stop 'consistency check #1 failed' chname = 'MSCI' nr = 1 chmoth = 'MSCP' x = 0.4 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 43 natmed = 'MSCL Air' nmat = 8 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'MSCL' chshap = 'BOX ' nmed = 43 npar = 3 par(1) = 0.1 par(2) = 0.5 par(3) = 1 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.43) stop 'consistency check #1 failed' irot = 24 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 1 chmoth = 'MSCI' x = 7.2 y = 0 z = -54.5 irot = 24 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'YD01' chmoth = 'MSCL' ndiv = 5 iaxis = 2 step = 0.2 c0 = -0.5 numed = 0 ndvmax = 0 call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax) imate = 13 namate = 'Scintillator' wmat(1) = 9 call gfmate(12,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf) wmat(2) = 10 call gfmate(2,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf) dens = 1.032 nlmat = -2 call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat) itmed = 44 natmed = 'MSFI Scintillator' nmat = 13 isvol = 1 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'MSFI' chshap = 'BOX ' nmed = 44 npar = 3 par(1) = 0.1 par(2) = 0.1 par(3) = 1 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.45) stop 'consistency check #1 failed' chname = 'MSFI' nr = 1 chmoth = 'YD01' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 25 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 2 chmoth = 'MSCI' x = 7.2 y = 0 z = -53.403 irot = 25 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 26 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 3 chmoth = 'MSCI' x = 7.2 y = 0 z = -52.306 irot = 26 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 27 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 4 chmoth = 'MSCI' x = 7.2 y = 0 z = -51.209 irot = 27 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 28 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 5 chmoth = 'MSCI' x = 7.2 y = 0 z = -50.112 irot = 28 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 29 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 6 chmoth = 'MSCI' x = 7.2 y = 0 z = -49.015 irot = 29 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 30 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 7 chmoth = 'MSCI' x = 7.2 y = 0 z = -47.918 irot = 30 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 31 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 8 chmoth = 'MSCI' x = 7.2 y = 0 z = -46.821 irot = 31 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 32 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 9 chmoth = 'MSCI' x = 7.2 y = 0 z = -45.724 irot = 32 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 33 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 10 chmoth = 'MSCI' x = 7.2 y = 0 z = -44.627 irot = 33 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 34 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 11 chmoth = 'MSCI' x = 7.2 y = 0 z = -43.53 irot = 34 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 35 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 12 chmoth = 'MSCI' x = 7.2 y = 0 z = -42.433 irot = 35 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 36 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 13 chmoth = 'MSCI' x = 7.2 y = 0 z = -41.336 irot = 36 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 37 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 14 chmoth = 'MSCI' x = 7.2 y = 0 z = -40.239 irot = 37 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 38 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 15 chmoth = 'MSCI' x = 7.2 y = 0 z = -39.142 irot = 38 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 39 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 16 chmoth = 'MSCI' x = 7.2 y = 0 z = -38.045 irot = 39 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 40 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 17 chmoth = 'MSCI' x = 7.2 y = 0 z = -36.948 irot = 40 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 41 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 18 chmoth = 'MSCI' x = 7.2 y = 0 z = -35.851 irot = 41 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 42 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 19 chmoth = 'MSCI' x = 7.2 y = 0 z = -34.754 irot = 42 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 43 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 20 chmoth = 'MSCI' x = 7.2 y = 0 z = -33.657 irot = 43 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 44 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 21 chmoth = 'MSCI' x = 7.2 y = 0 z = -32.56 irot = 44 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 45 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 22 chmoth = 'MSCI' x = 7.2 y = 0 z = -31.463 irot = 45 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 46 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 23 chmoth = 'MSCI' x = 7.2 y = 0 z = -30.366 irot = 46 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 47 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 24 chmoth = 'MSCI' x = 7.2 y = 0 z = -29.269 irot = 47 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 48 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 25 chmoth = 'MSCI' x = 7.2 y = 0 z = -28.172 irot = 48 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 49 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 26 chmoth = 'MSCI' x = 7.2 y = 0 z = -27.075 irot = 49 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 50 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 27 chmoth = 'MSCI' x = 7.2 y = 0 z = -25.978 irot = 50 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 51 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 28 chmoth = 'MSCI' x = 7.2 y = 0 z = -24.881 irot = 51 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 52 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 29 chmoth = 'MSCI' x = 7.2 y = 0 z = -23.784 irot = 52 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 53 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 30 chmoth = 'MSCI' x = 7.2 y = 0 z = -22.687 irot = 53 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 54 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 31 chmoth = 'MSCI' x = 7.2 y = 0 z = -21.59 irot = 54 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 55 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 32 chmoth = 'MSCI' x = 7.2 y = 0 z = -20.493 irot = 55 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 56 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 33 chmoth = 'MSCI' x = 7.2 y = 0 z = -19.396 irot = 56 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 57 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 34 chmoth = 'MSCI' x = 7.2 y = 0 z = -18.299 irot = 57 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 58 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 35 chmoth = 'MSCI' x = 7.2 y = 0 z = -17.202 irot = 58 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 59 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 36 chmoth = 'MSCI' x = 7.2 y = 0 z = -16.105 irot = 59 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 60 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 37 chmoth = 'MSCI' x = 7.2 y = 0 z = -15.008 irot = 60 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 61 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 38 chmoth = 'MSCI' x = 7.2 y = 0 z = -13.911 irot = 61 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 62 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 39 chmoth = 'MSCI' x = 7.2 y = 0 z = -12.814 irot = 62 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 63 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 40 chmoth = 'MSCI' x = 7.2 y = 0 z = -11.717 irot = 63 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 64 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 41 chmoth = 'MSCI' x = 7.2 y = 0 z = -10.62 irot = 64 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 65 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 42 chmoth = 'MSCI' x = 7.2 y = 0 z = -9.523 irot = 65 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 66 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 43 chmoth = 'MSCI' x = 7.2 y = 0 z = -8.426 irot = 66 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 67 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 44 chmoth = 'MSCI' x = 7.2 y = 0 z = -7.329 irot = 67 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 68 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 45 chmoth = 'MSCI' x = 7.2 y = 0 z = -6.232 irot = 68 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 69 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 46 chmoth = 'MSCI' x = 7.2 y = 0 z = -5.135 irot = 69 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 70 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 47 chmoth = 'MSCI' x = 7.2 y = 0 z = -4.038 irot = 70 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 71 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 48 chmoth = 'MSCI' x = 7.2 y = 0 z = -2.941 irot = 71 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 72 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 49 chmoth = 'MSCI' x = 7.2 y = 0 z = -1.844 irot = 72 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 73 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 50 chmoth = 'MSCI' x = 7.2 y = 0 z = -0.747 irot = 73 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 74 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 51 chmoth = 'MSCI' x = 7.2 y = 0 z = 0.35 irot = 74 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 75 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 52 chmoth = 'MSCI' x = 7.2 y = 0 z = 1.447 irot = 75 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 76 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 53 chmoth = 'MSCI' x = 7.2 y = 0 z = 2.544 irot = 76 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 77 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 54 chmoth = 'MSCI' x = 7.2 y = 0 z = 3.641 irot = 77 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 78 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 55 chmoth = 'MSCI' x = 7.2 y = 0 z = 4.738 irot = 78 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 79 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 56 chmoth = 'MSCI' x = 7.2 y = 0 z = 5.835 irot = 79 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 80 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 57 chmoth = 'MSCI' x = 7.2 y = 0 z = 6.932 irot = 80 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 81 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 58 chmoth = 'MSCI' x = 7.2 y = 0 z = 8.029 irot = 81 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 82 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 59 chmoth = 'MSCI' x = 7.2 y = 0 z = 9.126 irot = 82 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 83 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 60 chmoth = 'MSCI' x = 7.2 y = 0 z = 10.223 irot = 83 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 84 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 61 chmoth = 'MSCI' x = 7.2 y = 0 z = 11.32 irot = 84 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 85 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 62 chmoth = 'MSCI' x = 7.2 y = 0 z = 12.417 irot = 85 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 86 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 63 chmoth = 'MSCI' x = 7.2 y = 0 z = 13.514 irot = 86 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 87 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 64 chmoth = 'MSCI' x = 7.2 y = 0 z = 14.611 irot = 87 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 88 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 65 chmoth = 'MSCI' x = 7.2 y = 0 z = 15.708 irot = 88 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 89 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 66 chmoth = 'MSCI' x = 7.2 y = 0 z = 16.805 irot = 89 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 90 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 67 chmoth = 'MSCI' x = 7.2 y = 0 z = 17.902 irot = 90 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 91 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 68 chmoth = 'MSCI' x = 7.2 y = 0 z = 18.999 irot = 91 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 92 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 69 chmoth = 'MSCI' x = 7.2 y = 0 z = 20.096 irot = 92 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 93 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 70 chmoth = 'MSCI' x = 7.2 y = 0 z = 21.193 irot = 93 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 94 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 71 chmoth = 'MSCI' x = 7.2 y = 0 z = 22.29 irot = 94 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 95 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 72 chmoth = 'MSCI' x = 7.2 y = 0 z = 23.387 irot = 95 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 96 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 73 chmoth = 'MSCI' x = 7.2 y = 0 z = 24.484 irot = 96 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 97 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 74 chmoth = 'MSCI' x = 7.2 y = 0 z = 25.581 irot = 97 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 98 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 75 chmoth = 'MSCI' x = 7.2 y = 0 z = 26.678 irot = 98 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 99 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 76 chmoth = 'MSCI' x = 7.2 y = 0 z = 27.775 irot = 99 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 100 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 77 chmoth = 'MSCI' x = 7.2 y = 0 z = 28.872 irot = 100 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 101 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 78 chmoth = 'MSCI' x = 7.2 y = 0 z = 29.969 irot = 101 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 102 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 79 chmoth = 'MSCI' x = 7.2 y = 0 z = 31.066 irot = 102 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 103 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 80 chmoth = 'MSCI' x = 7.2 y = 0 z = 32.163 irot = 103 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 104 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 81 chmoth = 'MSCI' x = 7.2 y = 0 z = 33.26 irot = 104 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 105 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 82 chmoth = 'MSCI' x = 7.2 y = 0 z = 34.357 irot = 105 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 106 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 83 chmoth = 'MSCI' x = 7.2 y = 0 z = 35.454 irot = 106 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 107 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 84 chmoth = 'MSCI' x = 7.2 y = 0 z = 36.551 irot = 107 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 108 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 85 chmoth = 'MSCI' x = 7.2 y = 0 z = 37.648 irot = 108 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 109 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 86 chmoth = 'MSCI' x = 7.2 y = 0 z = 38.745 irot = 109 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 110 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 87 chmoth = 'MSCI' x = 7.2 y = 0 z = 39.842 irot = 110 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 111 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 88 chmoth = 'MSCI' x = 7.2 y = 0 z = 40.939 irot = 111 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 112 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 89 chmoth = 'MSCI' x = 7.2 y = 0 z = 42.036 irot = 112 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 113 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 90 chmoth = 'MSCI' x = 7.2 y = 0 z = 43.133 irot = 113 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 114 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 91 chmoth = 'MSCI' x = 7.2 y = 0 z = 44.23 irot = 114 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 115 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 92 chmoth = 'MSCI' x = 7.2 y = 0 z = 45.327 irot = 115 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 116 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 93 chmoth = 'MSCI' x = 7.2 y = 0 z = 46.424 irot = 116 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 117 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 94 chmoth = 'MSCI' x = 7.2 y = 0 z = 47.521 irot = 117 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 118 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 95 chmoth = 'MSCI' x = 7.2 y = 0 z = 48.618 irot = 118 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 119 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 96 chmoth = 'MSCI' x = 7.2 y = 0 z = 49.715 irot = 119 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 120 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 97 chmoth = 'MSCI' x = 7.2 y = 0 z = 50.812 irot = 120 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 121 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 98 chmoth = 'MSCI' x = 7.2 y = 0 z = 51.909 irot = 121 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 122 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 99 chmoth = 'MSCI' x = 7.2 y = 0 z = 53.006 irot = 122 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) irot = 123 theta1 = 79.5 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 10.5 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'MSCL' nr = 100 chmoth = 'MSCI' x = 7.2 y = 0 z = 54.103 irot = 123 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 45 natmed = 'FXSC Air' nmat = 8 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'FXSC' chshap = 'BOX ' nmed = 45 npar = 3 par(1) = 0.25 par(2) = 0.5 par(3) = 489.75 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.46) stop 'consistency check #1 failed' irot = 124 theta1 = 80.056 phi1 = 0 theta2 = 90 phi2 = 90 theta3 = 9.944 phi3 = 180 call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3) chname = 'FXSC' nr = 1 chmoth = 'AREA' x = -48.8734 y = -50 z = 1.95354 irot = 124 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'ZD01' chmoth = 'FXSC' ndiv = 141 iaxis = 3 step = 6.947 c0 = -489.274 numed = 0 ndvmax = 0 call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax) itmed = 46 natmed = 'FXCH Scintillator' nmat = 13 isvol = 1 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'FXCH' chshap = 'BOX ' nmed = 46 npar = 3 par(1) = 0.25 par(2) = 0.5 par(3) = 3.45 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.48) stop 'consistency check #1 failed' chname = 'FXCH' nr = 1 chmoth = 'ZD01' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 47 natmed = 'PITU Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'PITU' chshap = 'TUBE' nmed = 47 npar = 3 par(1) = 0 par(2) = 2 par(3) = 28.1 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.49) stop 'consistency check #1 failed' chname = 'PITU' nr = 1 chmoth = 'AREA' x = 130 y = -50 z = -272.085 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 48 natmed = 'PITW Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'PITW' chshap = 'TUBE' nmed = 48 npar = 3 par(1) = 1.8 par(2) = 2 par(3) = 28.1 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.50) stop 'consistency check #1 failed' chname = 'PITW' nr = 1 chmoth = 'PITU' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 49 natmed = 'PXTU Vacuum' nmat = 11 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'PXTU' chshap = 'TUBE' nmed = 49 npar = 3 par(1) = 0 par(2) = 2 par(3) = 1100 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.51) stop 'consistency check #1 failed' chname = 'PXTU' nr = 1 chmoth = 'AREA' x = 130 y = -50 z = 1164.45 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 50 natmed = 'PXTW Aluminum' nmat = 4 isvol = 0 ifield = 1 fieldm = 18 tmaxfd = 1 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'PXTW' chshap = 'TUBE' nmed = 50 npar = 3 par(1) = 1.8 par(2) = 2 par(3) = 1100 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.52) stop 'consistency check #1 failed' chname = 'PXTW' nr = 1 chmoth = 'PXTU' x = 0 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) itmed = 51 natmed = 'TUNL Air' nmat = 8 isvol = 0 ifield = 0 fieldm = 0 tmaxfd = 0 stemax = 0 deemax = 0 epsil = 0.001 stmin = 0 nwbuf = 0 call gstmed(itmed,natmed,nmat,isvol,ifield,fieldm,tmaxfd, + stemax,deemax,epsil,stmin,ubuf,nwbuf) chname = 'TUNL' chshap = 'BOX ' nmed = 51 npar = 3 par(1) = 280 par(2) = 225 par(3) = 250 call gsvolu(chname,chshap,nmed,par,npar,ivolu) if (ivolu.ne.53) stop 'consistency check #1 failed' chname = 'TUNL' nr = 1 chmoth = 'HILL' x = 0 y = 50 z = 2485 - (610-250) irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) C Trick to make the volume stick outside AREA chname = 'EDUV' chshap = 'BOX ' nmed = 39 npar = 3 par(1) = 5 par(2) = 2.5 par(3) = 600 call gsvolu(chname,chshap,nmed,par,npar,ivolu) chname = 'EDUV' nr = 1 chmoth = 'AREA' x = -227.401 y = -50 z = 1227 irot = 22 chonly = 'MANY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) chname = 'EDIV' chshap = 'BOX ' nmed = 40 npar = 3 par(1) = 4 par(2) = 2 par(3) = 600 call gsvolu(chname,chshap,nmed,par,npar,ivolu) chname = 'EDIV' nr = 1 chmoth = 'EDUV' x = -0.5 y = 0 z = 0 irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) C C -- Control Volume for Tagger Rates (1m^2 by 1 cm) chname = 'CVTG' chshap = 'BOX ' nmed = 40 npar = 3 par(1) = 0.5 par(2) = 50. par(3) = 50. call gsvolu(chname,chshap,nmed,par,npar,ivolu) chname = 'CVTG' nr = 1 chmoth = 'AREA' x = -145. y = -50. z = 550. irot = 0 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) C C -- Control Volume for Tagger Rates (11m^2 by 1 cm) chname = 'CVT1' chshap = 'BOX ' nmed = 40 npar = 3 par(1) = 0.5 par(2) = 50. par(3) = 550. call gsvolu(chname,chshap,nmed,par,npar,ivolu) chname = 'CVT1' nr = 1 chmoth = 'AREA' x = -50. y = -50. z = -42. irot = 13 chonly = 'ONLY' call gspos(chname,nr,chmoth,x,y,z,irot,chonly) end function getColumn() implicit none integer getColumn integer nlevel,names,number,lvolum common /gcvolu/nlevel,names(15),number(15),lvolum(15) integer i,istart(53) data (istart(i),i=1,53) / + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 1, 0, 0, 0, 101, 0, 0, 0, + 0, 0, 0/ integer lookup(241) data (lookup(i),i=1,100) / + 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, + 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, + 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, + 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, + 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, + 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, + 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, + 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, + 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, + 91, 92, 93, 94, 95, 96, 97, 98, 99, 100/ data (lookup(i),i=101,200) / + 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, + 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, + 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, + 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, + 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, + 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, + 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, + 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, + 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, + 91, 92, 93, 94, 95, 96, 97, 98, 99, 100/ data (lookup(i),i=201,241) / + 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, + 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, + 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, + 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, + 141/ integer level,index integer column getColumn = 0 do level=1,nlevel index = istart(lvolum(level)) if (index.gt.0) then column = lookup(index + number(level) - 1) if (column.gt.0) then getColumn = column endif endif enddo end function getMap() implicit none integer getMap integer nlevel,names,number,lvolum common /gcvolu/nlevel,names(15),number(15),lvolum(15) integer i,istart(53) data (istart(i),i=1,53) / + 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0/ integer lookup(2) data (lookup(i),i=1,2) / + 1, 0/ integer level,index integer map getMap = 0 do level=1,nlevel index = istart(lvolum(level)) if (index.gt.0) then map = lookup(index + number(level) - 1) if (map.gt.0) then getMap = map endif endif enddo end function getRow() implicit none integer getRow integer nlevel,names,number,lvolum common /gcvolu/nlevel,names(15),number(15),lvolum(15) integer i,istart(53) data (istart(i),i=1,53) / + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, + 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, + 0, 0, 0/ integer lookup(5) data (lookup(i),i=1,5) / + 1, 2, 3, 4, 5/ integer level,index integer row getRow = 0 do level=1,nlevel index = istart(lvolum(level)) if (index.gt.0) then row = lookup(index + number(level) - 1) if (row.gt.0) then getRow = row endif endif enddo end subroutine gufld_orig(r,B) implicit none real r(3),B(3) real rr(3),rs(3),BB(3) integer iregion integer getMap external getMap real orig1(3),rmat1(3,3) data orig1/0,0,0/ data rmat1/1,0,0, + 0,1,0, + 0,0,1/ real orig2(3),rmat2(3,3) data orig2/0,0,166.015/ data rmat2/1,0,0, + 0,1,0, + 0,0,1/ iregion = getMap() if (iregion.eq.0) then B(1) = 0 B(2) = 0 B(3) = 0 else if (iregion.eq.1) then rs(1) = r(1)-orig1(1) rs(2) = r(2)-orig1(2) rs(3) = r(3)-orig1(3) rr(1) = rs(1)*rmat1(1,1)+rs(2)*rmat1(1,2)+rs(3)*rmat1(1,3) rr(2) = rs(1)*rmat1(2,1)+rs(2)*rmat1(2,2)+rs(3)*rmat1(2,3) rr(3) = rs(1)*rmat1(3,1)+rs(2)*rmat1(3,2)+rs(3)*rmat1(3,3) call gufld1(rr,BB) B(1) = BB(1)*rmat1(1,1)+BB(2)*rmat1(2,1)+BB(3)*rmat1(3,1) B(2) = BB(1)*rmat1(1,2)+BB(2)*rmat1(2,2)+BB(3)*rmat1(3,2) B(3) = BB(1)*rmat1(1,3)+BB(2)*rmat1(2,3)+BB(3)*rmat1(3,3) B(1) = B(1)*0.001 B(2) = B(2)*0.001 B(3) = B(3)*0.001 endif end subroutine gufld1(r,B) implicit none real r(3),B(3) real rho,phi,alpha real u(3) real twopi parameter (twopi=6.28318530717959) real bound0(3,2) data bound0/10,0,-10, + -340,-1.5,1600/ integer reverse0(3) data reverse0/-1,1,-1/ real bound1(3,2) data bound1/10,0,-10, + -340,1.5,1600/ integer reverse1(3) data reverse1/1,1,1/ real Bmap(3,351,16,1611) integer nsites(3) data nsites/351,16,1611/ logical loaded data loaded/.false./ save Bmap,nsites,loaded integer i,i1,i2,i3 if (.not.loaded) then open(unit=78,file='taggerBfield-quad.map',status='old',err=7) read(unit=78,fmt=*,err=5,end=6) + ((((Bmap(i,i1,i2,i3),i=1,3), + i2=1,16), + i3=1,1611), + i1=1,351) go to 8 5 stop 'error reading magnetic field map, stop' 6 stop 'EOF encountered reading magnetic field map, stop' 7 stop 'error opening magnetic field map, stop' 8 loaded=.true. endif u(1) = (r(1)-bound0(1,1))/(bound0(1,2)-bound0(1,1)) u(2) = (r(2)-bound0(2,1))/(bound0(2,2)-bound0(2,1)) u(3) = (r(3)-bound0(3,1))/(bound0(3,2)-bound0(3,1)) if ((u(1).ge.0.and.u(1).le.1).and. + (u(2).ge.0.and.u(2).le.1).and. + (u(3).ge.0.and.u(3).le.1)) then call interpol3(Bmap,nsites,u,B) B(1)=B(1)*reverse0(1) B(2)=B(2)*reverse0(2) B(3)=B(3)*reverse0(3) return endif u(1) = (r(1)-bound1(1,1))/(bound1(1,2)-bound1(1,1)) u(2) = (r(2)-bound1(2,1))/(bound1(2,2)-bound1(2,1)) u(3) = (r(3)-bound1(3,1))/(bound1(3,2)-bound1(3,1)) if ((u(1).ge.0.and.u(1).le.1).and. + (u(2).ge.0.and.u(2).le.1).and. + (u(3).ge.0.and.u(3).le.1)) then call interpol3(Bmap,nsites,u,B) B(1)=B(1)*reverse1(1) B(2)=B(2)*reverse1(2) B(3)=B(3)*reverse1(3) return endif B(1) = 0 B(2) = 0 B(3) = 0 return end subroutine interpol3(Bmap,nsites,u,B) implicit none integer nsites(3) real Bmap(3,nsites(1),nsites(2),nsites(3)) real u(3),B(3) integer ir(3),ir0(3),ir1(3) real ur(3),dur(3),ugrad(3,3) integer i do i=1,3 ur(i)=u(i)*(nsites(i)-1)+1 ir(i)=nint(ur(i)) ir0(i)=max(ir(i)-1,1) ir1(i)=min(ir(i)+1,nsites(i)) dur(i)=(ur(i)-ir(i))/(ir1(i)-ir0(i)+1e-20) enddo ugrad(1,1)=(Bmap(1,ir1(1),ir(2),ir(3))-Bmap(1,ir0(1),ir(2),ir(3))) ugrad(2,1)=(Bmap(2,ir1(1),ir(2),ir(3))-Bmap(2,ir0(1),ir(2),ir(3))) ugrad(3,1)=(Bmap(3,ir1(1),ir(2),ir(3))-Bmap(3,ir0(1),ir(2),ir(3))) ugrad(1,2)=(Bmap(1,ir(1),ir1(2),ir(3))-Bmap(1,ir(1),ir0(2),ir(3))) ugrad(2,2)=(Bmap(2,ir(1),ir1(2),ir(3))-Bmap(2,ir(1),ir0(2),ir(3))) ugrad(3,2)=(Bmap(3,ir(1),ir1(2),ir(3))-Bmap(3,ir(1),ir0(2),ir(3))) ugrad(1,3)=(Bmap(1,ir(1),ir(2),ir1(3))-Bmap(1,ir(1),ir(2),ir0(3))) ugrad(2,3)=(Bmap(2,ir(1),ir(2),ir1(3))-Bmap(2,ir(1),ir(2),ir0(3))) ugrad(3,3)=(Bmap(3,ir(1),ir(2),ir1(3))-Bmap(3,ir(1),ir(2),ir0(3))) B(1)=Bmap(1,ir(1),ir(2),ir(3)) + +ugrad(1,1)*dur(1)+ugrad(1,2)*dur(2)+ugrad(1,3)*dur(3) B(2)=Bmap(2,ir(1),ir(2),ir(3)) + +ugrad(2,1)*dur(1)+ugrad(2,2)*dur(2)+ugrad(2,3)*dur(3) B(3)=Bmap(3,ir(1),ir(2),ir(3)) + +ugrad(3,1)*dur(1)+ugrad(3,2)*dur(2)+ugrad(3,3)*dur(3) end