Example 1: Use of Oh04

1. Requirements

• Installation of SenSE

2. Oh04 for different incidence angles

import numpy as np
import matplotlib.pyplot as plt
from sense.surface import Oh92, Oh04
from sense.util import f2lam

def db(x):
    return 10.*np.log10(x)

f  = 5.3 # GHz
lam = f2lam(f)  # m
k = 2.*np.pi/lam
mv = np.linspace(0.01,0.4)
theta = np.deg2rad(40)

f = plt.figure()
ax1 = f.add_subplot(121)
ax2 = f.add_subplot(122)

s = 0.3/100.
ks = k * s
Oh = Oh04(mv, ks, theta)
ax1.plot(db(Oh.hv), mv, label=['s=0.3 cm'])
ax2.plot(db(Oh.p), mv, label=['s=0.3 cm'])
s = 1.2/100.
ks = k * s
Oh = Oh04(mv, ks, theta)
ax1.plot(db(Oh.hv), mv, label=['s=1.2 cm'])
ax2.plot(db(Oh.p), mv, label=['s=1.2 cm'])
s = 4.8/100.
ks = k * s
Oh = Oh04(mv, ks, theta)
ax1.plot(db(Oh.hv), mv, label=['s=4.8 cm'])
ax2.plot(db(Oh.p), mv, label=['s=4.8 cm'])

ax1.grid()
ax1.set_xlim(-34.,-14.)
ax1.set_ylim(0.01,0.4)
ax1.set_xticks(np.arange(-34, -12, 2))
ax2.grid()
ax2.set_xlim(-4.,1.)
ax2.set_ylim(0.01,0.4)
ax2.set_xticks(np.arange(-4, 1, 0.5))

ax1.set_xlabel('Sigma VH [dB]')
ax2.set_xlabel('p (cross-pol ratio) [dB]')

ax1.set_ylabel('Mv [cm³/cm³]')
ax2.set_ylabel('p (cross-pol ratio) [dB]')
ax1.legend()
ax2.legend()
plt.show()

print(theta)

0.6981317007977318

type(mv)

numpy.ndarray