Adding cluster with newcorr

[unres.git] / PARAM / sc2scext.f

diff --git a/PARAM/sc2scext.f b/PARAM/sc2scext.f
new file mode 100644 (file)
index 0000000..1d23d3a
--- /dev/null
+++ b/PARAM/sc2scext.f
@@ -0,0 +1,37 @@
+      double precision eps(24,24),sigma0(24),sigii(24),chip(24),alp(24)
+      integer ind(4)
+      read(*,*) ipot,iexpon
+      read (*,'(4f20.10)')((eps(i,j),j=i,20),i=1,20)
+      read (*,'(4f20.10)')(sigma0(i),i=1,20),(sigii(i),i=1,20),
+     &  (chip(i),i=1,20),(alp(i),i=1,20)
+      ind(1)=2
+      ind(2)=7
+      ind(3)=9
+      ind(4)=9
+      do k=1,4
+        do j=1,k+20
+          if (j.gt.ind(k)) then
+            eps(j,k+20)=eps(ind(k),j)
+          else
+            eps(j,k+20)=eps(j,ind(k))
+          endif
+        enddo
+        do j=k+21,24
+          eps(k+20,j)=eps(ind(k),j)
+        enddo
+      enddo
+      do i=1,4
+        sigma0(i+20)=sigma0(ind(i)) 
+        sigii(i+20)=sigii(ind(i)) 
+        chip(i+20)=chip(ind(i)) 
+        alp(i+20)=alp(ind(i)) 
+      enddo
+      print *,ipot,iexpon
+      do i=1,24
+        print '(4f20.10)',(eps(i,j),j=i,24)
+      enddo
+      print '(4f20.10)',(sigma0(i),i=1,24)
+      print '(4f20.10)',(sigii(i),i=1,24)
+      print '(4f20.10)',(chip(i),i=1,24)
+      print '(4f20.10)',(alp(i),i=1,24)
+      end