refactor and clean up

Author: ChrisZYJ

examples/1D_exp_tube_phasechange/case.py

@@ -84,7 +84,7 @@
 # save frequency = SF + 1 (because the initial state, 0.dat, is also saved)
 SF = 200
 
-# making Nt divisible By SF
+# making Nt divisible by SF
 tendA = (xe - xb) / ss * 0.25
 
 # 1 - ensure NtA is sufficient to go a little beyond tendA

examples/2D_axisym_shockwatercavity/case.py

@@ -115,7 +115,7 @@
 
 # I calculating tend twice; first is an estimate, second is
 # the actual value used. This is because I am getting errors in the
-# post process part every time I approximate the actual Nt By an integer
+# post process part every time I approximate the actual Nt by an integer
 # number (think of a smarter way).
 
 # dimensionless time
@@ -135,7 +135,7 @@
 # Save Frequency. Note that the number of autosaves will be SF + 1, as th IC (0.dat) is also saved
 SF = 400
 
-## making Nt divisible By SF
+## making Nt divisible by SF
 # 1 - ensure NtA goes slightly beyond tendA
 NtA = int(tendA // dt + 1)
 

examples/2D_laplace_pressure_jump/case.py

@@ -1,6 +1,6 @@
 #!/usr/bin/env python3
 # This case file demonstrates the Laplace pressure jump of a water droplet in air. The laplace pressure jump
-# in 2D is given By delta = sigma / r where delta is the pressure jump, sigma is the surface tension coefficient,
+# in 2D is given by delta = sigma / r where delta is the pressure jump, sigma is the surface tension coefficient,
 # and r is the radius of the droplet. The results of this simulation agree with theory to well within 1%
 # relative error.
 

examples/2D_phasechange_bubble/case.py

@@ -150,7 +150,7 @@
 # Critical time-step
 tc = 0.915 * R0 * math.sqrt(rho0wl1 / p01)
 
-# making Nt divisible By SF
+# making Nt divisible by SF
 # tendA = 1.5 * tc
 tend = 1.2 * tc
 

examples/2D_phasechange_bubble/casefile.py

@@ -150,7 +150,7 @@
 # Critical time-step
 tc = 0.915 * R0 * math.sqrt(rho0wl1 / p01)
 
-# making Nt divisible By SF
+# making Nt divisible by SF
 # tendA = 1.5 * tc
 tend = 1.2 * tc
 

examples/2D_shocktube_phasechange/case.py

@@ -138,7 +138,7 @@
 # save frequency = SF + 1 (because the initial state, 0.dat, is also saved)
 SF = 40
 
-# making Nt divisible By SF
+# making Nt divisible by SF
 tendA = (xe - xb) / ss * 5 / 24
 
 # 1 - ensure NtA is sufficient to go a little beyond tendA

examples/3D_phasechange_bubble/case.py

@@ -150,7 +150,7 @@
 # Critical time-step
 tc = 0.915 * R0 * math.sqrt(rho0wl1 / p01)
 
-# making Nt divisible By SF
+# making Nt divisible by SF
 # tendA = 1.5 * tc
 tend = 1.2 * tc
 

examples/3D_shockdroplet/case.py

@@ -125,7 +125,7 @@
 
 # I calculating tend twice; first is an estimate, second is
 # the actual value used. This is because I am getting errors in the
-# post process part every time I approximate the actual Nt By an integer
+# post process part every time I approximate the actual Nt by an integer
 # number (think of a smarter way).
 
 # dimensionless time
@@ -145,7 +145,7 @@
 # Save Frequency. Note that the number of autosaves will be SF + 1, as th IC (0.dat) is also saved
 SF = 400
 
-## making Nt divisible By SF
+## making Nt divisible by SF
 # 1 - ensure NtA goes slightly beyond tendA
 NtA = int(tendA // dt + 1)
 

examples/scaling/case.py

@@ -179,7 +179,7 @@ def nxyz_from_ncells(ncells: float) -> typing.Tuple[int, int, int]:
 
 # I calculate tend twice; first is an estimate, second is
 # the actual value used. This is because I am getting errors in the
-# post process part every time I approximate the actual Nt By an integer
+# post process part every time I approximate the actual Nt by an integer
 # number (think of a smarter way).
 
 # dimensionless time
@@ -197,7 +197,7 @@ def nxyz_from_ncells(ncells: float) -> typing.Tuple[int, int, int]:
 # Save Frequency. Note that the number of autosaves will be SF + 1, as th IC (0.dat) is also saved
 SF = 400
 
-## making Nt divisible By SF
+## making Nt divisible by SF
 # 1 - ensure NtA goes slightly beyond tendA
 NtA = int(tendA // dt + 1)
 

src/common/m_variables_conversion.fpp

@@ -1589,7 +1589,7 @@ contains
     end subroutine s_finalize_variables_conversion_module
 
 #ifndef MFC_PRE_PROCESS
-    subroutine s_compute_speed_of_sound(pres, rho, gamma, pi_inf, H, adv, vel_sum, c_c, c)
+    pure subroutine s_compute_speed_of_sound(pres, rho, gamma, pi_inf, H, adv, vel_sum, c_c, c)
 #ifdef _CRAYFTN
         !DIR$ INLINEALWAYS s_compute_speed_of_sound
 #else