*
* 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 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 = '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 = 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 = '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 = 5
      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 = 1
      namate = 'Air'
      wmat(1) = 0.7494
      call gfmate(2,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(4,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      wmat(4) = 0.0008
      call gfmate(5,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 = 1
      natmed = 'SITE Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'SITE'
      chshap = 'BOX '
      nmed = 1
      npar = 3
      par(1) = 2500
      par(2) = 2500
      par(3) = 2500
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.1) stop 'consistency check #1 failed'

      imate = 7
      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 = 8
      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 = 9
      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 = 10
      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 = 6
      namate = 'Concrete'
      wmat(1) = 0.004
      call gfmate(5,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(7,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      wmat(4) = 0.07
      call gfmate(8,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf)
      wmat(5) = 0.345
      call gfmate(9,chnama,amat(5),zmat(5),dens,radl,absl,ubuf,nwbuf)
      wmat(6) = 0.038
      call gfmate(10,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 = 2
      natmed = 'SHLD Concrete'
      nmat = 6
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'SHLD'
      chshap = 'BOX '
      nmed = 2
      npar = 3
      par(1) = 275
      par(2) = 275
      par(3) = 625
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.2) stop 'consistency check #1 failed'

      chname = 'SHLD'
      nr = 1
      chmoth = 'SITE'
      x = 0
      y = 0
      z = -1625
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 3
      natmed = 'CAVE Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'CAVE'
      chshap = 'BOX '
      nmed = 3
      npar = 3
      par(1) = 225
      par(2) = 135
      par(3) = 600
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.3) stop 'consistency check #1 failed'

      chname = 'CAVE'
      nr = 1
      chmoth = 'SHLD'
      x = 0
      y = 35
      z = 25
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 4
      natmed = 'OSHO Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OSHO'
      chshap = 'TUBE'
      nmed = 4
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 17.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.4) stop 'consistency check #1 failed'

      chname = 'OSHO'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 567.5
      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 = 5
      natmed = 'ISHO Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ISHO'
      chshap = 'TUBE'
      nmed = 5
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 17.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.5) stop 'consistency check #1 failed'

      chname = 'ISHO'
      nr = 1
      chmoth = 'OSHO'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 13
      namate = 'Boron'
      a = 10.8
      z = 5
      dens = 2.34
      nlmat = 1
      wmat(1) = 1
      call gsmixt(imate,namate,a,z,dens,nlmat,wmat)

      imate = 12
      namate = 'BoronNitride'
      wmat(1) = 1
      call gfmate(13,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 1
      call gfmate(2,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      dens = 1.9
      nlmat = -2
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 6
      natmed = 'INSU BoronNitride'
      nmat = 12
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'INSU'
      chshap = 'TUBE'
      nmed = 6
      npar = 3
      par(1) = 0.25
      par(2) = 7.36
      par(3) = 2.1
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.6) stop 'consistency check #1 failed'

      chname = 'INSU'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -452.2
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 7
      natmed = 'HOUS Aluminum'
      nmat = 7
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'HOUS'
      chshap = 'TUBE'
      nmed = 7
      npar = 3
      par(1) = 0.25
      par(2) = 6.7
      par(3) = 1.6
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.7) stop 'consistency check #1 failed'

      chname = 'HOUS'
      nr = 1
      chmoth = 'INSU'
      x = 0
      y = 0
      z = 0.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 8
      natmed = 'AIRH Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIRH'
      chshap = 'TUBE'
      nmed = 8
      npar = 3
      par(1) = 0.25
      par(2) = 6.5
      par(3) = 1.45
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.8) stop 'consistency check #1 failed'

      chname = 'AIRH'
      nr = 1
      chmoth = 'HOUS'
      x = 0
      y = 0
      z = -0.15
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 9
      natmed = 'DIV1 Aluminum'
      nmat = 7
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'DIV1'
      chshap = 'BOX '
      nmed = 9
      npar = 3
      par(1) = 3.125
      par(2) = 0.05
      par(3) = 1.45
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.9) stop 'consistency check #1 failed'

      chname = 'DIV1'
      nr = 1
      chmoth = 'AIRH'
      x = 3.375
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 1
      theta1 = 90
      phi1 = 90
      theta2 = 90
      phi2 = 180
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'DIV1'
      nr = 2
      chmoth = 'AIRH'
      x = 2.06652e-16
      y = 3.375
      z = 0
      irot = 1
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 2
      theta1 = 90
      phi1 = 180
      theta2 = 90
      phi2 = -90
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'DIV1'
      nr = 3
      chmoth = 'AIRH'
      x = -3.375
      y = 4.13305e-16
      z = 0
      irot = 2
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 3
      theta1 = 90
      phi1 = -90
      theta2 = 90
      phi2 = -1.05247e-14
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'DIV1'
      nr = 4
      chmoth = 'AIRH'
      x = -6.19957e-16
      y = -3.375
      z = 0
      irot = 3
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 10
      natmed = 'DIV2 Aluminum'
      nmat = 7
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'DIV2'
      chshap = 'TUBE'
      nmed = 10
      npar = 3
      par(1) = 2.7
      par(2) = 2.8
      par(3) = 1.45
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.10) stop 'consistency check #1 failed'

      chname = 'DIV2'
      nr = 1
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 11
      natmed = 'ACWI Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ACWI'
      chshap = 'TUBS'
      nmed = 11
      npar = 5
      par(1) = 0.25
      par(2) = 2.5
      par(3) = 1.25
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.11) stop 'consistency check #1 failed'

      irot = 4
      theta1 = 90
      phi1 = 45
      theta2 = 90
      phi2 = 135
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWI'
      nr = 1
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 4
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 15
      chnama = 'Tungsten'
      a = 183.85
      z = 74
      dens = 19.3
      radl = 0.350259
      absl = 9.58549
      nwbuf = 0
      call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf)

      imate = 16
      chnama = 'Copper'
      a = 63.546
      z = 29
      dens = 8.96
      radl = 1.43527
      absl = 15.0558
      nwbuf = 0
      call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf)

      imate = 14
      namate = 'SoftTungsten'
      wmat(1) = 0.95
      call gfmate(15,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 0.05
      call gfmate(16,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      dens = 18.2471
      nlmat = 2
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 12
      natmed = 'ACBI SoftTungsten'
      nmat = 14
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ACBI'
      chshap = 'TUBS'
      nmed = 12
      npar = 5
      par(1) = 0.25
      par(2) = 2.5
      par(3) = 0.4
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.12) stop 'consistency check #1 failed'

      chname = 'ACBI'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = -0.85
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 13
      natmed = 'AIR1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR1'
      chshap = 'TUBS'
      nmed = 13
      npar = 5
      par(1) = 0.25
      par(2) = 0.33
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.13) stop 'consistency check #1 failed'

      chname = 'AIR1'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD01'
      chmoth = 'AIR1'
      ndiv = 3
      iaxis = 2
      step = 15
      c0 = -25
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 14
      natmed = 'AIS1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS1'
      chshap = 'TUBS'
      nmed = 14
      npar = 5
      par(1) = 0.25
      par(2) = 0.33
      par(3) = 0.85
      par(4) = -10
      par(5) = 10
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.15) stop 'consistency check #1 failed'

      irot = 5
      theta1 = 90
      phi1 = 2.5
      theta2 = 90
      phi2 = 92.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS1'
      nr = 1
      chmoth = 'SD01'
      x = 0
      y = 0
      z = 0
      irot = 5
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 15
      natmed = 'PIN1 SoftTungsten'
      nmat = 14
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'PIN1'
      chshap = 'BOX '
      nmed = 15
      npar = 3
      par(1) = 0.02
      par(2) = 0.02
      par(3) = 0.85
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.16) stop 'consistency check #1 failed'

      chname = 'PIN1'
      nr = 1
      chmoth = 'AIS1'
      x = 0.3
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 16
      natmed = 'AIR2 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR2'
      chshap = 'TUBS'
      nmed = 16
      npar = 5
      par(1) = 0.33
      par(2) = 0.42
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.17) stop 'consistency check #1 failed'

      chname = 'AIR2'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD02'
      chmoth = 'AIR2'
      ndiv = 4
      iaxis = 2
      step = 12
      c0 = -25.5
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 17
      natmed = 'AIS2 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS2'
      chshap = 'TUBS'
      nmed = 17
      npar = 5
      par(1) = 0.33
      par(2) = 0.42
      par(3) = 0.85
      par(4) = -7.5
      par(5) = 7.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.19) stop 'consistency check #1 failed'

      irot = 6
      theta1 = 90
      phi1 = 1.5
      theta2 = 90
      phi2 = 91.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS2'
      nr = 1
      chmoth = 'SD02'
      x = 0
      y = 0
      z = 0
      irot = 6
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 2
      chmoth = 'AIS2'
      x = 0.375
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 18
      natmed = 'AIR3 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR3'
      chshap = 'TUBS'
      nmed = 18
      npar = 5
      par(1) = 0.42
      par(2) = 0.49
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.20) stop 'consistency check #1 failed'

      chname = 'AIR3'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD03'
      chmoth = 'AIR3'
      ndiv = 5
      iaxis = 2
      step = 10
      c0 = -26
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 19
      natmed = 'AIS3 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS3'
      chshap = 'TUBS'
      nmed = 19
      npar = 5
      par(1) = 0.42
      par(2) = 0.49
      par(3) = 0.85
      par(4) = -6
      par(5) = 6
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.22) stop 'consistency check #1 failed'

      irot = 7
      theta1 = 90
      phi1 = 1
      theta2 = 90
      phi2 = 91
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS3'
      nr = 1
      chmoth = 'SD03'
      x = 0
      y = 0
      z = 0
      irot = 7
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 3
      chmoth = 'AIS3'
      x = 0.45
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 20
      natmed = 'AIR4 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR4'
      chshap = 'TUBS'
      nmed = 20
      npar = 5
      par(1) = 0.49
      par(2) = 0.56
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.23) stop 'consistency check #1 failed'

      chname = 'AIR4'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD04'
      chmoth = 'AIR4'
      ndiv = 6
      iaxis = 2
      step = 8.57
      c0 = -26.43
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 21
      natmed = 'AIS4 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS4'
      chshap = 'TUBS'
      nmed = 21
      npar = 5
      par(1) = 0.49
      par(2) = 0.56
      par(3) = 0.85
      par(4) = -5
      par(5) = 5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.25) stop 'consistency check #1 failed'

      irot = 8
      theta1 = 90
      phi1 = 0.715
      theta2 = 90
      phi2 = 90.715
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS4'
      nr = 1
      chmoth = 'SD04'
      x = 0
      y = 0
      z = 0
      irot = 8
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 4
      chmoth = 'AIS4'
      x = 0.525
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 22
      natmed = 'AIR5 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR5'
      chshap = 'TUBS'
      nmed = 22
      npar = 5
      par(1) = 0.56
      par(2) = 0.64
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.26) stop 'consistency check #1 failed'

      chname = 'AIR5'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD05'
      chmoth = 'AIR5'
      ndiv = 7
      iaxis = 2
      step = 7.5
      c0 = -26.79
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 23
      natmed = 'AIS5 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS5'
      chshap = 'TUBS'
      nmed = 23
      npar = 5
      par(1) = 0.56
      par(2) = 0.64
      par(3) = 0.85
      par(4) = -4.29
      par(5) = 4.28
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.28) stop 'consistency check #1 failed'

      irot = 9
      theta1 = 90
      phi1 = 0.54
      theta2 = 90
      phi2 = 90.54
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS5'
      nr = 1
      chmoth = 'SD05'
      x = 0
      y = 0
      z = 0
      irot = 9
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 5
      chmoth = 'AIS5'
      x = 0.6
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 24
      natmed = 'AIR6 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR6'
      chshap = 'TUBS'
      nmed = 24
      npar = 5
      par(1) = 0.64
      par(2) = 0.7
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.29) stop 'consistency check #1 failed'

      chname = 'AIR6'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD06'
      chmoth = 'AIR6'
      ndiv = 8
      iaxis = 2
      step = 6.67
      c0 = -27.08
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 25
      natmed = 'AIS6 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS6'
      chshap = 'TUBS'
      nmed = 25
      npar = 5
      par(1) = 0.64
      par(2) = 0.7
      par(3) = 0.85
      par(4) = -3.75
      par(5) = 3.75
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.31) stop 'consistency check #1 failed'

      irot = 10
      theta1 = 90
      phi1 = 0.415
      theta2 = 90
      phi2 = 90.415
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS6'
      nr = 1
      chmoth = 'SD06'
      x = 0
      y = 0
      z = 0
      irot = 10
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 6
      chmoth = 'AIS6'
      x = 0.675
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 26
      natmed = 'AIR7 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR7'
      chshap = 'TUBS'
      nmed = 26
      npar = 5
      par(1) = 0.7
      par(2) = 0.78
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.32) stop 'consistency check #1 failed'

      chname = 'AIR7'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD07'
      chmoth = 'AIR7'
      ndiv = 9
      iaxis = 2
      step = 6
      c0 = -27.33
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 27
      natmed = 'AIS7 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS7'
      chshap = 'TUBS'
      nmed = 27
      npar = 5
      par(1) = 0.7
      par(2) = 0.78
      par(3) = 0.85
      par(4) = -3.33
      par(5) = 3.34
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.34) stop 'consistency check #1 failed'

      irot = 11
      theta1 = 90
      phi1 = 0.33
      theta2 = 90
      phi2 = 90.33
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS7'
      nr = 1
      chmoth = 'SD07'
      x = 0
      y = 0
      z = 0
      irot = 11
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 7
      chmoth = 'AIS7'
      x = 0.75
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 28
      natmed = 'AIR8 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIR8'
      chshap = 'TUBS'
      nmed = 28
      npar = 5
      par(1) = 0.78
      par(2) = 0.85
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.35) stop 'consistency check #1 failed'

      chname = 'AIR8'
      nr = 1
      chmoth = 'ACWI'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD08'
      chmoth = 'AIR8'
      ndiv = 10
      iaxis = 2
      step = 5.45
      c0 = -27.55
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 29
      natmed = 'AIS8 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AIS8'
      chshap = 'TUBS'
      nmed = 29
      npar = 5
      par(1) = 0.78
      par(2) = 0.85
      par(3) = 0.85
      par(4) = -3
      par(5) = 3
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.37) stop 'consistency check #1 failed'

      irot = 12
      theta1 = 90
      phi1 = 0.275
      theta2 = 90
      phi2 = 90.275
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'AIS8'
      nr = 1
      chmoth = 'SD08'
      x = 0
      y = 0
      z = 0
      irot = 12
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 8
      chmoth = 'AIS8'
      x = 0.825
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 13
      theta1 = 90
      phi1 = 135
      theta2 = 90
      phi2 = -135
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWI'
      nr = 2
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 13
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 14
      theta1 = 90
      phi1 = -135
      theta2 = 90
      phi2 = -45
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWI'
      nr = 3
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 14
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 15
      theta1 = 90
      phi1 = -45
      theta2 = 90
      phi2 = 45
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWI'
      nr = 4
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 15
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 30
      natmed = 'ACWO Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ACWO'
      chshap = 'TUBS'
      nmed = 30
      npar = 5
      par(1) = 3
      par(2) = 6
      par(3) = 1.25
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.38) stop 'consistency check #1 failed'

      irot = 16
      theta1 = 90
      phi1 = 45
      theta2 = 90
      phi2 = 135
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWO'
      nr = 1
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 16
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 31
      natmed = 'ACBO SoftTungsten'
      nmat = 14
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ACBO'
      chshap = 'TUBS'
      nmed = 31
      npar = 5
      par(1) = 3
      par(2) = 6
      par(3) = 0.4
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.39) stop 'consistency check #1 failed'

      chname = 'ACBO'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = -0.85
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 32
      natmed = 'AOR1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR1'
      chshap = 'TUBS'
      nmed = 32
      npar = 5
      par(1) = 3
      par(2) = 3.08
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.40) stop 'consistency check #1 failed'

      chname = 'AOR1'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD09'
      chmoth = 'AOR1'
      ndiv = 40
      iaxis = 2
      step = 1.463
      c0 = -29.392
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 33
      natmed = 'AOS1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS1'
      chshap = 'TUBS'
      nmed = 33
      npar = 5
      par(1) = 3
      par(2) = 3.08
      par(3) = 0.85
      par(4) = -0.732
      par(5) = 0.731
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.42) stop 'consistency check #1 failed'

      chname = 'AOS1'
      nr = 1
      chmoth = 'SD09'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 9
      chmoth = 'AOS1'
      x = 3.05
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 34
      natmed = 'AOR2 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR2'
      chshap = 'TUBS'
      nmed = 34
      npar = 5
      par(1) = 3.08
      par(2) = 3.16
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.43) stop 'consistency check #1 failed'

      chname = 'AOR2'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0A'
      chmoth = 'AOR2'
      ndiv = 41
      iaxis = 2
      step = 1.429
      c0 = -29.383
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 35
      natmed = 'AOS2 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS2'
      chshap = 'TUBS'
      nmed = 35
      npar = 5
      par(1) = 3.08
      par(2) = 3.16
      par(3) = 0.85
      par(4) = -0.715
      par(5) = 0.714
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.45) stop 'consistency check #1 failed'

      chname = 'AOS2'
      nr = 1
      chmoth = 'SD0A'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 10
      chmoth = 'AOS2'
      x = 3.125
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 36
      natmed = 'AOR3 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR3'
      chshap = 'TUBS'
      nmed = 36
      npar = 5
      par(1) = 3.16
      par(2) = 3.24
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.46) stop 'consistency check #1 failed'

      chname = 'AOR3'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0B'
      chmoth = 'AOR3'
      ndiv = 42
      iaxis = 2
      step = 1.395
      c0 = -29.384
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 37
      natmed = 'AOS3 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS3'
      chshap = 'TUBS'
      nmed = 37
      npar = 5
      par(1) = 3.16
      par(2) = 3.24
      par(3) = 0.85
      par(4) = -0.698
      par(5) = 0.697
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.48) stop 'consistency check #1 failed'

      chname = 'AOS3'
      nr = 1
      chmoth = 'SD0B'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 11
      chmoth = 'AOS3'
      x = 3.2
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 38
      natmed = 'AOR4 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR4'
      chshap = 'TUBS'
      nmed = 38
      npar = 5
      par(1) = 3.24
      par(2) = 3.32
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.49) stop 'consistency check #1 failed'

      chname = 'AOR4'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0C'
      chmoth = 'AOR4'
      ndiv = 43
      iaxis = 2
      step = 1.364
      c0 = -29.385
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 39
      natmed = 'AOS4 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS4'
      chshap = 'TUBS'
      nmed = 39
      npar = 5
      par(1) = 3.24
      par(2) = 3.32
      par(3) = 0.85
      par(4) = -0.682
      par(5) = 0.682
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.51) stop 'consistency check #1 failed'

      chname = 'AOS4'
      nr = 1
      chmoth = 'SD0C'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 12
      chmoth = 'AOS4'
      x = 3.275
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 40
      natmed = 'AOR5 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR5'
      chshap = 'TUBS'
      nmed = 40
      npar = 5
      par(1) = 3.32
      par(2) = 3.38
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.52) stop 'consistency check #1 failed'

      chname = 'AOR5'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0D'
      chmoth = 'AOR5'
      ndiv = 44
      iaxis = 2
      step = 1.333
      c0 = -29.385
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 41
      natmed = 'AOS5 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS5'
      chshap = 'TUBS'
      nmed = 41
      npar = 5
      par(1) = 3.32
      par(2) = 3.38
      par(3) = 0.85
      par(4) = -0.667
      par(5) = 0.666
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.54) stop 'consistency check #1 failed'

      chname = 'AOS5'
      nr = 1
      chmoth = 'SD0D'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 13
      chmoth = 'AOS5'
      x = 3.35
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 42
      natmed = 'AOR6 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR6'
      chshap = 'TUBS'
      nmed = 42
      npar = 5
      par(1) = 3.38
      par(2) = 3.46
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.55) stop 'consistency check #1 failed'

      chname = 'AOR6'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0E'
      chmoth = 'AOR6'
      ndiv = 45
      iaxis = 2
      step = 1.304
      c0 = -29.385
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 43
      natmed = 'AOS6 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS6'
      chshap = 'TUBS'
      nmed = 43
      npar = 5
      par(1) = 3.38
      par(2) = 3.46
      par(3) = 0.85
      par(4) = -0.652
      par(5) = 0.652
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.57) stop 'consistency check #1 failed'

      chname = 'AOS6'
      nr = 1
      chmoth = 'SD0E'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 14
      chmoth = 'AOS6'
      x = 3.425
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 44
      natmed = 'AOR7 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR7'
      chshap = 'TUBS'
      nmed = 44
      npar = 5
      par(1) = 3.46
      par(2) = 3.54
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.58) stop 'consistency check #1 failed'

      chname = 'AOR7'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD0F'
      chmoth = 'AOR7'
      ndiv = 46
      iaxis = 2
      step = 1.277
      c0 = -29.387
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 45
      natmed = 'AOS7 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS7'
      chshap = 'TUBS'
      nmed = 45
      npar = 5
      par(1) = 3.46
      par(2) = 3.54
      par(3) = 0.85
      par(4) = -0.639
      par(5) = 0.638
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.60) stop 'consistency check #1 failed'

      chname = 'AOS7'
      nr = 1
      chmoth = 'SD0F'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 15
      chmoth = 'AOS7'
      x = 3.5
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 46
      natmed = 'AOR8 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR8'
      chshap = 'TUBS'
      nmed = 46
      npar = 5
      par(1) = 3.54
      par(2) = 3.62
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.61) stop 'consistency check #1 failed'

      chname = 'AOR8'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD10'
      chmoth = 'AOR8'
      ndiv = 47
      iaxis = 2
      step = 1.25
      c0 = -29.387
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 47
      natmed = 'AOS8 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS8'
      chshap = 'TUBS'
      nmed = 47
      npar = 5
      par(1) = 3.54
      par(2) = 3.62
      par(3) = 0.85
      par(4) = -0.625
      par(5) = 0.625
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.63) stop 'consistency check #1 failed'

      chname = 'AOS8'
      nr = 1
      chmoth = 'SD10'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 16
      chmoth = 'AOS8'
      x = 3.575
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 48
      natmed = 'AOR9 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOR9'
      chshap = 'TUBS'
      nmed = 48
      npar = 5
      par(1) = 3.62
      par(2) = 3.7
      par(3) = 0.85
      par(4) = -30
      par(5) = 30
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.64) stop 'consistency check #1 failed'

      chname = 'AOR9'
      nr = 1
      chmoth = 'ACWO'
      x = 0
      y = 0
      z = 0.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD11'
      chmoth = 'AOR9'
      ndiv = 48
      iaxis = 2
      step = 1.224
      c0 = -29.387
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 49
      natmed = 'AOS9 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'AOS9'
      chshap = 'TUBS'
      nmed = 49
      npar = 5
      par(1) = 3.62
      par(2) = 3.7
      par(3) = 0.85
      par(4) = -0.612
      par(5) = 0.612
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.66) stop 'consistency check #1 failed'

      chname = 'AOS9'
      nr = 1
      chmoth = 'SD11'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'PIN1'
      nr = 17
      chmoth = 'AOS9'
      x = 3.65
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 17
      theta1 = 90
      phi1 = 135
      theta2 = 90
      phi2 = -135
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWO'
      nr = 2
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 17
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 18
      theta1 = 90
      phi1 = -135
      theta2 = 90
      phi2 = -45
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWO'
      nr = 3
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 18
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 19
      theta1 = 90
      phi1 = -45
      theta2 = 90
      phi2 = 45
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'ACWO'
      nr = 4
      chmoth = 'AIRH'
      x = 0
      y = 0
      z = 0
      irot = 19
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 50
      natmed = 'HOUF Aluminum'
      nmat = 7
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'HOUF'
      chshap = 'TUBE'
      nmed = 50
      npar = 3
      par(1) = 0.25
      par(2) = 6.85
      par(3) = 0.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.67) stop 'consistency check #1 failed'

      chname = 'HOUF'
      nr = 1
      chmoth = 'INSU'
      x = 0
      y = 0
      z = -1.85
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 51
      natmed = 'DET1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'DET1'
      chshap = 'TUBE'
      nmed = 51
      npar = 3
      par(1) = 0
      par(2) = 50
      par(3) = 1
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.68) stop 'consistency check #1 failed'

      chname = 'DET1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -455.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 52
      natmed = 'COL1 Tungsten'
      nmat = 15
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'COL1'
      chshap = 'TUBE'
      nmed = 52
      npar = 3
      par(1) = 0.17
      par(2) = 10
      par(3) = 15
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.69) stop 'consistency check #1 failed'

      chname = 'COL1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -435
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 53
      natmed = 'MAG1 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'MAG1'
      chshap = 'BOX '
      nmed = 53
      npar = 3
      par(1) = 10
      par(2) = 25
      par(3) = 15
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.70) stop 'consistency check #1 failed'

      irot = 20
      theta1 = 90
      phi1 = 0
      theta2 = 3.50824e-15
      phi2 = 90
      theta3 = 90
      phi3 = -90
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'MAG1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -375
      irot = 20
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 54
      natmed = 'POL1 Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'POL1'
      chshap = 'BOX '
      nmed = 54
      npar = 3
      par(1) = 10
      par(2) = 25
      par(3) = 7
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.71) stop 'consistency check #1 failed'

      chname = 'POL1'
      nr = 1
      chmoth = 'MAG1'
      x = 0
      y = 0
      z = -8
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'POL1'
      nr = 2
      chmoth = 'MAG1'
      x = 0
      y = 0
      z = 8
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 17
      chnama = 'Lead'
      a = 207.2
      z = 82
      dens = 11.35
      radl = 0.561233
      absl = 17.0925
      nwbuf = 0
      call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf)

      itmed = 55
      natmed = 'GAP1 Lead'
      nmat = 17
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'GAP1'
      chshap = 'BOX '
      nmed = 55
      npar = 3
      par(1) = 4.5
      par(2) = 25
      par(3) = 1
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.72) stop 'consistency check #1 failed'

      chname = 'GAP1'
      nr = 1
      chmoth = 'MAG1'
      x = -5.5
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'GAP1'
      nr = 2
      chmoth = 'MAG1'
      x = 5.5
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 56
      natmed = 'ABS1 Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ABS1'
      chshap = 'BOX '
      nmed = 56
      npar = 3
      par(1) = 25
      par(2) = 25
      par(3) = 100
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.73) stop 'consistency check #1 failed'

      chname = 'ABS1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -240
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 57
      natmed = 'HOL1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'HOL1'
      chshap = 'BOX '
      nmed = 57
      npar = 3
      par(1) = 1
      par(2) = 1
      par(3) = 100
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.74) stop 'consistency check #1 failed'

      chname = 'HOL1'
      nr = 1
      chmoth = 'ABS1'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 58
      natmed = 'BLC1 Concrete'
      nmat = 6
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'BLC1'
      chshap = 'BOX '
      nmed = 58
      npar = 3
      par(1) = 100
      par(2) = 100
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.75) stop 'consistency check #1 failed'

      chname = 'BLC1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = -50
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 59
      natmed = 'BHO1 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'BHO1'
      chshap = 'BOX '
      nmed = 59
      npar = 3
      par(1) = 1
      par(2) = 1
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.76) stop 'consistency check #1 failed'

      chname = 'BHO1'
      nr = 1
      chmoth = 'BLC1'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 18
      namate = 'Nickel'
      a = 58.7
      z = 28
      dens = 8.876
      nlmat = 1
      wmat(1) = 1
      call gsmixt(imate,namate,a,z,dens,nlmat,wmat)

      itmed = 60
      natmed = 'COL2 Nickel'
      nmat = 18
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'COL2'
      chshap = 'TUBE'
      nmed = 60
      npar = 3
      par(1) = 0.5
      par(2) = 10
      par(3) = 25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.77) stop 'consistency check #1 failed'

      chname = 'COL2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 195
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 61
      natmed = 'MAG2 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'MAG2'
      chshap = 'BOX '
      nmed = 61
      npar = 3
      par(1) = 10
      par(2) = 25
      par(3) = 15
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.78) stop 'consistency check #1 failed'

      irot = 21
      theta1 = 90
      phi1 = 0
      theta2 = 3.50824e-15
      phi2 = 90
      theta3 = 90
      phi3 = -90
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'MAG2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 255
      irot = 21
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 62
      natmed = 'POL2 Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'POL2'
      chshap = 'BOX '
      nmed = 62
      npar = 3
      par(1) = 10
      par(2) = 25
      par(3) = 6.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.79) stop 'consistency check #1 failed'

      chname = 'POL2'
      nr = 1
      chmoth = 'MAG2'
      x = 0
      y = 0
      z = -8.75
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'POL2'
      nr = 2
      chmoth = 'MAG2'
      x = 0
      y = 0
      z = 8.75
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 63
      natmed = 'GAP2 Lead'
      nmat = 17
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'GAP2'
      chshap = 'BOX '
      nmed = 63
      npar = 3
      par(1) = 3.75
      par(2) = 25
      par(3) = 2.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.80) stop 'consistency check #1 failed'

      chname = 'GAP2'
      nr = 1
      chmoth = 'MAG2'
      x = -6.25
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'GAP2'
      nr = 2
      chmoth = 'MAG2'
      x = 6.25
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 64
      natmed = 'OMHO Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'OMHO'
      chshap = 'TUBE'
      nmed = 64
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.81) stop 'consistency check #1 failed'

      irot = 22
      theta1 = 90
      phi1 = 0
      theta2 = 3.50824e-15
      phi2 = 90
      theta3 = 90
      phi3 = -90
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'OMHO'
      nr = 1
      chmoth = 'MAG2'
      x = 0
      y = 0
      z = 0
      irot = 22
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 65
      natmed = 'IMHO Vacuum'
      nmat = 11
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'IMHO'
      chshap = 'TUBE'
      nmed = 65
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.82) stop 'consistency check #1 failed'

      chname = 'IMHO'
      nr = 1
      chmoth = 'OMHO'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 20
      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)

      imate = 19
      namate = 'Mylar'
      wmat(1) = 5
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 4
      call gfmate(5,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 2
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.39
      nlmat = -3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 66
      natmed = 'CAP1 Mylar'
      nmat = 19
      isvol = 0
      ifield = 2
      fieldm = 2
      tmaxfd = 1
      stemax = 1
      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 = 'CAP1'
      chshap = 'TUBE'
      nmed = 66
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 0.025
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.83) stop 'consistency check #1 failed'

      chname = 'CAP1'
      nr = 1
      chmoth = 'IMHO'
      x = 0
      y = 0
      z = 24.975
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 67
      natmed = 'OCH1 Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OCH1'
      chshap = 'TUBE'
      nmed = 67
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.84) stop 'consistency check #1 failed'

      chname = 'OCH1'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 285
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 68
      natmed = 'ICH1 Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ICH1'
      chshap = 'TUBE'
      nmed = 68
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.85) stop 'consistency check #1 failed'

      chname = 'ICH1'
      nr = 1
      chmoth = 'OCH1'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 69
      natmed = 'ABS2 Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ABS2'
      chshap = 'BOX '
      nmed = 69
      npar = 3
      par(1) = 50
      par(2) = 50
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.86) stop 'consistency check #1 failed'

      chname = 'ABS2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 340
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 70
      natmed = 'OAHO Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OAHO'
      chshap = 'TUBE'
      nmed = 70
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.87) stop 'consistency check #1 failed'

      chname = 'OAHO'
      nr = 1
      chmoth = 'ABS2'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 71
      natmed = 'IAHO Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'IAHO'
      chshap = 'TUBE'
      nmed = 71
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.88) stop 'consistency check #1 failed'

      chname = 'IAHO'
      nr = 1
      chmoth = 'OAHO'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 72
      natmed = 'OCH2 Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OCH2'
      chshap = 'TUBE'
      nmed = 72
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 20
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.89) stop 'consistency check #1 failed'

      chname = 'OCH2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 410
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 73
      natmed = 'ICH2 Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ICH2'
      chshap = 'TUBE'
      nmed = 73
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 20
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.90) stop 'consistency check #1 failed'

      chname = 'ICH2'
      nr = 1
      chmoth = 'OCH2'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 74
      natmed = 'BLC2 Concrete'
      nmat = 6
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'BLC2'
      chshap = 'BOX '
      nmed = 74
      npar = 3
      par(1) = 200
      par(2) = 125
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.91) stop 'consistency check #1 failed'

      chname = 'BLC2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -10
      z = 480
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 75
      natmed = 'OBHO Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OBHO'
      chshap = 'TUBE'
      nmed = 75
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.92) stop 'consistency check #1 failed'

      chname = 'OBHO'
      nr = 1
      chmoth = 'BLC2'
      x = 0
      y = -25
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 76
      natmed = 'IBHO Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'IBHO'
      chshap = 'TUBE'
      nmed = 76
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 50
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.93) stop 'consistency check #1 failed'

      chname = 'IBHO'
      nr = 1
      chmoth = 'OBHO'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 77
      natmed = 'OCH3 Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OCH3'
      chshap = 'TUBE'
      nmed = 77
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 55
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.94) stop 'consistency check #1 failed'

      chname = 'OCH3'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = -35
      z = 585
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 78
      natmed = 'ICH3 Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'ICH3'
      chshap = 'TUBE'
      nmed = 78
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 55
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.95) stop 'consistency check #1 failed'

      chname = 'ICH3'
      nr = 1
      chmoth = 'OCH3'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 79
      natmed = 'WALL Lead'
      nmat = 17
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'WALL'
      chshap = 'BOX '
      nmed = 79
      npar = 3
      par(1) = 235
      par(2) = 135
      par(3) = 7.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.96) stop 'consistency check #1 failed'

      chname = 'WALL'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = 0
      z = 592.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 80
      natmed = 'OWHO Iron'
      nmat = 10
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'OWHO'
      chshap = 'TUBE'
      nmed = 80
      npar = 3
      par(1) = 0
      par(2) = 2.5
      par(3) = 7.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.97) stop 'consistency check #1 failed'

      chname = 'OWHO'
      nr = 1
      chmoth = 'WALL'
      x = 0
      y = -35
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 81
      natmed = 'IWHO Vacuum'
      nmat = 11
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'IWHO'
      chshap = 'TUBE'
      nmed = 81
      npar = 3
      par(1) = 0
      par(2) = 2.2
      par(3) = 7.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.98) stop 'consistency check #1 failed'

      chname = 'IWHO'
      nr = 1
      chmoth = 'OWHO'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 82
      natmed = 'DET2 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'DET2'
      chshap = 'BOX '
      nmed = 82
      npar = 3
      par(1) = 225
      par(2) = 135
      par(3) = 1
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.99) stop 'consistency check #1 failed'

      chname = 'DET2'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = 0
      z = 75
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 83
      natmed = 'DET3 Air'
      nmat = 1
      isvol = 0
      ifield = 0
      fieldm = 0
      tmaxfd = 0
      stemax = 1
      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 = 'DET3'
      chshap = 'BOX '
      nmed = 83
      npar = 3
      par(1) = 225
      par(2) = 1
      par(3) = 600
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.100) stop 'consistency check #1 failed'

      chname = 'DET3'
      nr = 1
      chmoth = 'CAVE'
      x = 0
      y = 132
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 84
      natmed = 'HALL Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'HALL'
      chshap = 'BOX '
      nmed = 84
      npar = 3
      par(1) = 850
      par(2) = 750
      par(3) = 1500
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.101) stop 'consistency check #1 failed'

      chname = 'HALL'
      nr = 1
      chmoth = 'SITE'
      x = -150
      y = 350
      z = 500
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 85
      natmed = 'TARG Vacuum'
      nmat = 11
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'TARG'
      chshap = 'TUBE'
      nmed = 85
      npar = 3
      par(1) = 0
      par(2) = 4
      par(3) = 27.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.102) stop 'consistency check #1 failed'

      chname = 'TARG'
      nr = 1
      chmoth = 'HALL'
      x = 150
      y = -350
      z = -437.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 21
      namate = 'HighDensityROHACELL'
      wmat(1) = 11
      call gfmate(5,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 8
      call gfmate(20,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 1
      call gfmate(2,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      wmat(4) = 2
      call gfmate(3,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf)
      dens = 0.11
      nlmat = -4
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 86
      natmed = 'CYLW HighDensityROHACELL'
      nmat = 21
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'CYLW'
      chshap = 'TUBE'
      nmed = 86
      npar = 3
      par(1) = 3
      par(2) = 4
      par(3) = 27.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.103) stop 'consistency check #1 failed'

      chname = 'CYLW'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 87
      natmed = 'UWIT Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'UWIT'
      chshap = 'TUBE'
      nmed = 87
      npar = 3
      par(1) = 0
      par(2) = 2
      par(3) = 0.0035
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.104) stop 'consistency check #1 failed'

      chname = 'UWIT'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = -27.495
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 88
      natmed = 'DWIT Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DWIT'
      chshap = 'TUBE'
      nmed = 88
      npar = 3
      par(1) = 0
      par(2) = 2
      par(3) = 0.0035
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.105) stop 'consistency check #1 failed'

      chname = 'DWIT'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = 27.495
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 89
      natmed = 'UFLT Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'UFLT'
      chshap = 'TUBE'
      nmed = 89
      npar = 3
      par(1) = 2
      par(2) = 3
      par(3) = 0.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.106) stop 'consistency check #1 failed'

      chname = 'UFLT'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = -27.25
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 90
      natmed = 'DFLT HighDensityROHACELL'
      nmat = 21
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DFLT'
      chshap = 'TUBE'
      nmed = 90
      npar = 3
      par(1) = 2
      par(2) = 3
      par(3) = 0.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.107) stop 'consistency check #1 failed'

      chname = 'DFLT'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = 27
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 22
      namate = 'Kapton'
      wmat(1) = 22
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 10
      call gfmate(5,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 5
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      wmat(4) = 2
      call gfmate(2,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf)
      dens = 1.42
      nlmat = -4
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 91
      natmed = 'TGTV Kapton'
      nmat = 22
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'TGTV'
      chshap = 'TUBE'
      nmed = 91
      npar = 3
      par(1) = 0
      par(2) = 1.515
      par(3) = 15.015
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.108) stop 'consistency check #1 failed'

      chname = 'TGTV'
      nr = 1
      chmoth = 'TARG'
      x = 0
      y = 0
      z = 2.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 23
      chnama = 'LiqHydrogen'
      a = 1.00797
      z = 1
      dens = 0.0708
      radl = 865.537
      absl = 717.514
      nwbuf = 0
      call gsmate(imate,chnama,a,z,dens,radl,absl,ubuf,nwbuf)

      itmed = 92
      natmed = 'LIH2 LiqHydrogen'
      nmat = 23
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'LIH2'
      chshap = 'TUBE'
      nmed = 92
      npar = 3
      par(1) = 0
      par(2) = 1.5
      par(3) = 15
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.109) stop 'consistency check #1 failed'

      chname = 'LIH2'
      nr = 1
      chmoth = 'TGTV'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 93
      natmed = 'LASS Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'LASS'
      chshap = 'PCON'
      nmed = 93
      npar = 15
      par(1) = 0
      par(2) = 360
      par(3) = 4
      par(4) = -60.8
      par(5) = 4.5
      par(6) = 188
      par(7) = 92
      par(8) = 4.5
      par(9) = 188
      par(10) = 92
      par(11) = 0
      par(12) = 188
      par(13) = 460
      par(14) = 0
      par(15) = 188
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.110) stop 'consistency check #1 failed'

      chname = 'LASS'
      nr = 1
      chmoth = 'HALL'
      x = 150
      y = -350
      z = -500
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 94
      natmed = 'IYUP Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'IYUP'
      chshap = 'TUBE'
      nmed = 94
      npar = 3
      par(1) = 92.71
      par(2) = 187.96
      par(3) = 25.4
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.111) stop 'consistency check #1 failed'

      chname = 'IYUP'
      nr = 1
      chmoth = 'LASS'
      x = 0
      y = 0
      z = -25.4
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 95
      natmed = 'IYOK Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'IYOK'
      chshap = 'TUBE'
      nmed = 95
      npar = 3
      par(1) = 95
      par(2) = 187.96
      par(3) = 188.9
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.112) stop 'consistency check #1 failed'

      chname = 'IYOK'
      nr = 1
      chmoth = 'LASS'
      x = 0
      y = 0
      z = 188.9
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 96
      natmed = 'IYDN Iron'
      nmat = 10
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'IYDN'
      chshap = 'TUBE'
      nmed = 96
      npar = 3
      par(1) = 92.71
      par(2) = 187.96
      par(3) = 33
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.113) stop 'consistency check #1 failed'

      chname = 'IYDN'
      nr = 1
      chmoth = 'LASS'
      x = 0
      y = 0
      z = 410.8
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 24
      namate = 'LowDensityROHACELL'
      wmat(1) = 11
      call gfmate(5,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 8
      call gfmate(20,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 1
      call gfmate(2,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      wmat(4) = 2
      call gfmate(3,chnama,amat(4),zmat(4),dens,radl,absl,ubuf,nwbuf)
      dens = 0.032
      nlmat = -4
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 97
      natmed = 'STRT LowDensityROHACELL'
      nmat = 24
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STRT'
      chshap = 'PGON'
      nmed = 97
      npar = 25
      par(1) = 0
      par(2) = 360
      par(3) = 24
      par(4) = 7
      par(5) = -140
      par(6) = 6.951
      par(7) = 8.675
      par(8) = 51.383
      par(9) = 6.951
      par(10) = 8.675
      par(11) = 51.726
      par(12) = 6.809
      par(13) = 8.675
      par(14) = 53.938
      par(15) = 5.893
      par(16) = 7.759
      par(17) = 54.914
      par(18) = 4.916
      par(19) = 7.354
      par(20) = 57.915
      par(21) = 1.915
      par(22) = 4.353
      par(23) = 59.134
      par(24) = 3.134
      par(25) = 3.134
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.114) stop 'consistency check #1 failed'

      chname = 'STRT'
      nr = 1
      chmoth = 'LASS'
      x = 0
      y = 0
      z = 37.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 25
      namate = 'Scintillator'
      wmat(1) = 8
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 9
      call gfmate(5,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 = 98
      natmed = 'STRC Scintillator'
      nmat = 25
      isvol = 1
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STRC'
      chshap = 'PGON'
      nmed = 98
      npar = 25
      par(1) = -7.5
      par(2) = 15
      par(3) = 1
      par(4) = 7
      par(5) = 0
      par(6) = 7.481
      par(7) = 7.696
      par(8) = 51.488
      par(9) = 7.481
      par(10) = 7.696
      par(11) = 51.531
      par(12) = 7.463
      par(13) = 7.696
      par(14) = 54.238
      par(15) = 6.342
      par(16) = 6.575
      par(17) = 54.36
      par(18) = 6.22
      par(19) = 6.524
      par(20) = 58.29
      par(21) = 2.29
      par(22) = 2.594
      par(23) = 58.442
      par(24) = 2.442
      par(25) = 2.442
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.115) stop 'consistency check #1 failed'

      irot = 23
      theta1 = 90
      phi1 = 7.5
      theta2 = 90
      phi2 = 97.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 1
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 23
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 24
      theta1 = 90
      phi1 = 22.5
      theta2 = 90
      phi2 = 112.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 2
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 24
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 25
      theta1 = 90
      phi1 = 37.5
      theta2 = 90
      phi2 = 127.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 3
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 25
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 26
      theta1 = 90
      phi1 = 52.5
      theta2 = 90
      phi2 = 142.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 4
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 26
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 27
      theta1 = 90
      phi1 = 67.5
      theta2 = 90
      phi2 = 157.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 5
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 27
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 28
      theta1 = 90
      phi1 = 82.5
      theta2 = 90
      phi2 = 172.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 6
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 28
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 29
      theta1 = 90
      phi1 = 97.5
      theta2 = 90
      phi2 = -172.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 7
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 29
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 30
      theta1 = 90
      phi1 = 112.5
      theta2 = 90
      phi2 = -157.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 8
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 30
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 31
      theta1 = 90
      phi1 = 127.5
      theta2 = 90
      phi2 = -142.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 9
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 31
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 32
      theta1 = 90
      phi1 = 142.5
      theta2 = 90
      phi2 = -127.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 10
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 32
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 33
      theta1 = 90
      phi1 = 157.5
      theta2 = 90
      phi2 = -112.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 11
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 33
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 34
      theta1 = 90
      phi1 = 172.5
      theta2 = 90
      phi2 = -97.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 12
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 34
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 35
      theta1 = 90
      phi1 = -172.5
      theta2 = 90
      phi2 = -82.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 13
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 35
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 36
      theta1 = 90
      phi1 = -157.5
      theta2 = 90
      phi2 = -67.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 14
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 36
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 37
      theta1 = 90
      phi1 = -142.5
      theta2 = 90
      phi2 = -52.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 15
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 37
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 38
      theta1 = 90
      phi1 = -127.5
      theta2 = 90
      phi2 = -37.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 16
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 38
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 39
      theta1 = 90
      phi1 = -112.5
      theta2 = 90
      phi2 = -22.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 17
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 39
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 40
      theta1 = 90
      phi1 = -97.5
      theta2 = 90
      phi2 = -7.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 18
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 40
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 41
      theta1 = 90
      phi1 = -82.5
      theta2 = 90
      phi2 = 7.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 19
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 41
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 42
      theta1 = 90
      phi1 = -67.5
      theta2 = 90
      phi2 = 22.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 20
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 42
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 43
      theta1 = 90
      phi1 = -52.5
      theta2 = 90
      phi2 = 37.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 21
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 43
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 44
      theta1 = 90
      phi1 = -37.5
      theta2 = 90
      phi2 = 52.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 22
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 44
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 45
      theta1 = 90
      phi1 = -22.5
      theta2 = 90
      phi2 = 67.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 23
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 45
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      irot = 46
      theta1 = 90
      phi1 = -7.5
      theta2 = 90
      phi2 = 82.5
      theta3 = 0
      phi3 = 0
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STRC'
      nr = 24
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 46
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 99
      natmed = 'STSU HighDensityROHACELL'
      nmat = 21
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STSU'
      chshap = 'PGON'
      nmed = 99
      npar = 25
      par(1) = 0
      par(2) = 360
      par(3) = 24
      par(4) = 7
      par(5) = 0
      par(6) = 6.951
      par(7) = 7.481
      par(8) = 51.383
      par(9) = 6.951
      par(10) = 7.481
      par(11) = 51.488
      par(12) = 6.907
      par(13) = 7.481
      par(14) = 53.938
      par(15) = 5.893
      par(16) = 6.467
      par(17) = 54.238
      par(18) = 5.593
      par(19) = 6.342
      par(20) = 57.915
      par(21) = 1.915
      par(22) = 2.665
      par(23) = 58.29
      par(24) = 2.29
      par(25) = 2.29
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.116) stop 'consistency check #1 failed'

      chname = 'STSU'
      nr = 1
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 26
      namate = 'Cellulose'
      wmat(1) = 10
      call gfmate(5,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 6
      call gfmate(20,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 5
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.1
      nlmat = -3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 100
      natmed = 'STWR Cellulose'
      nmat = 26
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STWR'
      chshap = 'PGON'
      nmed = 100
      npar = 25
      par(1) = 0
      par(2) = 360
      par(3) = 24
      par(4) = 7
      par(5) = 0
      par(6) = 8.675
      par(7) = 8.725
      par(8) = 51.726
      par(9) = 8.675
      par(10) = 8.725
      par(11) = 51.736
      par(12) = 8.671
      par(13) = 8.725
      par(14) = 54.914
      par(15) = 7.354
      par(16) = 7.409
      par(17) = 54.942
      par(18) = 7.326
      par(19) = 7.397
      par(20) = 59.134
      par(21) = 3.134
      par(22) = 3.205
      par(23) = 59.17
      par(24) = 3.17
      par(25) = 3.17
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.117) stop 'consistency check #1 failed'

      chname = 'STWR'
      nr = 1
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'STWR'
      nr = 2
      chmoth = 'STRT'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 101
      natmed = 'TARM Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'TARM'
      chshap = 'TUBE'
      nmed = 101
      npar = 3
      par(1) = 6.951
      par(2) = 8.675
      par(3) = 70
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.118) stop 'consistency check #1 failed'

      chname = 'TARM'
      nr = 1
      chmoth = 'STRT'
      x = 0
      y = 0
      z = -70
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 27
      namate = 'Plexiglas'
      wmat(1) = 5
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 8
      call gfmate(5,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 2
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.2
      nlmat = -3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 102
      natmed = 'STLG Plexiglas'
      nmat = 27
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STLG'
      chshap = 'TUBE'
      nmed = 102
      npar = 3
      par(1) = 7.551
      par(2) = 8.675
      par(3) = 70
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.119) stop 'consistency check #1 failed'

      chname = 'STLG'
      nr = 1
      chmoth = 'TARM'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 103
      natmed = 'CDC Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'CDC'
      chshap = 'TUBE'
      nmed = 103
      npar = 3
      par(1) = 9
      par(2) = 59
      par(3) = 87
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.120) stop 'consistency check #1 failed'

      chname = 'CDC'
      nr = 1
      chmoth = 'LASS'
      x = 0
      y = 0
      z = 92
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 104
      natmed = 'DC1 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DC1'
      chshap = 'TUBE'
      nmed = 104
      npar = 3
      par(1) = 10.12
      par(2) = 11.8
      par(3) = 76.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.121) stop 'consistency check #1 failed'

      chname = 'DC1'
      nr = 1
      chmoth = 'CDC'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD12'
      chmoth = 'DC1'
      ndiv = 43
      iaxis = 2
      step = 8.37209
      c0 = -4.186
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 105
      natmed = 'DS1 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DS1'
      chshap = 'TUBS'
      nmed = 105
      npar = 5
      par(1) = 10.12
      par(2) = 11.8
      par(3) = 76.25
      par(4) = -4.186
      par(5) = 4.186
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.123) stop 'consistency check #1 failed'

      chname = 'DS1'
      nr = 1
      chmoth = 'SD12'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 29
      namate = 'CO2'
      wmat(1) = 1
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 2
      call gfmate(3,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      dens = 0.001977
      nlmat = -2
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      imate = 28
      namate = 'CDchamberGas'
      wmat(1) = 0.85
      call gfmate(4,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 0.15
      call gfmate(29,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      dens = 0.00170091
      nlmat = 2
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 106
      natmed = 'STRA CDchamberGas'
      nmat = 28
      isvol = 1
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STRA'
      chshap = 'TUBE'
      nmed = 106
      npar = 3
      par(1) = 0
      par(2) = 0.8
      par(3) = 76
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.124) stop 'consistency check #1 failed'

      chname = 'STRA'
      nr = 1
      chmoth = 'DS1'
      x = 10.96
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 30
      namate = 'AluminizedKapton'
      wmat(1) = 0.087
      call gfmate(7,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 0.913
      call gfmate(22,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      dens = 1.48109
      nlmat = 2
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 107
      natmed = 'STRW AluminizedKapton'
      nmat = 30
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STRW'
      chshap = 'TUBE'
      nmed = 107
      npar = 3
      par(1) = 0.7895
      par(2) = 0.8
      par(3) = 75
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.125) stop 'consistency check #1 failed'

      chname = 'STRW'
      nr = 1
      chmoth = 'STRA'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 108
      natmed = 'WIRE Tungsten'
      nmat = 15
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'WIRE'
      chshap = 'TUBE'
      nmed = 108
      npar = 3
      par(1) = 0
      par(2) = 0.0015
      par(3) = 76
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.126) stop 'consistency check #1 failed'

      chname = 'WIRE'
      nr = 1
      chmoth = 'STRA'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 109
      natmed = 'UPLG Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'UPLG'
      chshap = 'TUBE'
      nmed = 109
      npar = 3
      par(1) = 0
      par(2) = 0.8
      par(3) = 0.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.127) stop 'consistency check #1 failed'

      chname = 'UPLG'
      nr = 1
      chmoth = 'STRA'
      x = 0
      y = 0
      z = -75.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 31
      namate = 'Delrin'
      wmat(1) = 1
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 2
      call gfmate(5,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 1
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.42
      nlmat = -3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 110
      natmed = 'DPLG Delrin'
      nmat = 31
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DPLG'
      chshap = 'TUBE'
      nmed = 110
      npar = 3
      par(1) = 0
      par(2) = 0.8
      par(3) = 0.5
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.128) stop 'consistency check #1 failed'

      chname = 'DPLG'
      nr = 1
      chmoth = 'STRA'
      x = 0
      y = 0
      z = 75.5
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 111
      natmed = 'CDPU Aluminum'
      nmat = 7
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'CDPU'
      chshap = 'TUBE'
      nmed = 111
      npar = 3
      par(1) = 9
      par(2) = 59
      par(3) = 0.45
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.129) stop 'consistency check #1 failed'

      chname = 'CDPU'
      nr = 1
      chmoth = 'DS1'
      x = 0
      y = 0
      z = -75.45
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      imate = 33
      namate = 'CarbonFiber'
      wmat(1) = 1
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      dens = 2.265
      nlmat = -1
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      imate = 34
      namate = 'EpoxyResin'
      wmat(1) = 18
      call gfmate(20,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 20
      call gfmate(5,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 3
      call gfmate(3,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.3
      nlmat = -3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      imate = 32
      namate = 'CarbonFiberEpoxy'
      wmat(1) = 0.375
      call gfmate(33,chnama,amat(1),zmat(1),dens,radl,absl,ubuf,nwbuf)
      wmat(2) = 0.25
      call gfmate(34,chnama,amat(2),zmat(2),dens,radl,absl,ubuf,nwbuf)
      wmat(3) = 0.375
      call gfmate(33,chnama,amat(3),zmat(3),dens,radl,absl,ubuf,nwbuf)
      dens = 1.91046
      nlmat = 3
      call gsmixt(imate,namate,amat,zmat,dens,nlmat,wmat)

      itmed = 112
      natmed = 'CDPD CarbonFiberEpoxy'
      nmat = 32
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'CDPD'
      chshap = 'TUBE'
      nmed = 112
      npar = 3
      par(1) = 9
      par(2) = 59
      par(3) = 0.2
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.130) stop 'consistency check #1 failed'

      chname = 'CDPD'
      nr = 1
      chmoth = 'DS1'
      x = 0
      y = 0
      z = 75.2
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 113
      natmed = 'DC2 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DC2'
      chshap = 'TUBE'
      nmed = 113
      npar = 3
      par(1) = 11.901
      par(2) = 13.581
      par(3) = 76.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.131) stop 'consistency check #1 failed'

      chname = 'DC2'
      nr = 1
      chmoth = 'CDC'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD13'
      chmoth = 'DC2'
      ndiv = 50
      iaxis = 2
      step = 7.2
      c0 = -3.6
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 114
      natmed = 'DS2 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DS2'
      chshap = 'TUBS'
      nmed = 114
      npar = 5
      par(1) = 11.9
      par(2) = 13.58
      par(3) = 76.25
      par(4) = -3.6
      par(5) = 3.6
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.133) stop 'consistency check #1 failed'

      chname = 'DS2'
      nr = 1
      chmoth = 'SD13'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'STRA'
      nr = 2
      chmoth = 'DS2'
      x = 12.741
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'CDPU'
      nr = 2
      chmoth = 'DS2'
      x = 0
      y = 0
      z = -75.45
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'CDPD'
      nr = 2
      chmoth = 'DS2'
      x = 0
      y = 0
      z = 75.2
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 115
      natmed = 'DC3 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DC3'
      chshap = 'TUBE'
      nmed = 115
      npar = 3
      par(1) = 13.682
      par(2) = 15.362
      par(3) = 76.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.134) stop 'consistency check #1 failed'

      chname = 'DC3'
      nr = 1
      chmoth = 'CDC'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'SD14'
      chmoth = 'DC3'
      ndiv = 57
      iaxis = 2
      step = 6.31579
      c0 = -3.158
      numed = 0
      ndvmax = 0
      call gsdvx(chname,chmoth,ndiv,iaxis,step,c0,numed,ndvmax)

      itmed = 116
      natmed = 'DS3 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DS3'
      chshap = 'TUBS'
      nmed = 116
      npar = 5
      par(1) = 13.68
      par(2) = 15.36
      par(3) = 76.25
      par(4) = -3.158
      par(5) = 3.158
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.136) stop 'consistency check #1 failed'

      chname = 'DS3'
      nr = 1
      chmoth = 'SD14'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'STRA'
      nr = 3
      chmoth = 'DS3'
      x = 14.522
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'CDPU'
      nr = 3
      chmoth = 'DS3'
      x = 0
      y = 0
      z = -75.45
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      chname = 'CDPD'
      nr = 3
      chmoth = 'DS3'
      x = 0
      y = 0
      z = 75.2
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 117
      natmed = 'DC4 Air'
      nmat = 1
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'DC4'
      chshap = 'TUBE'
      nmed = 117
      npar = 3
      par(1) = 15.464
      par(2) = 23.88
      par(3) = 76.25
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.137) stop 'consistency check #1 failed'

      chname = 'DC4'
      nr = 1
      chmoth = 'CDC'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 118
      natmed = 'STLA CDchamberGas'
      nmat = 28
      isvol = 1
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STLA'
      chshap = 'TUBE'
      nmed = 118
      npar = 3
      par(1) = 0
      par(2) = 0.8
      par(3) = 76.42
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.138) stop 'consistency check #1 failed'

      irot = 47
      theta1 = 90
      phi1 = 0
      theta2 = 96
      phi2 = 90
      theta3 = 6
      phi3 = 90
      call gsrotm(irot,theta1,phi1,theta2,phi2,theta3,phi3)

      chname = 'STLA'
      nr = 1
      chmoth = 'DC4'
      x = 16.304
      y = 0
      z = 0
      irot = 47
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 119
      natmed = 'STLW AluminizedKapton'
      nmat = 30
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'STLW'
      chshap = 'TUBE'
      nmed = 119
      npar = 3
      par(1) = 0.7895
      par(2) = 0.8
      par(3) = 75.415
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.139) stop 'consistency check #1 failed'

      chname = 'STLW'
      nr = 1
      chmoth = 'STLA'
      x = 0
      y = 0
      z = 0
      irot = 0
      chonly = 'ONLY'
      call gspos(chname,nr,chmoth,x,y,z,irot,chonly)

      itmed = 120
      natmed = 'WILE Tungsten'
      nmat = 15
      isvol = 0
      ifield = 2
      fieldm = 22.4
      tmaxfd = 1
      stemax = 1
      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 = 'WILE'
      chshap = 'TUBE'
      nmed = 120
      npar = 3
      par(1) = 0
      par(2) = 0.0015
      par(3) = 76.415
      call gsvolu(chname,chshap,nmed,par,npar,ivolu)
      if (ivolu.ne.140) stop 'consistency check #1 failed'

      chname = 'WILE'
      nr = 1
      chmoth = 'STLA'
      x = 0
      y = 0
      z = 0