X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fcontrol.F90;h=2fd113e0e2b11296262b759f4034f978004c304c;hb=fa303bca64ee8c5549cdc075600e06fb60974c9c;hp=d70d14408556ec173468a96e10d3b2efa1336a37;hpb=3c2962eb07d557194f5368c06ed6a574a2f16e3b;p=unres4.git

diff --git a/source/unres/control.F90 b/source/unres/control.F90
index d70d144..2fd113e 100644
--- a/source/unres/control.F90
+++ b/source/unres/control.F90
@@ -237,11 +237,13 @@
 !      ielep_nucl= 131
       isidep_nucl=132
       iscpp_nucl=133
-
+      isidep_scbase=141
 
 
       iliptranpar=60
       itube=61
+!     IONS
+      iion=401
 #if defined(WHAM_RUN) || defined(CLUSTER)
 !
 ! setting the mpi variables for WHAM
@@ -696,7 +698,7 @@
         enddo
       enddo
       do i=nnt_molec(2),nct_molec(2)-1
-!        print*, "inloop2",i
+        print*, "inloop2",i
       call int_partition(ind_scint_nucl,my_sc_inds_nucl,my_sc_inde_nucl,i,&
            iatsc_s_nucl,iatsc_e_nucl,i+1,nct_molec(2),nint_gr_nucl(i), &
            istart_nucl(i,1),iend_nucl(i,1),*112)
@@ -728,7 +730,9 @@
       ispp=4 !?? wham ispp=2
 #ifdef MPI
 ! Now partition the electrostatic-interaction array
-      if (itype(nres_molec(1),1).eq.ntyp1_molec(1)) then
+      if (nres_molec(1).eq.0) then  
+       npept=0
+      elseif (itype(nres_molec(1),1).eq.ntyp1_molec(1)) then
       npept=nres_molec(1)-nnt-1
       else
       npept=nres_molec(1)-nnt
@@ -787,7 +791,7 @@
           ijunk,ielstart_nucl(i),ielend_nucl(i),*113)
       enddo ! i 
   113 continue
-      if (iatel_s.eq.0) iatel_s=1
+      if (iatel_s_nucl.eq.0) iatel_s_nucl=1
 
       nele_int_tot_vdw=(npept-2)*(npept-2+1)/2
 !      write (iout,*) "nele_int_tot_vdw",nele_int_tot_vdw
@@ -810,7 +814,7 @@
       if (iatel_s_vdw.eq.0) iatel_s_vdw=1
    15 continue
       if (iatel_s.eq.0) iatel_s=1
-
+      if (iatel_s_vdw.eq.0) iatel_s_vdw=1
       nele_int_tot_vdw_nucl=(npept_nucl-2)*(npept_nucl-2+1)/2
 !      write (iout,*) "nele_int_tot_vdw",nele_int_tot_vdw
       call int_bounds(nele_int_tot_vdw_nucl,my_ele_inds_vdw_nucl,&
@@ -831,7 +835,7 @@
 !        write (iout,*) i," ielstart_vdw",ielstart_vdw(i),
 !     &   " ielend_vdw",ielend_vdw(i)
       enddo ! i 
-      if (iatel_s_vdw.eq.0) iatel_s_vdw=1
+      if (iatel_s_vdw.eq.0) iatel_s_vdw_nucl=1
   115 continue
 
 #else
@@ -934,6 +938,7 @@
       enddo ! i
   114 continue
       print *, "after inloop3",iatscp_s_nucl,iatscp_e_nucl
+      if (iatscp_s_nucl.eq.0) iatscp_s_nucl=1
 #else
       iatscp_s=nnt
       iatscp_e=nct_molec(1)-1
@@ -1891,7 +1896,7 @@
 !-----------------------------------------------------------------------------
       subroutine setup_var
 
-      integer :: i
+      integer :: i,mnum
 !      implicit real*8 (a-h,o-z)
 !      include 'DIMENSIONS'
 !      include 'COMMON.IOUNITS'
@@ -1914,10 +1919,11 @@
       nvar=ntheta+nphi
       nside=0
       do i=2,nres-1
+      mnum=molnum(i)
 #ifdef WHAM_RUN
         if (itype(i,1).ne.10) then
 #else
-        if (itype(i,1).ne.10 .and. itype(i,1).ne.ntyp1) then
+        if (itype(i,1).ne.10 .and. itype(i,mnum).ne.ntyp1_molec(mnum) .and. mnum.ne.5) then
 #endif
 	  nside=nside+1
           ialph(i,1)=nvar+nside