FRCMFD-v2: TEST 1C-N — NARROW-BAND RESONANCE DIAGNOSTIC (0.28–0.32 v, +0.35)
""" FRCMFD-v2: TEST 1C-N — NARROW-BAND RESONANCE DIAGNOSTIC (0.28–0.32 v, +0.35) Focus: Clean COM, dv/dt, canonical momentum P_z, and momentum drift With periodic z-boundaries, consistent operators, and refined timestep. """ import numpy as np import scipy.sparse as sp from datetime import datetime import json import glob import os from scipy.signal import savgol_filter # ================================================================ # LOAD SOLITON # ================================================================ preferred = "/content/test_0A_soliton_20260523_195958.npz" if os.path.exists(preferred): soliton_file = preferred else: soliton_files = glob.glob("test_0A_*_soliton.npz") if not soliton_files: soliton_files = glob.glob("/content/drive/MyDrive/FRCMFD_v2_Backups/*/test_0A_*_soliton.npz") soliton_file = sorted(soliton_files)[-1] data = np.load(soliton_file) Psi_soliton = data["Psi_soli...