ndis=ncon_work*(ncon_work-1)/2
call work_partition(.true.,ndis)
#endif
ndis=ncon_work*(ncon_work-1)/2
call work_partition(.true.,ndis)
#endif
ind1=0
DO I=1,NCON_work-1
if (mod(i,100).eq.0) print *,'Calculating RMS i=',i
ind1=0
DO I=1,NCON_work-1
if (mod(i,100).eq.0) print *,'Calculating RMS i=',i