X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Fparmread.F;h=7550fd58b37d3c8bda754b010abba1d9bb7921b9;hb=f9783c27de185ff3e1d6e874a37a1a8c4d3c93a3;hp=2da8851ad4ca9b04f71c034176caccb5cb643b8a;hpb=c711143ad3fffb04d27b55aa823f399b8343c4c5;p=unres.git

diff --git a/source/unres/src-HCD-5D/parmread.F b/source/unres/src-HCD-5D/parmread.F
index 2da8851..7550fd5 100644
--- a/source/unres/src-HCD-5D/parmread.F
+++ b/source/unres/src-HCD-5D/parmread.F
@@ -69,6 +69,7 @@ C setenv LATEX YES
 C
       call getenv_loc("PRINT_PARM",lancuch)
       lprint = (ucase(lancuch).eq."YES" .or. ucase(lancuch).eq."Y")
+     & .and. (me.eq.king.or..not.out1file) .and. fg_rank.eq.0
       call getenv_loc("LATEX",lancuch)
       LaTeX = (ucase(lancuch).eq."YES" .or. ucase(lancuch).eq."Y")
 C
@@ -121,6 +122,7 @@ c
      &        vbldsc0(j,i),aksc(j,i),abond0(j,i)
           enddo
         enddo
+        call flush(iout)
       endif
 C reading lipid parameters
       if (lprint) then
@@ -132,6 +134,13 @@ C reading lipid parameters
        read(iliptranpar,*) liptranene(i)
        enddo
        close(iliptranpar)
+       if (lprint) then
+         write (iout,'(/a)') "Water-lipid transfer parameters"
+         write (iout,'(a3,3x,f10.5)') 'p',pepliptran
+         do i=1,ntyp
+           write (iout,'(a3,3x,f10.5)') restyp(i),liptranene(i)
+         enddo
+       endif
 #ifdef CRYST_THETA
 C
 C Read the parameters of the probability distribution/energy expression 
@@ -1452,7 +1461,7 @@ cc maxinter is maximum interaction sites
 
 #endif      
       if (lprint) then
-        write (iout,'(/a/)') 'Torsional constants:'
+        write (iout,'(/a/)') 'SCCor torsional constants:'
         do l=1,maxinter
         do i=1,nsccortyp
           do j=1,nsccortyp
@@ -1469,6 +1478,7 @@ cc maxinter is maximum interaction sites
           enddo
         enddo
         enddo
+        call flush(iout)
       endif
 
 C 
@@ -1951,6 +1961,7 @@ c
       call reada(weightcard,'WDFAB',wdfa_beta,0.0d0)
       call reada(weightcard,'SCAL14',scal14,0.4D0)
       call reada(weightcard,'SCALSCP',scalscp,1.0d0)
+      call reada(weightcard,'LIPSCALE',lipscale,1.0D0)
       call reada(weightcard,'CUTOFF',cutoff_corr,7.0d0)
       call reada(weightcard,'DELT_CORR',delt_corr,0.5d0)
       call reada(weightcard,'TEMP0',temp0,300.0d0)
@@ -2096,11 +2107,12 @@ C 12/1/95 Added weight for the multi-body term WCORR
       do i=1,maxres
         dyn_ss_mask(i)=.false.
       enddo
-      do i=1,maxres-1
-        do j=i+1,maxres
+      do i=1,max_cyst-1
+        do j=i+1,max_cyst
           dyn_ssbond_ij(i,j)=1.0d300
         enddo
       enddo
+      call flush(iout)
       call reada(weightcard,"HT",Ht,0.0D0)
       if (dyn_ss) then
         ss_depth=ebr/wsc-0.25*eps(1,1)
@@ -2131,7 +2143,7 @@ C 12/1/95 Added weight for the multi-body term WCORR
        write (iout,*) "BTRISS=", btriss
        write (iout,*) "CTRISS=", ctriss
        write (iout,*) "DTRISS=", dtriss
-       print *,'indpdb=',indpdb,' pdbref=',pdbref
+c       print *,'indpdb=',indpdb,' pdbref=',pdbref
       endif
       return
       end