diff --git a/examples/antora.yml b/examples/antora.yml
index cc051072..c84f1b4b 100644
--- a/examples/antora.yml
+++ b/examples/antora.yml
@@ -1,8 +1,9 @@
 name: cases
 title: Feel++ Toolbox Cases Guide
-version: '0.106'
+version: '0.106-case-10-opuseads-heatfluid'
 nav:
 - modules/ROOT/nav.adoc
+- modules/poisson/nav.adoc
 - modules/csm/nav.adoc
 - modules/cfd/nav.adoc
 - modules/heat/nav.adoc
diff --git a/examples/modules/heatfluid/assets/images/opus/eads_geometry.png b/examples/modules/heatfluid/assets/images/opus/eads_geometry.png
new file mode 100644
index 00000000..a72f9b0b
Binary files /dev/null and b/examples/modules/heatfluid/assets/images/opus/eads_geometry.png differ
diff --git a/examples/modules/heatfluid/examples/opus/eads.geo b/examples/modules/heatfluid/examples/opus/eads.geo
new file mode 100644
index 00000000..61717d89
--- /dev/null
+++ b/examples/modules/heatfluid/examples/opus/eads.geo
@@ -0,0 +1,107 @@
+h = 0.0001;
+// h     = DefineNumber[ 0.1, Name "h"     ];
+e_pcb = DefineNumber[ 2e-3, Name "e_pcb" ]; // 0.3
+e_a   = DefineNumber[ 4e-3, Name "e_a"   ]; // 0.7
+e_ic  = DefineNumber[ 2e-3, Name "e_ic"  ]; // 0.3
+h_1   = DefineNumber[ 2e-2, Name "h_1"   ]; // 0.2
+h_2   = DefineNumber[ 7e-2, Name "h_2"   ]; // 0.6
+h_ic  = DefineNumber[ 2e-2, Name "h_ic"  ]; // 0.2
+h_pcb = DefineNumber[13e-2, Name "h_pcb" ]; // 1
+
+// h_pad = DefineNumber[ 0.04, Name "h_pad" ];
+
+
+Point(1)  = {0, 0, 0, h};
+Point(2)  = {e_pcb, 0, 0, h};
+Point(3)  = {e_pcb + e_a, 0, 0, h};
+Point(4)  = {e_pcb, h_1, 0, h};
+Point(5)  = {e_pcb+e_ic, h_1, 0, h};
+Point(6)  = {e_pcb+e_ic, h_1+h_ic, 0, h};
+Point(7)  = {e_pcb, h_1+h_ic, 0, h};
+Point(8)  = {e_pcb, h_2, 0, h};
+Point(9)  = {e_pcb+e_ic, h_2, 0, h};
+Point(10) = {e_pcb+e_ic, h_2+h_ic, 0, h};
+Point(11) = {e_pcb, h_2+h_ic, 0, h};
+Point(12) = {0, h_pcb, 0, h};
+Point(13) = {e_pcb+e_a, h_pcb, 0, h};
+Point(14) = {e_pcb, h_pcb, 0, h};
+Point(15) = {e_pcb+e_ic, 0, 0, h};
+Point(16) = {e_pcb+e_ic, h_pcb, 0, h};
+
+// Point(17) = {e_pcb, h_pcb + h_pad, 0, h};
+// Point(18) = {e_pcb+e_a, h_pcb + h_pad, 0, h};
+
+
+Line(1)  = {4, 7};
+Line(2)  = {8, 11};
+Line(3)  = {1, 12};
+Line(4)  = {3, 13};
+Line(5)  = {2, 4};
+Line(6)  = {7, 8};
+Line(7)  = {11, 14};
+Line(8)  = {4, 5};
+Line(9)  = {5, 6};
+Line(10) = {6, 7};
+Line(11) = {8, 9};
+Line(12) = {9, 10};
+Line(13) = {10, 11};
+Line(14) = {15, 5};
+Line(15) = {6, 9};
+Line(16) = {10, 16};
+Line(17) = {1, 2};
+Line(18) = {12, 14};
+Line(19) = {2, 15};
+Line(20) = {14, 16};
+Line(21) = {15, 3};
+Line(22) = {16, 13};
+
+// Line(37) = {14, 17};
+// Line(38) = {17, 18};
+// Line(39) = {18, 13};
+
+
+Line Loop(23) = {17, 5, 1, 6, 2, 7, -18, -3};
+Plane Surface(24) = {23};
+Line Loop(25) = {13, -2, 11, 12};
+Plane Surface(26) = {25};
+Line Loop(27) = {10, -1, 8, 9};
+Plane Surface(28) = {27};
+Line Loop(29) = {20, -16, 13, 7};
+Plane Surface(30) = {29};
+Line Loop(31) = {11, -15, 10, 6};
+Plane Surface(32) = {31};
+Line Loop(33) = {8, -14, -19, 5};
+Plane Surface(34) = {33};
+Line Loop(35) = {21, 4, -22, -16, -12, -15, -9, -14};
+Plane Surface(36) = {35};
+
+// Line Loop(40) = {37, 38, 39, -22, -20};
+// Plane Surface(41) = {40};
+
+
+Physical Surface("Air")  = {30,32,34,36}; // 41
+Physical Surface("Pcb")  = {24};
+Physical Surface("IC1")  = {28};
+Physical Surface("IC2")  = {26};
+
+// Physical Line("Fixed")   = {21, 19, 17};
+
+Physical Line("Pcb/Input")  = {17};
+Physical Line("Pcb/Left")   = {3};
+Physical Line("Pcb/Output") = {18};
+
+// Physical Line("Upwards") = {18,20,22};
+
+Physical Line("IC1/Pcb") = {1};
+Physical Line("IC2/Pcb") = {2};
+
+// Physical Line("Air/Pcb") = {3,7,6,5,4};
+// Physical Line("Air/IC1") = {10,9,8};
+// Physical Line("Air/IC2") = {13,12,11};
+
+// Physical Line("Air/PostOutput") = {38};
+Physical Line("Air/Input")         = {21,19};
+Physical Line("Air/Walls")         = {7,6,5,10,9,8,13,12,11}; // 37,39
+Physical Line("Air/Right")         = {4};
+Physical Line("Air/Output")        = {22,20};
+
diff --git a/examples/modules/heatfluid/examples/opus/eads_normal.json b/examples/modules/heatfluid/examples/opus/eads_normal.json
new file mode 100644
index 00000000..e1449697
--- /dev/null
+++ b/examples/modules/heatfluid/examples/opus/eads_normal.json
@@ -0,0 +1,99 @@
+// -*- mode: javascript -*-
+{
+    "Name"      : "Eads",
+    "ShortName" : "Eads",
+
+    "Models" : {
+	"use-model-name" : 1,
+        "fluid" : { "equations" : "Navier-Stokes" }
+    },
+
+    "Materials" : {
+        "Pcb" : {
+	    "markers" : "Pcb",
+            "physics" : "heat",
+            "rho"     : "1.",
+            "k"       : "0.2",
+            "Cp"      : "2e6",
+	    "mu"      : "1e12"
+        },
+	"IC1" : {
+	    "markers" : "IC1",
+	    "physics" : "heat",
+	    "name"    : "IC1",
+	    "rho"     : "1.",
+	    "k"       : "2",
+	    "Cp"      : "1.4e6",
+	    "mu"      : "1e12"
+	},
+	"IC2" : {
+	    "markers" : "IC2",
+	    "physics" : "heat",
+	    "rho"     : "1.",
+	    "k"       : "2",
+	    "Cp"      : "1.4e6",
+	    "mu"      : "1e12"
+	},
+	"Air" : {
+	    "markers" : "Air",
+	    "physics" : ["heat","fluid"],
+	    "rho"     : "1.",
+	    "k"       : "3e-2",
+	    "Cp"      : "1100",
+	    "mu"      : "1.8e-5"
+	}
+    },
+
+    "BoundaryConditions" : {
+	"velocity" : {
+	    "Dirichlet" : {
+		"Air/Right"           : { "expr" : "{0,0}" },
+		"Air/Walls"           : { "expr" : "{0,0}" },
+		"Air/Input"           : { "expr" : "{0,(2e-3+2e-3<x)*(x<2e-3+4e-3)*7e-3*(3/2/(4e-3-2e-3))*(1-((x-(4e-3+2e-3)/2-2e-3)/((4e-3-2e-3)/2))^2)*(1-exp(-t/3))}:x:t" }
+	    }
+	},
+	"fluid" : {
+	    "outlet" : {
+		"Air/Output"          : { "expr" : "0." }
+	    }
+	},
+        "temperature" : {
+            "Dirichlet" : {
+                "Air/Input"           : { "expr" : "300" },
+		"Pcb/Input"           : { "expr" : "300" }
+            },
+            "Neumann" : {
+                "Air/Output"          : { "expr" : "0." },
+                "Pcb/Output"          : { "expr" : "0." },
+		"Air/Right"           : { "expr" : "0." },
+		"Pcb/Left"            : { "expr" : "0." }
+            },
+	    "VolumicForces" : {
+		"IC1"                 : { "expr" : "1e6*(1-exp(-t)):t" },
+		"IC2"                 : { "expr" : "1e6*(1-exp(-t)):t" }
+	    }
+	}
+    },
+
+    "InitialConditions" : {
+	"temperature" : {
+	    "" : {
+		"" : { "expr" : "300" }
+	    }
+	}
+    },
+
+    "PostProcess" : {
+	"use-model-name" : 1,
+	"heat-fluid" : {
+	    "Exports" : {
+		"fields" : [
+		    "fluid.velocity","fluid.pressure",
+		    "heat.temperature","fluid.pid"
+		]
+	    }
+	},
+	"fluid" : {}, "heat" : {}
+    }
+}
+
diff --git a/examples/modules/heatfluid/examples/opus/test1.cfg b/examples/modules/heatfluid/examples/opus/test1.cfg
new file mode 100644
index 00000000..661e5370
--- /dev/null
+++ b/examples/modules/heatfluid/examples/opus/test1.cfg
@@ -0,0 +1,28 @@
+directory=toolboxes/heatfluid/opus/test1
+
+case.dimension=2
+
+[heat-fluid]
+filename=$cfgdir/eads_normal.json
+mesh.filename=$cfgdir/eads.geo
+gmsh.hsize=5e-4
+
+[heat-fluid.heat]
+initial-solution.temperature=300
+stabilization-gls=1
+bdf.order=1
+
+# gmsh.partition=1
+
+[heat-fluid.fluid]
+bdf.order=1
+# solver=Newton
+pc-type=lu
+
+# gmsh.partition=1
+
+[ts]
+time-step=1
+time-final=1500
+# restart.at-last-save=true
+
diff --git a/examples/modules/heatfluid/nav.adoc b/examples/modules/heatfluid/nav.adoc
index edbbf238..9987c34e 100644
--- a/examples/modules/heatfluid/nav.adoc
+++ b/examples/modules/heatfluid/nav.adoc
@@ -1,6 +1,7 @@
 * xref:README.adoc[Heat & Fluid]
 ** xref:2Dbuilding_NS/README.adoc[2D Building with natural convection]
 ** xref:cabin/README.adoc[Aerothermal in an Airplane Cabin]
+** xref:opus/README.adoc[Cooling of electronic components]
 ** xref:TurekHron/README.adoc[Forced heat convection around a cylinder]
 ** xref:naturalconvection/README.adoc[Natural heat convection in a cavity]
 ** xref:pipestokesmass/README.adoc[Mass transport with a Stokes flow in a pipe]
diff --git a/examples/modules/heatfluid/pages/opus/README.adoc b/examples/modules/heatfluid/pages/opus/README.adoc
new file mode 100644
index 00000000..1bb17470
--- /dev/null
+++ b/examples/modules/heatfluid/pages/opus/README.adoc
@@ -0,0 +1,327 @@
+= Example : Electronic component cooling
+:stem: latexmath
+:page-plotly: true
+
+
+This test case has been proposed by mailto:Anabelle.Le-Hyaric@eads.net[Annabelle Le-Hyaric] and
+mailto:michel.fouquembergh@eads.net[Michel Fouquembergh] formerly at AIRBUS.
+
+////
+\label{sec:geometry-description}
+////
+
+We consider a 2D model representative of the neighboring of an electronic component submitted to a
+cooling air flow.  It is described by four geometrical domains in stem:[\mathbb{R}^2] named
+stem:[\Omega_i,i=1,2,3,4], see figure. We suppose the velocity stem:[\vec{v}] is known
+in each domain --- for instance in stem:[\Omega_4] it is the solution of previous Navier-Stokes
+computations.  --- The temperature stem:[T] of the domain stem:[\Omega = \cup_{i=1}^4 \Omega_i] is
+then solution of heat transfer equation :
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:1}
+  \rho C_i \Big( \frac{\partial T}{\partial t} + \vec{v} \cdot \nabla T \Big) 
+      - \nabla \cdot \left( k_i \nabla T \right) = Q_i,\quad i=1,2,3,4
+\end{equation}
+++++
+
+where stem:[t] is the time and in each sub-domain stem:[\Omega_i], stem:[\rho C_i] is
+the volumic thermal capacity, stem:[k_i] is thermal conductivity and stem:[Q_i]
+is a volumic heat dissipated.
+
+image::opus/eads_geometry.png[]
+
+One should notice that the convection term in heat transfer equation
+may lead to spatial oscillations which can be overcome by
+Petrov-Galerkin type or continuous interior penalty stabilization
+techniques.
+
+Integrated circuits (ICs) (domains stem:[\Omega_1] and stem:[\Omega_2] ) are respectively soldered
+on PCB at stem:[\mathbf{x1}=(e_{Pcb}, h_1)] and stem:[\mathbf{x_2}=(e_{Pcb}, h_2)]. They are
+considered as rectangles with width stem:[e_{IC}] and height stem:[h_{IC}]. The printed circuit
+board (PCB) is a rectangle stem:[\Omega_3] of width stem:[e_{PCB}] and height stem:[h_{PCB}]. The
+air(Air) is flowing along the PCB in domain stem:[\Omega_4]. Speed in the air channel
+stem:[\Omega_4] is supposed to have a parabolic profile function of stem:[x] coordinate. Its
+expression is simplified as follows (notice that stem:[\vec{v}] is normally solution of
+Navier-Stokes equations; the resulting temperature and velocity fields should be quite different
+from that simplified model), we have for all stem:[0 \leq y \leq h_{PCB}]
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:2}
+  \begin{array}[c]{rl}
+    e_{Pcb} + e_{Ic} \leq x \leq e_{Pcb} + e_a, & \displaystyle \vec{v} = \frac{3}{2 (e_a-e_{Ic})}\ D\ \Big( 1 - \Big( \frac{x - ( \frac{e_a+e_{Ic}}{2} + e_{Pcb})}{\frac{(e_a-e_{Ic})}{2}} \Big)^2 \Big)\ f(t)\ \vec{y}\\
+    e_{Pcb}  \leq x \leq e_{Pcb} + e_{Ic}, & \vec{v} = 0
+  \end{array}
+\end{equation}
+++++
+
+where stem:[f] is a function of time modelling the starting of the PCB ventilation, _i. e._
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:13}
+  f(t) = 1-\exp(-\frac{t}{3}),
+\end{equation}
+++++
+
+stem:[D] is the air flow rate, see table~\vref{tab:1} and stem:[\vec{y}=(0, 1)^T] is the unit
+vector along the stem:[y] axis. A quick verification shows that
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:18}
+  \int_{\Gamma_4 \cap \Omega_4} \vec{v} \cdot \vec{n} = \int_{\Gamma_4 \cap \Omega_4} v_y = D
+\end{equation}
+++++
+
+The medium velocity stem:[\vec{v}_i = \vec{0}, i=1,2,3] in the solid domains
+stem:[\Omega_i, i=1,2,3].
+
+ICs dissipate heat, we have respectively
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:12}
+  \begin{array}[c]{rl}
+    Q_1 &= Q \big( 1 - \exp(-t) \big)\quad \text{ in } \Omega_1\\
+    Q_2 &= Q \big( 1 - \exp(-t) \big)\quad \text{ in } \Omega_2\\
+  \end{array}
+\end{equation}
+++++
+
+where stem:[Q] is defined in table~\vref{tab:1}.
+
+We shall denote stem:[\vec{n}_{|{\Omega_i}} = \vec{n}_i] denotes the unit outward normal to
+stem:[\Omega_i] and stem:[\vec{n}_{|{\Omega_j}} = \vec{n}_j] denotes the unit outward normal to
+stem:[\Omega_j].
+
+== Boundary conditions
+////
+\label{sec:boundary-conditions}
+////
+
+We set
+
+ - on stem:[\Gamma_3 \cap \Omega_3], a zero flux (Neumann-like) condition
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:10}
+  -k_3\ \nabla T \cdot \vec{n}_3\ =\ 0;
+\end{equation}
+++++
+
+  - on stem:[\Gamma_3 \cap \Omega_4], a zero flux (Robin-like) condition
+
+[stem]
+++++
+  \begin{equation}
+    \label{eq:7}
+    (-k_4\ \nabla T + \rho C_4 T  \mathbf{v} ) \cdot \vec{n}_4\ =\ 0;
+  \end{equation}
+++++
+
+  - on stem:[\Gamma_4, (0 \leq x \leq e_{Pcb}+e_a, y=0)] the temperature is set (Dirichlet condition)
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:11}
+  T\ = T_0;
+\end{equation}
+++++
+
+  - between stem:[\Gamma_1] and stem:[\Gamma_2], periodic conditions
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:9}
+  \begin{array}{rl}
+  T_{\big|\mathbf{x} = 0} &= T_{\big|\mathbf{x} = e_{Pcb} + e_{a}} \\
+  -k_3\ \nabla T \cdot \vec{n}_{3_{\big|\mathbf{x} = 0}} &= k_4\ \nabla T \cdot \vec{n}_4{_{\big|\mathbf{x} = e_{Pcb} + e_{a}}};
+\end{array}
+\end{equation}
+++++
+
+  - at interfaces between the ICs and PCB, there is a thermal contact conductance:
+
+[stem]
+++++
+\begin{equation}
+\label{eq:5}
+\begin{array}{rl}
+  -k_1\ \nabla T \cdot \vec{n}_1 - k_3\ \nabla T \cdot \vec{n}_3 = r_{13}\big( T_{\partial \Omega_1} - T_{\partial \Omega_3} \big)\\
+  -k_2\ \nabla T \cdot \vec{n}_2 - k_3\ \nabla T \cdot \vec{n}_3 = r_{23}\big( T_{\partial \Omega_2} - T_{\partial \Omega_3} \big);\\
+\end{array}
+\end{equation}
+++++
+
+  - on other internal boundaries, the coontinuity of the heat flux and temperature, on
+    stem:[\Gamma_{ij} = \Omega_i \cap \Omega_j \neq \emptyset]
+
+[stem]
+++++
+\begin{equation}
+\label{eq:6}
+\begin{array}{rl}
+  T_i &= T_j \\
+  k_i\ \nabla T \cdot \vec{n}_i &= -k_j\ \nabla T \cdot \vec{n}_j.
+\end{array}
+\end{equation}
+++++
+
+
+== Initial condition
+////
+\label{sec:initial-condition}
+////
+
+At stem:[t=0s], we set stem:[T = T_0].
+
+=== Inputs
+////
+\label{sec:9999s}
+////
+
+The table displays the various fixed and variables parameters of this test-case.
+
+.Table of fixed and variable parameters
+[width="100%"]
+|=======================================================================
+| Name     | Description | Nominal Value | Range              | Units
+2+^| Parameters        3+|
+| stem:[t] | time        |               | stem:[[0, 1500]]   | stem:[s]
+| stem:[Q] | heat source | stem:[10^{6}] | stem:[[0, 10^{6}]] | stem:[W \cdot m^{-3}]
+5+|
+2+^| IC Parameters     3+| 
+| stem:[k_1 = k_2 =  k_{IC}] | thermal conductivity | stem:[2] | stem:[[0. 2,150 ]] | stem:[W \cdot m^{-1} \cdot K^{-1}]
+| stem:[r_{13} = r_{23} = r] | thermal conductance | stem:[100] | stem:[[10^{-1},10^{2} ]] | stem:[W \cdot m^{-2} \cdot K^{-1}]
+| stem:[\rho C_{IC}] | heat capacity | stem:[1. 4 \cdot 10^{6}]| |stem:[J \cdot m^{-3} \cdot K^{-1}]
+| stem:[e_{IC}] | thickness |  stem:[2\cdot 10^{-3}] | |stem:[m]
+| stem:[h_{IC} = L_{IC}] | height | stem:[2\cdot 10^{-2}] | |stem:[m]
+| stem:[h_{1}] | height | stem:[2\cdot 10^{-2}] | |stem:[m]
+| stem:[h_{2}] | height | stem:[7\cdot 10^{-2}] | |stem:[m]
+5+|
+2+^| PCB Parameters 3+|
+| stem:[k_3 =  k_{Pcb}] | thermal conductivity | stem:[0. 2] | | stem:[W \cdot m^{-1} \cdot K^{-1}]
+| stem:[\rho C_{3}] | heat capacity | stem:[2 \cdot 10^{6}]| |stem:[J \cdot m^{-3} \cdot K^{-1}]
+| stem:[e_{Pcb}] | thickness |  stem:[2\cdot 10^{-3}] | |stem:[m]
+| stem:[h_{Pcb}] | height | stem:[13\cdot 10^{-2}] | |stem:[m]
+5+|
+2+^| Air Parameters 3+|
+| stem:[T_0]        | Inflow temperature   | stem:[300]             | | stem:[K] 
+| stem:[D]          | Inflow rate          | stem:[7\cdot 10^{-3}]  | stem:[[5 \cdot 10^{-4} ,10^{-2}]] | stem:[m^2 \cdot s^{-1}]
+| stem:[k_4 ]       | thermal conductivity | stem:[3 \cdot 10^{-2}] | | stem:[W \cdot m^{-1} \cdot K^{-1}]
+| stem:[\rho C_{4}] | heat capacity        | stem:[1100]|           | stem:[J \cdot m^{-3} \cdot K^{-1}]
+| stem:[e_{a}]      | thickness            |  stem:[4\cdot 10^{-3}] | stem:[[2. 5 \cdot 10^{-3} , 5 \cdot 10^{-2}]] | stem:[m]
+|=======================================================================
+
+=== Outputs
+////
+\label{sec:outputs}
+////
+
+The outputs are _(i)_ the mean temperature stem:[s_1(\mu)] of the hottest IC
+
+[stem]
+++++
+\begin{equation}
+  \label{eq:3}
+  s_1(\mu) = \frac{1}{e_{IC} h_{IC}} \int_{\Omega_2} T
+\end{equation}
+++++
+
+and _(ii)_ mean temperature stem:[s_2(\mu)] of the air at the outlet
+
+[stem]
+++++
+  \begin{equation}
+    \label{eq:4}
+    s_2(\mu) = \frac{1}{e_{a}} \int_{\Omega_4 \cap \Gamma_3 } T
+  \end{equation}
+++++
+
+both depends on the solution of \eqref{eq:1} and are dependent on the parameter set stem:[\mu].
+
+We need to monitor stem:[s_1(\mu)] and stem:[s_2(\mu)] because stem:[s_1(\mu)] is the hottest part
+of the model and the IC can't have a temperature above stem:[340K]. stem:[s_2(\mu)] is the outlet
+of the air and in an industrial system we can have others components behind this outlet. So the
+temperature of the air doesn't have to be high to not interfere the proper functioning of these.
+
+= Testcases
+
+We build some example tests to compute for this problem.
+
+== Test 1
+
+This test is the base of all of the other tests, we bold the changes in the next tests. Only these
+parameters varies in these tests :
+
+[width="50%",cols="<2,>3"]
+|==================================================================================
+| Flow rate               | stem:[7 \cdot 10^{-3}]
+| Characteristic length   | stem:[5 \cdot 10^{-4}]
+| Triangle family         | P1
+| Temporal scheme         | BDF 1, stem:[\Delta t = 1s]
+| Stabilisation           | GLS
+|==================================================================================
+
+As the periodicity helps to dissipate the heat in both directions of the PCB, the results
+are attended to be higher than the true results.
+
+++++
+<div id="plotly_test1a"></div>
+<script type="text/javascript">
+  Plotly.d3.csv("https://girder.math.unistra.fr/api/v1/file/5afea86fb0e957402704804a/download",
+    function(err,rows) {
+      var data = [{
+        name: 'IC2 results of test 1',
+        type: 'scatter',
+        x: plotly_unpack(rows,'time'),
+        y: plotly_unpack(rows,'test1-s1'),
+        showlegend: true,
+        line: { color: '#CC3333' }
+      }{
+        name: 'Output results of test 1',
+        type: 'scatter',
+        x: plotly_unpack(rows,'time'),
+        y: plotly_unpack(rows,'test1-s2'),
+        showlegend: true,
+        line: { color: '#3333CC' }
+      }];
+
+      var layout = {
+        title: 'Temperature'
+      };
+      Plotly.plot(plotly_test1a,data,layout,{ showLink: false });
+
+      plotly_add_redimensionable_plot(plotly_test1a);
+    }
+  );
+</script>
+++++
+
+The Feel++ toolboxes are configured with two complementary files, the _.json_ file configure the
+specific materials related to the toolbox for the geometry. The _.cfg_ file configure the generic
+options, for instance the mesh file, the output directory or the solver. The _.json_ file is common
+to all examples, except the stationary problem.
+
+// highlight javascript better
+[source,json] 
+----
+include::{examplesdir}/opus/eads_normal.json[]
+----
+
+// highlight ini better
+[source,cfg]
+----
+include::{examplesdir}/opus/test1.cfg[]
+----