Compressible laminar boundary layer
 

I wrote this code which gives the disturbances of velocity and temperature on top of the steady flow, which is a steady laminar boundary layer with no pressure gradient (basically Blasius but with energy equation coupled). The code is parabolic and I am marching downstream from initial profiles. Variables are complex (also the BC) in the free-stream.

Everything is fine for Mach number smaller than 2.5, for higher ones the code blows up. Do you please have any suggestions for parabolic system with COMPLEX BC?

Thanx a lot!