pdbread-mult

[unres.git] / source / unres / src-HCD-5D / parmread.F

diff --git a/source/unres/src-HCD-5D/parmread.F b/source/unres/src-HCD-5D/parmread.F
index 2da8851..721d05b 100644 (file)
--- 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 
@@ -453,6 +462,7 @@ C here will be the apropriate recalibrating for D-aminoacid
 c      write (2,*) "Start reading THETA_PDB",ithep_pdb
       do i=1,ntyp
 c      write (2,*) 'i=',i
+      call flush(iout)
         read (ithep_pdb,*,err=111,end=111)
      &     a0thet(i),(athet(j,i,1,1),j=1,2),
      &    (bthet(j,i,1,1),j=1,2)
@@ -1452,7 +1462,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 +1479,7 @@ cc maxinter is maximum interaction sites
           enddo
         enddo
         enddo
+        call flush(iout)
       endif
 
 C 
@@ -1951,6 +1962,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 +2108,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 +2144,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