30 
\label{example4} 
\label{example4} 
31 
\sslist{example04a.py} 
\sslist{example04a.py} 
32 


33 
We modify the example in Chapter~\ref{CHAP HEAT 2a} in two ways. Firstly, we look at the steady state 
We modify the example in Chapter~\ref{CHAP HEAT 2a} in two ways: we look the steady state 
34 
case with slightly modified boundary conditions and then we use a more flexible tool 
case with slightly modified boundary conditions and use a more flexible tool 
35 
to generate the geometry. Lets look at the geometry first. 
to generate to generate the geometry. Lets look at the geometry first. 
36 


37 
We want to define a rectangular domain of width $5 km$ and depth $6 km$ below the surface of the Earth. The domain is subject to a few conditions. The temperature is known at the surface and the basement has a known heat flux. Each side of the domain is insulated and the aim is to calculate the final temperature distribution. 
We want to define a rectangular domain of width $5 km$ and depth $6 km$ below the surface of the Earth. The domain is subject to a few conditions. The temperature is known at the surface and the basement has a known heat flux. Each side of the domain is insulated and the aim is to calculate the final temperature distribution. 
38 


238 


239 
It is now possible to start defining our domain and boundaries. 
It is now possible to start defining our domain and boundaries. 
240 


241 
The curve defining our clinal structure is approximately located in the middle of the domain and has a sinusoidal shape. We define the curve by generating points at discrete intervals; $51$ in this case, and then create a smooth curve through the points using the \verb Spline() function. 
The curve defining our clinal structure is located approximately in the middle of the domain and has a sinusoidal shape. We define the curve by generating points at discrete intervals; $51$ in this case, and then create a smooth curve through the points using the \verb Spline() function. 
242 
\begin{python} 
\begin{python} 
243 
# Material Boundary 
# Material Boundary 
244 
x=[ Point(i*dsp\ 
x=[ Point(i*dsp\ 