iti1=itortyp(itype(i+1))
iti2=itortyp(itype(i+2))
iti3=itortyp(itype(i+3))
- write(iout,*) "iti1",iti1," iti2",iti2," iti3",iti3
+C write(iout,*) i,"iti1",iti1," iti2",iti2," iti3",iti3,itype(i+3)
call transpose2(EUg(1,1,i+1),e1t(1,1))
call transpose2(Eug(1,1,i+2),e2t(1,1))
call transpose2(Eug(1,1,i+3),e3t(1,1))