def test_plot2dcontourf(self): xc, yc = self.mesh2d.centers() f = field.fdata(self.model, self.mesh2d, [ 1., euler.datavector(0., 0.), self.fn(np.sqrt(np.square(xc - 8) + np.square(yc - 4))) ]) fig, ax = plt.subplots(1, 1) f.contourf('pressure', axes=ax, style='r-') return fig
def fuv(x,y): vmag = .01 ; k = 10. return euler.datavector(0.*x+.4, 0.*x+.2)
def fuv(x, y): return euler.datavector(0. * x + .4, 0. * x + .2)
def fuv(x, y): r = np.sqrt((x - lx / 2)**2 + (y - ly / 2)**2) return euler.datavector( 0.5 - b / (2 * np.pi) * np.exp(0.5 * (1 - r**2)) * (y - ly / 2), 0.5 + b / (2 * np.pi) * np.exp(0.5 * (1 - r**2)) * (x - lx / 2))