example_path = @__DIR__; # 我们得到包含当前脚本的目录的绝对路径
cd(example_path); # 转到示例目录

model_name = "ground_directional_overcurrent_relay.engee"
engee.open(model_name); # 打开模型
model = engee.load(model_name); # 下载模型以快速访问

path = engee.gcm().name # 模型名称
# 设置短路块的参数
t_sc = 0.5 # 我们设置短路的时间并将其发送到短路块
engee.set_param!(path * "/K2", "temporal" => true,
                "type" => "Single-phase to ground (a-g)")
duration = 10 # 将持续时间参数传递给块K3
# 启用平行线。 在AUV中，短路开始的参数被设置为明显长的模拟持续时间。
engee.set_param!(path * "/AUV V3", "t_start" => 10.0)
engee.set_param!(path * "/AUV B4", "t_start" => 10.0)
results = engee.run(model);

# 我们会得到模拟结果
I0_1 = results["|3I0_1|"].value
Iabc = results["I_1"].value
t = results["|3I0_1|"].time;

# 我们从C1侧获得第一和第二保护级l1的设定点值
I_trip_first = parse(Float64, engee.get_param(path * "/C1：TNZNP1st L1", "I_trip"))
I_trip_second = parse(Float64, engee.get_param(path * "/C1：TNZNP2st L1", "I_trip"));

# 计算步骤在模型中设置为100e-6。
# 瞬态过程在大约0.07秒内衰减，确定短路电流的值
index_sc = Int(floor((t_sc + 0.071) * 10^4)) # 我们得到了确定短路电流的指标
# 或者我们以这种方式获取索引。 无论计算步骤如何，此实现都允许您获取索引。
index_sc = findfirst(t .>= t_sc + 0.071)
# 我们保持短路电流
I_sc = I0_1[index_sc];

# 我们将从第一和第二阶段获得触发状态
signal_second_stage = results["S1：TNZNP1st L1.高级行政主任"].value[index_sc] # 当短路消失时，第一阶段已经被触发。
# 第二阶段稍后触发
# 接收第二级触发信号
t_second_trip = parse(Float64, engee.get_param(path * "/C1：TNZNP2st L1",
                    "t_trip")) + t_sc + 0.02
index_second_trip = findfirst(t .>= t_second_trip)
signal_second_stage = results["S1：TNZNP2st L1。高级行政主任"].value[index_second_trip];
# 如果触发任何阶段，我们认为保护已经起作用。
 if results["S1：TNZNP1st L1.高级行政主任"].value[index_sc] || results["S1：TNZNP2st L1。高级行政主任"].value[index_second_trip]
    trip_status = true
 else 
    trip_status = false
 end;

# 让我们将短路电流与设定点进行比较
if I_sc >= I_trip_first || I_sc >= I_trip_second
    trip_correct = true
else
    trip_correct = false
end;

# 输出结果
println("""
第一阶段的设定点是$I_trip_first和
第二阶段的设定点是$I_trip_second和
短路电流值：△(round(i_sc,digits=2))A
超过当前设定值：trip trip_correct
保护激活：trip trip_status
""")

# 绘制零序电流
# 我们构建主电流图
plot(
    t, I0_1,
    label="3I0_1",
    xlabel="t, c",
    ylabel="我，和",
    title="C1侧的电流曲线图3I0",
    color=:blue
)

# 我们添加水平线以图形方式显示当前设置。
hline!([I_trip_first], label="第一阶段", linestyle=:dash, color=:red)
hline!([I_trip_second], label="第二阶段", linestyle=:dash, color=:green)

I = reduce(hcat, Iabc)' # Транспонируем, чтобы получилась матрица (N×3)

# 构建图表
plot(t, I, lw=2, label=["A阶段" "B阶段" "C期"], color=[:orange :green :red])
title!("C1侧相电流曲线图")
xlabel!("时间，c")
ylabel!("电流，A")