enddo
enddo
endif
+
+
+
if (.not.allocated(ichargecat)) allocate (ichargecat(ntyp_molec(5)))
if (oldion.eq.1) then
do i=1,ntyp_molec(5)
enddo
print *, catprm
endif
+ allocate(catnuclprm(14,ntyp_molec(2),ntyp_molec(5)))
+ do i=1,ntyp_molec(5)
+ do j=1,ntyp_molec(2)
+ write(iout,*) i,j
+ read(iionnucl,*) (catnuclprm(k,j,i),k=1,14)
+ enddo
+ enddo
+ write(*,'(3(5x,a6)11(7x,a6))') "w1 ","w2 ","epslj ","pis1 ", &
+ "sigma0","epsi0 ","chi1 ","chip1 ","sig ","b1 ","b2 ", &
+ "b3 ","b4 ","chis1 "
+ do i=1,ntyp_molec(5)
+ do j=1,ntyp_molec(2)
+ write(*,'(3(f10.3,x),11(f12.6,x),a3,2a)') (catnuclprm(k,j,i),k=1,14), &
+ restyp(i,5),"-",restyp(j,2)
+ enddo
+ enddo
! read (iion,*) (vcatprm(k),k=1,ncatprotpram)
!----------------------------------------------------
allocate(a0thet(-ntyp:ntyp),theta0(-ntyp:ntyp))
open (itube,file=tubename,status='old')
call getenv_loc('IONPAR',ionname)
open (iion,file=ionname,status='old')
+ call getenv_loc('IONPAR_NUCL',ionnuclname)
+ open (iionnucl,file=ionnuclname,status='old')
#else
open(1,file=pref_orig(:ilen(pref_orig))//'.inp',status='old',&
readonly)
open (itube,file=tubename,status='old',action='read')
call getenv_loc('IONPAR',ionname)
open (iion,file=ionname,status='old',action='read')
+ call getenv_loc('IONPAR_NUCL',ionnuclname)
+ open (iionnucl,file=ionnuclname,status='old',action='read')
#ifndef CRYST_SC
call getenv_loc('ROTPARPDB',rotname_pdb)