c CODE TO INTEGRATE G.SPC OVER A GIVEN ENERGYRANGE

      open(unit=10,file='g.spc',status='unknown')
      write(*,*)'input elow,ehigh in MeV'
      read(*,*)elow,ehigh

      sum=0.
      do i=1,2300
      read(10,*)x,de,y
      if(x.ge.elow.and.x.le.ehigh)sum=sum+y*de

      enddo

100   format(e12.3)
      write(*,*)'lower and upper bounds on energy range'
      write(*,101)elow,ehigh
101   format(2f10.3)
      write(*,*)'integral='
      write(*,100)sum

      stop
	end