rminroz=100d2
c print *,"nz_start",nz_start," nz_end",nz_end
c if (symetr.le.1) then
+C print *,nperm,"nperm"
do kkk=1,nperm
c do i=nz_start,nz_end
c iatom=iatom+1
c endif
c ----- diagnostics
-c do kkk=1,nperm
-c write (iout,*) 'Ccopy and CREFcopy'
-c print '(i5,3f10.5,5x,3f10.5)',(k,(ccopy(j,k),j=1,3),
-c & (crefcopy(j,k),j=1,3),k=1,iatom)
-c write (iout,'(i5,3f10.5,5x,3f10.5)') (k,(ccopy(j,k),j=1,3),
-c & (crefcopy(j,k),j=1,3),k=1,iatom)
-c enddo
-c ----- end diagnostics
+C do kkk=1,nperm
+C write (iout,*) 'Ccopy and CREFcopy adasko',iatom
+C print '(i5,3f10.5,5x,3f10.5)',(k,(ccopy(j,k),j=1,3),
+C & (crefcopy(j,k),j=1,3),k=1,iatom)
+C write (iout,'(i5,3f10.5,5x,3f10.5)') (k,(ccopy(j,k),j=1,3),
+C & (crefcopy(j,k),j=1,3),k=1,iatom)
+C enddo
+Cc ----- end diagnostics
c do kkk=1,nperm
call fitsq(roznica,ccopy(1,1),crefcopy(1,1),iatom,
& przes,obrot,non_conv)
stop
#endif
endif
-c write (iout,*) "roznica", roznica,kkk
+C write (iout,*) "roznica", roznica,kkk
if (roznica.le.rminroz) rminroz=roznica
enddo
drms=dsqrt(dabs(rminroz))