X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fio_config.F90;h=61cb7b474b80e9ffbc68a6808efb912211efd6dc;hb=94d31d7bcffc5412d9a88a0aef46b68349b60cf3;hp=131c7894f4e71fb7024119012488a2607442d989;hpb=c6f305ff4a0b1fdc111572dc9692d866b9ffbb15;p=unres4.git

diff --git a/source/unres/io_config.F90 b/source/unres/io_config.F90
index 131c789..61cb7b4 100644
--- a/source/unres/io_config.F90
+++ b/source/unres/io_config.F90
@@ -867,9 +867,13 @@
           enddo
         enddo
       endif
+
+
+
+      if (.not.allocated(ichargecat)) allocate (ichargecat(ntyp_molec(5)))
        if (oldion.eq.1) then
             do i=1,ntyp_molec(5)
-             read(iion,*) msc(i,5),restok(i,5)
+             read(iion,*) msc(i,5),restok(i,5),ichargecat(i)
              print *,msc(i,5),restok(i,5)
             enddo
             ip(5)=0.2
@@ -881,6 +885,22 @@
             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))
@@ -3257,7 +3277,7 @@
       if (.not.allocated(chi2cat)) allocate(chi2cat(ntyp1,ntyp1)) !(ntyp,ntyp)
 
 
-      allocate (ichargecat(ntyp_molec(5)))
+            if (.not.allocated(ichargecat)) allocate (ichargecat(ntyp_molec(5)))
 ! i to SC, j to jon, isideocat - nazwa pliku z ktorego czytam parametry
        if (oldion.eq.0) then
             if (.not.allocated(icharge)) then ! this mean you are oprating in old sc-sc mode
@@ -3272,7 +3292,7 @@
              print *,msc(i,5),restok(i,5)
             enddo
             ip(5)=0.2
-
+!DIR$ NOUNROLL 
       do j=1,ntyp_molec(5)
        do i=1,ntyp
 !       do j=1,ntyp_molec(5)
@@ -3345,7 +3365,7 @@
 
        do i=1,ntyp
        do j=1,ntyp_molec(5)
-      if (i.eq.1) then
+      if (i.eq.10) then
       write (iout,*) 'i= ', i, ' j= ', j
       write (iout,*) 'epsi0= ', epscat(i,j)
       write (iout,*) 'sigma0= ', sigmacat(i,j)
@@ -3370,8 +3390,8 @@
       write (iout,*) 'epsin= ', epsintabcat(1,j), epsintabcat(j,1)
       write (iout,*) 'alphapol1= ',  alphapolcat(1,j)
       write (iout,*) 'alphapol2= ',  alphapolcat(j,1)
-      write (iout,*) 'w1= ', wqdipcat(1,1,j)
-      write (iout,*) 'w2= ', wqdipcat(2,1,j)
+      write (iout,*) 'w1= ', wqdipcat(1,i,j)
+      write (iout,*) 'w2= ', wqdipcat(2,i,j)
       write (iout,*) 'debaykapcat(i,j)= ',  debaykapcat(1,j)
       endif
 
@@ -4074,6 +4094,18 @@
          istype(i)=istype_temp(i)
         enddo
        enddo
+       if ((itype(1,1).eq.ntyp1).and.itype(2,5).ne.0) then
+! I have only ions now dummy atoms in the system        
+       molnum(1)=5
+       itype(1,5)=itype(2,5)
+       itype(1,1)=0
+       do i=2,nres
+         itype(i,5)=itype(i+1,5)
+       enddo
+       itype(nres,5)=0
+       nres=nres-1
+       nres_molec(1)=nres_molec(1)-1
+      endif
 !      if (itype(1,1).eq.ntyp1) then
 !        nsup=nsup-1
 !        nstart_sup=2
@@ -4432,7 +4464,7 @@
       endif
 
 ! CUTOFFF ON ELECTROSTATICS
-      call reada(controlcard,"R_CUT_ELE",r_cut_ele,15.0d0)
+      call reada(controlcard,"R_CUT_ELE",r_cut_ele,25.0d0)
       call reada(controlcard,"LAMBDA_ELE",rlamb_ele,0.3d0)
       write(iout,*) "R_CUT_ELE=",r_cut_ele
 ! Lipidic parameters
@@ -4573,6 +4605,7 @@
 !      enddo
       buff_shield=1.0d0
       endif
+      itime_mat=0
       return
       end subroutine read_control
 !-----------------------------------------------------------------------------
@@ -5323,6 +5356,8 @@
       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)
@@ -5382,6 +5417,8 @@
       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)