constants
p_om=(-1537:1536)/1537/p_om0(1)*pi/.6e-6';
