Plot an Internal Wave

This notebook plots Figure 2 in  "Anelastic Equations for the Atmosphere"

The Functions

?VectorPlot

uwplot[α_]:=VectorPlot[{-Cos[z] Sin[α x],Sin[z] Cos[α x]},{x,0,((2 π)/α)},{z,0,π},VectorPoints->{13,7},VectorStyle->Black]

phiplot[α_]:=ContourPlot[-((Cos[z] Sin[α x])/(1+(α^2))),{x,0,((2 π)/α)},{z,0,π},ColorFunction->MyTable,ContourStyle->None,AspectRatio->(α/2),Frame->False]

bplot[α_]:=ContourPlot[z+.1Sin[z]*Sin[α x],
		{x,0,2Pi/α},{z,0,Pi},ContourStyle->MyLine,ContourShading->False,AspectRatio->α/2,Frame->False]

BlueGrayRed[z_]:=RGBColor[z,2*z*(1-z),1-z]

MyGray[z_]:=RGBColor[.3+.6*z,.3+.6*z,.3+.6*z]

The Plots

MyLine=Green;

MyTable=BlueGrayRed;

uw1=uwplot[1];

ph1=phiplot[1];

b1=bplot[1];

super=Show[{ph1,b1,uw1}]

Export["colorwave.eps",super,"EPS"];