import numpy as np
[docs]
def get_reference_stats(ref_sys_id: int) -> dict:
"""
Returns the statistical parameters for a specified reference system.
Ported from 'sg_refsys_stat.m' (SynGrid).
Args:
ref_sys_id: 1, 2, or 3.
Returns:
Dictionary containing statistical tables and parameters.
"""
stats = {}
if ref_sys_id == 1:
stats['system_name'] = 'Reference System 1'
# --- Generation & Load Tables ---
stats['Tab_2D_Pgmax'] = np.array([
[0,0,0.0826359832635983,0.0962343096234310,0.0334728033472803,0.0271966527196653,0.0251046025104603,0.00941422594142259,0.00313807531380753,0.00627615062761506,0,0,0,0],
[0,0,0.109832635983264,0.0261506276150628,0.0135983263598326,0.0240585774058577,0.0104602510460251,0.0125523012552301,0.00627615062761506,0.00104602510460251,0.00104602510460251,0,0,0],
[0,0,0.0334728033472803,0.00627615062761506,0.00418410041841004,0.00732217573221757,0.0167364016736402,0.0115062761506276,0.00523012552301255,0.00104602510460251,0,0,0,0.00104602510460251],
[0,0,0.0742677824267782,0.0125523012552301,0.00313807531380753,0.0104602510460251,0.0167364016736402,0.0167364016736402,0.0125523012552301,0,0,0.00104602510460251,0,0],
[0,0,0.0491631799163180,0.00836820083682008,0.00313807531380753,0.0104602510460251,0.0135983263598326,0.0240585774058577,0.0177824267782427,0.00941422594142259,0.00418410041841004,0,0.00104602510460251,0],
[0,0,0.0125523012552301,0.00104602510460251,0.00523012552301255,0.00523012552301255,0.00941422594142259,0.0135983263598326,0.0135983263598326,0.00836820083682008,0.00836820083682008,0,0,0],
[0,0,0.00104602510460251,0,0,0.00418410041841004,0.00209205020920502,0.00418410041841004,0.00732217573221757,0.00313807531380753,0,0,0,0.00104602510460251],
[0,0,0,0.00104602510460251,0.00104602510460251,0,0.00209205020920502,0.00313807531380753,0.00313807531380753,0.00313807531380753,0.00209205020920502,0.00104602510460251,0.00104602510460251,0],
[0,0,0,0,0,0,0,0.00104602510460251,0.00209205020920502,0.00209205020920502,0.00104602510460251,0.00104602510460251,0.00104602510460251,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0.00104602510460251,0],
[0,0,0,0,0,0,0,0,0.00104602510460251,0,0.00104602510460251,0,0,0.00209205020920502],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0.00104602510460251,0,0,0],
[0,0,0,0,0,0,0.00209205020920502,0.00104602510460251,0,0,0,0,0,0]
])
stats['Tab_2D_load'] = np.array([
[0,0,0.0143769968051118,0.116613418530351,0.102236421725240,0.0559105431309904,0.0199680511182109,0.00798722044728435,0.00239616613418530,0.00159744408945687,0,0,0,0],
[0,0,0.00159744408945687,0.0383386581469649,0.0790734824281150,0.0439297124600639,0.0151757188498403,0.00559105431309904,0.00319488817891374,0.00239616613418530,0,0,0,0],
[0,0,0.000798722044728435,0.0207667731629393,0.0375399361022364,0.0303514376996805,0.0151757188498403,0.00399361022364217,0.00239616613418530,0.00159744408945687,0,0,0,0],
[0,0,0.000798722044728435,0.0135782747603834,0.0351437699680511,0.0535143769968051,0.0295527156549521,0.0255591054313099,0.00479233226837061,0.000798722044728435,0.000798722044728435,0,0,0],
[0,0,0,0.00479233226837061,0.00399361022364217,0.0151757188498403,0.0215654952076677,0.0295527156549521,0.0143769968051118,0.00798722044728435,0.00159744408945687,0,0,0],
[0,0,0,0,0.00479233226837061,0.00479233226837061,0.00479233226837061,0.0175718849840256,0.0119808306709265,0.00638977635782748,0.00159744408945687,0.000798722044728435,0,0],
[0,0,0,0,0,0.00239616613418530,0.00718849840255591,0.00958466453674121,0.00559105431309904,0.00479233226837061,0.000798722044728435,0.00159744408945687,0,0],
[0,0,0,0,0,0,0.000798722044728435,0.00399361022364217,0.00159744408945687,0.00399361022364217,0,0.00159744408945687,0,0],
[0,0,0,0,0,0,0.00159744408945687,0.000798722044728435,0,0.00319488817891374,0,0.00239616613418530,0,0.000798722044728435],
[0,0,0,0,0,0.000798722044728435,0,0,0,0,0.000798722044728435,0,0,0],
[0,0,0,0,0,0,0,0.000798722044728435,0,0,0.00159744408945687,0.00159744408945687,0.000798722044728435,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0.000798722044728435],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0.000798722044728435,0,0,0]
])
stats['Tab_2D_Pg'] = np.array([
[0,0,0,0,0,0.452586206896552,0.0517241379310345,0.0129310344827586,0.00862068965517241,0],
[0,0,0,0,0,0.0474137931034483,0.0129310344827586,0.0172413793103448,0,0],
[0,0,0,0,0,0.0172413793103448,0.00431034482758621,0.00431034482758621,0.00431034482758621,0.00431034482758621],
[0,0,0,0,0,0.0215517241379310,0.0301724137931035,0.0129310344827586,0.0301724137931035,0.0129310344827586],
[0,0,0,0,0,0.0129310344827586,0.0258620689655172,0,0.00431034482758621,0.00431034482758621],
[0,0,0,0,0,0.0129310344827586,0.0129310344827586,0.0129310344827586,0.00431034482758621,0.0129310344827586],
[0,0,0,0,0,0,0,0.0172413793103448,0.0129310344827586,0.0129310344827586],
[0,0,0,0,0,0,0.00431034482758621,0.00862068965517241,0.0129310344827586,0.0172413793103448],
[0,0,0,0,0,0,0.00431034482758621,0.00862068965517241,0.00862068965517241,0.00431034482758621],
[0,0,0,0,0,0.00431034482758621,0.00431034482758621,0,0,0.00431034482758621],
[0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0.00862068965517241],
[0,0,0,0,0,0,0,0,0.00431034482758621,0.00431034482758621],
[0,0,0,0,0,0,0,0.00862068965517241,0,0]
])
# --- Transmission Parameters ---
# Rows: Flow Class, Cols: Beta Class
stats['Tab_2D_FlBeta'] = np.array([
[0.0278086763070078,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0528364849833148,0,0.0116796440489433,0.0350389321468298,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0127919911012236,0.0233592880978865,0.0233592880978865,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0177975528364850,0.0344827586206897,0.00667408231368187,0.0300333704115684,0,0.00278086763070078,0,0,0,0,0,0,0,0,0,0],
[0,0.0139043381535039,0.00500556173526140,0.00945494994438265,0.0139043381535039,0.0133481646273637,0,0,0,0,0,0,0,0,0,0],
[0,0.0216907675194661,0.00389321468298109,0.0361512791991101,0.0200222469410456,0.0200222469410456,0.0111234705228031,0,0,0.00444938820912125,0,0,0,0,0,0],
[0,0.0233592880978865,0.00166852057842047,0.0177975528364850,0.0250278086763070,0,0.0211345939933259,0.00333704115684093,0,0,0,0,0,0,0,0],
[0,0.0105672969966630,0.0200222469410456,0.0200222469410456,0.0228031145717464,0.0250278086763070,0.00389321468298109,0.00611790878754171,0.00500556173526140,0,0,0,0,0,0,0],
[0,0,0.0417130144605117,0.0177975528364850,0.0200222469410456,0.0155728587319244,0.00889877641824249,0,0.00389321468298109,0.00333704115684093,0,0,0,0.00166852057842047,0,0],
[0,0,0,0.00667408231368187,0.0205784204671858,0.0200222469410456,0.0100111234705228,0.0127919911012236,0.00166852057842047,0.00389321468298109,0.00444938820912125,0.00333704115684093,0,0,0,0],
[0,0,0,0.00389321468298109,0.0127919911012236,0,0.00556173526140156,0.0111234705228031,0.00611790878754171,0.00166852057842047,0.00444938820912125,0,0,0,0,0],
[0,0,0,0.00166852057842047,0.0122358175750834,0.00945494994438265,0,0.0194660734149055,0,0.00333704115684093,0.00444938820912125,0,0,0,0,0],
[0,0,0,0,0.00834260289210234,0,0.00166852057842047,0,0.00723025583982202,0.00444938820912125,0,0,0,0,0,0],
[0,0,0,0,0.00389321468298109,0.00889877641824249,0.00444938820912125,0.00556173526140156,0,0.00222469410456062,0.00166852057842047,0,0.00166852057842047,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
])
stats['stab'] = [1.374, -0.838, 2.965, 85.801]
stats['Overload_b'] = 0.00083
stats['mu_beta'] = 0.27
stats['Zpr_pars'] = [-2.38, 1.99, 1.99] # LogN params (Clip?)
stats['mu_committed'] = 0.073
stats['Alpha_mod'] = 0
elif ref_sys_id == 2:
stats['system_name'] = 'Reference System 2'
stats['Tab_2D_Pgmax'] = np.array([
[0,0,0,0,0.123046092184369,0,0.0196392785571142,0.00480961923847695,0.00120240480961924,0.000801603206412826,0,0,0,0.000801603206412826],
[0,0,0,0,0.129458917835671,0,0.00721442885771543,0.00280561122244489,0.00120240480961924,0.000400801603206413,0.000400801603206413,0,0,0],
[0,0,0,0,0.0953907815631263,0,0.00480961923847695,0.00120240480961924,0,0.000400801603206413,0,0,0,0],
[0,0,0,0,0.228456913827655,0,0.0132264529058116,0.00240480961923848,0.000400801603206413,0,0,0.000801603206412826,0.00120240480961924,0],
[0,0,0,0,0.127454909819639,0,0.00320641282565130,0.00120240480961924,0.000400801603206413,0.000400801603206413,0,0,0,0],
[0,0,0,0,0.118637274549098,0,0.00200400801603206,0,0,0.000400801603206413,0,0,0,0],
[0,0,0,0,0.0392785571142285,0,0.00120240480961924,0,0,0,0,0,0,0],
[0,0,0,0,0.0288577154308617,0,0.00280561122244489,0,0,0,0,0,0,0],
[0,0,0,0,0.0128256513026052,0,0.00200400801603206,0,0,0,0,0,0,0],
[0,0,0,0,0.00561122244488978,0,0.000400801603206413,0,0,0,0,0,0,0],
[0,0,0,0,0.00440881763527054,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0.00400801603206413,0,0,0.000801603206412826,0,0,0,0,0,0],
[0,0,0,0,0.00120240480961924,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0.00280561122244489,0,0,0,0,0,0,0,0,0]
])
stats['Tab_2D_load'] = np.array([
[0,0,0,0.166854384553500,0,0.216572807723250,0.0328238133547868,0.0111021721641191,0.00772325020112631,0.00241351568785197,0.00193081255028158,0.000965406275140788,0.000160901045856798,0.000160901045856798],
[0,0,0,0.0872083668543846,0,0.125181013676589,0.0226870474658085,0.00981496379726468,0.00627514078841513,0.000643604183427192,0.00128720836685438,0.000482703137570394,0.000160901045856798,0.000160901045856798],
[0,0,0,0.0265486725663717,0,0.0571198712791633,0.0107803700724055,0.00450522928399035,0.00386162510056315,0.00128720836685438,0.000482703137570394,0.000160901045856798,0.000321802091713596,0.000160901045856798],
[0,0,0,0.0210780370072405,0,0.0619469026548673,0.0175382139983910,0.00756234915526951,0.00740144810941271,0.00160901045856798,0.00128720836685438,0.000804505229283990,0.000643604183427192,0.000321802091713596],
[0,0,0,0.00563153660498793,0,0.0210780370072405,0.00852775543041030,0.00418342719227675,0.00305711987127916,0.00112630732099759,0.000965406275140788,0.000321802091713596,0.000482703137570394,0.000160901045856798],
[0,0,0,0.00144810941271118,0,0.00611423974255833,0.00353982300884956,0.00209171359613838,0.00160901045856798,0,0,0,0.000160901045856798,0.000321802091713596],
[0,0,0,0.00209171359613838,0,0.000804505229283990,0.00193081255028158,0.000965406275140788,0.000160901045856798,0,0.000321802091713596,0,0,0],
[0,0,0,0.000321802091713596,0,0.000482703137570394,0.00128720836685438,0.000321802091713596,0.000482703137570394,0,0,0,0,0],
[0,0,0,0.000160901045856798,0,0.000160901045856798,0.000160901045856798,0.000321802091713596,0.000482703137570394,0,0,0.000160901045856798,0,0.000160901045856798],
[0,0,0,0,0,0.000965406275140788,0.000160901045856798,0.000482703137570394,0.000482703137570394,0.000321802091713596,0.000160901045856798,0,0,0],
[0,0,0,0,0,0.000321802091713596,0.000321802091713596,0.000643604183427192,0.000160901045856798,0,0.000160901045856798,0,0,0],
[0,0,0,0,0,0,0.000643604183427192,0.000482703137570394,0.000321802091713596,0.000160901045856798,0,0.000160901045856798,0,0],
[0,0,0,0,0,0,0.000160901045856798,0.000482703137570394,0.000643604183427192,0.000321802091713596,0.000160901045856798,0.000160901045856798,0,0],
[0,0,0,0,0,0.000160901045856798,0,0.000160901045856798,0,0,0,0.000160901045856798,0,0]
])
stats['Tab_2D_Pg'] = np.array([
[0,0,0,0,0,0.0227129337539432,0.0309148264984227,0.0309148264984227,0.0479495268138801,0.0479495268138801],
[0,0,0,0,0,0.00567823343848580,0.0107255520504732,0.0208201892744480,0.0479495268138801,0.0687697160883281],
[0,0,0,0,0,0,0.00946372239747634,0.0119873817034700,0.0271293375394322,0.0744479495268139],
[0,0,0,0,0,0,0.00189274447949527,0.0201892744479495,0.0460567823343849,0.133753943217666],
[0,0,0,0,0,0,0.000630914826498423,0.00567823343848580,0.0290220820189274,0.0940063091482650],
[0,0,0,0,0,0.000630914826498423,0,0.00378548895899054,0.0151419558359621,0.0876971608832808],
[0,0,0,0,0,0,0,0,0.00757097791798107,0.0283911671924290],
[0,0,0,0,0,0,0,0,0.00126182965299685,0.0214511041009464],
[0,0,0,0,0,0,0,0,0,0.0107255520504732],
[0,0,0,0,0,0,0,0,0.00126182965299685,0.00567823343848580],
[0,0,0,0,0,0,0,0,0,0.00441640378548896],
[0,0,0,0,0,0,0,0,0,0.00189274447949527],
[0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0.00126182965299685]
])
# Transmission Parameters
stats['Tab_2D_FlBeta'] = np.array([
[0.117410287897394,6.82221312593806e-05,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0165097557647701,0,0.000136444262518761,0,0,0,0,0,0,0,0,0,0,0],
[0.0180106426524765,0.000409332787556283,0.000204666393778142,0.000204666393778142,0,0,0,0,0,0,0,0,0,0],
[0.0484377131941602,0.00279710738163460,0.00102333196889071,0.000409332787556283,0,0,0,0,0,0,0,0,0,0],
[0.0535543730386137,0.00893709919497885,0.00191021967526266,0.00129622049392823,0.000682221312593805,0.000750443443853186,6.82221312593806e-05,0,0,0,0,0,0,0],
[0.0541683722199482,0.0141902033019512,0.00661754673215991,0.00429799426934097,0.000682221312593805,0.000545777050075044,0.000545777050075044,0,6.82221312593806e-05,0,0,0,0,0],
[0.0500068222131259,0.0266066311911584,0.00818665575112566,0.00893709919497885,0.00191021967526266,0.000545777050075044,0.00102333196889071,0.000818665575112567,0,6.82221312593806e-05,0,0,0,0],
[0.0325419566107245,0.0346568426797653,0.0119388729703916,0.0143266475644699,0.00491199345067540,0.00259244098785646,0.000682221312593805,0.00143266475644699,0.000682221312593805,0.000545777050075044,0,0,0,0],
[0.0221721926592987,0.0332924000545777,0.0140537590394324,0.0198526401964797,0.00866421066994133,0.00313821803793151,0.000818665575112567,0.000341110656296903,0.000886887706371947,0.000409332787556283,0,0,0,0],
[0.0161686451084732,0.0256515213535271,0.0185564197025515,0.0196479738027016,0.0146677582207668,0.00300177377541274,0.00238777459407832,0.000409332787556283,0.000682221312593805,0.000955109837631328,0.000477554918815664,0.000136444262518761,0,0],
[0.0127575385455042,0.0171919770773639,0.0166462000272889,0.0201255287215173,0.0112566516577978,0.00463910492563788,0.00211488606904080,0.00136444262518761,0.00109155410015009,0.000409332787556283,0.000341110656296903,0.000613999181334425,0.000272888525037522,0],
[0.00941465411379452,0.00948287624505390,0.0105744303452040,0.0171237549461045,0.0108473188702415,0.00252421885659708,0.00102333196889071,0.00184199754400327,0.000341110656296903,0.000409332787556283,0.000272888525037522,6.82221312593806e-05,0.00136444262518761,0.000272888525037522],
[0.00880065493246009,0.00689043525719744,0.00491199345067540,0.0121435393641697,0.00920998772001637,0.00334288443170965,0.00150088688770637,0.00115977623140947,0.000204666393778142,0.000272888525037522,0.000477554918815664,6.82221312593806e-05,0.000204666393778142,0.000136444262518761],
[0.00470732705689726,0.00532132623823168,0.00286532951289398,0.00682221312593805,0.00579888115704735,0.00266066311911584,0.00170555328148451,0.000545777050075044,6.82221312593806e-05,6.82221312593806e-05,6.82221312593806e-05,0.000409332787556283,0.000204666393778142,0.000682221312593805],
[0,6.82221312593806e-05,6.82221312593806e-05,0.000204666393778142,6.82221312593806e-05,0.000204666393778142,0,0,0,0,0.000204666393778142,0,6.82221312593806e-05,0.000613999181334425],
[0,0,0,0,0,0,0,0,0.000136444262518761,0,0,0,0,0]
])
stats['stab'] = [1.133,-0.988,3.886,85.713] # levy distribution parameters for transmission line gauge ratio generation
stats['Alpha_mod'] = 1
stats['Overload_b'] = 0.00111
stats['mu_beta'] = 0.28 # beta parameters for transmission line gauge ration generation ~ exponential distribution
stats['Zpr_pars'] = [-3.32, 1.85, 1.99]
stats['mu_committed'] = 0.069
elif ref_sys_id == 3:
stats['system_name'] = 'Reference System 3'
stats['Tab_2D_Pgmax'] = np.array([
[0,0,0,0,0,0.0297029702970297,0,0,0.00165016501650165,0,0,0,0,0],
[0,0,0,0,0,0.0561056105610561,0,0.00495049504950495,0,0,0,0,0,0],
[0,0,0,0,0,0.0544554455445545,0,0,0.00165016501650165,0.00165016501650165,0,0,0,0],
[0,0,0,0,0,0.244224422442244,0,0.0148514851485149,0,0,0,0,0,0],
[0,0,0,0,0,0.141914191419142,0,0.00495049504950495,0,0.00330033003300330,0,0,0,0],
[0,0,0,0,0,0.219471947194719,0,0.00660066006600660,0,0,0,0,0,0],
[0.00165016501650165,0,0,0,0,0.0627062706270627,0,0.00660066006600660,0.00165016501650165,0,0,0,0,0],
[0,0,0,0,0,0.0165016501650165,0,0,0,0,0,0.00165016501650165,0,0],
[0,0,0,0,0,0.0363036303630363,0,0,0,0,0,0,0,0],
[0.00165016501650165,0,0,0,0,0.0280528052805281,0,0.00165016501650165,0.00330033003300330,0.00330033003300330,0,0.00165016501650165,0.00825082508250825,0.00330033003300330],
[0,0,0,0,0,0.0115511551155116,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0.0181518151815182,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0.00660066006600660,0,0,0,0,0,0]
])
stats['Tab_2D_load'] = np.array([
[0,0,0,0,0.0711462450592885,0.145915678524374,0,0.0184453227931489,0.00724637681159420,0.00164690382081686,0,0.000988142292490119,0,0.000329380764163373],
[0,0,0,0,0.0306324110671937,0.108366271409750,0,0.0174571805006588,0.00494071146245059,0.000658761528326746,0,0,0,0],
[0,0,0,0,0.0158102766798419,0.0764163372859025,0,0.0108695652173913,0.00856389986824769,0.000658761528326746,0.000988142292490119,0.000329380764163373,0,0],
[0,0,0,0,0.0332674571805007,0.142951251646904,0,0.0210803689064559,0.0154808959156785,0.00197628458498024,0.00131752305665349,0.000658761528326746,0,0],
[0,0,0,0,0.0125164690382082,0.0800395256916996,0,0.0151515151515152,0.00955204216073781,0.00362318840579710,0.000658761528326746,0.000329380764163373,0,0.000658761528326746],
[0,0,0,0,0.00362318840579710,0.0497364953886693,0,0.0135046113306983,0.00625823451910408,0.00197628458498024,0.000329380764163373,0,0.000329380764163373,0,0.000329380764163373],
[0,0,0,0,0.00164690382081686,0.0217391304347826,0,0.00559947299077734,0.00362318840579710,0.00164690382081686,0,0.000658761528326746,0,0.000658761528326746],
[0,0,0,0,0.000988142292490119,0.00790513833992095,0,0.00197628458498024,0.000329380764163373,0.000988142292490119,0.000329380764163373,0.000658761528326746,0,0.000329380764163373],
[0,0,0,0,0.000329380764163373,0.00592885375494071,0,0.000658761528326746,0.000988142292490119,0,0.000329380764163373,0.000658761528326746,0,0],
[0,0,0,0,0,0.00131752305665349,0,0.000329380764163373,0.000658761528326746,0,0,0.000329380764163373,0,0],
[0,0,0,0,0,0.000329380764163373,0,0,0.000329380764163373,0.00131752305665349,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0.000329380764163373],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0.000329380764163373]
])
stats['Tab_2D_Pg'] = np.array([
[0,0,0,0,0,0.208530805687204,0.0236966824644550,0.0284360189573460,0.0142180094786730,0],
[0,0,0,0,0,0.123222748815166,0.00473933649289100,0,0.00947867298578199,0.0189573459715640],
[0,0,0,0,0,0.0142180094786730,0.0236966824644550,0.0284360189573460,0.0473933649289100,0.0189573459715640],
[0,0,0,0,0,0,0,0.0521327014218009,0.137440758293839,0.0947867298578199],
[0,0,0,0,0,0,0.00473933649289100,0.0142180094786730,0.0521327014218009,0.00947867298578199],
[0,0,0,0,0,0,0,0,0.00473933649289100,0.0284360189573460],
[0,0,0,0,0,0,0,0,0,0.00473933649289100],
[0,0,0,0,0,0,0,0,0.00947867298578199,0],
[0,0,0,0,0,0,0,0,0.00473933649289100,0.00947867298578199],
[0,0,0,0,0,0,0,0,0.00473933649289100,0.00473933649289100],
[0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0]
])
# Transmission Parameters
stats['Tab_2D_FlBeta'] = np.array([
[0.0568835768038955,0,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0163789287295263,0.00486941124391324,0,0,0,0,0,0,0,0,0,0,0,0],
[0.0177069499778663,0.00376272687029659,0.00177069499778663,0,0,0,0,0,0,0,0,0,0,0],
[0.0358565737051793,0.0208056662239929,0.00309871624612661,0.00641876936697654,0.00199203187250996,0,0,0,0,0,0,0,0,0],
[0.0196989818503763,0.0398406374501992,0.00907481186365649,0.0101814962372731,0.00376272687029659,0.000885347498893316,0.000442673749446658,0.000885347498893316,0,0,0,0,0,0],
[0.00885347498893316,0.0296591412129261,0.0329791943337760,0.0234617087206729,0.00730411686586985,0.00332005312084993,0.00243470562195662,0,0.000221336874723329,0,0,0,0,0],
[0.00686144311642320,0.0119521912350598,0.0239043824701195,0.0449313855688358,0.0148295706064630,0.00420540061974325,0.00464807436918991,0.00332005312084993,0.000221336874723329,0,0,0,0,0],
[0.00332005312084993,0.00752545374059318,0.0123948649845064,0.0548915449313856,0.0305444887118194,0.00818946436476317,0.00265604249667995,0.00398406374501992,0.00309871624612661,0.00110668437361664,0.000221336874723329,0,0,0],
[0.00154935812306330,0.00531208499335989,0.00486941124391324,0.0402833111996459,0.0276671093404161,0.0208056662239929,0.00708277999114653,0.00376272687029659,0.00376272687029659,0.00221336874723329,0.000442673749446658,0,0,0],
[0.000885347498893316,0.00442673749446658,0.00464807436918991,0.0135015493581231,0.0351925630810093,0.0263390880920761,0.0148295706064630,0.00442673749446658,0.00332005312084993,0.00420540061974325,0.00177069499778663,0.00154935812306330,0,0],
[0.000664010624169987,0.00309871624612661,0.00199203187250996,0.00951748561310314,0.0230190349712262,0.0166002656042497,0.0185922974767596,0.00464807436918991,0.00287737937140328,0.00221336874723329,0.000885347498893316,0.000664010624169987,0.000885347498893316,0],
[0,0.000442673749446658,0.00177069499778663,0.00287737937140328,0.00664010624169987,0.00907481186365649,0.00619743249225321,0.00818946436476317,0.00597609561752988,0.000442673749446658,0.00287737937140328,0,0.000442673749446658,0.000664010624169987],
[0,0,0.00154935812306330,0.00177069499778663,0.00265604249667995,0.00553342186808322,0.000664010624169987,0.00420540061974325,0.00177069499778663,0.00177069499778663,0.000221336874723329,0.000221336874723329,0.000442673749446658,0.000221336874723329],
[0.000221336874723329,0.000221336874723329,0,0.000664010624169987,0.000664010624169987,0.00177069499778663,0.000885347498893316,0.000664010624169987,0.00110668437361664,0.000664010624169987,0.000442673749446658,0.000442673749446658,0.000221336874723329,0],
[0,0,0.000221336874723329,0.000442673749446658,0.000221336874723329,0,0,0.000221336874723329,0.000442673749446658,0,0,0,0,0],
[0,0,0,0,0,0,0,0,0,0,0,0,0,0]
])
stats['stab'] = [1.073, -0.998, 4.197, 84.083]
stats['Alpha_mod'] = 0
stats['Overload_b'] = 0.00099
stats['mu_beta'] = 0.25
stats['Zpr_pars'] = [-3.72, 1.73, 1.99]
stats['mu_committed'] = 0.085
else:
# Fallback if invalid ID
raise ValueError("Invalid Reference System ID. Choose 1, 2, or 3.")
return stats