double precision msq_cs(0:2,-nf:nf,-nf:nf) common/msq_cs/msq_cs