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