add potential

.github/workflows/style/check_tabs.sh

@@ -0,0 +1,37 @@
+#!/usr/bin/env bash
+
+set -eu -o pipefail
+
+find . -type d \( -name .git \
+                  -o -path ./paper \
+                  -o -name build -o -name install \
+                  -o -name tmp_build_dir -o -name tmp_install_dir \
+                  -o -name Submodules \
+               \) -prune -o \
+       -type f \( \( -name "*.H" -o -name "*.h" -o -name "*.hh" -o -name "*.hpp" \
+                  -o -name "*.c" -o -name "*.cc" -o -name "*.cpp" -o -name "*.cxx" \
+                  -o -name "*.f" -o -name "*.F" -o -name "*.f90" -o -name "*.F90" \
+                  -o -name "*.py" \
+                  -o -name "*.md" -o -name "*.rst" \
+                  -o -name "*.sh" \
+                  -o -name "*.tex" \
+                  -o -name "*.txt" \
+                  -o -name "*.yml" \) \
+                 -a \( ! -name "*.tab.h" -a ! -name "*.tab.cpp" \
+                    -a ! -name "*.lex.h" -a ! -name "*.lex.cpp" \) \
+               \) \
+    -exec grep -Iq . {} \; \
+    -exec perl -i -pe's/\t/\ \ \ \ /g' {} +
+
+gitdiff=`git diff`
+
+if [ -z "$gitdiff" ]
+then
+    exit 0
+else
+    echo -e "\nTabs are not allowed. Changes suggested by"
+    echo -e "  ${0}\n"
+    git --no-pager diff
+    echo ""
+    exit 1
+fi

.github/workflows/style/check_trailing_whitespaces.sh

@@ -0,0 +1,37 @@
+#!/usr/bin/env bash
+
+set -eu -o pipefail
+
+find . -type d \( -name .git \
+                  -o -path ./paper \
+                  -o -name build -o -name install \
+                  -o -name tmp_build_dir -o -name tmp_install_dir \
+                  -o -name Submodules \
+               \) -prune -o \
+       -type f \( \( -name "*.H" -o -name "*.h" -o -name "*.hh" -o -name "*.hpp" \
+                  -o -name "*.c" -o -name "*.cc" -o -name "*.cpp" -o -name "*.cxx" \
+                  -o -name "*.f" -o -name "*.F" -o -name "*.f90" -o -name "*.F90" \
+                  -o -name "*.py" \
+                  -o -name "*.rst" \
+                  -o -name "*.sh" \
+                  -o -name "*.tex" \
+                  -o -name "*.txt" \
+                  -o -name "*.yml" \) \
+                 -a \( ! -name "*.tab.h" -a ! -name "*.tab.cpp" \
+                    -a ! -name "*.lex.h" -a ! -name "*.lex.cpp" \) \
+               \) \
+    -exec grep -Iq . {} \; \
+    -exec perl -i -pe's/[[:blank:]]+$//g' {} +
+
+gitdiff=`git diff`
+
+if [ -z "$gitdiff" ]
+then
+    exit 0
+else
+    echo -e "\nTrailing whitespaces at the end of a line are not allowed. Changes suggested by"
+    echo -e "  ${0}\n"
+    git --no-pager diff
+    echo ""
+    exit 1
+fi

exec/dean_kow/multilevel/amrvis.defaults

@@ -0,0 +1,10 @@
+palette               /home/almgren/bin/Palette
+initialderived        phi
+initialscale          6
+expansion             1
+numberformat          %8.6e
+maxpixmapsize         1000000
+reservesystemcolors   35
+showboxes             TRUE
+windowheight          650
+windowwidth           1100

src_compressible/main_driver.cpp

@@ -192,9 +192,13 @@ void main_driver(const char* argv)
 
     MultiFab surfcov;
     MultiFab dNadsdes;
+    MultiFab dNads;
+    MultiFab dNdes;
     if (n_ads_spec>0) {
         surfcov.define(ba,dmap,n_ads_spec,ngc);
         dNadsdes.define(ba,dmap,n_ads_spec,ngc);
+        dNads.define(ba,dmap,n_ads_spec,ngc);
+        dNdes.define(ba,dmap,n_ads_spec,ngc);
     }
 
     //statistics    
@@ -564,7 +568,7 @@ void main_driver(const char* argv)
 #ifdef MUI
         mui_push(cu, prim, dx, uniface, step);
 #endif
-        if (n_ads_spec>0) sample_MFsurfchem(cu, prim, surfcov, dNadsdes, geom, dt);
+        if (n_ads_spec>0) sample_MFsurfchem(cu, prim, surfcov, dNadsdes, dNads, dNdes, geom, dt);
 
         // FHD
         RK3step(cu, cup, cup2, cup3, prim, source, eta, zeta, kappa, chi, D, flux,
@@ -583,7 +587,7 @@ void main_driver(const char* argv)
 #endif
         if (n_ads_spec>0) {
 
-            update_MFsurfchem(cu, prim, surfcov, dNadsdes, geom);
+            update_MFsurfchem(cu, prim, surfcov, dNadsdes, dNads, dNdes, geom);
 
             conservedToPrimitive(prim, cu);