diff --git a/Magnetometer/Jacobian_Lib.pro b/Magnetometer/Jacobian_Lib.pro
index 4cd3cb13cb10190f2343d23b2446c1de72fcc239..3880822ab00d4556110af3c44423ef73d5f0e67c 100644
--- a/Magnetometer/Jacobian_Lib.pro
+++ b/Magnetometer/Jacobian_Lib.pro
@@ -42,6 +42,10 @@ Jacobian {
            { Region All ; Jacobian Vol ; }
     }
   }
+  { Name Sur ;
+    Case { { Region All ; Jacobian Sur ; }
+    }
+  }
 }
 
 
diff --git a/Magnetometer/Thermal.pro b/Magnetometer/Thermal.pro
index d22aa14049588815a2297d1d230d75491688c010..4bfab1a63d4028f5610b230cd0dcb560d7bd32ce 100644
--- a/Magnetometer/Thermal.pro
+++ b/Magnetometer/Thermal.pro
@@ -10,6 +10,14 @@ FunctionSpace {
   }
 }
 
+Group {
+  DefineGroup [ SurfaceConv_The, SurfaceRad_The ];
+}
+
+Function {
+  DefineFunction [ h, T0, rad ];
+}
+
 Formulation {
   { Name Thermal_T ; Type FemEquation ;
     Quantity {
@@ -26,18 +34,16 @@ Formulation {
 
       Galerkin { [ - /* 0.5 * */ 1/sigma[] * SquNorm[sigma[]*{d v}] , {T} ] ;
                  In Domain_The ; Integration GradGrad ; Jacobian Vol ; }
-      /*
+
       Galerkin { [ h[] * Dof{T} , {T} ] ; // Convection boundary condition
                  In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; }
-      Galerkin { [ -h[] * TemperatureConv[] , {T} ] ;
-                 In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; }
-      */
-      /*
-      Galerkin { [ rad[] * {T}^4 , {T} ] ; // Convection boundary condition
+      Galerkin { [ -h[] * T0[] , {T} ] ;
                  In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; }
-      Galerkin { [ -rad[] * TemperatureConv[]^4 , {T} ] ;
-                 In SurfaceConv_The ; Integration GradGrad ; Jacobian Sur ; }
-      */
+
+      Galerkin { [ rad[] * {T}^4 , {T} ] ; // Radiation boundary condition
+                 In SurfaceRad_The ; Integration GradGrad ; Jacobian Sur ; }
+      Galerkin { [ -rad[] * T0[]^4 , {T} ] ;
+                 In SurfaceRad_The ; Integration GradGrad ; Jacobian Sur ; }
     }
   }
 }