enddo
enddo
endif
+C reading lipid parameters
+ read(iliptranpar,*) pepliptran
+ do i=1,ntyp
+ read(iliptranpar,*) liptranene(i)
+ enddo
+ close(iliptranpar)
#ifdef CRYST_THETA
C
C Read the parameters of the probability distribution/energy expression
read (itorp,*,end=113,err=113) ntortyp,nterm_old
if (lprint)write (iout,*) 'ntortyp,nterm',ntortyp,nterm_old
read (itorp,*,end=113,err=113) (itortyp(i),i=1,ntyp)
- itortyp(ntyp1)=0
do i=1,ntortyp
do j=1,ntortyp
read (itorp,'(a)')