modernize the f77 eigenvalue module, prettify

.github/workflows/pretty.yml

@@ -0,0 +1,19 @@
+name: Pretty
+
+on:
+  push:
+
+  pull_request:
+
+  workflow_dispatch:
+
+jobs:
+  docs:
+    name:    Code formatting
+    runs-on: ubuntu-latest
+
+    steps:
+    - uses: actions/checkout@v3
+
+    - name: Count
+      run: ./mfc.sh format diff

mfc.sh

@@ -159,10 +159,21 @@ elif [ "$1" == "format" ]; then
         pip3 install --upgrade fprettify
     fi
 
-    fprettify ${@:-src} --exclude "src/*/autogen" --recursive --silent \
-        --indent 4 --c-relations --enable-replacements --enable-decl \
-        --whitespace-comma 1 --whitespace-multdiv 1 --whitespace-plusminus 1 \
-        --case 1 1 1 1 --strict-indent
+    opts="--recursive --silent --indent 4 --c-relations --enable-replacements --enable-decl --whitespace-comma 0 --whitespace-multdiv 0 --whitespace-plusminus 1 --case 1 1 1 1 --strict-indent --line-length 1000"
+    if [ "$1" == "diff" ]; then 
+        shift
+        out=$(fprettify ${@:-src} --exclude "src/*/autogen" $opts -s -d)
+        if [ -z "${out}" ]; then
+            echo "Already pretty!"
+            exit 0
+        else
+            error "Not pretty, run ./mfc.sh format"
+
+            exit 1
+        fi
+    fi
+
+    fprettify ${@:-src} --exclude "src/*/autogen" $opts 
     ret="$?"
 
     if [ "$ret" != '0' ]; then

src/common/m_derived_types.fpp

@@ -152,6 +152,9 @@ module m_derived_types
         real(kind(0d0)), dimension(num_fluids_max) :: alpha
         real(kind(0d0)) :: gamma
         real(kind(0d0)) :: pi_inf !<
+        real(kind(0d0)) :: cv !<
+        real(kind(0d0)) :: qv !<
+        real(kind(0d0)) :: qvp !<
 
 
         !! Primitive variables associated with the patch. In order, these include
@@ -178,6 +181,9 @@ module m_derived_types
         real(kind(0d0)) :: gamma   !< Sp. heat ratio
         real(kind(0d0)) :: pi_inf  !< Liquid stiffness
         real(kind(0d0)), dimension(2) :: Re      !< Reynolds number
+        REAL(KIND(0d0)) :: cv      !< heat capacity
+        REAL(KIND(0d0)) :: qv      !< reference energy per unit mass for SGEOS, q (see Le Metayer (2004))
+        REAL(KIND(0d0)) :: qvp     !< reference entropy per unit mass for SGEOS, q' (see Le Metayer (2004))
         real(kind(0d0)) :: mul0    !< Bubble viscosity
         real(kind(0d0)) :: ss      !< Bubble surface tension
         real(kind(0d0)) :: pv      !< Bubble vapour pressure