call reada(controlcard,"HOMOL_SCD",waga_d,1.0d0) ! new
call reada(controlcard,'DIST_CUT',dist_cut,5.0d0) ! for diff ways of calc sigma
call reada(controlcard,'DIST2_CUT',dist2_cut,9999.0d0)
- call readi(controlcard,"HOMOL_SET",homol_nset,1)
+ call readi(controlcard,"HOMOL_NSET",homol_nset,1)
read2sigma=(index(controlcard,'READ2SIGMA').gt.0)
call readi(controlcard,"IHSET",ihset,1)
if (homol_nset.gt.1)then