change to rp2_em_fwave speed s(i,k)=1.0