add regression test with init_wave_profile

test/CMakeLists.txt

@@ -246,6 +246,7 @@ add_test_re(zalesak_disk_scalar_vel)
 add_test_re(rain_drop)
 add_test_re(inertial_drop)
 add_test_re(ow_linear)
+add_test_re(ow_linear_init_waves)
 add_test_re(ow_stokes)
 add_test_re(scalar_advection_uniform)
 add_test_re(scalar_advection_refined)

test/test_files/ow_linear_init_waves/ow_linear_init_waves.inp

@@ -0,0 +1,70 @@
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            SIMULATION STOP            #
+#.......................................#
+time.stop_time               =   1     # Max (simulated) time to evolve
+time.max_step                =   20          # Max number of time steps
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#         TIME STEP COMPUTATION         #
+#.......................................#
+time.fixed_dt         =   0.01        # Use this constant dt if > 0
+time.cfl              =   0.95         # CFL factor
+time.use_force_cfl= false
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            INPUT AND OUTPUT           #
+#.......................................#
+time.plot_interval            =  10       # Steps between plot files
+time.checkpoint_interval      =  -1       # Steps between checkpoint files
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#               PHYSICS                 #
+#.......................................#
+transport.model = TwoPhaseTransport
+transport.viscosity_fluid1=0.0
+transport.viscosity_fluid2=0.0
+transport.turbulent_prandtl = 0.3333
+turbulence.model = Laminar 
+
+incflo.physics = MultiPhase OceanWaves 
+OceanWaves.label = Wave1
+OceanWaves.Wave1.type = LinearWaves
+OceanWaves.Wave1.initialize_wave_field = true
+OceanWaves.Wave1.wave_height=0.03
+OceanWaves.Wave1.wave_length=2.0
+OceanWaves.Wave1.water_depth=1
+OceanWaves.Wave1.relax_zone_gen_length=1.5
+OceanWaves.Wave1.numerical_beach_length=1.5
+OceanWaves.Wave1.numerical_beach_length_factor=2.0
+MultiPhase.density_fluid1=1000.
+MultiPhase.density_fluid2=1.
+ICNS.source_terms = GravityForcing
+ICNS.use_perturb_pressure = true
+ICNS.reconstruct_true_pressure = true
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#        ADAPTIVE MESH REFINEMENT       #
+#.......................................#
+amr.n_cell              = 192 16 32    # Grid cells at coarsest AMRlevel
+amr.max_level = 1
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#              GEOMETRY                 #
+#.......................................#
+geometry.prob_lo        =     0.0  -0.3  -1   # Lo corner coordinates
+geometry.prob_hi        =     7.2   0.3  0.2  # Hi corner coordinates
+geometry.is_periodic    =     0     1     0   # Periodicity x y z (0/1)
+
+xlo.type =     "wave_generation"
+xhi.type =     "slip_wall"
+
+zlo.type =     "slip_wall"
+zhi.type =     "slip_wall"
+
+tagging.labels = refine0
+tagging.refine0.type = GeometryRefinement
+tagging.refine0.shapes = b0
+tagging.refine0.level  = 0
+
+tagging.refine0.b0.type = box                 
+tagging.refine0.b0.origin = 0.0 -0.3 -0.1
+tagging.refine0.b0.xaxis = 7.2 0.0 0.0
+tagging.refine0.b0.yaxis = 0.0 0.6 0.0
+tagging.refine0.b0.zaxis = 0.0 0.0 0.2