if (iz_sc.lt.2) then
do ichain=1,nchain
indchain=tabpermchain(ichain,iperm)
+#define DEBUG
#ifdef DEBUG
write (iout,*) "ichain",ichain," indchain",indchain
write (iout,*) "chain_border",chain_border(1,ichain),
enddo
rmsside=dsqrt(rmsside/nnnn)
#ifdef DEBUG
- write (iout,*) iii,iref," nnnn",nnnn," rmsside",rmsside
+ write (iout,*) "nnnn",nnnn," rmsside",rmsside
#endif
+#undef DEBUG
rmscalc_side=rmsside
return
end