New Example CI + Fixes to some broken example!

.github/workflows/example.yml

@@ -0,0 +1,49 @@
+name: 'Example Suite'
+
+on: [push, pull_request, workflow_dispatch]
+
+jobs:
+  file-changes:
+    name: Detect File Changes
+    runs-on: 'ubuntu-latest'
+    outputs: 
+      checkall: ${{ steps.changes.outputs.checkall }}
+    steps:
+      - name: Clone
+        uses: actions/checkout@v4
+
+      - name: Detect Changes
+        uses: dorny/paths-filter@v3
+        id: changes
+        with: 
+          filters: ".github/file-filter.yml"
+
+  github:
+    name: Github
+    if: needs.file-changes.outputs.checkall == 'true'
+    needs: file-changes
+    continue-on-error: true
+    runs-on: macos-latest
+
+    steps:
+      - name: Clone
+        uses: actions/checkout@v4
+
+      - name: Setup MacOS
+        run: |
+          brew install wget make python make cmake coreutils gcc@14
+          echo "CC=gcc-14"      >> $GITHUB_ENV
+          echo "CXX=g++-14"     >> $GITHUB_ENV
+          echo "FC=gfortran-14" >> $GITHUB_ENV
+
+      - name: (MacOS) Build OpenMPI
+        run: |
+          brew install mpich
+
+      - name: Build
+        run:  |
+          /bin/bash mfc.sh build -j $(nproc) --no-debug --mpi
+
+      - name: Test
+        run:  |
+          /bin/bash mfc.sh test -o Example -j $(nproc) -r --generate

examples/1D_bubblescreen/case.py

@@ -193,6 +193,7 @@
     # Acoustic source ==========================================
     'Monopole'                     : 'T',
     'num_mono'                     : 1,
+    'Mono(1)%support'              : 1,
     'Mono(1)%loc(1)'               : -5.E-03/x0,
     'Mono(1)%npulse'               : 1,
     'Mono(1)%dir'                  : 1.,

examples/1D_exp_bubscreen/case.py

@@ -241,6 +241,7 @@
     # Acoustic source ==========================================
     'Monopole'                  : 'T',
     'num_mono'                  : 1,
+    'Mono(1)%support'           : 1,
     'Mono(1)%loc(1)'            : -0.05/x0,
     'Mono(1)%npulse'            : 1,
     'Mono(1)%dir'               : 1.,

examples/1D_impact/case.py

@@ -33,18 +33,18 @@
                     'model_eqns'                   : 2,                        
                     'alt_soundspeed'               : 'F',                      
                     'num_fluids'                   : 1,                        
-		    'adv_alphan'                   : 'T',                      
-		    'mpp_lim'                      : 'F',                      
-		    'mixture_err'                  : 'F',                      
-		    'time_stepper'                 : 3,                                               
+                    'adv_alphan'                   : 'T',                      
+                    'mpp_lim'                      : 'F',                      
+                    'mixture_err'                  : 'F',                      
+                    'time_stepper'                 : 3,                                               
                     'weno_order'                   : 3,                        
                     'weno_eps'                     : 1.E-16,
-		    'weno_Re_flux'                 : 'F',  
-    	 	    'weno_avg'                     : 'F',
+                    'weno_Re_flux'                 : 'F',  
+                    'weno_avg'                     : 'F',
                     'mapped_weno'                  : 'F',                      
                     'null_weights'                 : 'F',                      
                     'mp_weno'                      : 'F',                      
-		    'riemann_solver'               : 1,                        
+                    'riemann_solver'               : 1,                        
                     'wave_speeds'                  : 1,                        
                     'avg_state'                    : 2,                        
                     'bc_x%beg'                     : -3,                       
@@ -59,7 +59,7 @@
                     'format'                       : 1,                        
                     'precision'                    : 2,                        
                     'prim_vars_wrt'                :'T',                       
-		    'parallel_io'                  :'F',                       
+                    'parallel_io'                  :'F',                       
 		    # ==========================================================
 
 		    # Patch 1 L ================================================

examples/1D_qbmm/case.py

@@ -122,7 +122,7 @@
     'format'                       : 1,
     'precision'                    : 2,
     'prim_vars_wrt'                :'T',
-    'parallel_io'                  :'F',
+    'parallel_io'                  :'T',
     'fd_order'                     : 1,
     #'schlieren_wrt'                :'T',
     'probe_wrt'                    :'T',

examples/1D_shuosher_teno/case.py

@@ -51,7 +51,7 @@
     # ==========================================================================
 
     # Formatted Database Files Structure Parameters ============================
-    'format'                       : 2,
+    'format'                       : 1,
     'precision'                    : 2,
     'prim_vars_wrt'                :'T',
     'rho_wrt'                      :'T',

examples/1D_shuosher_wenojs/case.py

@@ -50,7 +50,7 @@
     # ==========================================================================
 
     # Formatted Database Files Structure Parameters ============================
-    'format'                       : 2,
+    'format'                       : 1,
     'precision'                    : 2,
     'prim_vars_wrt'                :'T',
     'rho_wrt'                      :'T',

examples/1D_shuosher_wenom/case.py

@@ -50,7 +50,7 @@
     # ==========================================================================
 
     # Formatted Database Files Structure Parameters ============================
-    'format'                       : 2,
+    'format'                       : 1,
     'precision'                    : 2,
     'prim_vars_wrt'                :'T',
     'rho_wrt'                      :'T',

examples/1D_shuosher_wenoz/case.py

@@ -50,7 +50,7 @@
     # ==========================================================================
 
     # Formatted Database Files Structure Parameters ============================
-    'format'                       : 2,
+    'format'                       : 1,
     'precision'                    : 2,
     'prim_vars_wrt'                :'T',
     'rho_wrt'                      :'T',

examples/1D_titarevtorro/case.py

@@ -67,7 +67,7 @@
     # ==========================================================================
 
     # Patch 2 R (-4.5 < x < 5) =================================================
-    'patch_icpp(2)%geometry'        : 15,
+    'patch_icpp(2)%geometry'        : 1,
     'patch_icpp(2)%x_centroid'      : 0.25,
     'patch_icpp(2)%length_x'        : 9.5,
     'patch_icpp(2)%vel(1)'          : 0.,

examples/2D_isentropicvortex/case.py

@@ -49,8 +49,8 @@
     'x_domain%end'                 : 3,
     'y_domain%beg'                 : -3,
     'y_domain%end'                 : 3,
-    'stretch_x'                    : True,
-    'stretch_y'                    : True,
+    'stretch_x'                    : 'T',
+    'stretch_y'                    : 'T',
     'loops_x'                      : 2,
     'loops_y'                      : 2,
     'a_x'                          : 1.03,

examples/2D_whale_bubble_annulus/case.py

@@ -186,6 +186,7 @@
 
     'Monopole'                     : 'T',
     'num_mono'                     : 1,
+    'Mono(1)%support'              : 1,
     'Mono(1)%loc(1)'               : -1.5,
     'Mono(1)%loc(2)'               : 0.0,
     'Mono(1)%npulse'               : 4,

src/pre_process/m_patches.fpp

@@ -969,9 +969,6 @@ contains
 
                         @:analytical()
 
-                        call s_assign_patch_primitive_variables(patch_id, i, j, 0, &
-                                                                eta, q_prim_vf, patch_id_fp)
-
                         if ((q_prim_vf(1)%sf(i, j, 0) < 1.e-10) .and. (model_eqns == 4)) then
                             !zero density, reassign according to Tait EOS
                             q_prim_vf(1)%sf(i, j, 0) = &
@@ -2067,6 +2064,9 @@ contains
                     ! if defining them analytically.
                     @:analytical()
 
+                    call s_assign_patch_primitive_variables(patch_id, i, j, k, &
+                                                            eta, q_prim_vf, patch_id_fp)
+
                 end do; end do; end do
 
         if (proc_rank == 0) then