diff --git a/mgrid_2/getmatrix.cpp b/mgrid_2/getmatrix.cpp
index 741b5e5..45dd276 100644
--- a/mgrid_2/getmatrix.cpp
+++ b/mgrid_2/getmatrix.cpp
@@ -438,22 +438,22 @@ double FEM_Matrix::Thermal_coefficient(const int subdomain)
     {
         // outside
         case 0:
-            lambda = 1.0;
+            lambda = 0.026;
             break;
 
         // ceramic
         case 1:
-            lambda = 1.0;
+            lambda = 3.0;           // anything from 1 to 4
             break;
 
         // water
         case 2:
-            lambda = 1.0;
+            lambda = 0.6;
             break;
 
         // air
         case 3:
-            lambda = 1.0;
+            lambda = 0.026;         // depends on temperature actually
             break;
 
         default: