X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;ds=sidebyside;f=source%2Funres%2Fio_config.f90;h=0c2d645a05b9409c40fc9dd34fc420e58a816889;hb=refs%2Fheads%2Fdevel;hp=d3c036849558e2d4da734ec10d2da5459926f511;hpb=1a2fa7f9c1e5cb97c4cd9adef9b745dfd01e93df;p=unres4.git diff --git a/source/unres/io_config.f90 b/source/unres/io_config.f90 index d3c0368..0c2d645 100644 --- a/source/unres/io_config.f90 +++ b/source/unres/io_config.f90 @@ -30,7 +30,7 @@ ! !----------------------------------------------------------------------------- contains -#if .not. defined(WHAM_RUN) && .not. defined(CLUSTER) +#if !defined(WHAM_RUN) && !defined(CLUSTER) !----------------------------------------------------------------------------- ! bank.F io_csa !----------------------------------------------------------------------------- @@ -735,7 +735,7 @@ character(len=1) :: toronelet(-2:2) = (/"p","a","G","A","P"/) logical :: lprint,LaTeX real(kind=8),dimension(3,3,maxlob) :: blower !(3,3,maxlob) - real(kind=8),dimension(13) :: b + real(kind=8),dimension(13) :: bN character(len=3) :: lancuch !,ucase !el local variables integer :: m,n,l,i,j,k,iblock,lll,llll,ll,nlobi,mm @@ -1782,88 +1782,88 @@ do i=0,nloctyp-1 read (ifourier,*,end=115,err=115) - read (ifourier,*,end=115,err=115) (b(ii),ii=1,13) + read (ifourier,*,end=115,err=115) (bN(ii),ii=1,13) if (lprint) then write (iout,*) 'Type',i - write (iout,'(a,i2,a,f10.5)') ('b(',ii,')=',b(ii),ii=1,13) + write (iout,'(a,i2,a,f10.5)') ('bN(',ii,')=',bN(ii),ii=1,13) endif - B1(1,i) = b(3) - B1(2,i) = b(5) - B1(1,-i) = b(3) - B1(2,-i) = -b(5) + B1(1,i) = bN(3) + B1(2,i) = bN(5) + B1(1,-i) = bN(3) + B1(2,-i) = -bN(5) ! b1(1,i)=0.0d0 ! b1(2,i)=0.0d0 - B1tilde(1,i) = b(3) - B1tilde(2,i) =-b(5) - B1tilde(1,-i) =-b(3) - B1tilde(2,-i) =b(5) + B1tilde(1,i) = bN(3) + B1tilde(2,i) =-bN(5) + B1tilde(1,-i) =-bN(3) + B1tilde(2,-i) =bN(5) ! b1tilde(1,i)=0.0d0 ! b1tilde(2,i)=0.0d0 - B2(1,i) = b(2) - B2(2,i) = b(4) - B2(1,-i) =b(2) - B2(2,-i) =-b(4) + B2(1,i) = bN(2) + B2(2,i) = bN(4) + B2(1,-i) =bN(2) + B2(2,-i) =-bN(4) ! b2(1,i)=0.0d0 ! b2(2,i)=0.0d0 - CC(1,1,i)= b(7) - CC(2,2,i)=-b(7) - CC(2,1,i)= b(9) - CC(1,2,i)= b(9) - CC(1,1,-i)= b(7) - CC(2,2,-i)=-b(7) - CC(2,1,-i)=-b(9) - CC(1,2,-i)=-b(9) + CC(1,1,i)= bN(7) + CC(2,2,i)=-bN(7) + CC(2,1,i)= bN(9) + CC(1,2,i)= bN(9) + CC(1,1,-i)= bN(7) + CC(2,2,-i)=-bN(7) + CC(2,1,-i)=-bN(9) + CC(1,2,-i)=-bN(9) ! CC(1,1,i)=0.0d0 ! CC(2,2,i)=0.0d0 ! CC(2,1,i)=0.0d0 ! CC(1,2,i)=0.0d0 - Ctilde(1,1,i)=b(7) - Ctilde(1,2,i)=b(9) - Ctilde(2,1,i)=-b(9) - Ctilde(2,2,i)=b(7) - Ctilde(1,1,-i)=b(7) - Ctilde(1,2,-i)=-b(9) - Ctilde(2,1,-i)=b(9) - Ctilde(2,2,-i)=b(7) + Ctilde(1,1,i)=bN(7) + Ctilde(1,2,i)=bN(9) + Ctilde(2,1,i)=-bN(9) + Ctilde(2,2,i)=bN(7) + Ctilde(1,1,-i)=bN(7) + Ctilde(1,2,-i)=-bN(9) + Ctilde(2,1,-i)=bN(9) + Ctilde(2,2,-i)=bN(7) ! Ctilde(1,1,i)=0.0d0 ! Ctilde(1,2,i)=0.0d0 ! Ctilde(2,1,i)=0.0d0 ! Ctilde(2,2,i)=0.0d0 - DD(1,1,i)= b(6) - DD(2,2,i)=-b(6) - DD(2,1,i)= b(8) - DD(1,2,i)= b(8) - DD(1,1,-i)= b(6) - DD(2,2,-i)=-b(6) - DD(2,1,-i)=-b(8) - DD(1,2,-i)=-b(8) + DD(1,1,i)= bN(6) + DD(2,2,i)=-bN(6) + DD(2,1,i)= bN(8) + DD(1,2,i)= bN(8) + DD(1,1,-i)= bN(6) + DD(2,2,-i)=-bN(6) + DD(2,1,-i)=-bN(8) + DD(1,2,-i)=-bN(8) ! DD(1,1,i)=0.0d0 ! DD(2,2,i)=0.0d0 ! DD(2,1,i)=0.0d0 ! DD(1,2,i)=0.0d0 - Dtilde(1,1,i)=b(6) - Dtilde(1,2,i)=b(8) - Dtilde(2,1,i)=-b(8) - Dtilde(2,2,i)=b(6) - Dtilde(1,1,-i)=b(6) - Dtilde(1,2,-i)=-b(8) - Dtilde(2,1,-i)=b(8) - Dtilde(2,2,-i)=b(6) + Dtilde(1,1,i)=bN(6) + Dtilde(1,2,i)=bN(8) + Dtilde(2,1,i)=-bN(8) + Dtilde(2,2,i)=bN(6) + Dtilde(1,1,-i)=bN(6) + Dtilde(1,2,-i)=-bN(8) + Dtilde(2,1,-i)=bN(8) + Dtilde(2,2,-i)=bN(6) ! Dtilde(1,1,i)=0.0d0 ! Dtilde(1,2,i)=0.0d0 ! Dtilde(2,1,i)=0.0d0 ! Dtilde(2,2,i)=0.0d0 - EE(1,1,i)= b(10)+b(11) - EE(2,2,i)=-b(10)+b(11) - EE(2,1,i)= b(12)-b(13) - EE(1,2,i)= b(12)+b(13) - EE(1,1,-i)= b(10)+b(11) - EE(2,2,-i)=-b(10)+b(11) - EE(2,1,-i)=-b(12)+b(13) - EE(1,2,-i)=-b(12)-b(13) + EE(1,1,i)= bN(10)+bN(11) + EE(2,2,i)=-bN(10)+bN(11) + EE(2,1,i)= bN(12)-bN(13) + EE(1,2,i)= bN(12)+bN(13) + EE(1,1,-i)= bN(10)+bN(11) + EE(2,2,-i)=-bN(10)+bN(11) + EE(2,1,-i)=-bN(12)+bN(13) + EE(1,2,-i)=-bN(12)-bN(13) ! ee(1,1,i)=1.0d0 ! ee(2,2,i)=1.0d0 @@ -2872,7 +2872,7 @@ return end subroutine readpdb -#if .not. defined(WHAM_RUN) && .not. defined(CLUSTER) +#if !defined(WHAM_RUN) && !defined(CLUSTER) !----------------------------------------------------------------------------- ! readrtns_CSA.F !-----------------------------------------------------------------------------