diff --git a/components/eamxx/src/physics/shoc/impl/shoc_main_impl.hpp b/components/eamxx/src/physics/shoc/impl/shoc_main_impl.hpp
index 1e9b4746701..11383e77245 100644
--- a/components/eamxx/src/physics/shoc/impl/shoc_main_impl.hpp
+++ b/components/eamxx/src/physics/shoc/impl/shoc_main_impl.hpp
@@ -686,6 +686,8 @@ Int Functions<S,D>::shoc_main(
                        w3_s, wqls_sec_s, brunt_s, isotropy_s);                // Diagnostic Output Variables
 
     shoc_output.pblh(i) = pblh_s;
+    shoc_output.ustar(i) = ustar_s;
+    shoc_output.obklen(i) = obklen_s;
   });
   Kokkos::fence();
 #else