In [ ]:

    

import numpy as np
from anemoi import SimpleSource, MiniZephyr25D, AnalyticalHelmholtz


    

In [ ]:

    

systemConfig = {
    'dx':       1.,     # m
    'dz':       1.,     # m
    'c':        2500.,  # m/s
    'rho':      1.,     # kg/m^3
    'nx':       100,    # count
    'nz':       200,    # count
    'freq':     2e2,  # Hz
    'nky':      80,
}


    

In [ ]:

    

sp = MiniZephyr25D(systemConfig)
ss = SimpleSource(systemConfig)


    

In [ ]:

    

q = ss(25, 25)
u = sp*q


    

In [ ]:

    

%pylab inline


    

In [ ]:

    

clip = 1e-1
imshow(u.reshape((200,100)).real, cmap=cm.bwr, vmin=-clip, vmax=clip)