call readi(controlcard,"HOMOL_NSET",homol_nset,1)
if (homol_nset.gt.1)then
- call readi(controlcard,"ISET",iset,homol_nset)
+ call readi(controlcard,"ISET",iset,1)
call card_concat(controlcard)
read(controlcard,*) (waga_homology(i),i=1,homol_nset)
else