from pylayers.antprop.slab import *

S = SlabDB()
fGHz = np.linspace(0.1,100,100)
theta = 0
S.add(name='5mm-thickGlass',lmatname=['ITU_GLASS'],lthick=[0.005])
S.add(name='5cm-thickGlass',lmatname=['ITU_GLASS'],lthick=[0.05])
S.add(name='5mm-Concrete',lmatname=['ITU_CONCRETE'],lthick=[0.005])
S.add(name='5cm-Concrete',lmatname=['ITU_CONCRETE'],lthick=[0.05])
S.add(name='0d5mum-Metal',lmatname=['ITU_METAL'],lthick=[0.5e-9])
S['5mm-thickGlass'].eval(fGHz,theta)
S['5cm-thickGlass'].eval(fGHz,theta)
S['5mm-Concrete'].eval(fGHz,theta)
S['5cm-Concrete'].eval(fGHz,theta)
S['0d5mum-Metal'].eval(fGHz,theta)
f,a = S['5mm-thickGlass'].plotwrt('f')
f,a = S['5cm-thickGlass'].plotwrt('f',fig=f,ax=a)
f,a = S['5mm-Concrete'].plotwrt('f',fig=f,ax=a)
f,a = S['5cm-Concrete'].plotwrt('f',fig=f,ax=a)
f,a = S['0d5mum-Metal'].plotwrt('f',fig=f,ax=a)

plt.show()