(alphasur(k,i,j),k=1,4),sigmap1(i,j),sigmap2(i,j), & !6 w tej linii
chis(i,j),chis(j,i) !2 w tej linii
endif
-! print *,eps(i,j),sigma(i,j),"SIGMAP",i,j,sigmap1(i,j),sigmap2(j,i)
+ print *,eps(i,j),sigma(i,j),"SIGMAP",i,j,sigmap1(i,j),sigmap2(j,i), wqdip(1,i,j)
END DO
END DO
do i=1,ntyp
enddo
ip(5)=0.2
- do i=1,ntyp
- do j=1,ntyp_molec(5)
+ do j=1,ntyp_molec(5)
+ do i=1,ntyp
+! do j=1,ntyp_molec(5)
! write (*,*) "Im in ALAB", i, " ", j
read(iion,*) &
epscat(i,j),sigmacat(i,j), &
do j=1,ntyp_molec(5)
epsij=epscat(i,j)
rrij=sigmacat(i,j)
+ rrij=rrij**expon
sigeps=dsign(1.0D0,epsij)
epsij=dabs(epsij)
aa_aq_cat(i,j)=epsij*rrij*rrij