Rainfall Distribution

import matplotlib.pyplot as plt
import numpy as np

fig, ax = plt.subplots(figsize=(10, 6), facecolor='#0a0a0f')
ax.set_facecolor('#0a0a0f')

np.random.seed(42)
rainfall = np.random.gamma(2, 15, 365)

ax.ecdf(rainfall, color='#38bdf8', linewidth=2.5)

dry_day = 1
heavy = 50
ax.axvline(x=dry_day, color='#fbbf24', linestyle='--', linewidth=2, label=f'Dry day (<{dry_day}mm)')
ax.axvline(x=heavy, color='#ef4444', linestyle='--', linewidth=2, label=f'Heavy rain (>{heavy}mm)')

ax.set_xlabel('Daily Rainfall (mm)', fontsize=12, color='white')
ax.set_ylabel('Cumulative Probability', fontsize=12, color='white')
ax.set_title('Daily Rainfall Distribution', fontsize=16, fontweight='bold', color='white', pad=20)
ax.tick_params(colors='white')
ax.legend(facecolor='#1a1a2f', edgecolor='#333', labelcolor='white')
ax.grid(True, alpha=0.2, color='white')

for spine in ax.spines.values():
    spine.set_color('#333')

plt.tight_layout()