! Diagonal elements of the dC part of A and the respective friction coefficients
ind=1
ind1=0
- print *,"TUTUTUT?!",nnt,nct-1
+! print *,"TUTUTUT?!",nnt,nct-1
do i=nnt,nct-1
mnum=molnum(i)
ind=ind+1
print *,i,coeff,mp(mnum)
massvec(ind1)=mp(mnum)
Gmat(ind,ind)=coeff
- print *,"i",mp(mnum)
+! print *,"i",mp(mnum)
A(ind1,ind)=0.5d0
enddo