subroutine promienie(*)
implicit none
include 'DIMENSIONS'
+ include 'DIMENSIONS.FREE'
include 'COMMON.CONTROL'
include 'COMMON.INTERACT'
include 'COMMON.IOUNITS'
enddo
close (isidep1)
do i=1,ntyp1
- if (i.eq.10 .or. i.eq.21) then
+ if (i.eq.10 .or. i.eq.ntyp1) then
dsc_inv(i)=0.0d0
else
dsc_inv(i)=1.0d0/dsc(i)