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

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

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

%%TOP BOUNDARY CONDITIONS       PHI=1
phi(IO,JM)=1;

%%BOTTOM BOUNDARY CONDITION     dq