(wqdipcat(k,i,j),k=1,2), &
alphapolcat(i,j),alphapolcat(j,i), &
(alphisocat(k,i,j),k=1,4),sigiso1cat(i,j),sigiso2cat(i,j),epsintabcat(i,j),debaykapcat(i,j)
+
+ if (chi1cat(i,j).gt.0.9) write (*,*) "WTF ANISO", i,j, chi1cat(i,j)
! print *,eps(i,j),sigma(i,j),"SIGMAP",i,j,sigmap1(i,j),sigmap2(j,i)
! if (i.eq.1) then
! write (iout,*) 'i= ', i, ' j= ', j
print_compon = index(controlcard,"PRINT_COMPON").gt.0
rattle = index(controlcard,"RATTLE").gt.0
preminim=(index(controlcard,'PREMINIM').gt.0)
+ forceminim=(index(controlcard,'FORCEMINIM').gt.0)
write (iout,*) "PREMINIM ",preminim
dccart=(index(controlcard,'CART').gt.0)
if (preminim) call read_minim