#include "trs5-inc.mm"

Func void set_parameters() {

     // set_parameters()

//     systype= WithVMdelay;
     systype= FirstOrderAVR;

     switch(systype){
     case FirstOrderAVR:	comp_dim= 4; break;
     case WithVMdelay: // == default
     default:			comp_dim= 5; break;
     }
// 状態変数 x の添字
     { Ndlt1, Nw1, NeqD1, NgV1, NVvm1 } = { 1, 2, 3, 4, 5 };

     // 計算打ち切り時刻を入力
     t_fin= 20.0;   //120.0
     // 数値積分の時間刻みを入力
     h= 0.02;       //0.02
     // 初期状態を入力
     //                     dlt[rad]    w[pu]    eqD[pu]   gV[pu]   Vvm[pu]
//     x_ini(Ndlt1:NVvm1) = [1.35        0.0       1.2       2        0.5 ]';
     x_ini(Ndlt1:NVvm1) = [1.66        0.0       1.18       2.23      0.5 ]';

     // 事故を起こす場合のパラメータ ------------------------------
     // 三相短絡事故が生じる時刻を入力
     fault_time= Inf;			// 事故なしと同じ
     // 事故継続時間を入力
     fault_duration= 0.0666666665;		// 4 cycles
     // 事故除去から再閉路までの時間を入力
     reclose_interval= 0.8333333332;	// 50 cycles
     // 一回線と二回線ある内、二回線の比率 a を入力
     twoline_rate= 0.5;
     // 二回線上での三相短絡事故が生じる位置 (一回線との接続点からの) b を入力
     fault_position= 0.1;

     // 距離を計る時の基準値を入力
     // これを適切な値にしておかないと
     // bifudiag_output, 2pbifu_func での周期や分岐点の判定がおかしくなる!
     xbase(Ndlt1:NVvm1) = [2.0   0.02   0.4  20  1]';

     // ---------------- パラメータ ----------------
     wB=	120*PI;
     Xd1 =	1.03;
     Xq1 =	0.618;
     XdD1 =	0.326;
     Vinf1 =	1.0;
     M1 =	10.59;
     D1 =	3.0;
     Re1 =	0.073;
     // Xe= 0.47 を XtransとXline に分ける.        
     Xtrans1 =	0.1;
     Xline1 =	0.37;
     Vts1 =	1.05;
     // 励磁制御なしの時の ef はf()の中で直接指定してある.
     efs1 =	1.05;
     UlV1 =	3.08;
     LlV1 =	-2.01;
     // TdoD= Xf/(wB*rf)  ただしXf[Ohm], rf[Ohm]
     TdoD1 =	7.61;
     Pms1  =	1.4;
     GV1   =	108.7;
     TV1   =	0.05;
     // 
     Tvm1  =	0.15;        
}