Avoid division by zero

ElmerCSC · Aug 26, 2024 · ec82dfe · ec82dfe
1 parent a178076
commit ec82dfe
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/fem/src/modules/VectorHelmholtz.F90 b/fem/src/modules/VectorHelmholtz.F90
@@ -263,13 +263,15 @@ SUBROUTINE VectorHelmholtzSolver( Model,Solver,dt,Transient )
 
   Found = .FALSE.
   IF( ASSOCIATED( Model % Constants ) ) THEN
-    mu0inv = 1.0_dp / GetConstReal( Model % Constants,  'Permeability of Vacuum', Found )
+    IF (ListCheckPresent(Model % Constants, 'Permeability of Vacuum')) &
+        mu0inv = 1.0_dp / GetConstReal( Model % Constants, 'Permeability of Vacuum', Found )
   END IF
   IF(.NOT. Found ) mu0inv = 1.0_dp / ( PI * 4.0d-7 )
 
   Found = .FALSE.
   IF( ASSOCIATED( Model % Constants ) ) THEN
-    eps0 = GetConstReal ( Model % Constants, 'Permittivity of Vacuum', Found )
+    IF (ListCheckPresent(Model % Constants, 'Permittivity of Vacuum')) &
+        eps0 = GetConstReal ( Model % Constants, 'Permittivity of Vacuum', Found )
   END IF
   IF(.NOT. Found ) eps0 = 8.854187817d-12