""" NSSCII - FEM. Input to be parsed through. SI-units to be used: + T in K + L in m + k in W/(mK) + q in W/m^2 - ad Neumann + P in W - ad nodal forces """ # Group number. groupnr = 1 # Length in x- and y-direction. L = 0.01 # Thickness (z-direction). hz = 0.0005 # Thermal conductivity (k=k_xx=k_yy, k_xy = 0.). k = 429. # Factor c for modifying thermal conductivity k for # elements in elements_to_be_modified. c = 10. # Elements to be modified. elements_to_be_modified = [ 41-47, 59-63, 77-79, 95 ] # Boundary conditions. q(y=0) = 2000000. T(y=L) = 293.