function w=Phys_to_Riemann(p,v);
%
%	Converts (p,v) to Riemann coordinates
%
w(1,:)=v-0.25*log(p);
w(2,:)=v+0.25*log(p);
