+ print *,"before inloop3",iatscp_s,iatscp_e,iscp_nucl
+ nscp_int_tot_nucl=(npept_nucl-iscp_nucl+1)*(npept_nucl-iscp_nucl+1)
+ call int_bounds(nscp_int_tot_nucl,my_scp_inds_nucl,my_scp_inde_nucl)
+ if (lprint) write (iout,*) 'Processor',fg_rank,' CG group',kolor,&
+ ' absolute rank',myrank,&
+ ' nscp_int_tot',nscp_int_tot_nucl,' my_scp_inds=',my_scp_inds_nucl,&
+ ' my_scp_inde',my_scp_inde_nucl
+ print *,"nscp_int_tot_nucl",nscp_int_tot_nucl,my_scp_inds_nucl,my_scp_inde_nucl
+ iatscp_s_nucl=0
+ iatscp_e_nucl=0
+ ind_scpint_nucl=0
+ ind_scpint_old_nucl=0
+ do i=nnt_molec(2),nct_molec(2)-1
+ print *,"inloop3",i,nnt_molec(2)+iscp,nct_molec(2)-iscp
+ if (i.lt.nnt_molec(2)+iscp) then
+!d write (iout,*) 'i.le.nnt+iscp'
+ call int_partition(ind_scpint_nucl,my_scp_inds_nucl,&
+ my_scp_inde_nucl,i,iatscp_s_nucl,iatscp_e_nucl,i+iscp,&
+ nct_molec(2),nscp_gr_nucl(i),iscpstart_nucl(i,1),&
+ iscpend_nucl(i,1),*114)
+ else if (i.gt.nct_molec(2)-iscp) then
+!d write (iout,*) 'i.gt.nct-iscp'
+ call int_partition(ind_scpint_nucl,my_scp_inds_nucl,&
+ my_scp_inde_nucl,i,&
+ iatscp_s_nucl,iatscp_e_nucl,nnt_molec(2),i-iscp,nscp_gr_nucl(i),&
+ iscpstart_nucl(i,1),&
+ iscpend_nucl(i,1),*114)
+ else
+ call int_partition(ind_scpint_nucl,my_scp_inds_nucl,&
+ my_scp_inde_nucl,i,iatscp_s_nucl,iatscp_e_nucl,nnt_molec(2),&
+ i-iscp,nscp_gr_nucl(i),iscpstart_nucl(i,1),&
+ iscpend_nucl(i,1),*114)
+ ii=nscp_gr_nucl(i)+1
+ call int_partition(ind_scpint_nucl,my_scp_inds_nucl,&
+ my_scp_inde_nucl,i,iatscp_s_nucl,iatscp_e_nucl,i+iscp,&
+ nct_molec(2),nscp_gr_nucl(i),iscpstart_nucl(i,ii),&
+ iscpend_nucl(i,ii),*114)
+ endif
+ enddo ! i
+ 114 continue
+ print *, "after inloop3",iatscp_s_nucl,iatscp_e_nucl
+ if (iatscp_s_nucl.eq.0) iatscp_s_nucl=1