restore original b4gic

test/data/matpower/b4gic.m

@@ -0,0 +1,143 @@
+%% MATPOWER Case Format : Version 2
+function mpc = b4gic
+mpc.version = '2';
+
+
+%%-----  Power Flow Data  -----%%
+
+%% system MVA base
+mpc.baseMVA = 100;
+
+
+%% bus data
+%    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    zone    Vmax    Vmin
+mpc.bus = [
+	1	1	0	   0	  0	0	1	1.100000	0.000000	765	1	1.15	0.85
+	2	1	0	   0	  0	0	1	1.100000	0.000000	765	1	1.15	0.85
+	3	1	1000 200	0	0	1	1.100000	0.000000	20	1	1.15	0.85
+	4	3	0	   0	  0	0	1	1.000000	0.000000	20	1	1.15	0.85
+];
+
+
+%% generator data
+%    bus    Pg    Qg    Qmax    Qmin    Vg    mBase    status    Pmax    Pmin    Pc1    Pc2    Qc1min    Qc1max    Qc2min    Qc2max    ramp_agc    ramp_10    ramp_30    ramp_q    apf
+mpc.gen = [
+	4	0.0	0.0	2000.0	-2000.0	1.1	100	1	2000.0	0.0	0	0	0	0	0	0	0	0	0	-1000	0
+];
+
+
+%% branch data
+%    fbus    tbus    r    x    b    rateA    rateB    rateC    ratio    angle    status    angmin    angmax
+mpc.branch = [
+	1	3	0.0001	  0.004	0	2000.0	0.0	0.0	1	0.0	1	-30.0	30.0
+	1	2	0.000513	0.01	0	2000.0	0.0	0.0	1	0.0	1	-30.0	30.0
+	2	4	0.0001	  0.004	0	2000.0	0.0	0.0	1	0.0	1	-30.0	30.0
+];
+
+
+%%-----  OPF Data  -----%%
+
+%% generator cost data
+%    1    startup    shutdown    n    x1    y1    ...    xn    yn
+%    2    startup    shutdown    n    c(n-1)    ...    c0
+mpc.gencost = [
+	2	0	0	3	0.11	5.0	0
+];
+
+
+%%-----  GMD - Thermal Data  -----%%
+
+%% gmd_bus data
+%column_names% parent_index status g_gnd name
+mpc.gmd_bus = {
+	1	1	5	'dc_sub1'
+	2	1	5	'dc_sub2'
+	1	1	0	'dc_bus1'
+	2	1	0	'dc_bus2'
+	3	1	0	'dc_bus3'
+	4	1	0	'dc_bus4'
+};
+
+
+%% gmd_branch data
+%column_names% f_bus t_bus parent_index br_status br_r br_v len_km name
+mpc.gmd_branch = {
+	3	1	1	1	0.1	0	0	'dc_xf1_hi'
+	3	4	2	1	1.00073475	170.78806587354	170.78806587354	'dc_br1'
+	4	2	3	1	0.1	0	0	'dc_xf2_hi'
+};
+
+
+%% branch_gmd data
+%column_names% hi_bus lo_bus gmd_br_hi gmd_br_lo gmd_k gmd_br_series gmd_br_common baseMVA type config
+mpc.branch_gmd = {
+	1	3	1	-1	1.793	-1	-1	100	'xfmr'	'gwye-delta'
+	1	2	-1	-1	-1	-1	-1	100	'line'	'none'
+	2	4	3	-1	1.793	-1	-1	100	'xfmr'	'gwye-delta'
+};
+
+
+%% branch_thermal data
+%column_names% xfmr temperature_ambient hotspot_instant_limit hotspot_avg_limit hotspot_rated topoil_time_const topoil_rated topoil_init topoil_initialized hotspot_coeff
+mpc.branch_thermal = {
+	1	25	280	240	150	71	75	0	1	0.63
+	0	-1	-1	-1	-1	-1	-1	-1	-1	-1
+	1	25	280	240	150	71	75	0	1	0.63
+};
+
+
+%% bus_gmd data
+%column_names% lat lon
+mpc.bus_gmd = {
+	40	-89
+	40	-87
+	40	-89
+	40	-87
+};
+
+
+%% time_elapsed
+%column_names% seconds
+mpc.time_elapsed = 10.0;
+
+
+%% thermal caps
+% thermal_cap_x0 ([per unit])
+%column_names% A B C D E F G H I J K
+mpc.thermal_cap_x0 = [
+	0.23033 0.25000 0.26438 0.27960 0.30000 0.31967 0.33942 0.36153 0.38444 0.40000 0.43894
+];
+% thermal_cap_y0 ([percent per unit])
+%column_names% A B C D E F G H I J K
+mpc.thermal_cap_y0 = [
+	100.0 93.94 90.0 85.42 80.0 74.73 70.0 64.94 59.97 56.92 50.0
+];
+% Values are from Fig.2. of https://arxiv.org/pdf/1701.01469.pdf paper
+
+
+%%-----  SourceID Data  -----%%
+
+%% bus_sourceid data
+%column_names% bus_sid
+mpc.bus_sourceid = [
+	1;
+	2;
+	3;
+	4;
+];
+
+
+%% gen_sourceid data
+%column_names% bus_i gen_sid
+mpc.gen_sourceid = [
+	4 '1 ';
+];
+
+
+%% branch_sourceid data
+%column_names% fbus tbus branch_sid
+mpc.branch_sourceid = [
+	1 3 '1 ';
+	1 2 '1 ';
+	2 4 '1 ';
+];

test/test_cases.jl

@@ -19,7 +19,7 @@ data_epri21_verification = "../test/data/pw_csv/verification/EpriCase_Mar2012_GI
 
 data_b4gic_vm = "../test/data/pw_csv/verification/B4GIC_Slack_Only_Bus_1Vkm.csv"
 data_epri21_vm = "../test/data/pw_csv/verification/EpriCase_Mar2012_Bus_1Vkm.csv"
-b4gic_default = "../test/data/matpower/b4gic.m"
+b4gic_default = "../test/data/matpower/b4gic_two_gens.m"
 b4gic_offbase = "../test/data/matpower/b4gic_offbase.m"
 autotransformer = "../test/data/matpower/autotransformer.m"
 epricase = "../test/data/matpower/epricase_aug2022_v22_fix.m"