revert tau output & add example

Author: ChrisZYJ

examples/2D_axisym_hypoelasticity/case.py

@@ -0,0 +1,101 @@
+#!/usr/bin/env python3
+import json
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0,
+            "x_domain%end": 0.001,
+            "y_domain%beg": 0.0,
+            "y_domain%end": 0.0005,
+            "cyl_coord": "T",
+            "m": 49,
+            "n": 24,
+            "p": 0,
+            "dt": 2e-9,
+            "t_step_start": 0,
+            "t_step_stop": 400,
+            "t_step_save": 20,
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "model_eqns": 2,
+            "alt_soundspeed": "F",
+            "num_fluids": 2,
+            "mpp_lim": "F",
+            "mixture_err": "F",
+            "time_stepper": 3,
+            "weno_order": 5,
+            "weno_eps": 1.0e-16,
+            "teno": "T",
+            "teno_CT": 1e-8,
+            "null_weights": "F",
+            "mp_weno": "F",
+            "riemann_solver": 1,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -6,
+            "bc_x%end": -6,
+            "bc_y%beg": -2,
+            "bc_y%end": -6,
+            # Hypoelasticity
+            "hypoelasticity": "T",
+            "fd_order": 4,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "T",
+            # Patch 1 Liquid
+            "patch_icpp(1)%geometry": 3,
+            "patch_icpp(1)%x_centroid": 0.0005,
+            "patch_icpp(1)%y_centroid": 0.00025,
+            "patch_icpp(1)%length_x": 0.001,
+            "patch_icpp(1)%length_y": 0.0005,
+            "patch_icpp(1)%vel(1)": 0.0,
+            "patch_icpp(1)%vel(2)": 0.0,
+            "patch_icpp(1)%pres": 1e05,
+            "patch_icpp(1)%alpha_rho(1)": 1100 * (1.0 - 1e-6),
+            "patch_icpp(1)%alpha(1)": 1.0 - 1e-6,
+            "patch_icpp(1)%alpha_rho(2)": 1100 * 1e-6,
+            "patch_icpp(1)%alpha(2)": 1e-6,
+            # Patch 2 Solid
+            "patch_icpp(2)%alter_patch(1)": "T",
+            "patch_icpp(2)%geometry": 3,
+            "patch_icpp(2)%x_centroid": 0.0005,
+            "patch_icpp(2)%y_centroid": 0.000125,
+            "patch_icpp(2)%length_x": 0.0005,
+            "patch_icpp(2)%length_y": 0.00025,
+            "patch_icpp(2)%vel(1)": 0.0,
+            "patch_icpp(2)%vel(2)": 0.0,
+            "patch_icpp(2)%pres": 1e05,
+            "patch_icpp(2)%alpha_rho(1)": 1100 * 1e-6,
+            "patch_icpp(2)%alpha(1)": 1e-6,
+            "patch_icpp(2)%alpha_rho(2)": 1100 * (1.0 - 1e-6),
+            "patch_icpp(2)%alpha(2)": 1.0 - 1e-6,
+            # Acoustic source
+            "acoustic_source": "T",
+            "num_source": 1,
+            "acoustic(1)%support": 6,
+            "acoustic(1)%loc(1)": 0.00006,
+            "acoustic(1)%loc(2)": 0.0,
+            "acoustic(1)%pulse": 2,
+            "acoustic(1)%npulse": 1,
+            "acoustic(1)%mag": 1.0,
+            "acoustic(1)%gauss_sigma_time": 2e-8,
+            "acoustic(1)%foc_length": 0.00054,
+            "acoustic(1)%aperture": 0.0008,
+            "acoustic(1)%delay": 1e-7,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 4.4e00 * 5.57e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%G": 0.0,
+            "fluid_pp(2)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(2)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(2)%G": 1.0e09,
+        }
+    )
+)

src/post_process/m_data_output.fpp

@@ -344,7 +344,7 @@ contains
             if (pres_wrt .or. prim_vars_wrt) dbvars = dbvars + 1
 
             ! Elastic stresses
-            if (elasticity .and. (tau_wrt .or. prim_vars_wrt)) dbvars = dbvars + (num_dims*(num_dims + 1))/2
+            if (hypoelasticity) dbvars = dbvars + (num_dims*(num_dims + 1))/2
 
             ! Magnetic field
             if (mhd) then