do k=1,21
enetb(k,iii+1,iparm)=energia(k)
enddo
-#define DEBUG
#ifdef DEBUG
write (iout,'(2i5,f10.1,3e15.5)') i,iii,
& 1.0d0/(beta_h(ib,ipar)*1.987D-3),energia(0),eini,efree
call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev)
close(ipdb)
#endif
-#undef DEBUG
endif
enddo ! iparm