call reada(weightcard,'TEMP0',temp0,300.0d0)
if (index(weightcard,'SOFT').gt.0) ipot=6
call reada(weightcard,'WBOND_NUCL',wbond_nucl,1.0D0)
+ call reada(weightcard,'WCATCAT',wcatcat,0.0d0)
+ call reada(weightcard,'WCATPROT',wcatprot,0.0d0)
! 12/1/95 Added weight for the multi-body term WCORR
call reada(weightcard,'WCORRH',wcorr,1.0D0)
if (wcorr4.gt.0.0d0) wcorr=wcorr4