sumQ_p(j,k,iparm)=sumQ_p(j,k,iparm)+q(j,t)*weight
sumQsq_p(j,k,iparm)=sumQsq_p(j,k,iparm)+q(j,t)**2*weight
sumEQ_p(j,k,iparm)=sumEQ_p(j,k,iparm)
sumQ_p(j,k,iparm)=sumQ_p(j,k,iparm)+q(j,t)*weight
sumQsq_p(j,k,iparm)=sumQsq_p(j,k,iparm)+q(j,t)**2*weight
sumEQ_p(j,k,iparm)=sumEQ_p(j,k,iparm)