$calculatekp
$calculatekp all θ ϕ kx_nmax ky_nmax dk A1 A2 A3 A4 A5 A6 A7 A1 A2 A3 A4 A5 A6 A7 A1 A2 A3 A4 A5 A6 A7 A1 A2 A3 A4 A5 A6 A7 ..... (n layers) Incom D1 D2 D3 a c c11 c12 c13 c33 c44 c66