|
Size: 410
Comment:
|
Size: 1430
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 13: | Line 13: |
{{{#!python ef eps_water_func(T,S,f): """Complex dielectric constant of seawater by Klein and Swift (1977)""" T=T-273.15 omega=2*pi*f eps_0=8.854*10**(-12) eps_inf=4.9 eps_s_T=87.134-1.949*10**(-1)*T-1.276*10**(-2)*T**2+2.491*10**(-4)*T**3 a_ST=1.+1.613*10**(-5)*S*T-3.656*10**(-3)*S+3.210*10**(-5)*S**2-4.232*10**(-7)*S**3 eps_s=eps_s_T*a_ST tau_T0=1.768*10**(-11)-6.086*10**(-13)*T+1.104*10**(-14)*T**2-8.111*10**(-17)*T**3 b_ST=1.+2.282*10**(-5)*S*T-7.638*10**(-4)*S-7.760*10**(-6)*S**2+1.105*10**(-8)*S**3 tau=tau_T0*b_ST delta=25-T beta=2.0333*10**(-2)+1.266*10**(-4)*delta+2.464*10**(-6)*delta**2-S*(1.849*10**(-5)-2.551*10**(-7)*delta+2.551*10**(-8)*delta**2) sigma_25S=S*(0.182521-1.46192*10**(-3)*S+2.09324*10**(-5)*S**2-1.28205*10**(-7)*S**3) sigma=sigma_25S*exp(-delta*beta) eps_water=eps_inf+(eps_s-eps_inf)/(1+1j*omega*tau)-1j*sigma/(omega*eps_0) return eps_water }}} |
Microwave remote sensing of the ocean
Exercises
PDF Brightness temperature of ocean surface.pdf
IPython notebook Brightness temperature of ocean surface.ipynb (see Python quick start)
1 ef eps_water_func(T,S,f):
2
3 """Complex dielectric constant of seawater by Klein and Swift (1977)"""
4
5 T=T-273.15
6
7 omega=2*pi*f
8
9 eps_0=8.854*10**(-12)
10
11 eps_inf=4.9
12
13
14
15 eps_s_T=87.134-1.949*10**(-1)*T-1.276*10**(-2)*T**2+2.491*10**(-4)*T**3
16
17 a_ST=1.+1.613*10**(-5)*S*T-3.656*10**(-3)*S+3.210*10**(-5)*S**2-4.232*10**(-7)*S**3
18
19 eps_s=eps_s_T*a_ST
20
21
22
23 tau_T0=1.768*10**(-11)-6.086*10**(-13)*T+1.104*10**(-14)*T**2-8.111*10**(-17)*T**3
24
25 b_ST=1.+2.282*10**(-5)*S*T-7.638*10**(-4)*S-7.760*10**(-6)*S**2+1.105*10**(-8)*S**3
26
27 tau=tau_T0*b_ST
28
29
30
31 delta=25-T
32
33 beta=2.0333*10**(-2)+1.266*10**(-4)*delta+2.464*10**(-6)*delta**2-S*(1.849*10**(-5)-2.551*10**(-7)*delta+2.551*10**(-8)*delta**2)
34
35 sigma_25S=S*(0.182521-1.46192*10**(-3)*S+2.09324*10**(-5)*S**2-1.28205*10**(-7)*S**3)
36
37 sigma=sigma_25S*exp(-delta*beta)
38
39 eps_water=eps_inf+(eps_s-eps_inf)/(1+1j*omega*tau)-1j*sigma/(omega*eps_0)
40
41 return eps_water