2 layer config

swash/config.ini

@@ -4,14 +4,12 @@ out=../data/out
 [swash]
 input=sws/SPEC_buoy.sws
 path=/data/code/swash
-out=out/spec_ot
+out=out/spec_ot_2lay
 mpi=8
 
 [post]
-inp=inp_post/spec_ot
-compare=inp_post_nb
-out=out_post
-#nperseg=1024
-dt=0.25
+inp=inp_post/spec_ot_2lay
+compare=inp_post/spec_ot_2lay
+out=out_post/spec_ot_2lay
 x0=-1250
 t0=180

swash/processing/post.py

@@ -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(
-        (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]))
+       (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_))
-    # )
 
 
 # Plotting

swash/processing/sws_npz.py

@@ -30,14 +30,15 @@ np.save(inp.joinpath("tsec"), rsws.read_time(sws_out.joinpath("tsec.dat")))
 np.save(inp.joinpath("xp"), rsws.read_x(sws_out.joinpath("xp.dat")))
 
 var = {
-    "dep": rsws.read_scalar,
-    "botl": rsws.read_const,
-    "watl": rsws.read_scalar,
-    "pressk": rsws.read_scalar_lay,
-    "nhprsk": rsws.read_scalar_lay,
+    #"dep": rsws.read_scalar,
+    #"botl": rsws.read_const,
+    #"watl": rsws.read_scalar,
+    #"pressk": rsws.read_scalar_lay,
+    #"nhprsk": rsws.read_scalar_lay,
     "zk": rsws.read_scalar_lay,
-    "velk": rsws.read_vector_lay,
-    "vz": rsws.read_scalar_lay,
+    #"velk": rsws.read_vector_lay,
+    #"vz": rsws.read_scalar_lay,
+    #"vel": rsws.read_vector,
 }
 
 with ThreadPool() as pool: