call multreadi(controlcard,"ISAMPL",isampl,nparmset,1)
write (iout,*) "MaxSlice",MaxSlice
call readi(controlcard,"NSLICE",nslice,1)
+ call reada(controlcard,"FRAC_CUTOFF",frac_cutoff,10.0d0)
+ write (iout,*) "RMSD cutoff for fraction calculations",frac_cutoff
call flush(iout)
if (nslice.gt.MaxSlice) then
write (iout,*) "Error: parameter out of range: NSLICE",nslice,