diff --git a/analysis/flume_tank/02_trawl_geometry.R b/analysis/flume_tank/02_trawl_geometry.R
index 7433e9d..17a07a0 100644
--- a/analysis/flume_tank/02_trawl_geometry.R
+++ b/analysis/flume_tank/02_trawl_geometry.R
@@ -1,14 +1,44 @@
-# Flume tank analysis
+# Visualize flume tank results
 
 library(trawlmetrics)
 library(xlsx)
+library(ggthemes)
+
+# ggthemes::colorblind_pal()(8)
+#colors
+# "#000000" "#E69F00" "#56B4E9" "#009E73" "#F0E442" "#0072B2" "#D55E00" "#CC79A7"
 
 # Load flume tank data
 
-path_flume_xlsx <- here::here("analysis", "flume_tank", "data", "test_data.xlsx")
+path_flume_xlsx <- here::here("analysis", "flume_tank", "data", "flume_tank_data.xlsx")
+# sheet_flume_xlsx <- "data"
+sheet_flume_xlsx <- "example"
 
 flume_data <- xlsx::read.xlsx(file = path_flume_xlsx,
-                              sheetName = 'data')                                     
+                              sheetName = sheet_flume_xlsx) |>
+  dplyr::mutate(GEAR_NAME = "83-112",
+                type = "Flume tank")
+
+# Plot field vs. flume height/spread ----
+ggplot() +
+  geom_point(data = trawlmetrics::bts_geom |>
+               dplyr::filter(SURVEY_DEFINITION_ID != 78) |>
+               dplyr::mutate(type = "Survey"),
+             mapping = aes(x = NET_WIDTH_M,
+                           y = NET_HEIGHT_M,
+                           color = type),
+             alpha = 0.5,
+             size = 0.1) +
+  geom_point(data = flume_data,
+             mapping = aes(x = spread_u_wing_m,
+                           y = height_headline_m,
+                           color = type)) +
+  facet_wrap(~GEAR_NAME, scales = "free") +
+  scale_x_continuous(name = "Net width (m)") +
+  scale_y_continuous(name = "Net height (m)") +
+  scale_color_manual(name = "Type",
+                     values = c("Survey" = "#000000", 
+                                "Flume tank" = "#E69F00"))
     
 # 
 
@@ -16,7 +46,7 @@ ggplot() +
   geom_point(data = flume_data,
              mapping = aes(x = towing_speed_kn,
                            y = spread_mean_we_m)) +
-  facet_wrap(~round(door_m))
+  facet_wrap(~round(spread_door_m))
 
 ggplot(data = flume_data,
        mapping = aes(x = bridle_angle_deg,
@@ -32,10 +62,3 @@ ggplot(data = flume_data,
   geom_point() +
   geom_smooth() +
   facet_wrap(~towing_speed_kn)
-
-ggplot(data = flume_data,
-       mapping = aes(x = toowing,
-                     y = spread_mean_we_m,
-                     color = factor(towing_speed_kn))) +
-  geom_point() +
-  geom_smooth()