call flush(iout)
enddo
+ write (iout,*) "HOMOL_NSET",homol_nset
+ if (homol_nset.gt.1) then
+ write (iout,*) "HOMOL_NSET: nT_h",nT_h(iparm)
+ do ib=1,nT_h(iparm)
+ nR(ib,iparm)=homol_nset
+ write (iout,*) "iparm",iparm," ib",ib," nR",nR(ib,iparm)
+ enddo
+ endif
+
enddo
if (hamil_rep) then