Merge pull request #101 from wilfonba/viscousFix

Author: sbryngelson

src/simulation/p_main.fpp

@@ -156,7 +156,7 @@ program p_main
     if (monopole) then
         call s_initialize_monopole_module()
     end if
-    if (any(Re_size > 1)) then
+    if (any(Re_size > 0)) then
         call s_initialize_viscous_module()
     end if
     call s_initialize_rhs_module()

tests/0A022883/golden.txt

@@ -2,7 +2,7 @@
 
 from ..    import common
 from .case import TestCase, create_case, CaseGeneratorStack
-
+from ..printer import cons
 
 def get_bc_mods(bc: int, dimInfo):
     params = {}
@@ -130,6 +130,15 @@ def alter_num_fluids(dimInfo, dimParams):
 
             alter_riemann_solvers(num_fluids)
 
+            if num_fluids == 1:
+                stack.push(f"Viscous", {
+                    'fluid_pp(1)%Re(1)' : 0.0001, 'dt' : 1e-11})
+
+                cases.append(create_case(stack, "",             {'weno_Re_flux': 'F'}))             
+                cases.append(create_case(stack, "weno_Re_flux", {'weno_Re_flux': 'T'}))
+
+                stack.pop()
+
             if num_fluids == 2:
                 stack.push(f"Viscous", {
                     'fluid_pp(1)%Re(1)' : 0.001, 'fluid_pp(1)%Re(2)' : 0.001,
@@ -141,8 +150,8 @@ def alter_num_fluids(dimInfo, dimParams):
                 stack.pop()
                 stack.pop()
 
-            stack.pop()        
-    
+            stack.pop() 
+
     def alter_2d(dimInfo, dimParams):
         stack.push("Axisymmetric", {
             'num_fluids' : 2, 'bc_y%beg': -2, 'cyl_coord': 'T',