Adding cluster with newcorr
[unres.git] / PARAM / sc2scext.f
1       double precision eps(24,24),sigma0(24),sigii(24),chip(24),alp(24)
2       integer ind(4)
3       read(*,*) ipot,iexpon
4       read (*,'(4f20.10)')((eps(i,j),j=i,20),i=1,20)
5       read (*,'(4f20.10)')(sigma0(i),i=1,20),(sigii(i),i=1,20),
6      &  (chip(i),i=1,20),(alp(i),i=1,20)
7       ind(1)=2
8       ind(2)=7
9       ind(3)=9
10       ind(4)=9
11       do k=1,4
12         do j=1,k+20
13           if (j.gt.ind(k)) then
14             eps(j,k+20)=eps(ind(k),j)
15           else
16             eps(j,k+20)=eps(j,ind(k))
17           endif
18         enddo
19         do j=k+21,24
20           eps(k+20,j)=eps(ind(k),j)
21         enddo
22       enddo
23       do i=1,4
24         sigma0(i+20)=sigma0(ind(i)) 
25         sigii(i+20)=sigii(ind(i)) 
26         chip(i+20)=chip(ind(i)) 
27         alp(i+20)=alp(ind(i)) 
28       enddo
29       print *,ipot,iexpon
30       do i=1,24
31         print '(4f20.10)',(eps(i,j),j=i,24)
32       enddo
33       print '(4f20.10)',(sigma0(i),i=1,24)
34       print '(4f20.10)',(sigii(i),i=1,24)
35       print '(4f20.10)',(chip(i),i=1,24)
36       print '(4f20.10)',(alp(i),i=1,24)
37       end