program iwr c implicit real*8 (a-h,o-z) implicit integer*4 (i-n) dimension bmm(100) dimension bmmn(1000) c write(6,6001) 6001 format(' This program rewrites integrals from BBCD.SCR',/, 1 'into RBCD.SCR using a cut-off',/,/, 1 'Input the limit:') read(5,*)OL write(6,*)'NMO:' read(5,*)nmo c ks=1 je=nmo p=0.0d0 pt=0.0d0 open(53,file='BBCD.SCR',form='unformatted',status='old') do 57 IDB=ks,je c c read from 100 word blocks: read(53)ntest i1=1 60001 i2=min(i1+99,ntest) read(53)(bmm(ibb-i1+1),ibb=i1,i2) c c look for big integrals do 1 i8=i1,i2 if(abs(bmm(i8-i1+1)).gt.OL)then p=p+1.0d0 endif 1 pt=pt+1.0d0 c if(i2.lt.ntest)then i1=i2+1 goto 60001 endif c 57 continue close(53) write(6,6002)p,pt,p/pt*100.0d0 6002 format(f20.0,' big',/,f20.0,' total',/,20(1h-),/,f20.1,' %') stop end