call reada(weightcard,'WCATPROT',wcatprot,0.0d0)
call reada(weightcard,'WPEPBASE',wpepbase,1.0d0)
call reada(weightcard,'WSCPHO',wscpho,1.0d0)
+ call reada(weightcard,'WPEPPHO',wpeppho,1.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