function [phi]=Murman_Cole_bc(phi,IM,JM,IL,JL,M_inf,dys,del_y_min)

%%LEFT BOUNDARY CONDITION       PHI=1
phi(JM,1)=1;

%%RIGHT BOUNDARY CONDITION      PHI=1
phi(JM,IL)=1;

%%TOP BOUNDARY CONDITIONS       PHI=1
phi(JL,IM)=1;

%%BOTTOM BOUNDARY CONDITION     dphi/dy=v_inf*dys
phi(1,IM)=-del_y_min*dys(IM)+phi(2,IM);