GRIM_GUI.m

% ************************************************************
%    GRIM(TM) - Generic Restoration Intelligent Milestone
%    Version 3.0.0 (2026)
%
%    A tool of the GriRes(TM) Suite
%    (Grid Resilience Suite)
%
% Copyright (c) 2017-2026 Francis Yunhe Hou
% SPDX-License-Identifier: GPL-3.0-or-later
%
% Author:   Dr. Francis Yunhe Hou
% Contact:  yunhe.hou@gmail.com
% Website:  https://github.com/GriRes/GRIM
%
% ============================================================
% LICENSE
% ============================================================
% This program is free software: you can redistribute it and/or
% modify it under the terms of the GNU General Public License as
% published by the Free Software Foundation, either version 3 of
% the License, or (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see <https://www.gnu.org/licenses/>.
%
% DUAL LICENSING:
%   This software is also available under a commercial license for
%   closed-source or proprietary use where GPL terms cannot be met.
%   Contact: yunhe.hou@gmail.com
%
% ============================================================
% AI-ASSISTED DEVELOPMENT NOTICE
% ============================================================
% Portions of this software were developed with the assistance of
% AI-based programming tools (including OpenAI ChatGPT, Anthropic
% Claude, and xAI Grok). All AI-generated code has been reviewed,
% tested, and validated by the author. The author assumes full
% responsibility for the correctness and quality of the final code.
%
% ============================================================
% HERITAGE & EVOLUTION
% ============================================================
% GRIM v3.0 is the open-source continuation and significant
% evolution of prior work on power system restoration decision
% support tools. It builds directly upon the Generic Restoration
% Milestone (GRM) methodology, originally developed and published in:
%
%   Y. Hou et al., "Computation of Milestones for Decision Support
%   during System Restoration," IEEE Trans. Power Syst., vol. 26,
%   no. 3, pp. 1399-1409, Aug. 2011.
%
% Previous implementations include:
%   HITS (Healing Touch for Power Systems) - open-source MATLAB
%   toolbox, versions 1.0 and 2.0, developed by Prof. Zhijun Qin,
%   released under the permissive 3-clause BSD license.
%
% ============================================================
% CITATION
% ============================================================
% If you use this software, please cite BOTH:
%   [1] Y. Hou, C.-C. Liu, K. Sun, P. Zhang, S. Liu, and
%       D. Mizumura, "Computation of Milestones for Decision
%       Support during System Restoration," IEEE Trans. Power
%       Syst., vol. 26, pp. 1399-1409, 2011.
%   [2] K. Sun, Y. Hou, W. Sun, and J. Qi, Power System Control
%       Under Cascading Failures, John Wiley & Sons, 2019.
%
% GRIM(TM) and GriRes(TM) are trademarks of Francis Yunhe Hou.
% ************************************************************


function GRIM_GUI()
% GRIM_GUI  Unified control centre for the GRIM restoration toolbox
%
% Copyright (c) 2017-2026 Dr. Francis Yunhe Hou
% SPDX-License-Identifier: GPL-3.0-or-later

    %% ===== Locate GRIM root =====
    grimRoot = fileparts(mfilename('fullpath'));
    if isempty(grimRoot), grimRoot = pwd; end
    gridDataDir = fullfile(grimRoot, 'GridData');
    monDir      = fullfile(grimRoot, 'Monitoring');

    %% ===== Application state =====
    S.systems    = {};
    S.sysIdx     = 1;
    S.sysName    = '';
    S.running    = false;
    S.Grid       = [];
    S.hasResults = false;
    S.modes      = {'backtrack','greedy'};
    S.modeNames  = {'Backtrack  (tries alternatives)', ...
                    'Greedy  (fast single-path)'};
    S.hCon       = [];
    S.animating  = false;
    S.animIdx    = 0;
    S.solveInfo  = '';
    S.childFigs  = [];
    S.solveTimer = [];
    S.progFig    = [];            % uifigure for progress dialog
    S.solveStartTic = [];
    addpath(fullfile(grimRoot,'Config'));       % ensure GRIM_Config is on path
    S.config     = GRIM_Config();  % solver configuration (passed as override)
    scanSystems();

    %% ===== Colour palette =====
    C.bg     = [0.067 0.075 0.102];
    C.side   = [0.086 0.098 0.133];
    C.panel  = [0.106 0.118 0.157];
    C.card   = [0.133 0.149 0.192];
    C.input  = [0.114 0.129 0.169];
    C.accent = [0.000 0.545 0.647];
    C.green  = [0.118 0.573 0.345];
    C.amber  = [0.839 0.514 0.075];
    C.red    = [0.745 0.235 0.188];
    C.blue   = [0.149 0.447 0.776];
    C.purple = [0.404 0.275 0.675];
    C.txt    = [0.843 0.867 0.910];
    C.dim    = [0.431 0.471 0.537];
    C.menuTxtIdle = [0.65 0.68 0.74];   % brighter idle colour for menu labels
    C.bdr    = [0.176 0.196 0.255];
    C.con    = [0.043 0.051 0.071];
    C.w      = [1 1 1];
    C.mile   = [0.910 0.682 0.137];
    C.load   = [0.850 0.280 0.450];
    resClrs  = {C.accent,C.accent,C.mile,C.mile,C.mile,...
                C.blue,C.blue,C.green,C.load};

    C.btnMain  = [0.235 0.537 0.635];
    C.btnMainT = [1 1 1];
    C.btnGray  = [0.141 0.173 0.220];
    C.btnTxt   = [0.902 0.929 0.965];
    C.btnViewT = [0.749 0.890 0.933];
    C.btnDataT = [0.839 0.635 0.227];
    C.btnExit  = [0.557 0.247 0.271];
    C.btnExitT = [0.750 0.420 0.440];

    C.menuBg   = [0.094 0.106 0.145];
    C.menuHov  = [0.133 0.149 0.192];
    C.menuTxt  = [0.843 0.867 0.910];
    C.menuDrop = [0.114 0.129 0.169];
    C.menuSep  = [0.176 0.196 0.255];

    %% ===== Main window =====
    scr = get(0,'ScreenSize');
    FW = min(1380, scr(3)-40);
    FH = min(840,  scr(4)-80);
    fig = figure( ...
        'Name',['GRIM' char(8482) ' v3.0.0 — Generic Restoration Intelligent Milestone'], ...
        'NumberTitle','off','MenuBar','none','ToolBar','none', ...
        'Color',C.bg,'Resize','off','DockControls','off', ...
        'Position',[(scr(3)-FW)/2 (scr(4)-FH)/2 FW FH], ...
        'CloseRequestFcn',@onCloseFig);

    %% ===== Custom menu bar (modern flat style with hover) =====
    MENU_H  = 28;
    menuY   = FH - MENU_H;
    uipanel(fig,'Units','pixels','Position',[0 menuY FW MENU_H], ...
        'BackgroundColor',C.menuBg,'BorderType','none');
    % Thin bottom accent line
    uipanel(fig,'Units','pixels','Position',[0 menuY FW 1], ...
        'BackgroundColor',[0.20 0.22 0.28],'BorderType','none');

    % Menu definitions with icons (char(0) = no icon placeholder)
    menuDefs = { ...
      'File', 50, {
        [char(9881) '  Generate Data...'],       @(~,~)dataGenDlg()
        '---',                                   []
        [char(9998) '  Edit System Data...'],     @onDataEdit
        '---',                                   []
        [char(10005) '  Exit                            Ctrl+Q'],  @onExitBtn
      }
      'Settings', 68, {
        [char(9881) '  Solver Settings...'],             @onSettings
        '---',                                           []
        [char(9851) '  Clear Monitoring Data...'],       @onClearMonitoring
      }
      'Analysis', 74, {
        [char(9654) '  Restoration Strategy          Ctrl+R'], @onRun
        '---',                                   []
        [char(9733) '  Power Trajectory'],       @onPowerPlot
        [char(9632) '  Step Viewer'],             @onStepViewer
        [char(9788) '  HTML Viewer...'],          @onGenViewer
        '---',                                   []
        [char(9776) '  View Report'],             @onOpenReport
      }
      'Help', 46, {
        [char(9432) '  User Guide'],             @(~,~)userGuideDlg()
        '---',                                   []
        [char(9679) '  About GRIM'],              @(~,~)aboutDlg()
        [char(10070) '  Citation Info'],          @(~,~)citeDlg()
      }
    };

    S.menuPanels  = {};
    S.menuBtns    = {};
    S.menuItemH   = {};     % cell of arrays: item handles per menu
    S.menuItemCB  = {};     % cell of arrays: callbacks per menu
    mx = 8;
    ITEM_H  = 26;
    ITEM_FS = 10;
    TOP_FS  = 10;
    SEP_H   = 5;

    for mi = 1:size(menuDefs,1)
        label  = menuDefs{mi,1};
        mw     = menuDefs{mi,2};
        items  = menuDefs{mi,3};

        % Top-level button: minimal text, no visual button appearance
        hBtn = uicontrol(fig,'Style','text','String',label, ...
            'Position',[mx menuY+1 mw MENU_H-2], ...
            'FontSize',TOP_FS,'FontName','Segoe UI','FontWeight','normal', ...
            'BackgroundColor',C.menuBg,'ForegroundColor',C.menuTxtIdle, ...
            'HorizontalAlignment','center','Enable','inactive', ...
            'ButtonDownFcn',@(src,evt)toggleMenu(mi));

        nItems  = size(items, 1);
        nSep    = sum(cellfun(@(x)strcmp(x,'---'), items(:,1)));
        nReal   = nItems - nSep;
        dropH   = nReal * ITEM_H + nSep * SEP_H + 8;
        dropW   = 310;
        dropX   = max(mx, 2);
        dropY   = menuY - dropH;

        hDrop = uipanel(fig,'Units','pixels', ...
            'Position',[dropX dropY dropW dropH], ...
            'BackgroundColor',C.menuDrop,'BorderType','line', ...
            'HighlightColor',[0.22 0.24 0.30],'Visible','off');

        itemH  = [];
        itemCB = {};
        iy = dropH - 4;
        for ii = 1:nItems
            iLabel = items{ii,1};
            iCB    = items{ii,2};
            if strcmp(iLabel,'---')
                iy = iy - SEP_H;
                uipanel(hDrop,'Units','pixels', ...
                    'Position',[12 iy+SEP_H/2-1 dropW-24 1], ...
                    'BackgroundColor',[0.20 0.22 0.28],'BorderType','none');
            else
                iy = iy - ITEM_H;
                hItem = uicontrol(hDrop,'Style','text', ...
                    'String',['  ' iLabel], ...
                    'Position',[4 iy dropW-8 ITEM_H], ...
                    'FontSize',ITEM_FS,'FontName','Segoe UI', ...
                    'BackgroundColor',C.menuDrop,'ForegroundColor',C.menuTxt, ...
                    'HorizontalAlignment','left','Enable','inactive', ...
                    'ButtonDownFcn',@(src,evt)menuItemClick(mi, iCB, src, evt));
                itemH(end+1)  = hItem; %#ok<AGROW>
                itemCB{end+1} = iCB;   %#ok<AGROW>
            end
        end

        S.menuBtns{mi}    = hBtn;
        S.menuPanels{mi}  = hDrop;
        S.menuItemH{mi}   = itemH;
        S.menuItemCB{mi}  = itemCB;
        mx = mx + mw + 2;
    end

    set(fig,'WindowButtonDownFcn',@onFigClick);
    set(fig,'WindowButtonMotionFcn',@onFigMotion);

    S.lastHoverItem = [];   % track last highlighted item to reduce flicker

    function toggleMenu(idx)
        vis = get(S.menuPanels{idx},'Visible');
        closeAllMenus();
        if strcmp(vis,'off')
            uistack(S.menuPanels{idx},'top');
            set(S.menuPanels{idx},'Visible','on');
            set(S.menuBtns{idx},'BackgroundColor',C.menuHov, ...
                'ForegroundColor',C.menuTxt,'FontWeight','bold');
        end
    end

    function menuItemClick(menuIdx, cb, src, evt) %#ok<INUSL>
        closeAllMenus();
        if ~isempty(cb), cb(src, evt); end
    end

    function closeAllMenus()
        for ci2 = 1:numel(S.menuPanels)
            set(S.menuPanels{ci2},'Visible','off');
            set(S.menuBtns{ci2},'BackgroundColor',C.menuBg, ...
                'ForegroundColor',C.menuTxtIdle,'FontWeight','normal');
        end
        % Reset hover state
        if ~isempty(S.lastHoverItem) && ishandle(S.lastHoverItem)
            set(S.lastHoverItem,'BackgroundColor',C.menuDrop, ...
                'ForegroundColor',C.menuTxt);
        end
        S.lastHoverItem = [];
    end

    function onFigMotion(~,~)
        % Hover highlighting for open dropdown items
        if ~ishandle(fig); return; end
        try
        cp = get(fig,'CurrentPoint');
        newHover = [];
        for cm = 1:numel(S.menuPanels)
            if ~ishandle(S.menuPanels{cm}); continue; end
            if strcmp(get(S.menuPanels{cm},'Visible'),'on')
                ppos = get(S.menuPanels{cm},'Position');
                for ci = 1:numel(S.menuItemH{cm})
                    h = S.menuItemH{cm}(ci);
                    if ~ishandle(h); continue; end
                    ip = get(h,'Position');
                    % Convert item coords to figure coords
                    ax = ppos(1) + ip(1);
                    ay = ppos(2) + ip(2);
                    if cp(1)>=ax && cp(1)<=ax+ip(3) && ...
                       cp(2)>=ay && cp(2)<=ay+ip(4)
                        newHover = h;
                        break;
                    end
                end
                break;  % only one menu open at a time
            end
        end
        % Update highlight
        changed = false;
        if isempty(newHover) && isempty(S.lastHoverItem)
            changed = false;
        elseif isempty(newHover) || isempty(S.lastHoverItem)
            changed = true;
        elseif newHover ~= S.lastHoverItem
            changed = true;
        end
        if changed
            if ~isempty(S.lastHoverItem) && ishandle(S.lastHoverItem)
                set(S.lastHoverItem,'BackgroundColor',C.menuDrop, ...
                    'ForegroundColor',C.menuTxt);
            end
            if ~isempty(newHover) && ishandle(newHover)
                set(newHover,'BackgroundColor',[0.18 0.20 0.26], ...
                    'ForegroundColor',[1 1 1]);
            end
            S.lastHoverItem = newHover;
        end
        catch
            % Silently ignore stale-handle errors during hover
        end
    end

    function onFigClick(~,~)
        cp = get(fig,'CurrentPoint');
        for ci3 = 1:numel(S.menuPanels)
            if strcmp(get(S.menuPanels{ci3},'Visible'),'on')
                pos = get(S.menuPanels{ci3},'Position');
                if cp(1)>=pos(1) && cp(1)<=pos(1)+pos(3) && ...
                   cp(2)>=pos(2) && cp(2)<=pos(2)+pos(4), return; end
                bpos = get(S.menuBtns{ci3},'Position');
                if cp(1)>=bpos(1) && cp(1)<=bpos(1)+bpos(3) && ...
                   cp(2)>=bpos(2) && cp(2)<=bpos(2)+bpos(4), return; end
            end
        end
        closeAllMenus();
    end

    %% ===== Keyboard shortcuts =====
    set(fig,'KeyPressFcn',@onKeyPress);
    function onKeyPress(~,evt)
        if ismember('control', evt.Modifier)
            switch lower(evt.Key)
                case 'q', onExitBtn([],[]);
                case 'r', onRun([],[]);
            end
        end
    end

    %% ===== Header =====
    HH = 72;  SBH = 32;  HDR_TOP = menuY;

    uipanel(fig,'Units','pixels','Position',[0 HDR_TOP-3 FW 3], ...
        'BackgroundColor',C.accent,'BorderType','none');
    uipanel(fig,'Units','pixels','Position',[20 HDR_TOP-HH-2 FW-40 1], ...
        'BackgroundColor',C.bdr,'BorderType','none');

    %% ===== Status bar =====
    uipanel(fig,'Units','pixels','Position',[0 0 FW SBH], ...
        'BackgroundColor',C.side,'BorderType','none');
    uipanel(fig,'Units','pixels','Position',[0 SBH FW 1], ...
        'BackgroundColor',C.bdr,'BorderType','none');
    S.hStat = uicontrol(fig,'Style','text','String','Ready', ...
        'Position',[14 5 400 20],'FontSize',9,'FontName','Segoe UI', ...
        'ForegroundColor',C.green,'BackgroundColor',C.side,'HorizontalAlignment','left');
    if exist('loadcase','file'), mpT='MATPOWER: OK'; mpC=C.green;
    else, mpT='MATPOWER: not found'; mpC=C.amber; end
    uicontrol(fig,'Style','text','String',mpT, ...
        'Position',[FW-200 5 180 20],'FontSize',9,'FontName','Segoe UI', ...
        'ForegroundColor',mpC,'BackgroundColor',C.side,'HorizontalAlignment','right');
    S.hNSys = uicontrol(fig,'Style','text', ...
        'String',sprintf('Systems: %d',numel(S.systems)), ...
        'Position',[FW-420 5 200 20],'FontSize',9,'FontName','Segoe UI', ...
        'ForegroundColor',C.dim,'BackgroundColor',C.side,'HorizontalAlignment','right');

    %% ===== Sidebar =====
    cTop = HDR_TOP-HH-6;  cBot = SBH+6;  cH = cTop-cBot;
    SX = 16;  SW = 264;
    uipanel(fig,'Units','pixels','Position',[SX cBot SW cH], ...
        'BackgroundColor',C.side,'BorderType','line','HighlightColor',C.bdr);
    bw = SW-32;  yy = cBot+cH-36;

    sideHdr(SX+16,yy,'SYSTEM');
    yy = yy-28;
    S.hSysPop = uicontrol(fig,'Style','popupmenu', ...
        'String',sysPopList(),'Value',max(1,S.sysIdx), ...
        'Position',[SX+16 yy bw 26],'FontSize',10,'FontName','Segoe UI', ...
        'BackgroundColor',C.input,'ForegroundColor',C.txt,'Callback',@onSysChange);
    yy = yy-32;
    uicontrol(fig,'Style','pushbutton','String','Refresh List', ...
        'Position',[SX+16 yy bw 28],'FontSize',9,'FontName','Segoe UI', ...
        'BackgroundColor',C.btnGray,'ForegroundColor',C.btnTxt,'Callback',@onRefresh);

    yy = yy-38;
    sideHdr(SX+16,yy,'SOLVER MODE');
    yy = yy-28;
    S.hModePop = uicontrol(fig,'Style','popupmenu', ...
        'String',S.modeNames,'Value',1, ...
        'Position',[SX+16 yy bw 26],'FontSize',9,'FontName','Segoe UI', ...
        'BackgroundColor',C.input,'ForegroundColor',C.txt);

    yy = yy-34;
    sideHdr(SX+16,yy,'CONFIG PROFILE');
    yy = yy-28;
    S.hCfgPop = uicontrol(fig,'Style','popupmenu', ...
        'String',{'Default'},'Value',1, ...
        'Position',[SX+16 yy bw 26],'FontSize',9,'FontName','Segoe UI', ...
        'BackgroundColor',C.input,'ForegroundColor',C.txt, ...
        'Callback',@onCfgProfileChange);
    refreshCfgProfiles();

    yy = yy-38;
    sideHdr(SX+16,yy,'ANALYSIS');
    btnsA = {
        [char(9654) '  Generate Restoration Strategy'], C.btnMain, @onRun,        44, C.btnMainT
        'Power Trajectory',  C.btnGray, @onPowerPlot,  34, C.btnViewT
        'Step Viewer',       C.btnGray, @onStepViewer, 34, C.btnViewT
        'HTML Viewer',       C.btnGray, @onGenViewer,  34, C.btnViewT
        'View Report',       C.btnGray, @onOpenReport, 34, C.btnViewT
    };
    S.hBtnSolve = gobjects(1);
    for k = 1:size(btnsA,1)
        bh = btnsA{k,4};  yy = yy-bh-6;
        hb = uicontrol(fig,'Style','pushbutton','String',btnsA{k,1}, ...
            'Position',[SX+16 yy bw bh], ...
            'FontSize',10,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',btnsA{k,2},'ForegroundColor',btnsA{k,5}, ...
            'Callback',btnsA{k,3});
        if k==1, S.hBtnSolve = hb; end
    end

    yy = yy-20;
    sideHdr(SX+16,yy,'DATA TOOLS');
    btnsD = {
        'Edit System Data',  C.btnGray, @onDataEdit,        34, C.btnDataT
        'Generate New Data', C.btnGray, @(~,~)dataGenDlg(), 34, C.btnDataT
    };
    for k = 1:size(btnsD,1)
        bh = btnsD{k,4};  yy = yy-bh-6;
        uicontrol(fig,'Style','pushbutton','String',btnsD{k,1}, ...
            'Position',[SX+16 yy bw bh], ...
            'FontSize',10,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',btnsD{k,2},'ForegroundColor',btnsD{k,5}, ...
            'Callback',btnsD{k,3});
    end

    uicontrol(fig,'Style','pushbutton','String',[char(10006) '  Exit GRIM'], ...
        'Position',[SX+16 cBot+30 bw 34], ...
        'FontSize',10,'FontWeight','bold','FontName','Segoe UI', ...
        'BackgroundColor',C.btnExit,'ForegroundColor',C.btnExitT, ...
        'Callback',@onExitBtn);

    uicontrol(fig,'Style','text','String',[char(169) ' 2017-2026 Dr. Francis Yunhe Hou'], ...
        'Position',[SX+16 cBot+8 bw 14],'FontSize',7,'FontName','Segoe UI', ...
        'ForegroundColor',C.dim,'BackgroundColor',C.side,'HorizontalAlignment','center');

    %% ===== Main content area =====
    MX = SX+SW+12;  MW = FW-MX-16;

    sumH = 108;  sumY = cBot+cH-sumH;
    uipanel(fig,'Units','pixels','Position',[MX sumY MW sumH], ...
        'Title','  SYSTEM OVERVIEW  ','FontSize',9,'FontWeight','bold','FontName','Segoe UI', ...
        'TitlePosition','lefttop','ForegroundColor',C.accent,'BackgroundColor',C.panel, ...
        'HighlightColor',C.bdr);
    cardL = {'Buses','Generators','Lines','Transformers','BS Gens','Critical Ld','Disp. Ld'};
    nc = numel(cardL);  cg = 6;  cw = floor((MW-20-(nc-1)*cg)/nc);  ch = sumH-42;
    S.hCard = gobjects(1,nc);
    for k = 1:nc
        cx = MX+10+(k-1)*(cw+cg);  cy = sumY+8;
        uipanel(fig,'Units','pixels','Position',[cx cy cw ch], ...
            'BackgroundColor',C.card,'BorderType','none');
        uicontrol(fig,'Style','text','String',cardL{k}, ...
            'Position',[cx+2 cy+ch-20 cw-4 14],'FontSize',7,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.card,'HorizontalAlignment','center');
        S.hCard(k) = uicontrol(fig,'Style','text','String','--', ...
            'Position',[cx+2 cy+4 cw-4 30],'FontSize',16,'FontWeight','bold','FontName','Consolas', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.card,'HorizontalAlignment','center');
    end

    resH = 100;  resY = sumY-8-resH;
    S.hResPanel = uipanel(fig,'Units','pixels','Position',[MX resY MW resH], ...
        'Title','  RESULTS SUMMARY  ','FontSize',9,'FontWeight','bold','FontName','Segoe UI', ...
        'TitlePosition','lefttop','ForegroundColor',C.mile,'BackgroundColor',C.panel, ...
        'HighlightColor',C.bdr);
    resL = {'Steps','Time(min)','M1: Start','M2: Crank','M3: MinOut', ...
            'Buses On','Gens On','Gen(pu)','Load(pu)'};
    nr = numel(resL);  rg = 5;  rw = floor((MW-20-(nr-1)*rg)/nr);  rh = resH-42;
    S.hRes = gobjects(1,nr);
    for k = 1:nr
        rx = MX+10+(k-1)*(rw+rg);  ry = resY+8;
        uipanel(fig,'Units','pixels','Position',[rx ry rw rh], ...
            'BackgroundColor',C.card,'BorderType','none');
        uicontrol(fig,'Style','text','String',resL{k}, ...
            'Position',[rx+2 ry+rh-18 rw-4 12],'FontSize',7,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.card,'HorizontalAlignment','center');
        S.hRes(k) = uicontrol(fig,'Style','text','String','--', ...
            'Position',[rx+2 ry+2 rw-4 26],'FontSize',14,'FontWeight','bold','FontName','Consolas', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.card,'HorizontalAlignment','center');
    end

    banH = 30;  banY = resY-6-banH;
    uipanel(fig,'Units','pixels','Position',[MX banY MW banH], ...
        'BackgroundColor',C.panel,'BorderType','line','HighlightColor',C.bdr);
    S.hOutcome = uicontrol(fig,'Style','text', ...
        'String','  No results yet — run solver to generate', ...
        'Position',[MX+6 banY+3 MW-12 banH-6], ...
        'FontSize',11,'FontWeight','bold','FontName','Segoe UI', ...
        'ForegroundColor',C.dim,'BackgroundColor',C.panel,'HorizontalAlignment','left');

    conY = cBot;  conH = banY-6-conY;
    uipanel(fig,'Units','pixels','Position',[MX conY MW conH], ...
        'Title','  CONSOLE  ','FontSize',9,'FontWeight','bold','FontName','Segoe UI', ...
        'TitlePosition','lefttop','ForegroundColor',C.accent,'BackgroundColor',C.panel, ...
        'HighlightColor',C.bdr);
    conIX = MX+6;  conIY = conY+6;  conIW = MW-12;  conIH = conH-32;
    S.hCon = uicontrol(fig,'Style','listbox', ...
        'Max',2,'Min',0,'Value',[], ...
        'String',welcomeText(), ...
        'Position',[conIX conIY conIW conIH], ...
        'FontSize',9,'FontName','Consolas', ...
        'ForegroundColor',C.green,'BackgroundColor',C.con,'Enable','on');

    %% ===== Header text with logo =====
    logoFile52 = fullfile(grimRoot,'Config','grim_logo_52.png');
    txtX = 28;  % default if no logo
    if exist(logoFile52,'file')
        logoErr = '';
        logoImg52 = [];
        try logoImg52 = imread(logoFile52);
        catch ME, logoErr = ['imread: ' ME.message]; end
        if isempty(logoErr) && ~isempty(logoImg52)
            try
                imgH = size(logoImg52,1);  imgW = size(logoImg52,2);
                % Position axes to match actual image aspect ratio
                axY = HDR_TOP - HH + round((56-imgH)/2) + 2;
                hLogoAx = axes('Parent',fig,'Units','pixels', ...
                    'Position',[22 axY imgW imgH], ...
                    'Visible','off');
                image('CData',logoImg52,'Parent',hLogoAx);
                set(hLogoAx,'XTick',[],'YTick',[], ...
                    'XLim',[0.5 imgW+0.5],'YLim',[0.5 imgH+0.5], ...
                    'YDir','reverse','Visible','off');
                uistack(hLogoAx,'top');
                txtX = 22 + imgW + 6;
            catch ME, logoErr = ['axes/image: ' ME.message];
            end
        end
        if ~isempty(logoErr)
            fprintf('GRIM_GUI logo: %s\n', logoErr);
        end
    else
        fprintf('GRIM_GUI logo: file not found: %s\n', logoFile52);
    end
    % GRIM (large) + TM (superscript)
    uicontrol(fig,'Style','text','String','GRIM', ...
        'Position',[txtX HDR_TOP-HH+14 150 52], ...
        'FontSize',32,'FontWeight','bold','FontName','Consolas', ...
        'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','left');
    uicontrol(fig,'Style','text','String',char(8482), ...
        'Position',[txtX+122 HDR_TOP-HH+40 30 22], ...
        'FontSize',14,'FontWeight','bold','FontName','Segoe UI', ...
        'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','left');
    uicontrol(fig,'Style','text','String','Generic Restoration Intelligent Milestone', ...
        'Position',[txtX+170 HDR_TOP-HH+28 430 22], ...
        'FontSize',13,'FontWeight','bold','FontName','Segoe UI', ...
        'ForegroundColor',C.txt,'BackgroundColor',C.bg,'HorizontalAlignment','left');
    uicontrol(fig,'Style','text', ...
        'String',['V3.0.0  |  A tool of the GriRes' char(8482) ' Suite (Grid Resilience Suite)'], ...
        'Position',[txtX+170 HDR_TOP-HH+8 500 18], ...
        'FontSize',9,'FontName','Segoe UI', ...
        'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','left');

    %% ===== Initial setup =====
    setupPaths();
    if ~isempty(S.systems)
        S.sysName = S.systems{S.sysIdx};
        loadSystemInfo();
    end

    %% ============================================================
    %%                         CALLBACKS
    %% ============================================================

    % ---- Settings menu callbacks ----
    function refreshCfgProfiles()
        cfgDir2 = fullfile(grimRoot, 'Config');
        items = {'Default'};
        if exist(cfgDir2,'dir')
            mats = dir(fullfile(cfgDir2, '*.mat'));
            for km = 1:length(mats)
                [~, fn] = fileparts(mats(km).name);
                items{end+1} = fn; %#ok<AGROW>
            end
        end
        set(S.hCfgPop, 'String', items, 'Value', 1);
    end

    function onCfgProfileChange(~,~)
        items = get(S.hCfgPop, 'String');
        idx   = get(S.hCfgPop, 'Value');
        sel   = items{idx};
        if strcmp(sel, 'Default')
            S.config = GRIM_Config();
            cLog('Config: reset to Default.');
        else
            cfgDir2 = fullfile(grimRoot, 'Config');
            fp = fullfile(cfgDir2, [sel '.mat']);
            if exist(fp, 'file')
                tmp = load(fp);
                if isfield(tmp, 'grimConfig')
                    S.config = tmp.grimConfig;
                    % Fill missing fields from default
                    def = GRIM_Config();
                    fns = fieldnames(def);
                    for fi = 1:length(fns)
                        if ~isfield(S.config, fns{fi})
                            S.config.(fns{fi}) = def.(fns{fi});
                        end
                    end
                    cLog(sprintf('Config: loaded "%s".', sel));
                else
                    cLog(sprintf('Config: "%s.mat" is not a valid config file.', sel));
                    set(S.hCfgPop, 'Value', 1);
                end
            end
        end
    end

    function onSettings(~,~)
        newCfg = GRIM_Settings(S.config);
        S.config = newCfg;
        refreshCfgProfiles();
        cLog('Settings updated.');
    end

    function onClearMonitoring(~,~)
        if ~exist(monDir,'dir')
            cLog('Monitoring folder does not exist.');
            return;
        end
        ch = questdlg( ...
            sprintf('Delete all files in:\n%s\n\nThis cannot be undone.', monDir), ...
            'Clear Monitoring Data', 'Delete All', 'Cancel', 'Cancel');
        if strcmp(ch, 'Delete All')
            fList = dir(fullfile(monDir, '*'));
            nDel = 0;
            for fi = 1:length(fList)
                if fList(fi).isdir; continue; end
                try delete(fullfile(monDir, fList(fi).name)); nDel = nDel+1;
                catch, end
            end
            cLog(sprintf('Cleared %d files from Monitoring/', nDel));
            setStat(sprintf('Monitoring data cleared (%d files)', nDel), C.green);
        end
    end

    function onExitBtn(~,~)
        if S.running, confirmExit('Solver is still running. Exit anyway?');
        else,         confirmExit('Exit GRIM?'); end
    end

    function confirmExit(msg)
        dW = 360;  dH = 150;
        df = figure('Name','Confirm Exit','NumberTitle','off', ...
            'MenuBar','none','ToolBar','none','Color',C.bg, ...
            'Resize','off','WindowStyle','modal', ...
            'Position',[(scr(3)-dW)/2 (scr(4)-dH)/2 dW dH]);
        uipanel(df,'Units','pixels','Position',[0 dH-3 dW 3], ...
            'BackgroundColor',C.red,'BorderType','none');
        uicontrol(df,'Style','text','String',msg, ...
            'Position',[20 70 dW-40 40], ...
            'FontSize',13,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        uicontrol(df,'Style','pushbutton','String','Exit', ...
            'Position',[dW/2-130 18 120 36], ...
            'FontSize',11,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',C.red,'ForegroundColor',C.w, ...
            'Callback',@(~,~)doExit(df));
        uicontrol(df,'Style','pushbutton','String','Cancel', ...
            'Position',[dW/2+10 18 120 36], ...
            'FontSize',11,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',C.card,'ForegroundColor',C.txt, ...
            'Callback',@(~,~)delete(df));
    end

    function doExit(dlgFig)
        try delete(dlgFig); catch, end
        cleanupAndClose();
    end

    function onCloseFig(~,~)
        if S.running, confirmExit('Solver is still running. Exit anyway?');
        else,         cleanupAndClose(); end
    end

    function cleanupAndClose()
        S.animating = false;
        S.running   = false;
        % Close progress dialog if open
        if ~isempty(S.progFig)
            try delete(S.progFig); catch, end
            S.progFig = [];
        end
        if ~isempty(S.solveTimer)
            try stop(S.solveTimer);   catch, end
            try delete(S.solveTimer); catch, end
            S.solveTimer = [];
        end
        try
            allTimers = timerfindall;
            for ti = 1:numel(allTimers)
                try
                    tName = get(allTimers(ti),'Tag');
                    if contains(tName,'GRIM') || isempty(tName)
                        stop(allTimers(ti)); delete(allTimers(ti));
                    end
                catch
                    try stop(allTimers(ti));  catch, end
                    try delete(allTimers(ti)); catch, end
                end
            end
        catch, end
        S.childFigs(~ishandle(S.childFigs)) = [];
        for ci4 = 1:numel(S.childFigs)
            try delete(S.childFigs(ci4)); catch, end
        end
        try
            allFigs = findall(0,'Type','figure');
            for fi = 1:numel(allFigs)
                nm = get(allFigs(fi),'Name');
                if ~isempty(nm) && startsWith(nm,'GRIM') && allFigs(fi)~=fig
                    try delete(allFigs(fi)); catch, end
                end
            end
        catch, end
        try delete(fig); catch, end
    end

    function onSysChange(~,~)
        v = get(S.hSysPop,'Value');
        if v <= numel(S.systems)
            S.sysIdx = v;  S.sysName = S.systems{v};
            loadSystemInfo();
            cLog(sprintf('System: %s', S.sysName));
        end
    end

    function onRefresh(~,~)
        scanSystems();
        set(S.hSysPop,'String',sysPopList(),'Value',max(1,S.sysIdx));
        set(S.hNSys,'String',sprintf('Systems: %d',numel(S.systems)));
        if ~isempty(S.systems)
            S.sysIdx  = min(S.sysIdx, numel(S.systems));
            S.sysName = S.systems{S.sysIdx};
            loadSystemInfo();
        else
            clearCards(); clearResults();
        end
        cLog('System list refreshed.');
    end

    %% ================================================================
    %%  SOLVE — uifigure + uihtml with JS elapsed timer
    %%  Progress window stays open until Results Summary is populated.
    %% ================================================================

    function onRun(~,~)
        if isempty(S.sysName),    cLog('ERROR: No system selected.'); return; end
        if S.running,             cLog('Solver already running.');     return; end
        if ~exist('GRIM_Run','file'), cLog('ERROR: GRIM_Run.m not found.'); return; end

        sName = S.sysName;
        sMode = S.modes{get(S.hModePop,'Value')};

        S.running = true;
        set(S.hBtnSolve,'Enable','off','String','  Solving...');
        set(fig,'Pointer','watch');
        setStat(sprintf('Solving %s (%s) ...', sName, sMode), C.amber);
        clearResults();
        S.solveInfo = sprintf('%s  (%s)', sName, sMode);
        set(S.hOutcome,'String', ...
            sprintf('  Solving  %s  %s ...   please wait', char(8212), S.solveInfo), ...
            'ForegroundColor', C.amber);
        drawnow;

        sep = repmat(char(9552),1,52);
        cLog(sep);
        cLog(sprintf('  GRIM_Run(''%s'', ''%s'', true, configOverride)', sName, sMode));
        cLog(sep);
        drawnow;

        % Config override as local variable (evalc runs in caller workspace)
        GRIM_GUI_cfgOvr__ = S.config; %#ok<NASGU>
        cmd = sprintf('GRIM_Run(''%s'',''%s'',true,GRIM_GUI_cfgOvr__)', sName, sMode);
        S.solveStartTic = tic;

        % ============================================================
        %  Progress window — uifigure + uihtml
        %  JS elapsed timer runs in Chromium (CEF), independent of
        %  MATLAB main thread blocking in evalc().
        %  Window stays open until Results Summary is fully populated.
        % ============================================================
        try
            scr2 = get(0,'ScreenSize');
            pw = 460;  ph = 240;

            spinHtml = [ ...
                '<html><body style="margin:0;padding:0;' ...
                'background:#111319;height:100vh;' ...
                'font-family:Segoe UI,Arial,sans-serif;overflow:hidden;">' ...
                '<div style="height:3px;width:100%;' ...
                'background:linear-gradient(90deg,#008b96,#1b7a85)"></div>' ...
                '<div style="display:flex;flex-direction:column;' ...
                'align-items:center;justify-content:center;' ...
                'height:calc(100vh - 3px)">' ...
                '<style>' ...
                '@keyframes spin{from{transform:rotate(0deg)}' ...
                'to{transform:rotate(360deg)}}' ...
                '@keyframes pulse{0%,100%{opacity:.35}50%{opacity:1}}' ...
                '</style>' ...
                '<div style="width:48px;height:48px;' ...
                'border:4px solid #1b2135;border-top-color:#008b96;' ...
                'border-radius:50%;animation:spin .8s linear infinite"></div>' ...
                '<div style="color:#d7dde8;font-size:15px;font-weight:bold;' ...
                'margin-top:18px">Generating Restoration Strategy</div>' ...
                '<div style="color:#6e7789;font-size:12px;margin-top:8px">' ...
                'Solving  ' sName '  (' sMode ')</div>' ...
                '<div id="elapsed" style="color:#008b96;font-size:26px;' ...
                'font-weight:bold;font-family:Consolas,monospace;' ...
                'margin-top:16px;letter-spacing:2px">00:00</div>' ...
                '<div style="color:#6e7789;font-size:10px;margin-top:10px;' ...
                'animation:pulse 2s ease-in-out infinite">' ...
                'Please wait ' char(8212) ' this may take a while</div>' ...
                '</div>' ...
                '<script>' ...
                'var t0=Date.now();' ...
                'setInterval(function(){' ...
                'var e=Math.floor((Date.now()-t0)/1000);' ...
                'var h=Math.floor(e/3600);' ...
                'var m=Math.floor((e%3600)/60);' ...
                'var s=e%60;' ...
                'var str;' ...
                'if(h>0){str=(h<10?"0":"")+h+":"+(m<10?"0":"")+m+":"+(s<10?"0":"")+s}' ...
                'else{str=(m<10?"0":"")+m+":"+(s<10?"0":"")+s}' ...
                'document.getElementById("elapsed").textContent=str;' ...
                '},1000);' ...
                '</script>' ...
                '</body></html>'];

            tmpHtml = fullfile(tempdir, 'grim_progress.html');
            fid3 = fopen(tmpHtml,'w','n','UTF-8');
            fwrite(fid3, spinHtml, 'char');
            fclose(fid3);

            S.progFig = uifigure( ...
                'Name','GRIM — Solving', ...
                'Position',[round((scr2(3)-pw)/2) round((scr2(4)-ph)/2) pw ph], ...
                'Color',[0.067 0.075 0.102], ...
                'Resize','off', ...
                'CloseRequestFcn',@(~,~)0);   % block accidental close

            try
                uihtml(S.progFig,'HTMLSource',tmpHtml, ...
                    'Position',[0 0 pw ph]);
            catch
                % pre-R2019b fallback: static label, no animation
                uilabel(S.progFig, ...
                    'Text',sprintf('Solving  %s  (%s) ...\n\nPlease wait', sName, sMode), ...
                    'Position',[30 30 pw-60 ph-60], ...
                    'FontSize',14,'FontColor',[0.843 0.867 0.910], ...
                    'BackgroundColor',[0.067 0.075 0.102], ...
                    'HorizontalAlignment','center','WordWrap','on');
            end

            drawnow;
            pause(0.6);    % let CEF fully render before blocking
        catch
            S.progFig = [];
        end

        runOK       = false;
        output      = '';
        runErr      = [];
        isInterrupt = false;

        try
            output = evalc(cmd);
            runOK  = true;
        catch ME
            runErr = ME;
            if strcmp(ME.identifier,'MATLAB:interrupt') || ...
               contains(ME.message,'terminated','IgnoreCase',true) || ...
               contains(ME.message,'interrupt','IgnoreCase',true)
                isInterrupt = true;
            end
        end

        % Clean up temp config variable from base workspace
        % (config override was local variable, no base cleanup needed)

        % NOTE: progress window is NOT closed here.
        % It stays open until Results Summary is fully populated below.

        % ---- Elapsed time ----
        if ~isempty(S.solveStartTic), solveElapsed = toc(S.solveStartTic);
        else,                         solveElapsed = 0; end
        if solveElapsed < 60
            etStr = sprintf('%.1fs', solveElapsed);
        elseif solveElapsed < 3600
            etStr = sprintf('%dm %02ds', floor(solveElapsed/60), round(mod(solveElapsed,60)));
        else
            etStr = sprintf('%.1fh', solveElapsed/3600);
        end

        % ---- Log solver output ----
        if ~isempty(output)
            outLines = strsplit(output, {char(10), char(13)});
            for li = 1:numel(outLines)
                t = outLines{li};
                if ~isempty(strtrim(t)), cLogRaw(t); end
            end
        end
        cLog(sep);
        cLog(sprintf('  Elapsed: %s', etStr));

        % ---- Process results ----
        if runOK
            resFile = fullfile(monDir, sprintf('GRIM_results_%s.mat', sName));
            if exist(resFile,'file')
                d = load(resFile);
                S.Grid = d.Grid;  S.hasResults = true;
                try populateResults();
                catch ME_pop, cLog(sprintf('WARNING: populateResults failed: %s', ME_pop.message));
                end

                G2 = S.Grid;
                ns = findLastValidStep(G2);

                if ns == 0
                    nGenAll = numel(G2.Generator.ID);
                    nBusAll = numel(G2.Bus.ID);
                    cLog('');
                    cLog([char(9888) '  NO SOLUTION  ' char(9888)]);
                    cLog('  Solver stopped without completing any step.');
                    setStat([char(9888) ' No solution found  (' etStr ')'], C.amber);
                    set(S.hOutcome,'String', ...
                        sprintf('  %s  NO SOLUTION %s  Solver could not complete any step (%d gens, %d buses)  [%s]', ...
                        char(9679), char(8212), nGenAll, nBusAll, etStr), ...
                        'ForegroundColor', C.amber);
                    set(S.hResPanel,'ForegroundColor', C.amber);
                else
                    fn = sprintf('Step_%d', ns);
                    nGenOn  = sum(G2.Generator.IsOn.(fn));
                    nGenAll = numel(G2.Generator.ID);
                    nBusOn  = sum(G2.Bus.IsOn.(fn));
                    nBusAll = numel(G2.Bus.ID);

                    cLog('');
                    if nGenOn==nGenAll && nBusOn==nBusAll
                        cLog([char(9733) char(9733) char(9733) ...
                              '  RESTORATION COMPLETE  ' ...
                              char(9733) char(9733) char(9733)]);
                        cLog(sprintf('  All %d generators online, all %d buses energized.', nGenAll, nBusAll));
                        setStat(sprintf('%s  Complete — %d/%d gens, %d/%d buses  [%s]', ...
                            char(10004), nGenOn, nGenAll, nBusOn, nBusAll, etStr), C.green);
                        setOutcome(nGenOn,nGenAll,nBusOn,nBusAll,'complete');
                    elseif nGenOn==nGenAll
                        cLog([char(9733) '  ALL GENERATORS STARTED  ' char(9733)]);
                        cLog(sprintf('  %d/%d buses energized.', nBusOn,nBusAll));
                        setStat(sprintf('%s  Done — %d/%d gens, %d/%d buses  [%s]', ...
                            char(10004), nGenOn,nGenAll,nBusOn,nBusAll, etStr), C.green);
                        setOutcome(nGenOn,nGenAll,nBusOn,nBusAll,'done');
                    else
                        cLog([char(9888) '  PARTIAL RESTORATION  ' char(9888)]);
                        cLog(sprintf('  %d/%d gens, %d/%d buses.', nGenOn,nGenAll,nBusOn,nBusAll));
                        setStat(sprintf('%s  Partial — %d/%d gens, %d/%d buses  [%s]', ...
                            char(9888), nGenOn,nGenAll,nBusOn,nBusAll, etStr), C.amber);
                        setOutcome(nGenOn,nGenAll,nBusOn,nBusAll,'partial');
                    end
                end
            else
                cLog('WARNING: Results file not found after run.');
                setStat(['Done (no results file)  [' etStr ']'], C.amber);
                set(S.hOutcome,'String', ...
                    sprintf('  %s  WARNING %s  Results file not found  [%s]', char(9679),char(8212),etStr), ...
                    'ForegroundColor',C.amber);
            end

        elseif isInterrupt
            cLog([char(9632) '  Solver interrupted by user (Ctrl+C).']);
            resFile = fullfile(monDir, sprintf('GRIM_results_%s.mat', sName));
            if exist(resFile,'file')
                try
                    d = load(resFile);
                    S.Grid = d.Grid;  S.hasResults = true;
                    try populateResults(); catch, end
                    ns = findLastValidStep(S.Grid);
                    if ns > 0
                        fn = sprintf('Step_%d', ns);
                        nGenOn  = sum(S.Grid.Generator.IsOn.(fn));
                        nGenAll = numel(S.Grid.Generator.ID);
                        nBusOn  = sum(S.Grid.Bus.IsOn.(fn));
                        nBusAll = numel(S.Grid.Bus.ID);
                        cLog(sprintf('  Partial results loaded: step %d, %d/%d gens, %d/%d buses.', ...
                            ns, nGenOn, nGenAll, nBusOn, nBusAll));
                        setStat(sprintf('Interrupted at step %d — %d/%d gens, %d/%d buses  [%s]', ...
                            ns, nGenOn, nGenAll, nBusOn, nBusAll, etStr), C.amber);
                        set(S.hOutcome,'String', ...
                            sprintf('  %s  INTERRUPTED %s  Stopped at step %d  (%d/%d gens, %d/%d buses)  [%s]', ...
                            char(9632), char(8212), ns, nGenOn, nGenAll, nBusOn, nBusAll, etStr), ...
                            'ForegroundColor', C.amber);
                        set(S.hResPanel,'ForegroundColor', C.amber);
                    else
                        cLog('  No completed steps found in partial results.');
                        setStat(['Interrupted — no completed steps  [' etStr ']'], C.amber);
                        set(S.hOutcome,'String', ...
                            sprintf('  %s  INTERRUPTED %s  No completed steps  [%s]', ...
                            char(9632), char(8212), etStr), 'ForegroundColor', C.amber);
                    end
                catch
                    cLog('  Could not load partial results.');
                    setStat(['Interrupted  [' etStr ']'], C.amber);
                    set(S.hOutcome,'String', ...
                        sprintf('  %s  INTERRUPTED %s  Solver stopped by user  [%s]', ...
                        char(9632), char(8212), etStr), 'ForegroundColor', C.amber);
                end
            else
                cLog('  No results file found.');
                setStat(['Interrupted — no results saved  [' etStr ']'], C.amber);
                set(S.hOutcome,'String', ...
                    sprintf('  %s  INTERRUPTED %s  No results saved yet  [%s]', ...
                    char(9632), char(8212), etStr), 'ForegroundColor', C.amber);
            end

        else
            cLog(sprintf('ERROR: %s', runErr.message));
            for si = 1:min(8,numel(runErr.stack))
                cLog(sprintf('  in %s (line %d)',runErr.stack(si).name,runErr.stack(si).line));
            end
            setStat([char(10008) ' Error — see console  [' etStr ']'], C.red);
            set(S.hOutcome,'String', ...
                sprintf('  %s  ERROR %s  %s  [%s]', char(9679),char(8212),runErr.message,etStr), ...
                'ForegroundColor',C.red);
            set(S.hResPanel,'ForegroundColor',C.red);
        end
        cLog(sep); cLog('');
        scrollBottom();

        % ---- Close progress window AFTER results are populated ----
        % Re-capture elapsed time so status bar matches the progress dialog
        if ~isempty(S.solveStartTic)
            solveElapsed = toc(S.solveStartTic);
            if solveElapsed < 60
                etStr = sprintf('%.1fs', solveElapsed);
            elseif solveElapsed < 3600
                etStr = sprintf('%dm %02ds', floor(solveElapsed/60), round(mod(solveElapsed,60)));
            else
                etStr = sprintf('%.1fh', solveElapsed/3600);
            end
        end
        drawnow;                                                  % <-- flush UI so Results Summary is visible
        try delete(S.progFig); catch, end
        S.progFig = [];

        % Refresh status bar with final elapsed time (matches progress dialog)
        try
            curStat = get(S.hStat,'String');
            curStat = regexprep(curStat, '\[[\d.]+[smh][\d ]*s?\]', ['[' etStr ']']);
            set(S.hStat,'String',curStat);
        catch
        end

        S.running = false;
        set(S.hBtnSolve,'Enable','on','String',[char(9654) '  Generate Restoration Strategy']);
        set(fig,'Pointer','arrow');
    end

    %% ---- Data Editor ----
    function onDataEdit(~,~)
        if exist('GRIM_DataEditor','file')
            cLog('Opening Data Editor...');
            try
                GRIM_DataEditor();
                hDE = findall(0,'Type','figure','Name','GRIM Data Editor');
                if ~isempty(hDE), trackChild(hDE(1)); end
                cLog('Data Editor opened.');
            catch ME
                cLog(sprintf('DataEditor error: %s', ME.message));
                for si = 1:min(5,numel(ME.stack))
                    cLog(sprintf('  in %s (line %d)',ME.stack(si).name,ME.stack(si).line));
                end
            end
        else
            cLog('ERROR: GRIM_DataEditor.m not found on path.');
        end
    end

    function onGenViewer(~,~)
        if isempty(S.sysName), cLog('No system selected.'); return; end
        resFile = fullfile(monDir,sprintf('GRIM_results_%s.mat',S.sysName));
        if ~exist(resFile,'file'), cLog('No results. Run solver first.'); return; end
        if ~exist('GRIM_GenViewer','file'), cLog('ERROR: GRIM_GenViewer.m not found.'); return; end
        cLog('Generating HTML Viewer...'); setStat('Generating HTML...',C.amber); drawnow;
        try
            owd = pwd; cd(grimRoot);
            GRIM_GenViewer(S.sysName); cd(owd);
            htmlFile = fullfile(grimRoot,'Monitoring',sprintf('GRIM_Viewer_%s.html',S.sysName));
            if exist(htmlFile,'file')
                web(htmlFile,'-browser');
                cLog(sprintf('HTML Viewer: %s',htmlFile));
            else, cLog('HTML generated (check Monitoring/).');
            end
            setStat('HTML generated',C.green);
        catch ME
            try cd(owd); catch, end
            cLog(sprintf('GenViewer: %s',ME.message)); setStat('Error',C.red);
        end
    end

    function onStepViewer(~,~)
        if isempty(S.sysName), cLog('No system selected.'); return; end
        resFile = fullfile(monDir,sprintf('GRIM_results_%s.mat',S.sysName));
        if ~exist(resFile,'file'), cLog('No results. Run solver first.'); return; end
        if exist('GRIM_StepViewer','file')
            cLog('Opening Step Viewer...');
            try GRIM_StepViewer(S.sysName);
            catch ME, cLog(sprintf('StepViewer: %s',ME.message)); end
        else, cLog('ERROR: GRIM_StepViewer.m not found.');
        end
    end

    function onPowerPlot(~,~)
        if isempty(S.sysName), cLog('No system selected.'); return; end
        resFile = fullfile(monDir,sprintf('GRIM_results_%s.mat',S.sysName));
        if ~exist(resFile,'file'), cLog('No results. Run solver first.'); return; end
        if exist('GRIM_PlotPower','file')
            cLog('Generating power plot...');
            try d = load(resFile); GRIM_PlotPower(d.Grid,S.sysName); cLog('Plot generated.');
            catch ME, cLog(sprintf('PlotPower: %s',ME.message)); end
        else, cLog('ERROR: GRIM_PlotPower.m not found.');
        end
    end

    function onOpenReport(~,~)
        if isempty(S.sysName), cLog('No system selected.'); return; end
        mdFile = fullfile(grimRoot,'Monitoring',sprintf('GRIM_Report_%s.md',S.sysName));
        if ~exist(mdFile,'file'), cLog('Report not found. Run solver first.'); return; end
        cLog(sprintf('Opening report: %s',mdFile));

        % Determine initial showInput state from Config or bus count
        resFile = fullfile(monDir,sprintf('GRIM_results_%s.mat',S.sysName));
        showInput = false;
        try
            busFile = fullfile(gridDataDir,[S.sysName '_bus.txt']);
            nBus2 = cntDataLines(busFile);
            showInput = (nBus2 <= 30);
        catch
        end

        % Generate report with current showInput setting
        rptSysName = S.sysName;
        generateAndShow();

        function generateAndShow()
            % Regenerate markdown report with showInput flag
            if exist(resFile,'file')
                try
                    d2 = load(resFile);
                    busF = fullfile(gridDataDir,[rptSysName '_bus.txt']);
                    nb = cntDataLines(busF);
                    GRIM_ReportExport(d2.Grid, rptSysName, 'backtrack', mdFile, nb, showInput);
                catch
                end
            end

            fid = fopen(mdFile,'r','n','UTF-8');
            if fid < 0, cLog('Cannot read file.'); return; end
            raw = fread(fid,Inf,'*char')'; fclose(fid);
            htmlStr = md2html(raw);

            tmpHtmlFile = fullfile(tempdir, sprintf('GRIM_Report_%s.html', rptSysName));
            fid3 = fopen(tmpHtmlFile,'w','n','UTF-8');
            fwrite(fid3, htmlStr, 'char'); fclose(fid3);

            rW = 860; rH = 700;
            rf = uifigure('Name',['GRIM Report: ' rptSysName], ...
                'Position',[(scr(3)-rW)/2 (scr(4)-rH)/2 rW rH]);
            trackChild(rf);

            hHtml = uihtml(rf, 'HTMLSource', tmpHtmlFile, ...
                'Position', [10 50 rW-20 rH-60]);

            % Toggle button: Show / Hide Input Data
            if showInput
                togLabel = 'Hide Input Data';
                togClr   = C.amber;
            else
                togLabel = 'Show Input Data';
                togClr   = C.green;
            end
            hToggle = uibutton(rf, 'Text', togLabel, ...
                'Position', [rW/2-310 10 140 32], ...
                'FontSize', 10, 'FontWeight', 'bold', ...
                'BackgroundColor', togClr, 'FontColor', C.w);
            % Set callback AFTER hToggle exists
            hToggle.ButtonPushedFcn = @(~,~)toggleInputData(hHtml, hToggle);

            uibutton(rf, 'Text', 'Open in Browser', ...
                'Position', [rW/2-155 10 140 32], ...
                'FontSize', 10, 'FontWeight', 'bold', ...
                'BackgroundColor', C.blue, 'FontColor', C.w, ...
                'ButtonPushedFcn', @(~,~)openInBrowser(htmlStr,rptSysName));
            uibutton(rf, 'Text', 'Open .md in Editor', ...
                'Position', [rW/2+0 10 140 32], ...
                'FontSize', 10, 'FontWeight', 'bold', ...
                'BackgroundColor', C.card, 'FontColor', C.w, ...
                'ButtonPushedFcn', @(~,~)tryEdit(mdFile));
            uibutton(rf, 'Text', 'Close', ...
                'Position', [rW/2+155 10 140 32], ...
                'FontSize', 10, 'FontWeight', 'bold', ...
                'BackgroundColor', C.card, 'FontColor', C.w, ...
                'ButtonPushedFcn', @(~,~)delete(rf));
        end

        function toggleInputData(hHtml2, hTog)
            showInput = ~showInput;
            % Regenerate report
            if exist(resFile,'file')
                try
                    d2 = load(resFile);
                    busF = fullfile(gridDataDir,[rptSysName '_bus.txt']);
                    nb = cntDataLines(busF);
                    GRIM_ReportExport(d2.Grid, rptSysName, 'backtrack', mdFile, nb, showInput);
                catch
                end
            end
            % Re-read and convert
            fid4 = fopen(mdFile,'r','n','UTF-8');
            if fid4 < 0; return; end
            raw2 = fread(fid4,Inf,'*char')'; fclose(fid4);
            html2 = md2html(raw2);
            % Write to UNIQUE temp file (bust uihtml cache)
            tmpFile2 = fullfile(tempdir, sprintf('GRIM_Report_%s_%d.html', rptSysName, round(rem(now,1)*1e8)));
            fid5 = fopen(tmpFile2,'w','n','UTF-8');
            fwrite(fid5, html2, 'char'); fclose(fid5);
            % Refresh viewer
            if isvalid(hHtml2)
                hHtml2.HTMLSource = tmpFile2;
            end
            % Update button
            if showInput
                set(hTog, 'Text', 'Hide Input Data', 'BackgroundColor', C.amber);
            else
                set(hTog, 'Text', 'Show Input Data', 'BackgroundColor', C.green);
            end
        end
    end

    function openInBrowser(htmlStr,sn)
        tmp = fullfile(tempdir,['GRIM_Report_' sn '.html']);
        fid2 = fopen(tmp,'w','n','UTF-8'); fwrite(fid2,htmlStr,'char'); fclose(fid2);
        web(tmp,'-browser');
    end

    function tryEdit(f)
        try edit(f); catch, web(['file:///' strrep(f,'\','/')],'-browser'); end
    end

    %% ============================================================
    %%                        DIALOGS
    %% ============================================================

    function dataGenDlg()
        if ~exist('GRIM_DataGen','file')
            cLog('ERROR: GRIM_DataGen.m not found on path.'); return;
        end
        if ~exist('loadcase','file')
            cLog('WARNING: MATPOWER not on path — Data Generator may not work.');
        end
        cLog('Opening GRIM Data Generator...');
        try
            GRIM_DataGen();
            cLog('Data Generator opened.');
            cLog('Click [Refresh List] after generating files.');
            setStat('Data Generator opened', C.green);
        catch ME
            cLog(sprintf('GRIM_DataGen error: %s', ME.message));
            setStat('Error opening Data Generator', C.red);
        end
    end

    function userGuideDlg()
        gW = 680;  gH = 640;
        gf2 = figure('Name','GRIM User Guide','NumberTitle','off', ...
            'MenuBar','none','ToolBar','none','Color',C.bg,'Resize','on', ...
            'Position',[(scr(3)-gW)/2 (scr(4)-gH)/2 gW gH]);
        trackChild(gf2);
        uipanel(gf2,'Units','pixels','Position',[0 gH-3 gW 3], ...
            'BackgroundColor',C.accent,'BorderType','none');
        uicontrol(gf2,'Style','text','String','GRIM User Guide', ...
            'Position',[20 gH-42 gW-40 28], ...
            'FontSize',15,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','left');
        uipanel(gf2,'Units','pixels','Position',[20 gH-48 gW-40 1], ...
            'BackgroundColor',C.bdr,'BorderType','none');

        ln = repmat(char(9472),1,38);
        guide = {
            ''
            '  GRIM (Generic Restoration Intelligent'
            '  Milestone) — decision-support tool for'
            '  power system restoration after blackout.'
            ''
            ['  ' ln]
            '  QUICK START'
            ['  ' ln]
            ''
            '  1. Select system from dropdown'
            '  2. Choose mode: Backtrack / Greedy'
            '  3. Click [Generate Restoration Strategy]'
            '     Console shows solver output after run.'
            '     Coloured banner shows outcome.'
            '  4. Inspect via viewers / report'
            ''
            '  To stop a long-running solve:'
            '     Press Ctrl+C in the Command Window.'
            '     Partial results will be loaded if available.'
            ''
            ['  ' ln]
            '  ANALYSIS BUTTONS'
            ['  ' ln]
            ''
            '  [Generate Restoration Strategy]'
            '    Runs GRIM_Run(). Output is captured'
            '    and replayed in the console. The'
            '    outcome banner shows:'
            '      Green dot  = full restoration'
            '      Amber dot  = partial restoration'
            '      Red dot    = solver error'
            '      Square     = interrupted by user'
            ''
            '  [Power Trajectory]'
            '    Generation/load plot over all steps.'
            ''
            '  [Step Viewer]'
            '    Step-by-step MATLAB visualisation.'
            ''
            '  [HTML Viewer]'
            '    Generates interactive HTML, opens'
            '    in your default browser.'
            ''
            '  [View Report]'
            '    Opens the Markdown report rendered'
            '    as styled HTML inside MATLAB.'
            ''
            ['  ' ln]
            '  DATA TOOLS'
            ['  ' ln]
            ''
            '  [Edit System Data]'
            '    Spreadsheet editor for buses, gens,'
            '    lines, transformers, and loads.'
            ''
            '  [Generate New Data]'
            '    Opens the GRIM Data Generator GUI.'
            '    Converts MATPOWER cases to GRIM'
            '    data files. Click Refresh List in the'
            '    sidebar after generating.'
            ''
            ['  ' ln]
            '  MILESTONES'
            ['  ' ln]
            ''
            '  M1  Last NBS gen receives crank power'
            '  M2  Last NBS gen finishes cranking'
            '  M3  Last NBS gen reaches min output'
            ''
            ['  ' ln]
            '  EXIT'
            ['  ' ln]
            ''
            '  Click [Exit GRIM] button (sidebar) or'
            '  File > Exit (Ctrl+Q). A confirmation'
            '  dialog will appear. All child windows,'
            '  timers, and resources are cleaned up.'
            ''
            ['  ' ln]
            '  SHORTCUTS'
            ['  ' ln]
            ''
            '  Ctrl+R   Solve restoration'
            '  Ctrl+Q   Quit GRIM GUI'
            '  Ctrl+C   Stop solver (in Command Window)'
            ''
        };
        uicontrol(gf2,'Style','edit','Max',2,'Min',0, ...
            'String',guide,'Position',[14 52 gW-28 gH-110], ...
            'FontSize',10,'FontName','Consolas', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.con, ...
            'HorizontalAlignment','left','Enable','inactive');
        uicontrol(gf2,'Style','pushbutton','String','Close', ...
            'Position',[(gW-140)/2 12 140 32], ...
            'FontSize',11,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',C.accent,'ForegroundColor',C.w, ...
            'Callback',@(~,~)delete(gf2));
    end

    function aboutDlg()
        TM = char(8482);
        aW = 560;  aH = 640;
        af = figure('Name',['About GRIM' TM],'NumberTitle','off', ...
            'MenuBar','none','ToolBar','none','Color',C.bg, ...
            'Resize','off','WindowStyle','modal', ...
            'Position',[(scr(3)-aW)/2 (scr(4)-aH)/2 aW aH]);
        uipanel(af,'Units','pixels','Position',[0 aH-3 aW 3], ...
            'BackgroundColor',C.accent,'BorderType','none');

        pH = aH - 54;

        % ---- Page 1: Main info ----
        pg1 = uipanel(af,'Units','pixels','Position',[0 48 aW pH], ...
            'BackgroundColor',C.bg,'BorderType','none');

        % Logo
        logoFile96 = fullfile(grimRoot,'Config','grim_logo_96.png');
        y1 = pH - 20;
        if exist(logoFile96,'file')
            logoErr2 = '';
            im96 = [];
            try im96 = imread(logoFile96);
            catch ME, logoErr2 = ['imread: ' ME.message]; end
            if isempty(logoErr2) && ~isempty(im96)
                try
                    imgH2 = size(im96,1);  imgW2 = size(im96,2);
                    y1 = y1 - imgH2;
                    hLogoAx2 = axes('Parent',pg1,'Units','pixels', ...
                        'Position',[(aW-imgW2)/2 y1 imgW2 imgH2], ...
                        'Visible','off');
                    image('CData',im96,'Parent',hLogoAx2);
                    set(hLogoAx2,'XTick',[],'YTick',[], ...
                        'XLim',[0.5 imgW2+0.5],'YLim',[0.5 imgH2+0.5], ...
                        'YDir','reverse','Visible','off');
                    y1 = y1 - 8;
                catch ME, logoErr2 = ['axes/image: ' ME.message];
                    y1 = pH - 50;
                end
            else
                y1 = y1 - 30;
            end
            if ~isempty(logoErr2)
                fprintf('About logo: %s\n', logoErr2);
                y1 = pH - 50;
            end
        else
            y1 = y1 - 30;
        end

        % GRIM (centered large) + TM (small superscript)
        y1 = y1 - 65;
        uicontrol(pg1,'Style','text','String','GRIM', ...
            'Position',[0 y1 aW 62], ...
            'FontSize',36,'FontWeight','bold','FontName','Consolas', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        % TM superscript: right of "GRIM" center
        uicontrol(pg1,'Style','text','String',TM, ...
            'Position',[aW/2+76 y1+48 28 18], ...
            'FontSize',12,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','left');
        y1 = y1 - 30;
        uicontrol(pg1,'Style','text', ...
            'String','Generic Restoration Intelligent Milestone', ...
            'Position',[0 y1 aW 28], ...
            'FontSize',13,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 26;
        uicontrol(pg1,'Style','text','String','Version 3.0.0', ...
            'Position',[0 y1 aW 22], ...
            'FontSize',11,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');

        % Separator
        y1 = y1 - 16;
        uipanel(pg1,'Units','pixels','Position',[80 y1 aW-160 1], ...
            'BackgroundColor',C.bdr,'BorderType','none');

        % Purpose + Suite
        y1 = y1 - 30;
        uicontrol(pg1,'Style','text', ...
            'String','Power System Restoration Decision Support', ...
            'Position',[0 y1 aW 26], ...
            'FontSize',12,'FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 30;
        uicontrol(pg1,'Style','text', ...
            'String',['A tool of the GriRes' TM ' Suite'], ...
            'Position',[0 y1 aW 24], ...
            'FontSize',11,'FontName','Segoe UI','FontWeight','bold', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 22;
        uicontrol(pg1,'Style','text','String','(Grid Resilience Suite)', ...
            'Position',[0 y1 aW 20], ...
            'FontSize',10,'FontName','Segoe UI','FontAngle','italic', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');

        % Separator
        y1 = y1 - 16;
        uipanel(pg1,'Units','pixels','Position',[80 y1 aW-160 1], ...
            'BackgroundColor',C.bdr,'BorderType','none');

        % Author
        y1 = y1 - 30;
        uicontrol(pg1,'Style','text','String','Dr. Francis Yunhe Hou', ...
            'Position',[0 y1 aW 24], ...
            'FontSize',12,'FontName','Segoe UI', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 24;
        uicontrol(pg1,'Style','text','String','yunhe.hou@gmail.com', ...
            'Position',[0 y1 aW 20], ...
            'FontSize',10,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');

        % Copyright + License
        y1 = y1 - 32;
        uicontrol(pg1,'Style','text', ...
            'String',[char(169) ' 2017-2026 Francis Yunhe Hou'], ...
            'Position',[0 y1 aW 20], ...
            'FontSize',9,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 20;
        uicontrol(pg1,'Style','text','String','Licensed under GPL-3.0-or-later', ...
            'Position',[0 y1 aW 18], ...
            'FontSize',9,'FontName','Segoe UI', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 18;
        uicontrol(pg1,'Style','text','String','Commercial licensing available', ...
            'Position',[0 y1 aW 18], ...
            'FontSize',9,'FontName','Segoe UI','FontAngle','italic', ...
            'ForegroundColor',C.dim,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        y1 = y1 - 18;
        uicontrol(pg1,'Style','text','String','github.com/GriRes/GRIM', ...
            'Position',[0 y1 aW 18], ...
            'FontSize',9,'FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','center');

        % ---- Page 2: Details ----
        pg2 = uipanel(af,'Units','pixels','Position',[0 48 aW pH], ...
            'BackgroundColor',C.bg,'BorderType','none','Visible','off');
        ln38 = repmat(char(9472),1,38);
        detBody = {
            ''
            ['  ' ln38]
            '  Heritage & Evolution'
            ['  ' ln38]
            ''
            ['  GRIM' TM ' v3.0.0 builds upon the Generic']
            '  Restoration Milestone (GRM) framework,'
            '  originally developed and published in:'
            ''
            '    Y. Hou, C.-C. Liu, K. Sun, P. Zhang,'
            '    S. Liu, and D. Mizumura,'
            '    "Computation of Milestones for Decision'
            '    Support during System Restoration,"'
            '    IEEE Trans. Power Syst., vol. 26,'
            '    pp. 1399-1409, 2011.'
            ''
            ['  ' ln38]
            '  Previous Implementations'
            ['  ' ln38]
            ''
            '  HITS (Healing Touch for Power Systems)'
            '    MATLAB toolbox v1.0 & v2.0'
            '    by Prof. Zhijun Qin'
            '    3-clause BSD license'
            ''
            '  Commercial versions developed in'
            '  parallel for industry applications.'
            ''
            ['  ' ln38]
            '  References (please cite BOTH)'
            ['  ' ln38]
            ''
            '  [1] Y. Hou, C.-C. Liu, K. Sun,'
            '      P. Zhang, S. Liu, and D. Mizumura,'
            '      "Computation of Milestones for'
            '      Decision Support during System'
            '      Restoration,"'
            '      IEEE Trans. Power Syst., vol. 26,'
            '      pp. 1399-1409, 2011.'
            ''
            '  [2] K. Sun, Y. Hou, W. Sun, and J. Qi,'
            '      Power System Control Under'
            '      Cascading Failures,'
            '      John Wiley & Sons, 2019.'
            ''
        };
        uicontrol(pg2,'Style','edit','Max',2,'Min',0, ...
            'String',detBody,'Position',[14 6 aW-28 pH-12], ...
            'FontSize',10,'FontName','Consolas', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.con, ...
            'HorizontalAlignment','left','Enable','inactive');

        aboutPage = 1;
        hPrev = uicontrol(af,'Style','pushbutton','String',[char(9664) ' Back'], ...
            'Position',[20 10 100 30], ...
            'FontSize',10,'FontName','Segoe UI', ...
            'BackgroundColor',C.card,'ForegroundColor',C.txt, ...
            'Enable','off','Callback',@(~,~)flipPage(-1));
        hNext = uicontrol(af,'Style','pushbutton','String',['Details ' char(9654)], ...
            'Position',[aW-120 10 100 30], ...
            'FontSize',10,'FontName','Segoe UI', ...
            'BackgroundColor',C.accent,'ForegroundColor',C.w, ...
            'Callback',@(~,~)flipPage(1));
        uicontrol(af,'Style','pushbutton','String','Close', ...
            'Position',[(aW-100)/2 10 100 30], ...
            'FontSize',10,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',C.card,'ForegroundColor',C.w, ...
            'Callback',@(~,~)delete(af));

        function flipPage(d)
            aboutPage = aboutPage + d;
            if aboutPage <= 1
                aboutPage = 1;
                set(pg1,'Visible','on');  set(pg2,'Visible','off');
                set(hPrev,'Enable','off');
                set(hNext,'Enable','on','String',['Details ' char(9654)]);
            else
                aboutPage = 2;
                set(pg1,'Visible','off'); set(pg2,'Visible','on');
                set(hPrev,'Enable','on');
                set(hNext,'Enable','off');
            end
        end
    end

    function citeDlg()
        cW2 = 580;  cH2 = 420;
        cf = figure('Name','GRIM Citation','NumberTitle','off', ...
            'MenuBar','none','ToolBar','none','Color',C.bg, ...
            'Resize','off','WindowStyle','modal', ...
            'Position',[(scr(3)-cW2)/2 (scr(4)-cH2)/2 cW2 cH2]);
        uipanel(cf,'Units','pixels','Position',[0 cH2-3 cW2 3], ...
            'BackgroundColor',C.accent,'BorderType','none');
        uicontrol(cf,'Style','text', ...
            'String','If you use GRIM, please cite BOTH:', ...
            'Position',[20 cH2-42 cW2-40 26], ...
            'FontSize',12,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.bg,'HorizontalAlignment','center');
        citB = {
            ''
            '  [1] Y. Hou, C.-C. Liu, K. Sun,'
            '      P. Zhang, S. Liu, and D. Mizumura,'
            '      "Computation of Milestones for'
            '      Decision Support during System'
            '      Restoration,"'
            '      IEEE Trans. Power Syst., vol. 26,'
            '      pp. 1399-1409, 2011.'
            ''
            '  [2] K. Sun, Y. Hou, W. Sun, and J. Qi,'
            '      Power System Control Under'
            '      Cascading Failures,'
            '      John Wiley & Sons, 2019.'
            ''
            ['  ' repmat(char(9472),1,44)]
            ''
            '  BibTeX:'
            ''
            '  @article{Hou2011Milestones,'
            '    author  = {Y. Hou and C.-C. Liu and'
            '               K. Sun and P. Zhang and'
            '               S. Liu and D. Mizumura},'
            '    title   = {Computation of Milestones'
            '               for Decision Support during'
            '               System Restoration},'
            '    journal = {IEEE Trans. Power Syst.},'
            '    volume  = {26},'
            '    pages   = {1399--1409},'
            '    year    = {2011}}'
            ''
            '  @book{Sun2019Cascading,'
            '    author    = {K. Sun and Y. Hou and'
            '                 W. Sun and J. Qi},'
            '    title     = {Power System Control'
            '                 Under Cascading Failures},'
            '    publisher = {John Wiley \& Sons},'
            '    year      = {2019}}'
        };
        uicontrol(cf,'Style','edit','Max',2,'Min',0, ...
            'String',citB,'Position',[20 52 cW2-40 cH2-108], ...
            'FontSize',9,'FontName','Consolas', ...
            'ForegroundColor',C.txt,'BackgroundColor',C.con, ...
            'HorizontalAlignment','left','Enable','inactive');
        uicontrol(cf,'Style','pushbutton','String','Close', ...
            'Position',[(cW2-130)/2 12 130 32], ...
            'FontSize',11,'FontWeight','bold','FontName','Segoe UI', ...
            'BackgroundColor',C.accent,'ForegroundColor',C.w, ...
            'Callback',@(~,~)delete(cf));
    end

    %% ============================================================
    %%                    UTILITY FUNCTIONS
    %% ============================================================

    function scanSystems()
        S.systems = {};
        if ~exist(gridDataDir,'dir'), return; end
        entries = dir(gridDataDir);
        for k2 = 1:numel(entries)
            if entries(k2).isdir && entries(k2).name(1) ~= '.'
                gf3 = fullfile(gridDataDir,entries(k2).name,[entries(k2).name '_gen.txt']);
                if exist(gf3,'file'), S.systems{end+1} = entries(k2).name; end %#ok
            end
        end
        S.systems = sort(S.systems);
        if ~isempty(S.systems)
            S.sysIdx = 1; S.sysName = S.systems{1};
        else
            S.sysIdx = 0; S.sysName = '';
        end
    end

    function lst = sysPopList()
        if isempty(S.systems), lst = {'(no systems)'}; else, lst = S.systems; end
    end

    function loadSystemInfo()
        if isempty(S.sysName), clearCards(); clearResults(); return; end
        sd = fullfile(gridDataDir,S.sysName);  pre = S.sysName;
        nb  = cntDataLines(fullfile(sd,[pre '_bus.txt']));
        ng  = cntDataLines(fullfile(sd,[pre '_gen.txt']));
        nl  = cntDataLines(fullfile(sd,[pre '_line.txt']));
        nt  = cntDataLines(fullfile(sd,[pre '_tran.txt']));
        ncl = cntDataLines(fullfile(sd,[pre '_criticalLoad.txt']));
        ndl = cntDataLines(fullfile(sd,[pre '_dispatchableLoad.txt']));
        nbs = 0;
        gf3 = fullfile(sd,[pre '_gen.txt']);
        if exist(gf3,'file')
            try
                fid = fopen(gf3);
                D = textscan(fid,'%d%s%d%f%d%f%f%f%f%f%f%d%*[^\n]', ...
                    'Delimiter',',','CommentStyle','%');
                fclose(fid);
                if numel(D) >= 5 && ~isempty(D{5}), nbs = sum(D{5}==1); end
            catch, end
        end
        vals = {nb, ng, nl, nt, nbs, ncl, ndl};
        for k2 = 1:min(numel(vals),numel(S.hCard))
            set(S.hCard(k2),'String',num2str(vals{k2}));
        end
        clearResults();
    end

    function clearCards()
        for k2 = 1:numel(S.hCard), set(S.hCard(k2),'String','--'); end
    end

    function clearResults()
        S.Grid = [];  S.hasResults = false;
        for k2 = 1:numel(S.hRes)
            set(S.hRes(k2),'String','--','ForegroundColor',C.dim);
        end
        set(S.hResPanel,'ForegroundColor',C.mile);
        set(S.hOutcome,'String', ...
            '  No results yet — run solver to generate', ...
            'ForegroundColor',C.dim);
    end

    function populateResults()
        G2 = S.Grid;
        ns = findLastValidStep(G2);
        if ns == 0
            for k2 = 1:numel(S.hRes)
                set(S.hRes(k2),'String','--','ForegroundColor',C.dim);
            end
            set(S.hOutcome,'String','  No completed steps in results','ForegroundColor',C.amber);
            set(S.hResPanel,'ForegroundColor',C.amber);
            return;
        end

        fn = sprintf('Step_%d', ns);
        tFinal   = safeField(G2.Time,           fn, 0);
        busIsOn  = safeField(G2.Bus.IsOn,       fn, zeros(numel(G2.Bus.ID),1));
        genIsOn  = safeField(G2.Generator.IsOn, fn, zeros(numel(G2.Generator.ID),1));
        pOut     = safeField(G2.Generator.P_Out,fn, zeros(numel(G2.Generator.ID),1));

        nBusOn   = sum(busIsOn);
        nGenOn   = sum(genIsOn);
        nBusAll  = numel(G2.Bus.ID);
        nGenAll  = numel(G2.Generator.ID);

        clPon = 0;
        try
            clIsOn = G2.Load.CriticalLoad.IsOn.(fn);
            clPon  = sum(G2.Load.CriticalLoad.P(clIsOn > 0));
        catch, end
        dlPon = 0;
        try
            dlIsOn = G2.Load.DispatchableLoad.IsOn.(fn);
            dlPout = G2.Load.DispatchableLoad.P_Out.(fn);
            dlPon  = sum(dlPout(dlIsOn > 0));
        catch, end
        totalLoad = clPon + dlPon;

        nbsIdx = find(~G2.Generator.BlackStartFlag);
        if ~isempty(nbsIdx)
            M1 = max(G2.Generator.StartTime(nbsIdx));
            M2 = max(G2.Generator.StartTime(nbsIdx) + ...
                     double(G2.Generator.TimeToCrank(nbsIdx)));
            M3 = 0;
            for gi = nbsIdx(:)'
                sr  = G2.Generator.StartRequirement(gi);
                Pmt = G2.Generator.Capacity(gi) * G2.Generator.MinOutputRate(gi);
                crankEnd = G2.Generator.StartTime(gi) + ...
                           double(G2.Generator.TimeToCrank(gi));
                rampMin  = G2.Generator.RampingRate(gi) / 60;
                if rampMin > 0, M3 = max(M3, crankEnd + (Pmt+sr)/rampMin); end
            end
            % M1 also requires all CLs connected
            nCL_g = length(G2.Load.CriticalLoad.ID);
            for ci = 1:nCL_g
                for s = 1:G2.step
                    clOn_s = GetStepValue(G2, 'Grid.Load.CriticalLoad.IsOn', s);
                    if ~isempty(clOn_s) && clOn_s(ci)
                        fn = sprintf('Step_%d', s);
                        if isfield(G2.Time, fn)
                            M1 = max(M1, G2.Time.(fn));
                        end
                        break;
                    end
                end
            end
        else
            M1 = 0;  M2 = 0;  M3 = 0;
        end

        vals = {num2str(ns), sprintf('%.1f',double(tFinal)), ...
                sprintf('%.1f',M1), sprintf('%.1f',M2), sprintf('%.1f',M3), ...
                sprintf('%d/%d',nBusOn,nBusAll), sprintf('%d/%d',nGenOn,nGenAll), ...
                sprintf('%.3f',sum(pOut(pOut>0))), sprintf('%.3f',totalLoad)};
        for k2 = 1:min(numel(vals),numel(S.hRes))
            set(S.hRes(k2),'String',vals{k2},'ForegroundColor',resClrs{k2});
        end
        updateOutcome(nGenOn,nGenAll,nBusOn,nBusAll);
    end

    function updateOutcome(ngo,nga,nbo,nba)
        if ngo==nga && nbo==nba,  setOutcome(ngo,nga,nbo,nba,'complete');
        elseif ngo==nga,          setOutcome(ngo,nga,nbo,nba,'done');
        else,                     setOutcome(ngo,nga,nbo,nba,'partial');
        end
    end

    function setOutcome(ngo,nga,nbo,nba,state)
        dot = char(9679);  dash = char(8212);
        switch state
            case 'complete'
                set(S.hOutcome,'String', ...
                    sprintf('  %s  COMPLETE %s  All %d generators online, all %d buses energized', ...
                    dot,dash,nga,nba), 'ForegroundColor',C.green);
                set(S.hResPanel,'ForegroundColor',C.green);
            case 'done'
                set(S.hOutcome,'String', ...
                    sprintf('  %s  DONE %s  All %d generators online, %d/%d buses energized', ...
                    dot,dash,nga,nbo,nba), 'ForegroundColor',C.green);
                set(S.hResPanel,'ForegroundColor',C.green);
            case 'partial'
                set(S.hOutcome,'String', ...
                    sprintf('  %s  PARTIAL %s  %d/%d generators, %d/%d buses', ...
                    dot,dash,ngo,nga,nbo,nba), 'ForegroundColor',C.amber);
                set(S.hResPanel,'ForegroundColor',C.amber);
        end
    end

    function n = cntDataLines(filepath)
        n = 0;
        if ~exist(filepath,'file'), return; end
        fid = fopen(filepath);  if fid < 0, return; end
        while ~feof(fid)
            ln2 = fgetl(fid);
            if ischar(ln2)
                t = strtrim(ln2);
                if ~isempty(t) && t(1) ~= '%', n = n+1; end
            end
        end
        fclose(fid);
    end

    function setupPaths()
        pa = fullfile(grimRoot,'GRMs_PathAdd.m');
        if exist(pa,'file')
            try
                owd2 = pwd; cd(grimRoot); GRMs_PathAdd; cd(owd2);
                cLog('GRIM paths configured.'); return;
            catch ME
                cLog(['Path warning: ' ME.message]);
                try cd(owd2); catch, end
            end
        end
        addpath(grimRoot);
        skip = {'GridData','Monitoring','.git','__pycache__'};
        dd = dir(grimRoot);  cnt2 = 0;
        for k2 = 1:numel(dd)
            if dd(k2).isdir && dd(k2).name(1) ~= '.' && ~ismember(dd(k2).name,skip)
                addpath(fullfile(grimRoot,dd(k2).name));  cnt2 = cnt2+1;
            end
        end
        cLog(sprintf('Added %d subdirectories to path.', cnt2));
    end

    function trackChild(hFig)
        if ~isempty(hFig) && ishandle(hFig)
            S.childFigs(end+1) = hFig;
            S.childFigs(~ishandle(S.childFigs)) = [];
        end
    end

    function cLog(msg)
        ts = sprintf('[%s]  %s', datestr(now,'HH:MM:SS'), msg); %#ok
        appendToConsole(ts);
    end

    function cLogRaw(msg)
        appendToConsole(['        ' msg]);
    end

    function appendToConsole(txt)
        if isempty(S.hCon) || ~ishandle(S.hCon), return; end
        old = get(S.hCon,'String');
        if ischar(old), old = {old}; end
        old{end+1} = txt;
        set(S.hCon,'String',old);
        drawnow;
    end

    function scrollBottom()
        if isempty(S.hCon) || ~ishandle(S.hCon), return; end
        n = numel(get(S.hCon,'String'));
        if n > 0, set(S.hCon,'ListboxTop',n); drawnow; end
    end

    function setStat(msg,clr)
        if isempty(S.hStat) || ~ishandle(S.hStat), return; end
        set(S.hStat,'String',msg,'ForegroundColor',clr);
    end

    function sideHdr(x,y,str)
        uicontrol(fig,'Style','text','String',str, ...
            'Position',[x y bw 16],'FontSize',9,'FontWeight','bold','FontName','Segoe UI', ...
            'ForegroundColor',C.accent,'BackgroundColor',C.side,'HorizontalAlignment','left');
    end

    function txt = welcomeText()
        TM = char(8482);
        txt = {
            ''
            ['  GRIM' TM '  v3.0.0']
            '  Generic Restoration Intelligent Milestone'
            ['  A tool of the GriRes' TM ' Suite (Grid Resilience Suite)']
            ''
            ['  ' char(169) ' 2017-2026  Dr. Francis Yunhe Hou']
            '  Licensed under GPL-3.0-or-later  |  Commercial license available'
            '  github.com/GriRes/GRIM'
            ''
            ['  Root:     ' grimRoot]
            ['  GridData: ' gridDataDir]
            ''
            sprintf('  %d system(s) available.', numel(S.systems))
            ''
            '  Workflow:'
            '    1. Select system from dropdown'
            '    2. Choose solver mode & config profile'
            '    3. Click [Generate Restoration Strategy]'
            '    4. View results (HTML / Step / Plot / Report)'
            ''
            '  Ctrl+C in Command Window to stop a running solve.'
            '  Help > User Guide for full instructions.'
            ''
        };
    end

    function ns = findLastValidStep(G2)
        ns = 0;
        if ~isstruct(G2),                                    return; end
        if ~isfield(G2,'Time') || ~isstruct(G2.Time),        return; end
        if isfield(G2,'step') && isnumeric(G2.step) && G2.step > 0
            candidate = sprintf('Step_%d', G2.step);
            if isfield(G2.Time, candidate), ns = G2.step; return; end
        end
        fns = fieldnames(G2.Time);
        for ii = 1:numel(fns)
            tok = regexp(fns{ii}, '^Step_(\d+)$', 'tokens');
            if ~isempty(tok)
                sn = str2double(tok{1}{1});
                if sn > ns, ns = sn; end
            end
        end
    end

    function val = safeField(strct, fieldName, default)
        if isstruct(strct) && isfield(strct, fieldName)
            val = strct.(fieldName);
        else
            val = default;
        end
    end

    function html = md2html(raw)
        lines = strsplit(raw, char(10));
        css = ['body{margin:16px;font-family:Segoe UI,Arial,sans-serif;' ...
               'font-size:13px;line-height:1.6;}' ...
               'h1{color:#00b4d8;font-size:18px;' ...
               'border-bottom:1px solid #2d3241;padding-bottom:6px;}' ...
               'h2{color:#e8ad23;font-size:15px;margin-top:18px;}' ...
               'h3{color:#1e72c6;font-size:13px;margin-top:14px;}' ...
               'table{border-collapse:collapse;margin:10px 0;}' ...
               'th,td{border:1px solid #2d3241;' ...
               'padding:4px 10px;text-align:left;font-size:12px;}' ...
               'th{background:#1a2030;color:#00b4d8;}' ...
               'code{background:#1a2030;padding:1px 4px;' ...
               'font-family:Consolas,monospace;}' ...
               'hr{border:none;border-top:1px solid #2d3241;margin:14px 0;}' ...
               'strong{color:#f0f0f0;}'];
        html = ['<html><head><style>' css '</style></head>' ...
                '<body bgcolor="#0b0d12" text="#d7dde8">'];
        inTable  = false;
        firstTR  = false;
        for ii = 1:numel(lines)
            ln2 = lines{ii};
            tr  = strtrim(ln2);
            if isempty(tr)
                if inTable, html = [html '</table>']; inTable = false; end %#ok
                continue;
            end
            if ~inTable && numel(tr) >= 3 && (all(tr=='-') || all(tr=='='))
                html = [html '<hr>']; continue; %#ok
            end
            if numel(tr) > 4 && strcmp(tr(1:4),'### ')
                if inTable, html = [html '</table>']; inTable = false; end %#ok
                html = [html '<h3>' inlineFmt(tr(5:end)) '</h3>']; continue; %#ok
            elseif numel(tr) > 3 && strcmp(tr(1:3),'## ')
                if inTable, html = [html '</table>']; inTable = false; end %#ok
                html = [html '<h2>' inlineFmt(tr(4:end)) '</h2>']; continue; %#ok
            elseif numel(tr) > 2 && strcmp(tr(1:2),'# ')
                if inTable, html = [html '</table>']; inTable = false; end %#ok
                html = [html '<h1>' inlineFmt(tr(3:end)) '</h1>']; continue; %#ok
            end
            if tr(1) == '|'
                if ~isempty(regexp(tr,'^\|[\s\-:|]+\|$','once')), continue; end
                if ~inTable
                    html = [html '<table border="1" cellpadding="4" cellspacing="0">']; %#ok
                    inTable = true;  firstTR = true;
                end
                cells = strsplit(tr,'|');
                if ~isempty(cells) && isempty(strtrim(cells{1})),   cells = cells(2:end);   end
                if ~isempty(cells) && isempty(strtrim(cells{end})), cells = cells(1:end-1); end
                if firstTR, tag = 'th'; firstTR = false; else, tag = 'td'; end
                html = [html '<tr>']; %#ok
                for ci5 = 1:numel(cells)
                    html = [html '<' tag '>' inlineFmt(strtrim(cells{ci5})) '</' tag '>']; %#ok
                end
                html = [html '</tr>']; continue; %#ok
            end
            if inTable, html = [html '</table>']; inTable = false; end %#ok
            html = [html '<p>' inlineFmt(tr) '</p>']; %#ok
        end
        if inTable, html = [html '</table>']; end %#ok
        html = [html '</body></html>'];
    end

    function out = inlineFmt(txt)
        out = regexprep(txt,'\*\*([^*]+)\*\*','<strong>$1</strong>');
        out = regexprep(out,'`([^`]+)`','<code>$1</code>');
    end

    function out = imresize_nn(img, tgt)
    %imresize_nn  Nearest-neighbour resize (no Image Processing Toolbox)
        [h, w, ch] = size(img);
        if h == w
            sz = tgt;
        else
            sz = round([tgt, tgt * w / h]);
        end
        rr = round(linspace(1, h, sz(1)));
        cc = round(linspace(1, w, sz(end)));
        out = img(rr, cc, :);
    end

end