Fusion Reactor Plasma

import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.colors import LinearSegmentedColormap

fig = plt.figure(figsize=(12, 8), facecolor='#ffffff')
ax = fig.add_subplot(111, projection='3d')
ax.set_facecolor('#ffffff')

sectors = 6
modes = 4
xpos = np.arange(sectors)
zpos = np.arange(modes)
xpos, zpos = np.meshgrid(xpos, zpos)
xpos = xpos.flatten()
zpos = zpos.flatten()
ypos = np.zeros_like(xpos)

dx = 0.55
dz = 0.6
np.random.seed(666)
dy = np.random.exponential(50, size=24) + 100

cmap = LinearSegmentedColormap.from_list('neon', ['#276CF5', '#27D3F5', '#F5D327', '#F54927', '#C82909'])
norm = plt.Normalize(min(dy), max(dy))
colors = [cmap(norm(v)) for v in dy]

ax.bar3d(xpos, ypos, zpos, dx, dy, dz, color=colors, alpha=0.9, edgecolor='#e5e7eb', linewidth=0.4)

ax.set_xlabel('Sector', fontsize=11, color='#1f2937', labelpad=12)
ax.set_ylabel('Temp (M °C)', fontsize=11, color='#1f2937', labelpad=12)
ax.set_zlabel('Mode', fontsize=11, color='#1f2937', labelpad=10)
ax.set_title('Tokamak Plasma Temperature Profile', fontsize=14, color='#1f2937', fontweight='bold', pad=20)

ax.set_xticks(range(6))
ax.set_xticklabels(['S1', 'S2', 'S3', 'S4', 'S5', 'S6'], fontsize=8, color='#374151')
ax.set_zticks(range(4))
ax.set_zticklabels(['L-mode', 'H-mode', 'I-mode', 'QH-mode'], fontsize=7, color='#374151')
ax.tick_params(colors='#374151', labelsize=9)

ax.xaxis.pane.fill = False
ax.yaxis.pane.fill = False
ax.zaxis.pane.fill = False
ax.xaxis.pane.set_edgecolor('#e5e7eb')
ax.yaxis.pane.set_edgecolor('#e5e7eb')
ax.zaxis.pane.set_edgecolor('#e5e7eb')
ax.grid(True, alpha=0.3, linewidth=0.5)

ax.view_init(elev=25, azim=50)
plt.tight_layout()
plt.show()