Merge branch 'jgfouca/more_gw_bridges' into next (PR #7461)

Add remaining cxx -> f90 bridges and bfb unit tests for GWD subroutines

GWD subroutines completed here:
* gw_convect_init
* gw_convect_project_winds
* gw_heating_depth
* gw_storm_speed
* gw_convect_gw_sources
* gw_beres_src
* gw_ediff
* gw_diff_tend
* gw_oro_src

With this PR, all the bridges we need to begin porting are in place (I think).

[BFB] except will require new bfb unit baselines

Author: jgfouca

components/eam/src/physics/cam/gw/gw_convect.F90

@@ -17,6 +17,11 @@ module gw_convect
 
 public :: gw_convect_init
 public :: gw_beres_src
+! Only public for bridging/testing
+public :: gw_convect_project_winds
+public :: gw_heating_depth
+public :: gw_storm_speed
+public :: gw_convect_gw_sources
 
 ! Dimension for heating depth.
 integer :: maxh
@@ -42,14 +47,19 @@ subroutine gw_convect_init( plev_src_wind, mfcc_in, errstring)
   real(r8), intent(in) :: mfcc_in(:,:,:)       ! Source spectra to keep as table
   character(len=*), intent(out) :: errstring   ! Report any errors from this routine
   integer :: ierr
+#ifndef SCREAM_CONFIG_IS_CMAKE
   integer :: k
+#endif
 
   errstring = ""
 
-#ifndef SCREAM_CONFIG_IS_CMAKE
+#ifdef SCREAM_CONFIG_IS_CMAKE
+  ! Just set k_src_wind to pver
+  k_src_wind = pver
+#else
   do k = 0, pver
-    if ( pref_edge(k+1) < plev_src_wind ) k_src_wind = k+1
- end do
+     if ( pref_edge(k+1) < plev_src_wind ) k_src_wind = k+1
+  end do
 #endif
 
 #ifndef SCREAM_CONFIG_IS_CMAKE

components/eamxx/scripts/gen_boiler.py

@@ -2494,10 +2494,14 @@ def gen_cxx_bfb_unit_impl(self, phys, sub, force_arg_data=None):
     }"""
 
         _, _, _, _, scalars, real_data, int_data, bool_data = group_data(arg_data, filter_out_intent="in")
-        check_scalars, check_arrays = "", ""
+        check_scalars, check_arrays, scalar_comments = "", "", ""
         for scalar in scalars:
             check_scalars += f"        REQUIRE(d_baseline.{scalar[0]} == d_test.{scalar[0]});\n"
 
+        _, _, _, all_dims, input_scalars, _, _, _ = group_data(arg_data, filter_out_intent="out")
+        all_scalar_inputs = all_dims + [scalar_name for scalar_name, _ in input_scalars]
+        scalar_comments = "// " + ", ".join(all_scalar_inputs)
+
         if has_array:
             c2f_transpose_code = "" if not need_transpose else \
 """
@@ -2531,6 +2535,8 @@ def gen_cxx_bfb_unit_impl(self, phys, sub, force_arg_data=None):
     // Set up inputs
     {data_struct} baseline_data[] = {{
       // TODO
+      {scalar_comments}
+      {data_struct}(),
     }};
 
     static constexpr Int num_runs = sizeof(baseline_data) / sizeof({data_struct});{gen_random}
@@ -2539,6 +2545,7 @@ def gen_cxx_bfb_unit_impl(self, phys, sub, force_arg_data=None):
     // inout data is in original state
     {data_struct} test_data[] = {{
       // TODO
+      {data_struct}(baseline_data[0]),
     }};
 
     // Read baseline data
@@ -2567,6 +2574,7 @@ def gen_cxx_bfb_unit_impl(self, phys, sub, force_arg_data=None):
       }}
     }}
   }} // run_bfb""".format(data_struct=data_struct,
+                          scalar_comments=scalar_comments,
                           sub=sub,
                           gen_random=gen_random,
                           c2f_transpose_code=c2f_transpose_code,

components/eamxx/src/physics/gw/CMakeLists.txt

@@ -23,6 +23,14 @@ if (NOT EAMXX_ENABLE_GPU OR Kokkos_ENABLE_CUDA_RELOCATABLE_DEVICE_CODE OR Kokkos
     eti/gw_gw_front_project_winds.cpp
     eti/gw_gw_front_gw_sources.cpp
     eti/gw_gw_cm_src.cpp
+    eti/gw_gw_convect_project_winds.cpp
+    eti/gw_gw_heating_depth.cpp
+    eti/gw_gw_storm_speed.cpp
+    eti/gw_gw_convect_gw_sources.cpp
+    eti/gw_gw_beres_src.cpp
+    eti/gw_gw_ediff.cpp
+    eti/gw_gw_diff_tend.cpp
+    eti/gw_gw_oro_src.cpp
   ) # GW ETI SRCS
 endif()
 

components/eamxx/src/physics/gw/eti/gw_gw_beres_src.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_beres_src_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_beres_src on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_convect_gw_sources.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_convect_gw_sources_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_convect_gw_sources on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_convect_project_winds.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_convect_project_winds_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_convect_project_winds on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_diff_tend.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_diff_tend_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_diff_tend on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_ediff.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_ediff_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_ediff on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_heating_depth.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_heating_depth_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_heating_depth on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream

components/eamxx/src/physics/gw/eti/gw_gw_oro_src.cpp

@@ -0,0 +1,14 @@
+#include "impl/gw_gw_oro_src_impl.hpp"
+
+namespace scream {
+namespace gw {
+
+/*
+ * Explicit instantiation for doing gw_oro_src on Reals using the
+ * default device.
+ */
+
+template struct Functions<Real,DefaultDevice>;
+
+} // namespace gw
+} // namespace scream