program ro c rewrite outputs for Lamb crated by 'grep 'kcal' *.o parameter (nr0=11) dimension hf(nr0,nr0),b3(nr0,nr0) character*80 s80 do 1 i=1,nr0 do 1 j=1,nr0 b3(i,j)=0.0 1 hf(i,j)=0.0 open(44,file='o') 2 read(44,4400,end=999,err=999)s80 4400 format(a80) ir=0 ib3=0 ln=len(s80) i1s=0 i2s=0 i2p=0 i2sr=0 i2pr=0 il=0 ic=0 inu=0 do i=1,ln if(s80(i:i).eq.':')ic=i if(s80(i:i).eq.'=')ir=1 if(s80(i:min(i+3,ln)).eq.'E nu')inu=1 if(s80(i:min(i+1,ln)).eq.'ss')il=1 if(s80(i:min(i+1,ln)).eq.'pp')il=2 if(s80(i:min(i+2,ln)).eq.'vdz')il=3 if(s80(i:min(i+2,ln)).eq.'vtz')il=4 if(s80(i:min(i+2,ln)).eq.'vqz')il=5 if(s80(i:min(i+2,ln)).eq.'v5z')il=6 if(s80(i:min(i+2,ln)).eq.'v6z')il=7 if(s80(i:min(i+2,ln)).eq.'b3l')ib3=1 if(s80(i:min(i+3,ln)).eq.'1s.o')i1s=1 if(s80(i:min(i+3,ln)).eq.'2s.o')i2s=1 if(s80(i:min(i+3,ln)).eq.'2p.o')i2p=1 if(s80(i:min(i+4,ln)).eq.'2pr.o')i2pr=1 if(s80(i:min(i+4,ln)).eq.'2sr.o')i2sr=1 enddo c write(6,4400)s80 c write(6,640)ic,inu,il,ib3,i1s,i2s,i2p,i2pr,i2sr 640 format(10i2) c read(5,*)ioo if(ib3.eq.1)then if(inu.eq.0.and.ir.eq.0)then if(i1s.eq.1)read(s80(ic+1:ic+12),*)b3(il,1) if(i2sr.eq.1)read(s80(ic+1:ic+12),*)b3(il,2) if(i2sr.eq.1)read(s80(ic+38:ic+49),*)b3(il,3) if(i2pr.eq.1)read(s80(ic+1:ic+12),*)b3(il,4) if(i2pr.eq.1)read(s80(ic+38:ic+49),*)b3(il,5) if(i2s.eq.1)read(s80(ic+1:ic+12),*)b3(il,7) if(i2s.eq.1)read(s80(ic+38:ic+49),*)b3(il,8) if(i2p.eq.1)read(s80(ic+1:ic+12),*)b3(il,9) if(i2p.eq.1)read(s80(ic+38:ic+49),*)b3(il,10) else if(inu.eq.1)then if(i2s.eq.1)read(s80(ic+8:ic+23),*)b3(il,11) if(i2sr.eq.1)read(s80(ic+8:ic+23),*)b3(il,6) endif endif else if(inu.eq.0.and.ir.eq.0)then if(i1s.eq.1)read(s80(ic+1:ic+12),*)hf(il,1) if(i2s.eq.1)read(s80(ic+1:ic+12),*)hf(il,7) if(i2s.eq.1)read(s80(ic+38:ic+49),*)hf(il,8) if(i2sr.eq.1)read(s80(ic+1:ic+12),*)hf(il,2) if(i2sr.eq.1)read(s80(ic+38:ic+49),*)hf(il,3) if(i2p.eq.1)read(s80(ic+1:ic+12),*)hf(il,9) if(i2p.eq.1)read(s80(ic+38:ic+49),*)hf(il,10) if(i2pr.eq.1)read(s80(ic+1:ic+12),*)hf(il,4) if(i2pr.eq.1)read(s80(ic+38:ic+49),*)hf(il,5) else if(inu.eq.1)then if(i2s.eq.1)read(s80(ic+8:ic+23),*)hf(il,11) if(i2sr.eq.1)read(s80(ic+8:ic+23),*)hf(il,6) endif endif endif goto 2 999 close(44) do i=1,7 write(6,600)(b3(i,j),j=1,11) enddo do i=1,7 write(6,600)(hf(i,j),j=1,11) enddo 600 format(11f13.7) stop end