- call reada(controlcard,"HOMOL_DIST",waga_dist,1.0d0)
- call reada(controlcard,"HOMOL_ANGLE",waga_angle,1.0d0)
+ call readi(controlcard,"HOMOL_NSET",homol_nset,1)
+ if (homol_nset.gt.1)then
+ call card_concat(controlcard)
+ read(controlcard,*) (waga_dist(i),i=1,homol_nset)
+ call card_concat(controlcard)
+ read(controlcard,*) (waga_angle(i),i=1,homol_nset)
+ write(iout,*) "iset distance_weight angle_weight"
+ do i=1,homol_nset
+ write(iout,*) i,waga_dist(i),waga_angle(i)
+ enddo
+ else
+ iset=1
+ call reada(controlcard,"HOMOL_DIST",waga_dist(1),1.0d0)
+ call reada(controlcard,"HOMOL_ANGLE",waga_angle(1),1.0d0)
+ endif