nntheterm=max0(ntheterm,ntheterm2,ntheterm3)
write(iout,*) "I am here",ntyp1
read (ithep,*,err=111,end=111) (ithetyp(i),i=1,ntyp1)
nntheterm=max0(ntheterm,ntheterm2,ntheterm3)
write(iout,*) "I am here",ntyp1
read (ithep,*,err=111,end=111) (ithetyp(i),i=1,ntyp1)