+ wcorr3_nucl=ww(38)
+ wcorr_nucl=ww(37)
+ wtor_d_nucl=ww(36)
+ wtor_nucl=ww(35)
+ wsbloc=ww(34)
+ wang_nucl=ww(33)
+ wbond_nucl=ww(32)
+ welsb=ww(31)
+ wvdwsb=ww(30)
+ welpsb=ww(29)
+ wvdwpsb=ww(28)
+ welpp=ww(27)
+ wvdwpp_nucl=ww(26)
+ wscbase=ww(46)
+ wpepbase=ww(47)
+ wscpho=ww(48)
+ wpeppho=ww(49)
+! print *,"KURWA",ww(48)
+! "WSCBASE ","WPEPBASE ","WSCPHO ","WPEPPHO "
+! "WVDWPP ","WELPP ","WVDWPSB ","WELPSB ","WVDWSB ",&
+! "WELSB ","WBOND ","WANG ","WSBLOC ","WTOR ",&
+! "WTORD ","WCORR ","WCORR3 ","WNULL ","WNULL ",&
+! "WCATPROT ","WCATCAT
+! +wbond_nucl*estr_nucl+wang_nucl*ebe_nucl&
+! +wvdwpp_nucl*evdwpp+welpp*eespp+wvdwpsb*evdwpsb+welpsb*eelpsb&
+! +wvdwsb*evdwsb+welsb*eelsb+wsbloc*esbloc+wtor_nucl*etors_nucl&
+! +wtor_d_nucl*etors_d_nucl+wcorr_nucl*ecorr_nucl+wcorr3_nucl*ecorr3_nucl&