!#endif
! include 'COMMON.IOUNITS'
! include 'COMMON.TIME1'
- logical :: lprn = .true.
+ logical :: lprn = .false.
logical :: osob
#ifndef FIVEDIAG
real(kind=8),dimension(2*nres,2*nres) :: Bmat,matmult
integer :: i,j,ind,ind1,k,l,ii,jj,m,m1,ii1,iti,nres2,ierr,nind,mark
real(kind=8) :: ip4
integer :: iz,mnum
+ print *,"just entered"
relfeh=1.0d-14
nres2=2*nres
-
+ print *,"FIVEDIAG"
write (iout,*) "before FIVEDIAG"
#ifndef FIVEDIAG
write (iout,*) "ALLOCATE"
+ print *,"ALLOCATE"
if(.not.allocated(Gcopy)) allocate(Gcopy(nres2,nres2)) !(maxres2,maxres2)
if(.not.allocated(Ghalf)) allocate(Ghalf(nres2*(nres2+1)/2)) !mmaxres2=(maxres2*(maxres+1)/2)
#endif