1
Fork 0

Compare nolay and 2lay

This commit is contained in:
Edgar P. Burkhart 2022-03-30 12:22:57 +02:00
parent 858928353a
commit dbced643ca
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227
2 changed files with 12 additions and 12 deletions

View file

@ -9,7 +9,7 @@ mpi=8
[post]
inp=inp_post/spec_2lay
compare=inp_post/spec_2lay
out=out_post/spec_2lay
compare=inp_post/spec_nolay
out=out_post/spec_2layvsnolay
x0=-1250
t0=180

View file

@ -59,10 +59,10 @@ R = np.sqrt(
(np.abs(phi_eta[1]) + np.abs(phi_u[1]) - 2 * np.abs(phi_eta_u[1]))
/ (np.abs(phi_eta[1]) + np.abs(phi_u[1]) + 2 * np.abs(phi_eta_u[1]))
)
# R = np.sqrt(
# (1 + G**2 - 2 * G * np.cos(th_eta_u))
# / (1 + G**2 + 2 * G * np.cos(th_eta_u))
# )
#R = np.sqrt(
# (1 + G**2 - 2 * G * np.cos(th_eta_u))
# / (1 + G**2 + 2 * G * np.cos(th_eta_u))
#)
if config.has_option("post", "compare"):
inp_comp = pathlib.Path(config.get("post", "compare"))
@ -88,14 +88,14 @@ if config.has_option("post", "compare"):
G_ = H_ / U_
th_eta_u_ = np.angle(phi_eta_u_[1])
#R_ = np.sqrt(
# (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) - 2 * np.abs(phi_eta_u_[1]))
# / (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) + 2 * np.abs(phi_eta_u_[1]))
#)
R_ = np.sqrt(
(1 + G_**2 - 2 * G_ * np.cos(th_eta_u_))
/ (1 + G_**2 + 2 * G_ * np.cos(th_eta_u_))
(np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) - 2 * np.abs(phi_eta_u_[1]))
/ (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) + 2 * np.abs(phi_eta_u_[1]))
)
#R_ = np.sqrt(
# (1 + G_**2 - 2 * G_ * np.cos(th_eta_u_))
# / (1 + G_**2 + 2 * G_ * np.cos(th_eta_u_))
#)
# Plotting