Dear Manato,
ADDA cannot currently simulate a single polarization of the incident field
(
https://github.com/adda-team/adda/issues/30). It always considers two orthogonal polarizations and computes the
amplitude or Mueller scattering matrix - see Section 11 "What Scattering Quantities Are Calculated" of the manual. From
these matrices you can compute any quantities corresponding to any incident polarization that you want.
Files with cross sections are always named X and Y, but it does not mean that they are indeed x and y polarizations (in
the laboratory reference frame). Log file provides the full vectors for two polarizations corresponding to x and y
ones. For oblique incidence, one of them would be exactly the p polarization that you want.
Surface mode is compatible with -prop ..., so you can specify any incident direction. However, direction strictly
parallel to the surface is not allowed, as it would be ambiguous. You can use
-prop 1 0 -0.0000001
this will be almost parallel, but still above the surface.
Let me know, if you need more details on any of the above.
Maxim.