diff --git a/code/steps/STEPS.cbp b/code/steps/STEPS.cbp
index 5fd914d5..670566a3 100644
--- a/code/steps/STEPS.cbp
+++ b/code/steps/STEPS.cbp
@@ -8,7 +8,7 @@
 		<Option compiler="gcc" />
 		<Build>
 			<Target title="release">
-				<Option output="../../build/bin/release/STEPS" imp_lib="$(TARGET_OUTPUT_DIR)$(TARGET_OUTPUT_BASENAME).a" def_file="$(TARGET_OUTPUT_DIR)$(TARGET_OUTPUT_BASENAME).def" prefix_auto="1" extension_auto="1" />
+				<Option output="../../build/bin/release/libSTEPS" imp_lib="$(TARGET_OUTPUT_DIR)$(TARGET_OUTPUT_BASENAME).a" def_file="$(TARGET_OUTPUT_DIR)$(TARGET_OUTPUT_BASENAME).def" prefix_auto="1" extension_auto="1" />
 				<Option object_output="../../build/obj/release/" />
 				<Option type="3" />
 				<Option compiler="gcc" />
diff --git a/code/steps/STEPS.depend b/code/steps/STEPS.depend
index 920c24a7..cda3b6ab 100644
--- a/code/steps/STEPS.depend
+++ b/code/steps/STEPS.depend
@@ -11,7 +11,7 @@
 	"header/prepare_for_tests/prepare_models_for_test.h"
 	"header/data_imexporter/equivalent_model_imexporter.h"
 
-1645580274 d:\steps\code\steps\header\power_system_database.h
+1649395482 d:\steps\code\steps\header\power_system_database.h
 	"header/device/bus.h"
 	"header/device/generator.h"
 	"header/device/wt_generator.h"
@@ -58,7 +58,7 @@
 	<vector>
 	<map>
 
-1645580274 d:\steps\code\steps\header\basic\device_id.h
+1649395458 d:\steps\code\steps\header\basic\device_id.h
 	"header/basic/constants.h"
 	"header/basic/steps_enum.h"
 	<vector>
@@ -66,7 +66,7 @@
 	<functional>
 	"terminal.h"
 
-1598590744 d:\steps\code\steps\header\basic\constants.h
+1649395458 d:\steps\code\steps\header\basic\constants.h
 
 1625794452 d:\steps\code\steps\header\basic\terminal.h
 	<vector>
@@ -76,7 +76,7 @@
 	"header/basic/steps_enum.h"
 	<vector>
 
-1645580274 d:\steps\code\steps\header\basic\steps_enum.h
+1649395482 d:\steps\code\steps\header\basic\steps_enum.h
 	<iso646.h>
 
 1645580274 d:\steps\code\steps\header\model\model.h
@@ -364,13 +364,13 @@
 1630827861 d:\steps\code\steps\header\basic\bus_index.h
 	<vector>
 
-1585314815 d:\steps\code\steps\header\basic\device_index_map.h
+1649395482 d:\steps\code\steps\header\basic\device_index_map.h
 	"header/basic/device_id.h"
 	<functional>
 	<map>
 	<unordered_map>
 
-1645580274 d:\steps\code\steps\header\steps_namespace.h
+1649395458 d:\steps\code\steps\header\steps_namespace.h
 	"header/STEPS.h"
 	"header/basic/steps_define.h"
 	"header/basic/constants.h"
@@ -463,7 +463,7 @@
 	"header/basic/sparse_matrix_define.h"
 	<ctime>
 
-1645580274 d:\steps\code\steps\header\toolkit\dynamic_simulator\dynamic_simulator.h
+1649400120 d:\steps\code\steps\header\toolkit\dynamic_simulator\dynamic_simulator.h
 	"header/basic/power_mismatch_struct.h"
 	"header/meter/meter.h"
 	"header/network/network_matrix.h"
@@ -472,20 +472,20 @@
 
 1583938871 d:\steps\code\steps\header\basic\power_mismatch_struct.h
 
-1583938871 d:\steps\code\steps\header\apis\steps_api_search_buffer.h
+1649395482 d:\steps\code\steps\header\apis\steps_api_search_buffer.h
 	"header/power_system_database.h"
 	"header/basic/constants.h"
 	<ctime>
 	<string>
 
-1645580274 d:\steps\code\steps\header\basic\utility.h
+1649307618 d:\steps\code\steps\header\basic\utility.h
 	"header/power_system_database.h"
 	"header/network/network_matrix.h"
 	"header/toolkit/powerflow_solver/powerflow_solver.h"
 	"header/toolkit/dynamic_simulator/dynamic_simulator.h"
 	<complex>
 
-1645580274 d:\steps\code\steps\header\data_imexporter\psse_imexporter.h
+1649395482 d:\steps\code\steps\header\data_imexporter\psse_imexporter.h
 	"header/data_imexporter/data_imexporter.h"
 	"header/device/source.h"
 	"header/device/transformer.h"
@@ -522,7 +522,7 @@
 	"header/toolkit/dynamic_simulator/dynamic_simulator.h"
 	<fstream>
 
-1645584219 source:d:\steps\code\steps\main_dynamics_simulation.cpp
+1649395482 source:d:\steps\code\steps\main_dynamics_simulation.cpp
 	<istream>
 	<iostream>
 	"header/power_system_database.h"
@@ -552,7 +552,7 @@
 1645580274 d:\steps\code\steps\header\meter\meter_setter.h
 	"header/meter/meter.h"
 
-1645580274 source:d:\steps\code\steps\main_powerflow.cpp
+1649395482 source:d:\steps\code\steps\main_powerflow.cpp
 	<istream>
 	<iostream>
 	"header/power_system_database.h"
@@ -1030,7 +1030,7 @@
 	"header/data_imexporter/steps_imexporter.h"
 	"header/STEPS.h"
 
-1645669255 d:\steps\code\steps\header\data_imexporter\steps_imexporter.h
+1649395482 d:\steps\code\steps\header\data_imexporter\steps_imexporter.h
 	"header/data_imexporter/data_imexporter.h"
 	"header/device/source.h"
 	"header/device/transformer.h"
@@ -2342,14 +2342,14 @@
 	"header/basic/continuous_buffer.h"
 	<complex>
 
-1583938871 source:d:\steps\code\steps\source\apis\steps_api_add_devices.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_add_devices.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
 	"header/data_imexporter/psse_imexporter.h"
 	"header/data_imexporter/bpa_imexporter.h"
 
-1645580274 d:\steps\code\steps\header\apis\steps_api.h
+1649403749 d:\steps\code\steps\header\apis\steps_api.h
 	<cstdlib>
 	"header/basic/constants.h"
 	"header/basic/dc_device_id.h"
@@ -2372,7 +2372,7 @@
 	"header/toolkit/powerflow_solver/powerflow_solver.h"
 	"header/toolkit/dynamic_simulator/dynamic_simulator.h"
 
-1634637238 source:d:\steps\code\steps\source\apis\steps_api_data_import_export.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_data_import_export.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -2380,14 +2380,14 @@
 	"header/data_imexporter/bpa_imexporter.h"
 	"header/data_imexporter/steps_imexporter.h"
 
-1645580274 source:d:\steps\code\steps\source\apis\steps_api_dynamic_simulation.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_dynamic_simulation.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/toolkit/cct_searcher/cct_searcher.h"
 	"header/steps_namespace.h"
 	<iostream>
 
-1583938871 source:d:\steps\code\steps\source\apis\steps_api_get_device_count.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_get_device_count.cpp
 	"header/apis/steps_api.h"
 	"header/apis/steps_api_common.h"
 	"header/basic/utility.h"
@@ -2395,7 +2395,7 @@
 	"header/data_imexporter/psse_imexporter.h"
 	"header/data_imexporter/bpa_imexporter.h"
 
-1583938871 source:d:\steps\code\steps\source\apis\steps_api_is_device_exist.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_is_device_exist.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -2406,7 +2406,7 @@
 	"header/steps_namespace.h"
 	<iostream>
 
-1583938871 source:d:\steps\code\steps\source\apis\steps_api_remove_devices.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_remove_devices.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -2417,7 +2417,7 @@
 	"header/apis/steps_api_search_buffer.h"
 	<iostream>
 
-1611548671 source:d:\steps\code\steps\source\apis\steps_api_search_devices.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_search_devices.cpp
 	<istream>
 	"header/apis/steps_api.h"
 	"header/apis/steps_api_common.h"
@@ -2559,7 +2559,7 @@
 	"header/data_imexporter/psse_imexporter.h"
 	"header/data_imexporter/bpa_imexporter.h"
 
-1583938871 source:d:\steps\code\steps\source\apis\steps_api_show_device_data.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_show_device_data.cpp
 	"header/apis/steps_api.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -2611,7 +2611,7 @@
 	<iostream>
 	<cstdio>
 
-1629270179 source:d:\steps\code\steps\source\apis\steps_api_utilities.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_utilities.cpp
 	"header/apis/steps_api.h"
 	"header/apis/steps_api_common.h"
 	"header/basic/utility.h"
@@ -2682,7 +2682,7 @@
 	<cstdio>
 	<cmath>
 
-1645580274 source:d:\steps\code\steps\source\basic\device_id.cpp
+1649395458 source:d:\steps\code\steps\source\basic\device_id.cpp
 	"header/basic/device_id.h"
 	"header/basic/exception.h"
 	"header/basic/utility.h"
@@ -2704,7 +2704,7 @@
 	<iostream>
 	<cstdio>
 
-1585314815 source:d:\steps\code\steps\source\basic\device_index_map.cpp
+1649395482 source:d:\steps\code\steps\source\basic\device_index_map.cpp
 	"header/basic/device_index_map.h"
 	"header/basic/utility.h"
 	"header/basic/constants.h"
@@ -3166,7 +3166,7 @@
 	<iostream>
 	<cstdio>
 
-1645580274 source:d:\steps\code\steps\source\basic\utility.cpp
+1649307624 source:d:\steps\code\steps\source\basic\utility.cpp
 	"header/basic/utility.h"
 	"header/basic/constants.h"
 	"header/steps_namespace.h"
@@ -3394,7 +3394,7 @@
 	<iostream>
 	<fstream>
 
-1634637238 source:d:\steps\code\steps\source\data_imexporter\bpa_imexporter_test.cpp
+1649395482 source:d:\steps\code\steps\source\data_imexporter\bpa_imexporter_test.cpp
 	<cmath>
 	"cpptest.h"
 	"header/basic/test_macro.h"
@@ -3480,7 +3480,7 @@
 	<iostream>
 	<cstdio>
 
-1645580274 source:d:\steps\code\steps\source\data_imexporter\psse_powerflow_imexporter.cpp
+1649395482 source:d:\steps\code\steps\source\data_imexporter\psse_powerflow_imexporter.cpp
 	"header/data_imexporter/psse_imexporter.h"
 	"header/data_imexporter/steps_imexporter.h"
 	"header/basic/utility.h"
@@ -3497,7 +3497,7 @@
 	<iostream>
 	<algorithm>
 
-1645669468 source:d:\steps\code\steps\source\data_imexporter\steps_dynamics_imexporter.cpp
+1649395482 source:d:\steps\code\steps\source\data_imexporter\steps_dynamics_imexporter.cpp
 	"header/data_imexporter/steps_imexporter.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -3805,7 +3805,7 @@
 	<iostream>
 	<cstdio>
 
-1645669416 source:d:\steps\code\steps\source\device\vsc_hvdc.cpp
+1649395482 source:d:\steps\code\steps\source\device\vsc_hvdc.cpp
 	"header/device/vsc_hvdc.h"
 	"header/basic/utility.h"
 	"header/basic/steps_define.h"
@@ -3817,7 +3817,7 @@
 	<cstdio>
 	<cmath>
 
-1645542341 d:\steps\code\steps\header\device\vsc_hvdc.h
+1649395482 d:\steps\code\steps\header\device\vsc_hvdc.h
 	"header/device/nonbus_device.h"
 	"header/device/vsc_hvdc_var.h"
 	"header/basic/device_id.h"
@@ -5657,7 +5657,7 @@
 	"header/steps_namespace.h"
 	"header/basic/utility.h"
 
-1645580274 source:d:\steps\code\steps\source\steps.cpp
+1649395458 source:d:\steps\code\steps\source\steps.cpp
 	"header/STEPS.h"
 	"header/basic/utility.h"
 	"header/steps_namespace.h"
@@ -5666,7 +5666,7 @@
 	<chrono>
 	<thread>
 
-1645580274 source:d:\steps\code\steps\source\steps_namespace.cpp
+1649307607 source:d:\steps\code\steps\source\steps_namespace.cpp
 	"header/steps_namespace.h"
 
 1634637239 source:d:\steps\code\steps\source\toolkit\cct_searcher\cct_searcher.cpp
@@ -5740,7 +5740,7 @@
 	"header/toolkit/contingency_screener/contingency_screener.h"
 	"header/STEPS.h"
 
-1645543051 source:d:\steps\code\steps\source\toolkit\dynamic_simulator\dynamic_simulator.cpp
+1649400112 source:d:\steps\code\steps\source\toolkit\dynamic_simulator\dynamic_simulator.cpp
 	"header/toolkit/dynamic_simulator/dynamic_simulator.h"
 	"header/basic/constants.h"
 	"header/basic/utility.h"
@@ -5799,7 +5799,7 @@
 	<cstdio>
 	<cmath>
 
-1645580274 source:d:\steps\code\steps\source\toolkit\powerflow_solver\powerflow_solver.cpp
+1646044699 source:d:\steps\code\steps\source\toolkit\powerflow_solver\powerflow_solver.cpp
 	"header/toolkit/powerflow_solver/powerflow_solver.h"
 	"header/basic/utility.h"
 	"header/basic/constants.h"
@@ -6075,7 +6075,7 @@
 	<map>
 	<functional>
 
-1634637238 source:d:\steps\code\steps\source\basic\string_int_map.cpp
+1649395458 source:d:\steps\code\steps\source\basic\string_int_map.cpp
 	"header/basic/string_int_map.h"
 	"header/basic/utility.h"
 	"header/basic/constants.h"
@@ -12254,7 +12254,7 @@
 	<istream>
 	<iostream>
 
-1645580274 source:d:\steps\code\steps\source\toolkit\dynamic_simulator\dynamic_simulator_prepare_meter.cpp
+1649395482 source:d:\steps\code\steps\source\toolkit\dynamic_simulator\dynamic_simulator_prepare_meter.cpp
 	"header/toolkit/dynamic_simulator/dynamic_simulator.h"
 	"header/basic/constants.h"
 	"header/basic/utility.h"
@@ -12336,7 +12336,7 @@
 	<ctime>
 	<chrono>
 
-1645669819 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsc_hvdc_converter_model.cpp
+1645673172 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsc_hvdc_converter_model.cpp
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/vsc_hvdc_converter_model.h"
 	"header/device/vsc_hvdc.h"
 	"header/basic/utility.h"
@@ -12355,11 +12355,12 @@
 1643805345 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_project_model\vsc_hvdc_project_model.h
 	"header/model/vsc_hvdc_model/vsc_hvdc_model.h"
 
-1645581445 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsc_hvdc_converter_models.h
+1649395482 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsc_hvdc_converter_models.h
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSCHVDCC0.h"
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSG0.h"
+	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSG1.h"
 
-1645546130 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vschvdcc0.h
+1649395482 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vschvdcc0.h
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/vsc_hvdc_converter_model.h"
 	"header/block/pid_block.h"
 	"header/block/first_order_block.h"
@@ -12367,7 +12368,7 @@
 1643800465 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_project_model\vsc_hvdc_project_model.cpp
 	"header/model/vsc_hvdc_model/vsc_hvdc_project_model/vsc_hvdc_project_model.h"
 
-1645546092 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vschvdcc0.cpp
+1649395482 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vschvdcc0.cpp
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSCHVDCC0.h"
 	"header/basic/utility.h"
 	"header/STEPS.h"
@@ -12378,7 +12379,7 @@
 	"header/basic/utility.h"
 	"header/STEPS.h"
 
-1645580274 source:d:\steps\code\steps\source\apis\steps_api_set_get_vsc_hvdc_related_model_parameters.cpp
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_set_get_vsc_hvdc_related_model_parameters.cpp
 	"header/apis/steps_api.h"
 	"header/apis/steps_api_common.h"
 	"header/basic/utility.h"
@@ -12481,7 +12482,7 @@
 
 1647153997 e:\users\zma\desktop\steps\code\steps\header\basic\steps_define.h
 
-1645582171 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg0.h
+1649395482 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg0.h
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/vsc_hvdc_converter_model.h"
 	"header/block/integral_block.h"
 	"header/block/pid_block.h"
@@ -12489,7 +12490,7 @@
 	"header/block/lead_lag_block.h"
 	"header/block/first_order_block.h"
 
-1645670114 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg0.cpp
+1649395482 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg0.cpp
 	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSG0.h"
 	"header/basic/utility.h"
 	"header/STEPS.h"
@@ -12531,3 +12532,25 @@
 	"header/data_imexporter/psse_imexporter.h"
 	"header/data_imexporter/bpa_imexporter.h"
 
+1649395482 d:\steps\code\steps\header\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg1.h
+	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/vsc_hvdc_converter_model.h"
+	"header/block/integral_block.h"
+	"header/block/pid_block.h"
+	"header/block/pi_block.h"
+	"header/block/lead_lag_block.h"
+	"header/block/first_order_block.h"
+
+1649395482 source:d:\steps\code\steps\source\apis\steps_api_set_get_vsc_hvdc_parameters.cpp
+	"header/apis/steps_api.h"
+	"header/apis/steps_api_common.h"
+	"header/basic/utility.h"
+	"header/steps_namespace.h"
+	"header/data_imexporter/psse_imexporter.h"
+	"header/data_imexporter/bpa_imexporter.h"
+
+1649395482 source:d:\steps\code\steps\source\model\vsc_hvdc_model\vsc_hvdc_converter_model\vsg1.cpp
+	"header/model/vsc_hvdc_model/vsc_hvdc_converter_model/VSG1.h"
+	"header/basic/utility.h"
+	"header/STEPS.h"
+	<iostream>
+
diff --git a/code/steps/header/apis/steps_api.h b/code/steps/header/apis/steps_api.h
index ecc813aa..97693898 100644
--- a/code/steps/header/apis/steps_api.h
+++ b/code/steps/header/apis/steps_api.h
@@ -56,10 +56,10 @@ EXPORT_STEPS_DLL void api_set_owner_capacity(unsigned int cap, unsigned int tool
 
 
 EXPORT_STEPS_DLL void api_load_powerflow_data_from_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
-EXPORT_STEPS_DLL void api_load_vsc_powerflow_data_from_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
+EXPORT_STEPS_DLL void api_load_vsc_hvdc_powerflow_data_from_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
 EXPORT_STEPS_DLL void api_load_powerflow_result_from_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
 EXPORT_STEPS_DLL void api_save_powerflow_data_to_file(char* file, char* file_type, bool export_zero_impedance_line=true, bool export_out_of_service_bus=true, bool export_internal_bus_number=false, unsigned int powerflow_data_save_mode=0, unsigned int toolkit_index=INDEX_NOT_EXIST);
-EXPORT_STEPS_DLL void api_save_vsc_powerflow_data_to_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
+EXPORT_STEPS_DLL void api_save_vsc_hvdc_powerflow_data_to_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
 EXPORT_STEPS_DLL void api_load_dynamic_data_from_file(char* file, char* file_type, unsigned int toolkit_index=INDEX_NOT_EXIST);
 EXPORT_STEPS_DLL void api_save_dynamic_data_to_file(char* file, char* file_type, bool export_internal_bus_number=false, unsigned int toolkit_index=INDEX_NOT_EXIST);
 
diff --git a/code/steps/header/toolkit/dynamic_simulator/dynamic_simulator.h b/code/steps/header/toolkit/dynamic_simulator/dynamic_simulator.h
index c8b0306a..48325251 100644
--- a/code/steps/header/toolkit/dynamic_simulator/dynamic_simulator.h
+++ b/code/steps/header/toolkit/dynamic_simulator/dynamic_simulator.h
@@ -182,7 +182,7 @@ class DYNAMICS_SIMULATOR
         double get_hvdc_power_order_in_MW(const DEVICE_ID& hvdc_id);
         void change_hvdc_power_order_in_MW(const DEVICE_ID& hvdc_id, double porder);
 
-        void set_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID line_did, unsigned int side_bus, double location, double fault_shunt);
+        void set_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID line_did, unsigned int side_bus, double location, double fault_r);
         void clear_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID line_did);
 
         void trip_vsc_hvdc_line(string vsc_name, DC_DEVICE_ID line_did);
diff --git a/code/steps/main_dynamics_simulation.cpp b/code/steps/main_dynamics_simulation.cpp
index a92c1e62..86c58c88 100644
--- a/code/steps/main_dynamics_simulation.cpp
+++ b/code/steps/main_dynamics_simulation.cpp
@@ -50,7 +50,7 @@ int main()
 
     //importer.load_dynamic_data("../../../bench/shandong_load/5_bus.dyr");
     importer.load_dynamic_data("../../../bench/39_bus/IEEE39.dyr");
-    STEPS_SHOW_FILE_FUNCTION_AND_LINE_INFO
+
     //vector<GENERATOR*> gens = psdb.get_all_generators();
     //unsigned int n = gens.size();
     //for(unsigned int i=0; i!=n; ++i)
diff --git a/code/steps/source/apis/steps_api_data_import_export.cpp b/code/steps/source/apis/steps_api_data_import_export.cpp
index e12c6d12..96b06bc3 100644
--- a/code/steps/source/apis/steps_api_data_import_export.cpp
+++ b/code/steps/source/apis/steps_api_data_import_export.cpp
@@ -24,7 +24,7 @@ void api_load_powerflow_data_from_file(char* file, char* file_type, unsigned int
     }
 }
 
-void api_load_vsc_powerflow_data_from_file(char* file, char* file_type, unsigned int toolkit_index)
+void api_load_vsc_hvdc_powerflow_data_from_file(char* file, char* file_type, unsigned int toolkit_index)
 {
     STEPS& toolkit = get_toolkit(toolkit_index);
     string string_file_type = string2upper(file_type);
@@ -107,7 +107,7 @@ void api_save_powerflow_data_to_file(char* file, char* file_type,
     }
 }
 
-void api_save_vsc_powerflow_data_to_file(char* file, char* file_type, unsigned int toolkit_index)
+void api_save_vsc_hvdc_powerflow_data_to_file(char* file, char* file_type, unsigned int toolkit_index)
 {
     STEPS& toolkit = get_toolkit(toolkit_index);
     string string_file_type = string2upper(file_type);
diff --git a/code/steps/source/device/vsc_hvdc.cpp b/code/steps/source/device/vsc_hvdc.cpp
index 16e3ff3a..56a0da5a 100644
--- a/code/steps/source/device/vsc_hvdc.cpp
+++ b/code/steps/source/device/vsc_hvdc.cpp
@@ -2239,7 +2239,6 @@ void VSC_HVDC::set_model(MODEL* model)
 {
     if(model != NULL and model->has_allowed_device_type(STEPS_VSC_HVDC))
     {
-        STEPS_SHOW_FILE_FUNCTION_AND_LINE_INFO
         model->set_device_id(get_device_id());
         if(model->get_model_type()=="VSC HVDC NETWORK")
         {
diff --git a/code/steps/source/power_system_database.cpp b/code/steps/source/power_system_database.cpp
index 8992c277..90e49195 100644
--- a/code/steps/source/power_system_database.cpp
+++ b/code/steps/source/power_system_database.cpp
@@ -65,7 +65,7 @@ void POWER_SYSTEM_DATABASE::set_database_capacity()
     unsigned int line_capacity = (unsigned int)(round(bus_capacity*3.0));
     unsigned int transformer_capacity = bus_capacity;
     unsigned int hvdc_capacity = (unsigned int)(round(bus_capacity*0.1));
-    unsigned int vsc_hvdc_capacity = (unsigned int)(round(bus_capacity*0.1));
+    unsigned int vsc_hvdc_capacity = (unsigned int)(round(bus_capacity*0.01));
     unsigned int equivalentdevice_capacity = (unsigned int)(round(bus_capacity*0.1));
     unsigned int energy_storage_capacity = (unsigned int)(round(bus_capacity*0.1));
     unsigned int lcc_hvdc_capacity = hvdc_capacity;
@@ -111,7 +111,7 @@ void POWER_SYSTEM_DATABASE::set_database_capacity()
             hvdc_capacity = (unsigned int)(document["hvdc"].GetInt());
 
         if(document.HasMember("vschvdc") and document["vschvdc"].IsInt())
-        hvdc_capacity = (unsigned int)(document["vschvdc"].GetInt());
+            vsc_hvdc_capacity = (unsigned int)(document["vschvdc"].GetInt());
 
         if(document.HasMember("equivalentDevice") and document["equivalentDevice"].IsInt())
             equivalentdevice_capacity = (unsigned int)(document["equivalentDevice"].GetInt());
diff --git a/code/steps/source/toolkit/dynamic_simulator/dynamic_simulator.cpp b/code/steps/source/toolkit/dynamic_simulator/dynamic_simulator.cpp
index 3d1be164..6e3a0abf 100644
--- a/code/steps/source/toolkit/dynamic_simulator/dynamic_simulator.cpp
+++ b/code/steps/source/toolkit/dynamic_simulator/dynamic_simulator.cpp
@@ -4264,7 +4264,7 @@ void DYNAMICS_SIMULATOR::change_hvdc_power_order_in_MW(const DEVICE_ID& hvdc_id,
 }
 
 
-void DYNAMICS_SIMULATOR::set_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID line_did, unsigned int side_bus, double location, double fault_shunt)
+void DYNAMICS_SIMULATOR::set_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID line_did, unsigned int side_bus, double location, double fault_r)
 {
     POWER_SYSTEM_DATABASE& psdb = toolkit->get_power_system_database();
     VSC_HVDC* vsc_hvdc = psdb.get_vsc_hvdc(vsc_name);
@@ -4272,7 +4272,7 @@ void DYNAMICS_SIMULATOR::set_vsc_hvdc_line_fault(string vsc_name, DC_DEVICE_ID l
     if(line_index!=INDEX_NOT_EXIST)
     {
         vsc_hvdc->set_dc_line_fault_location(line_index, side_bus, location);
-        vsc_hvdc->set_dc_line_fault_r_in_ohm(line_index, fault_shunt);
+        vsc_hvdc->set_dc_line_fault_r_in_ohm(line_index, fault_r);
     }
 }
 
diff --git a/code/steps/source/toolkit/powerflow_solver/powerflow_solver.cpp b/code/steps/source/toolkit/powerflow_solver/powerflow_solver.cpp
index 338d600b..541f40e5 100644
--- a/code/steps/source/toolkit/powerflow_solver/powerflow_solver.cpp
+++ b/code/steps/source/toolkit/powerflow_solver/powerflow_solver.cpp
@@ -1646,18 +1646,14 @@ void POWERFLOW_SOLVER::check_SLACK_bus_constraint_of_physical_bus(unsigned int p
                 {
                     if(vsc_hvdcs[i]->get_converter_active_power_operation_mode(index)==VSC_AC_VOLTAGE_ANGLE_CONTROL)
                     {
-                        //STEPS_SHOW_FILE_FUNCTION_AND_LINE_INFO
                         double P_loading_in_MW = vsc_hvdcs[i]->get_converter_Pmax_in_MW(index) - vsc_hvdcs[i]->get_converter_Pmin_in_MW(index);
                         P_loading_in_MW = P_loading_in_MW*P_loading_percentage + vsc_hvdcs[i]->get_converter_Pmin_in_MW(index);
-                        //cout<<"P_loading_in_MW: "<<P_loading_in_MW<<endl;
                         vsc_hvdcs[i]->set_converter_P_to_AC_bus_in_MW(index, P_loading_in_MW);
                     }
                     if(vsc_hvdcs[i]->get_converter_reactive_power_operation_mode(index)==VSC_AC_VOLTAGE_CONTROL)
                     {
-                        //STEPS_SHOW_FILE_FUNCTION_AND_LINE_INFO
                         double Q_loading_in_MVar = vsc_hvdcs[i]->get_converter_Qmax_in_MVar(index) - vsc_hvdcs[i]->get_converter_Qmin_in_MVar(index);
                         Q_loading_in_MVar = Q_loading_in_MVar*Q_loading_percentage + vsc_hvdcs[i]->get_converter_Qmin_in_MVar(index);
-                        //cout<<"Q_loading_in_MVar: "<<Q_loading_in_MVar<<endl;
                         vsc_hvdcs[i]->set_converter_Q_to_AC_bus_in_MVar(index, Q_loading_in_MVar);
                     }
                 }
diff --git a/python/stepspy-current/stepspy-dev/libsteps/pylibsteps.py b/python/stepspy-current/stepspy-dev/libsteps/pylibsteps.py
index a3e1587e..170cad2c 100644
--- a/python/stepspy-current/stepspy-dev/libsteps/pylibsteps.py
+++ b/python/stepspy-current/stepspy-dev/libsteps/pylibsteps.py
@@ -109,6 +109,8 @@ def load_library():
 
     libsteps.api_load_powerflow_data_from_file.restype = None
     libsteps.api_load_powerflow_data_from_file.argtypes = (c_char_p, c_char_p, c_uint)
+    libsteps.api_load_vsc_hvdc_powerflow_data_from_file.restype = None
+    libsteps.api_load_vsc_hvdc_powerflow_data_from_file.argtypes = (c_char_p, c_char_p, c_uint)   
     libsteps.api_load_powerflow_result_from_file.restype = None
     libsteps.api_load_powerflow_result_from_file.argtypes = (c_char_p, c_char_p, c_uint)
     libsteps.api_load_dynamic_data_from_file.restype = None
@@ -116,6 +118,8 @@ def load_library():
 
     libsteps.api_save_powerflow_data_to_file.restype = None
     libsteps.api_save_powerflow_data_to_file.argtypes = (c_char_p, c_char_p, c_bool, c_bool, c_bool, c_uint, c_uint)
+    libsteps.api_save_vsc_hvdc_powerflow_data_to_file.restype = None
+    libsteps.api_save_vsc_hvdc_powerflow_data_to_file.argtypes = (c_char_p, c_char_p, c_bool, c_bool, c_bool, c_uint, c_uint)
     libsteps.api_save_dynamic_data_to_file.restype = None
     libsteps.api_save_dynamic_data_to_file.argtypes = (c_char_p, c_char_p, c_bool, c_uint)
     
@@ -195,6 +199,8 @@ def load_library():
     libsteps.api_is_transformer_exist.argtypes = (c_uint, c_uint, c_uint, c_char_p, c_uint)
     libsteps.api_is_hvdc_exist.restype = c_bool
     libsteps.api_is_hvdc_exist.argtypes = (c_uint, c_uint, c_char_p, c_uint)
+    libsteps.api_is_vsc_hvdc_exist.restype = c_bool
+    libsteps.api_is_vsc_hvdc_exist.argtypes = (c_char_p, c_uint)    
     libsteps.api_is_equivalent_device_exist.restype = c_bool
     libsteps.api_is_equivalent_device_exist.argtypes = (c_uint, c_char_p, c_uint)
 
@@ -274,6 +280,8 @@ def load_library():
     libsteps.api_add_line.argtypes = (c_uint, c_uint, c_char_p, c_uint)
     libsteps.api_add_hvdc.restype = None
     libsteps.api_add_hvdc.argtypes = (c_uint, c_uint, c_char_p, c_uint)
+    libsteps.api_add_vsc_hvdc.restype = None
+    libsteps.api_add_vsc_hvdc.argtypes = (c_char_p, c_uint)   
     libsteps.api_add_transformer.restype = None
     libsteps.api_add_transformer.argtypes = (c_uint, c_uint, c_uint, c_char_p, c_uint)
     libsteps.api_add_equivalent_device.restype = None
@@ -303,6 +311,8 @@ def load_library():
     libsteps.api_remove_line.argtypes = (c_uint, c_uint, c_char_p, c_uint)
     libsteps.api_remove_hvdc.restype = None
     libsteps.api_remove_hvdc.argtypes = (c_uint, c_uint, c_char_p, c_uint)
+    libsteps.api_remove_vsc_hvdc.restype = None
+    libsteps.api_remove_vsc_hvdc.argtypes = (c_char_p, c_uint)  
     libsteps.api_remove_transformer.restype = None
     libsteps.api_remove_transformer.argtypes = (c_uint, c_uint, c_uint, c_char_p, c_uint)
     libsteps.api_remove_equivalent_device.restype = None
@@ -436,6 +446,23 @@ def load_library():
     libsteps.api_set_hvdc_boolean_data.restype = None
     libsteps.api_set_hvdc_boolean_data.argtypes = (c_uint, c_uint, c_char_p, c_char_p, c_char_p, c_bool, c_uint)
 
+    libsteps.api_get_vsc_hvdc_integer_data.restype = (c_int)
+    libsteps.api_get_vsc_hvdc_integer_data.argtypes = (c_char_p, c_char_p, c_char_p, c_uint, c_uint)
+    libsteps.api_get_vsc_hvdc_float_data.restype = (c_double)
+    libsteps.api_get_vsc_hvdc_float_data.argtypes = (c_char_p, c_char_p, c_char_p, c_uint, c_uint)
+    libsteps.api_get_vsc_hvdc_string_data.restype = (c_char_p)
+    libsteps.api_get_vsc_hvdc_string_data.argtypes = (c_char_p, c_char_p, c_char_p, c_uint, c_uint)
+    libsteps.api_get_vsc_hvdc_boolean_data.restype = (c_bool)
+    libsteps.api_get_vsc_hvdc_boolean_data.argtypes = (c_char_p, c_char_p, c_char_p, c_uint, c_uint)
+    libsteps.api_set_vsc_hvdc_integer_data.restype = None
+    libsteps.api_set_vsc_hvdc_integer_data.argtypes = (c_char_p, c_char_p, c_char_p, c_int, c_uint, c_uint)
+    libsteps.api_set_vsc_hvdc_float_data.restype = None
+    libsteps.api_set_vsc_hvdc_float_data.argtypes = (c_char_p, c_char_p, c_char_p, c_double, c_uint, c_uint)
+    libsteps.api_set_vsc_hvdc_string_data.restype = None
+    libsteps.api_set_vsc_hvdc_string_data.argtypes = (c_char_p, c_char_p, c_char_p, c_char_p, c_uint, c_uint)
+    libsteps.api_set_vsc_hvdc_boolean_data.restype = None
+    libsteps.api_set_vsc_hvdc_boolean_data.argtypes = (c_char_p, c_char_p, c_char_p, c_bool, c_uint, c_uint)
+
     libsteps.api_get_area_integer_data.restype = (c_int)
     libsteps.api_get_area_integer_data.argtypes = (c_uint, c_char_p, c_uint)
     libsteps.api_get_area_float_data.restype = (c_double)
@@ -561,7 +588,6 @@ def load_library():
     libsteps.api_get_vsc_hvdc_related_model_float_parameter_name.restype = c_char_p
     libsteps.api_get_vsc_hvdc_related_model_float_parameter_name.argtypes = (c_char_p, c_char_p, c_uint, c_uint, c_uint)
 
-
     libsteps.api_get_powerflow_solver_integer_parameter.restype = (c_uint)
     libsteps.api_get_powerflow_solver_integer_parameter.argtypes = (c_char_p, c_uint)
     libsteps.api_get_powerflow_solver_float_parameter.restype = (c_double)
@@ -665,6 +691,8 @@ def load_library():
     libsteps.api_prepare_transformer_related_meters.argtypes = (c_uint, )
     libsteps.api_prepare_hvdc_related_meters.restype = None
     libsteps.api_prepare_hvdc_related_meters.argtypes = (c_uint, )
+    libsteps.api_prepare_vsc_hvdc_related_meters.restype = None
+    libsteps.api_prepare_vsc_hvdc_related_meters.argtypes = (c_uint, )   
     libsteps.api_prepare_equivalent_device_related_meters.restype = None
     libsteps.api_prepare_equivalent_device_related_meters.argtypes = (c_uint, )
 
@@ -687,6 +715,8 @@ def load_library():
     libsteps.api_prepare_transformer_related_meter.argtypes = (c_uint, c_uint, c_uint, c_char_p, c_char_p, c_char_p, c_char_p, c_uint)
     libsteps.api_prepare_hvdc_related_meter.restype = None
     libsteps.api_prepare_hvdc_related_meter.argtypes = (c_uint, c_uint, c_char_p, c_char_p, c_char_p, c_char_p, c_uint)
+    libsteps.api_prepare_vsc_hvdc_related_meter.restype = None
+    libsteps.api_prepare_vsc_hvdc_related_meter.argtypes = (c_char_p, c_char_p, c_uint,c_char_p, c_uint)    
     libsteps.api_prepare_equivalent_device_related_meter.restype = None
     libsteps.api_prepare_equivalent_device_related_meter.argtypes = (c_uint, c_char_p, c_char_p, c_char_p, c_uint)
 
@@ -803,6 +833,15 @@ def load_library():
     libsteps.api_set_hvdc_power_order_in_MW.restype = None
     libsteps.api_set_hvdc_power_order_in_MW.argtypes = (c_uint, c_uint, c_char_p, c_double, c_uint)
 
+    libsteps.api_set_vsc_hvdc_dc_line_fault.restype = None
+    libsteps.api_set_vsc_hvdc_dc_line_fault.argtypes = (c_char_p, c_uint, c_uint, c_char_p, c_double,c_double, c_uint)
+    libsteps.api_clear_vsc_hvdc_dc_line_fault.restype = None
+    libsteps.api_clear_vsc_hvdc_dc_line_fault.argtypes = (c_char_p, c_uint, c_uint, c_char_p, c_uint)
+    libsteps.api_trip_vsc_hvdc_dc_line.restype = None
+    libsteps.api_trip_vsc_hvdc_dc_line.argtypes = (c_char_p, c_uint, c_uint, c_char_p, c_uint)
+    libsteps.api_close_vsc_hvdc_dc_line.restype = None
+    libsteps.api_close_vsc_hvdc_dc_line.argtypes = (c_char_p, c_uint, c_uint, c_char_p, c_uint)
+
     libsteps.api_search_cct.restype = (c_double)
     libsteps.api_search_cct.argtypes = (c_char_p, c_char_p, c_uint, c_uint, c_char_p, c_uint, c_uint, c_uint)
 
diff --git a/python/stepspy-current/stepspy-dev/stepspy.py b/python/stepspy-current/stepspy-dev/stepspy.py
index a42bcbd2..ecf05dc3 100644
--- a/python/stepspy-current/stepspy-dev/stepspy.py
+++ b/python/stepspy-current/stepspy-dev/stepspy.py
@@ -23,9 +23,8 @@ def __init__(self, is_default=False, log_file=""):
         self.__encoding = "gb18030"
         if 'STEPS_LIB' not in globals():
             STEPS_LIB = pylibsteps.load_library()
-        self.STEPS_INDEX_NOT_EXIST = STEPS_LIB.api_get_const_INDEX_NOT_EXIST()
         if is_default==True:
-            self.toolkit_index = self.STEPS_INDEX_NOT_EXIST
+            self.toolkit_index = STEPS_LIB.api_get_const_INDEX_NOT_EXIST()
             self.set_toolkit_log_file(log_file)
         else:
             log_file = self.__get_c_char_p_of_string(log_file)
@@ -535,6 +534,17 @@ def get_hvdc_capacity(self):
         device = self.__get_c_char_p_of_string("HVDC")
         return int(STEPS_LIB.api_get_device_capacity(device, self.toolkit_index))
         
+    def get_vsc_hvdc_capacity(self):
+        """
+        Get capacity for storing VSC HVDC projects in database.
+        Args: N/A
+        Rets:
+            (1) Maximum number of VSC HVDC projects that can be stored in the database.
+        """
+        global STEPS_LIB
+        device = self.__get_c_char_p_of_string("VSC HVDC")
+        return int(STEPS_LIB.api_get_device_capacity(device, self.toolkit_index))    
+
     def get_equivalent_device_capacity(self):
         """
         Get capacity for storing equivalent devices in database.
@@ -626,6 +636,7 @@ def set_wt_generator_capacity(self, capacity):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Wt Generator")
         return STEPS_LIB.api_set_device_capacity(device, capacity, self.toolkit_index)
+
     def set_pv_unit_capacity(self, capacity):
         """
         Set capacity for storing PV units in database.
@@ -724,6 +735,19 @@ def set_hvdc_capacity(self, capacity):
         device = self.__get_c_char_p_of_string("HVDC")
         return STEPS_LIB.api_set_device_capacity(device, capacity, self.toolkit_index)
         
+    def set_vsc_hvdc_capacity(self, capacity):
+        """
+        Set capacity for storing VSC HVDC projects in database.
+        Args:
+            (1) capacity: Maximum number of VSC HVDC projects that can be stored in the database.
+        Rets: N/A
+        Tips:
+            This function SHOULD be called before adding devices to the database.
+        """
+        global STEPS_LIB
+        device = self.__get_c_char_p_of_string("VSC HVDC")
+        return STEPS_LIB.api_set_device_capacity(device, capacity, self.toolkit_index)  
+
     def set_equivalent_device_capacity(self, capacity):
         """
         Set capacity for storing equivalent devices in database.
@@ -913,6 +937,19 @@ def add_hvdc(self, hvdc):
         STEPS_LIB.api_add_hvdc(ibus, jbus, ickt, self.toolkit_index)
         return
 
+    def add_vsc_hvdc(self, vschvdc):
+        """
+        Add new VSC HVDC project with project name.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+        Rets: N/A
+        Example: N/A
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        STEPS_LIB.api_add_vsc_hvdc(vschvdc, self.toolkit_index)
+        return
+
     def add_transformer(self, transformer):
         """
         Add new transformer with device id.
@@ -1109,6 +1146,19 @@ def remove_hvdc(self, hvdc):
         STEPS_LIB.api_remove_hvdc(ibus, jbus, ickt, self.toolkit_index)
         return
 
+    def remove_vsc_hvdc(self, vschvdc):
+        """
+        Remove VSC HVDC project with project name.
+        Args:
+            (1) vschvdc: VSC HVDC projecct name in format of string.
+        Rets: N/A
+        Example: N/A
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        STEPS_LIB.api_remove_vsc_hvdc(vschvdc, self.toolkit_index)
+        return
+
     def remove_transformer(self, transformer):
         """
         Remove transformer with device id.
@@ -1231,6 +1281,7 @@ def get_bus_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("BUS")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_generator_count(self):
         """
         Return number of generators.
@@ -1242,6 +1293,7 @@ def get_generator_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Generator")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_wt_generator_count(self):
         """
         Return number of wind turbine generators.
@@ -1253,6 +1305,7 @@ def get_wt_generator_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("WT Generator")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_pv_unit_count(self):
         """
         Return number of PV units.
@@ -1264,6 +1317,7 @@ def get_pv_unit_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("PV Unit")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_load_count(self):
         """
         Return number of loads.
@@ -1275,6 +1329,7 @@ def get_load_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Load")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_fixed_shunt_count(self):
         """
         Return number of fixed shunt.
@@ -1286,6 +1341,7 @@ def get_fixed_shunt_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Fixed Shunt")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_line_count(self):
         """
         Return number of transmission lines.
@@ -1297,6 +1353,7 @@ def get_line_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Line")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_transformer_count(self):
         """
         Return number of transformers.
@@ -1308,6 +1365,7 @@ def get_transformer_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Transformer")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_hvdc_count(self):
         """
         Return number of HVDC links.
@@ -1319,6 +1377,18 @@ def get_hvdc_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("HVDC")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
+    def get_vsc_hvdc_count(self):
+        """
+        Return number of VSC HVDC projects.
+        Args: N/A
+        Rets: 
+            (1) Number of VSC HVDC projects.
+        """
+        global STEPS_LIB
+        device = self.__get_c_char_p_of_string("VSC HVDC")
+        return STEPS_LIB.api_get_device_count(device, self.toolkit_index)    
+
     def get_equivalent_device_count(self):
         """
         Return number of equivalent devices.
@@ -1330,6 +1400,7 @@ def get_equivalent_device_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Equivalent Device")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_energy_storage_count(self):
         """
         Return number of energy storages.
@@ -1341,6 +1412,7 @@ def get_energy_storage_count(self):
         global STEPS_LIB
         device = self.__get_c_char_p_of_string("Energy Storage")
         return STEPS_LIB.api_get_device_count(device, self.toolkit_index)
+
     def get_area_count(self):
         """
         Return number of areas.
@@ -1351,6 +1423,7 @@ def get_area_count(self):
         """
         global STEPS_LIB
         return STEPS_LIB.api_get_area_count(self.toolkit_index)
+
     def get_zone_count(self):
         """
         Return number of zones.
@@ -1361,6 +1434,7 @@ def get_zone_count(self):
         """
         global STEPS_LIB
         return STEPS_LIB.api_get_zone_count(self.toolkit_index)
+
     def get_owner_count(self):
         """
         Return number of owners.
@@ -1371,6 +1445,7 @@ def get_owner_count(self):
         """
         global STEPS_LIB
         return STEPS_LIB.api_get_owner_count(self.toolkit_index)
+        
     def get_in_service_bus_count(self):
         """
         Return number of in-service buses.
@@ -1382,6 +1457,7 @@ def get_in_service_bus_count(self):
         global STEPS_LIB
         STEPS_LIB.api_get_in_service_bus_count(self.toolkit_index)
         return
+
     def update_overshadowed_buses(self):
         """
         Update overshowed buses.
@@ -1391,6 +1467,7 @@ def update_overshadowed_buses(self):
         """
         global STEPS_LIB
         return STEPS_LIB.api_update_overshadowed_buses(self.toolkit_index)
+
     def set_all_buses_un_overshadowed(self):
         """
         Set all buses as un-overshowed.
@@ -1401,6 +1478,7 @@ def set_all_buses_un_overshadowed(self):
         global STEPS_LIB
         STEPS_LIB.api_set_all_buses_un_overshadowed(self.toolkit_index)
         return
+
     def get_overshadowed_bus_count(self):
         """
         Get number of overshowed buses. If there are n buses directly connected by zero impedance line or lines, n-1 buses are overshadowed by one of them. 
@@ -1540,6 +1618,19 @@ def is_hvdc_exist(self, hvdc):
         ickt = self.__get_c_char_p_of_string(ickt)
         return STEPS_LIB.api_is_hvdc_exist(ibus, jbus, ickt, self.toolkit_index)
 
+    def is_vsc_hvdc_exist(self, vschvdc):
+        """
+        Check if given VSC HVDC project exists or not.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+        Rets:
+            (1) True if the VSC HVDC project exists, False otherwise.
+        Example: N/A
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        return STEPS_LIB.api_is_vsc_hvdc_exist(vschvdc, self.toolkit_index)
+
     def is_equivalent_device_exist(self, equivalent_device):
         """
         Check if given equivalent device exists or not.
@@ -2066,6 +2157,39 @@ def get_hvdcs_between_buses(self, ibus, jbus):
                     dcs.append(dc) 
             return tuple(dcs)
 
+    def get_all_vsc_hvdcs(self):
+        """
+        Get all VSC HVDC projects in the database.
+        Args: N/A
+        Rets:
+            (1) Tuple of all VSC HVDC projects. Empty tuple if no VSC HVDC projects in the database.
+        """
+        return self.get_vsc_hvdcs_at_bus(0)
+
+    def get_vsc_hvdcs_at_bus(self, bus):
+        """
+        Get all VSC HVDC projects at bus with input bus number.
+        Args:
+            (1) bus: bus number
+        Rets:
+            (1) Tuple of all VSC HVDC projects at given bus. Empty tuple if no VSC HVDC projects at given bus.
+        """
+        global STEPS_LIB
+        device = "VSC HVDC"
+        device = self.__get_c_char_p_of_string(device)
+        blank = self.__get_c_char_p_of_string("")
+        STEPS_LIB.api_initialize_device_search(device, bus, self.toolkit_index)
+        vscs = []
+        while True:
+            nbus = STEPS_LIB.api_get_current_device_bus_number(device, blank, self.toolkit_index)
+            if nbus==0:
+                break
+            id = STEPS_LIB.api_get_current_device_identifier(device, self.toolkit_index)
+            id = self.__get_string_from_c_char_p(id)
+            vscs.append((int(nbus), id))
+            STEPS_LIB.api_goto_next_device(device, self.toolkit_index)
+        return tuple(vscs)
+
     def get_generators_with_constraints(self, area=0, zone=0):
         """
         Get all generators satisfying area and zone constraints.
@@ -2592,6 +2716,40 @@ def get_hvdc_data(self, hvdc, par_type, side, par_name):
             return self.__get_string_from_c_char_p(STEPS_LIB.api_get_hvdc_string_data(ibus, jbus, ickt, side, par_name, self.toolkit_index))
         return None
 
+    def get_vsc_hvdc_data(self, vschvdc, dev_type, par_type, par_name, index):
+        """
+        Get VSC HVDC project data.
+        Args:
+            (1) vschvdc: VSC HVDC project name
+            (2) dev_type: String of device type. Choose one from("PROJECT", "CONVERTER", "DC BUS", "DC LINE")
+            (3) par_type: String of parameter type. Choose one from {"I", "F", "D", "S", "B"}.
+            (4) par_name: String of parameter name.
+            (5) index: Int of converter or dc bus or dc line index
+        Rets:
+            (1) Value of parameter.
+        Tips:
+            The par_type meaning: "I": integer number, "F" or "D": float number, "S": string, "B": boolean data.
+            The type of given parameter MUST be consistent with the given parameter type. Otherwise, 0, 0.0, "", or False will be returned.
+        """
+        global STEPS_LIB
+        par_type = par_type.upper()
+        if par_type not in ['I', 'INT', 'INTEGER', 'F', 'D', 'FLOAT', 'DOUBLE', 'B', 'BOOL', 'BOOLEAN', 'S', 'STRING']:
+            return None
+        dev_type = dev_type.upper()
+        par_name = par_name.upper()
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        dev_type = self.__get_c_char_p_of_string(dev_type)
+        par_name = self.__get_c_char_p_of_string(par_name)
+        if par_type in ['I', 'INT', 'INTEGER']:
+            return int(STEPS_LIB.api_get_vsc_hvdc_integer_data(vschvdc, dev_type, par_name, index, self.toolkit_index))
+        if par_type in ['F', 'D', 'FLOAT', 'DOUBLE']:
+            return STEPS_LIB.api_get_vsc_hvdc_float_data(vschvdc, dev_type, par_name, index, self.toolkit_index)
+        if par_type in ['B', 'BOOL', 'BOOLEAN']:
+            return STEPS_LIB.api_get_vsc_hvdc_boolean_data(vschvdc, dev_type, par_name, index, self.toolkit_index)
+        if par_type in ['S', 'STRING']:
+            return self.__get_string_from_c_char_p(STEPS_LIB.api_get_vsc_hvdc_string_data(vschvdc, dev_type, par_name, index, self.toolkit_index))
+        return None
+
     def get_area_data(self, area, par_type, par_name):
         """
         Get area data.
@@ -3104,6 +3262,7 @@ def set_hvdc_data(self, hvdc, par_type, side, par_name, value):
             value = self.__get_c_char_p_of_string(value)
             return STEPS_LIB.api_set_hvdc_string_data(ibus, jbus, ickt, side, par_name, value, self.toolkit_index)
         return
+
     def set_hvdc_power(self, hvdc, p):
         """
         Set HVDC link power command.
@@ -3118,6 +3277,41 @@ def set_hvdc_power(self, hvdc, p):
         self.set_hvdc_data(hvdc, "F", "HVDC", "PDCN_MW", p)
         return
         
+    def set_vsc_hvdc_data(self, vschvdc, dev_type, par_type, par_name, value, index):
+        """
+        Set VSC HVDC project data.
+        Args:
+            (1) vschvdc: VSC HVDC project name
+            (2) dev_type: String of device type. Choose one from("PROJECT", "CONVERTER", "DC BUS", "DC LINE")
+            (3) par_type: String of parameter type. Choose one from {"I", "F", "D", "S", "B"}.
+            (4) par_name: String of parameter name.
+            (5) value: Value of parameter.
+            (6) index: Int of converter or dc bus or dc line index
+        Rets: N/A
+        Tips:
+            The par_type meaning: "I": integer number, "F" or "D": float number, "S": string, "B": boolean data.
+            The type of given parameter MUST be consistent with the given parameter type. Otherwise, nothing will be done.
+        """
+        global STEPS_LIB
+        par_type = par_type.upper()
+        if par_type not in ['I', 'INT', 'INTEGER', 'F', 'D', 'FLOAT', 'DOUBLE', 'B', 'BOOL', 'BOOLEAN', 'S', 'STRING']:
+            return None
+        dev_type = dev_type.upper()
+        par_name = par_name.upper()
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        dev_type = self.__get_c_char_p_of_string(dev_type)
+        par_name = self.__get_c_char_p_of_string(par_name)
+        if par_type in ['I', 'INT', 'INTEGER']:
+            return int(STEPS_LIB.api_set_vsc_hvdc_integer_data(vschvdc, dev_type, par_name, value, index, self.toolkit_index))
+        if par_type in ['F', 'D', 'FLOAT', 'DOUBLE']:
+            return STEPS_LIB.api_set_vsc_hvdc_float_data(vschvdc, dev_type, par_name, value, index, self.toolkit_index)
+        if par_type in ['B', 'BOOL', 'BOOLEAN']:
+            return STEPS_LIB.api_set_vsc_hvdc_boolean_data(vschvdc, dev_type, par_name, value, index, self.toolkit_index)
+        if par_type in ['S', 'STRING']:
+            value = self.__get_c_char_p_of_string(value)
+            return self.__get_string_from_c_char_p(STEPS_LIB.api_set_vsc_hvdc_string_data(vschvdc, dev_type, par_name, value, index, self.toolkit_index))
+        return None
+
     def set_area_data(self, area, par_type, par_name, value):
         """
         Set area data.
@@ -3225,6 +3419,19 @@ def load_powerflow_data(self, file, ftype):
         ftype = self.__get_c_char_p_of_string(ftype)
         STEPS_LIB.api_load_powerflow_data_from_file(file, ftype, self.toolkit_index)
 
+    def load_vsc_hvdc_powerflow_data(self, file, ftype):
+        """
+        Load VSC HVDC powerflow data from file.
+        Args:
+            (1) file: string, source VSC HVDC powerflow file name.
+            (2) ftype: string, powerflow data format.
+        Rets: N/A
+        """
+        global STEPS_LIB
+        file = self.__get_c_char_p_of_string(file)
+        ftype = self.__get_c_char_p_of_string(ftype)
+        STEPS_LIB.api_load_vsc_hvdc_powerflow_data_from_file(file, ftype, self.toolkit_index)
+
     def check_powerflow_data(self):
         """
         Check powerflow data. If any inappropriate data is set, report will be sent to log file.
@@ -3261,6 +3468,20 @@ def save_powerflow_data(self, file, ftype, export_zero_line=True, export_out_of_
         else:
             print("parameter export_mode is invalid in save_powerflow_data()")
         
+    def save_vsc_hvdc_powerflow_data(self, file, ftype):
+        """
+        Save VSC HVDC powerflow data to file
+        Args:
+            (1) file: string, target powerflow file name.
+            (2) ftype: string, powerflow data format.
+        Rets: N/A
+        """
+
+        global STEPS_LIB
+        file = self.__get_c_char_p_of_string(file)
+        ftype = self.__get_c_char_p_of_string(ftype)
+        STEPS_LIB.api_save_vsc_hvdc_powerflow_data_to_file(file, ftype, self.toolkit_index)    
+
     def __save_powerflow_data_in_keep_mode(self, file, ftype, export_zero_line=True, export_out_of_service_bus=True, export_internal_bus_number=False):
         """
         Save powerflow data to file in keep as original mode
@@ -4314,7 +4535,7 @@ def get_hvdc_related_model_data(self, hvdc, model_type, par_name):
         
     def set_hvdc_related_model_data(self, hvdc, model_type, par_name, value):
         """
-        Set HVDC link related model data.
+        Set HVDC linke related model data.
         Args:
             (1) hvdc: HVDC link device id in format of (ibus, jbus, ickt).
             (2) model_type: String of model type.
@@ -4358,13 +4579,12 @@ def get_hvdc_related_model_parameter_pair(self, hvdc, model_type):
             parameters.append((par_name, par_value))
         return tuple(parameters)
         
-    def get_vsc_hvdc_related_model_name(self, vsc_hvdc_name, model_type, converter_index=self.STEPS_INDEX_NOT_EXIST):
+    def get_vsc_hvdc_related_model_name(self, vschvdc, model_type):
         """
-        Get VSC HVDC related model name.
+        Get VSC HVDC project related model name.
         Args:
-            (1) vsc_hvdc_name: String of VSC HVDC name.
+            (1) vschvdc: VSC HVDC project name in format of string.
             (2) model_type: String of model type.
-            (3) converter_index: Integer if model_type is "VSC HVDC CONVERTER"
         Rets:
             (1) String of model name.
         Tips:
@@ -4372,21 +4592,20 @@ def get_vsc_hvdc_related_model_name(self, vsc_hvdc_name, model_type, converter_i
         Example: N/A
         """
         global STEPS_LIB
-        vsc_hvdc_name = self.__get_c_char_p_of_string(vsc_hvdc_name)
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
         model_type = self.__get_c_char_p_of_string(model_type)
-        model_name = self.__get_string_from_c_char_p(STEPS_LIB.api_get_vsc_hvdc_related_model_name(vsc_hvdc_name, model_type, converter_index, self.toolkit_index))
+        model_name = self.__get_string_from_c_char_p(STEPS_LIB.api_get_vsc_hvdc_related_model_name(vschvdc, model_type, self.toolkit_index))
         if model_name=="":
             model_name = None
         return model_name
 
-    def get_vsc_hvdc_related_model_data(self, vsc_hvdc_name, model_type, par_name, converter_index=self.STEPS_INDEX_NOT_EXIST):
+    def get_vsc_hvdc_related_model_data(self, vschvdc, model_type, par_name):
         """
-        Get VSC HVDC related model data.
+        Get VSC HVDC project related model data.
         Args:
-            (1) vsc_hvdc_name: String of VSC HVDC name.
+            (1) vschvdc: VSC HVDC project name in format of string.
             (2) model_type: String of model type.
             (3) par_name: String of parameter name.
-            (4) converter_index: Integer if model_type is "VSC HVDC CONVERTER"
         Rets:
             (1) Value of parameter
         Tips:
@@ -4394,20 +4613,19 @@ def get_vsc_hvdc_related_model_data(self, vsc_hvdc_name, model_type, par_name, c
         Example: N/A
         """
         global STEPS_LIB
-        vsc_hvdc_name = self.__get_c_char_p_of_string(vsc_hvdc_name)
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
         model_type = self.__get_c_char_p_of_string(model_type)
         par_name = self.__get_c_char_p_of_string(par_name)
-        return STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter(vsc_hvdc_name, model_type, par_name, converter_index, self.toolkit_index)
+        return STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter(vschvdc, model_type, par_name, self.toolkit_index)
         
-    def set_vsc_hvdc_related_model_data(self, vsc_hvdc_name, model_type, par_name, value, converter_index=self.STEPS_INDEX_NOT_EXIST):
+    def set_vsc_hvdc_related_model_data(self, vschvdc, model_type, par_name, value):
         """
-        Set VSC HVDC related model data.
+        Set VSC HVDC project related model data.
         Args:
-            (1) vsc_hvdc_name: String of VSC HVDC name.
+            (1) vschvdc: VSC HVDC project name in format of string.
             (2) model_type: String of model type.
             (3) par_name: String of parameter name.
             (4) value: Value of parameter.
-            (5) converter_index: Integer if model_type is "VSC HVDC CONVERTER"
         Rets: N/A
         Tips:
             If model type or parameter name is not supported, nothing will be changed.
@@ -4415,18 +4633,17 @@ def set_vsc_hvdc_related_model_data(self, vsc_hvdc_name, model_type, par_name, v
         Example: N/A
         """
         global STEPS_LIB
-        vsc_hvdc_name = self.__get_c_char_p_of_string(vsc_hvdc_name)
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
         model_type = self.__get_c_char_p_of_string(model_type)
         par_name = self.__get_c_char_p_of_string(par_name)
-        return STEPS_LIB.api_set_vsc_hvdc_related_model_float_parameter(vsc_hvdc_name, model_type, par_name, value, converter_index, self.toolkit_index)
+        return STEPS_LIB.api_set_vsc_hvdc_related_model_float_parameter(vschvdc, model_type, par_name, value, self.toolkit_index)
         
-    def get_vsc_hvdc_related_model_parameter_pair(self, vsc_hvdc_name, model_type, converter_index=self.STEPS_INDEX_NOT_EXIST):
+    def get_vsc_hvdc_related_model_parameter_pair(self, vschvdc, model_type):
         """
-        Get VSC HVDC related model parameter pair.
+        Get VSC HVDC project related model parameter pair.
         Args:
-            (1) vsc_hvdc_name: String of VSC HVDC name.
+            (1) vschvdc: VSC HVDC project name in format of string.
             (2) model_type: String of model type.
-            (3) converter_index: Integer if model_type is "VSC HVDC CONVERTER"
         Rets:
             (1) Tuple of parameter name and value pairs. Each parameter name and value pair is in format of (string of parameter name, value of parameter).
         Tips:
@@ -4434,17 +4651,17 @@ def get_vsc_hvdc_related_model_parameter_pair(self, vsc_hvdc_name, model_type, c
         Example: N/A
         """
         global STEPS_LIB
-        vsc_hvdc_name = self.__get_c_char_p_of_string(vsc_hvdc_name)
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
         model_type = self.__get_c_char_p_of_string(model_type)
         parameters = []
-        n = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter_count(vsc_hvdc_name, model_type, converter_index, self.toolkit_index)
+        n = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter_count(vschvdc, model_type, self.toolkit_index)
         for parameter_index in range(n):
-            par_name = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter_name(vsc_hvdc_name, model_type, parameter_index, converter_index, self.toolkit_index)
-            par_value = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter(vsc_hvdc_name, model_type, par_name, converter_index, self.toolkit_index);
+            par_name = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter_name(vschvdc, model_type, parameter_index, self.toolkit_index)
+            par_value = STEPS_LIB.api_get_vsc_hvdc_related_model_float_parameter(vschvdc, model_type, par_name, self.toolkit_index);
             par_name = self.__get_string_from_c_char_p(par_name)
             parameters.append((par_name, par_value))
         return tuple(parameters)
-        
+
     def get_generator_governor_pmax(self, gen):
         """
         Get generator pmax in MW from turbine governor model.
@@ -4732,13 +4949,12 @@ def prepare_meters(self, device_type):
                 DYNAMIC_SIMULATOR::prepare_line_related_meters()
                 DYNAMIC_SIMULATOR::prepare_transformer_related_meters()
                 DYNAMIC_SIMULATOR::prepare_hvdc_related_meters()
+                DYNAMIC_SIMULATOR::prepare_vsc_hvdc_related_meters()
                 DYNAMIC_SIMULATOR::prepare_equivalent_device_related_meters()                
         Example: N/A
         """
         global STEPS_LIB
         device_type = device_type.upper()
-        if device_type not in ['ALL', 'BUS', 'GENERATOR', 'WT GENERATOR', 'PV UNIT', 'ENERGY STORAGE', 'LOAD', 'LINE', 'TRANSFORMER', 'HVDC', 'EQUIVALENT DEVICE']:
-            return
         if device_type in ['ALL']:
             STEPS_LIB.api_prepare_meters(self.toolkit_index)
             return
@@ -4769,6 +4985,9 @@ def prepare_meters(self, device_type):
         if device_type in ['HVDC']:
             STEPS_LIB.api_prepare_hvdc_related_meters(self.toolkit_index)
             return
+        if device_type in ['VSC HVDC']:
+            STEPS_LIB.api_prepare_vsc_hvdc_related_meters(self.toolkit_index)
+            return
         if device_type in ['EQUIVALENT DEVICE']:
             STEPS_LIB.api_prepare_equivalent_device_related_meters(self.toolkit_index)
             return
@@ -4966,6 +5185,27 @@ def prepare_hvdc_meter(self, hvdc, meter_type, side, var_name=""):
         STEPS_LIB.api_prepare_hvdc_related_meter(ibus, jbus, ickt, meter_type, side, var_name, self.toolkit_index)
         return
 
+    def prepare_vsc_hvdc_meter(self, vschvdc, meter_type, side, var_name=""):
+        """
+        Prepare specific VSC HVDC project meter.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+            (2) meter_type: String of meter type.
+            (3) side: index of Converter or DC line.
+            (4) var_name: String of internal variable name. Only valid when meter type is internal variable of dynamic models.
+        Rets: N/A
+        Tips:
+            For the supported meter types, see implementation of the following function of STEPS:
+                DYNAMIC_SIMULATOR::prepare_vsc_hvdc_related_meter()
+        Example: N/A
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        meter_type = self.__get_c_char_p_of_string(meter_type)
+        var_name = self.__get_c_char_p_of_string(var_name)
+        STEPS_LIB.api_prepare_vsc_hvdc_related_meter(vschvdc, meter_type, side, var_name, self.toolkit_index)
+        return
+
     def prepare_equivalent_device_meter(self, edevice, meter_type, var_name=""):
         """
         Prepare specific equivalent device meter.
@@ -5879,4 +6119,91 @@ def set_hvdc_power_order_in_MW(self, hvdc, value):
         ibus, jbus, ickt = self.__extract_double_bus_device_id(hvdc)
         ickt = self.__get_c_char_p_of_string(ickt)
         return STEPS_LIB.api_set_hvdc_power_order_in_MW(ibus, jbus, ickt, value, self.toolkit_index)
-        
\ No newline at end of file
+        
+    def set_vsc_hvdc_dc_line_fault(self, vschvdc, line, fault_location, fault_r):
+        """
+        Set VSC HVDC DC line fault.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+            (2) line: DC line device id in format of (ibus, jbus, ickt).
+            (3) fault_location: Relative fault location to ibus.
+            (4) fault_r: Fault resistance in Ohm.
+        Rets: N/A
+        Tips:
+            The fault location should be in the range of [0, 1.0], including 0 and 1.0. It represent the relative location of the fault on the line to the ibus.
+            For example, 0.5 means the fault is set at the middle of the line. 0 means the fault is set at exactly ibus. 1.0 means the fault is set at exactly jbus.
+            The absolute value of the fault shunt should not be too great. Otherwise, network solution may fail to converge. Typically, |b|<1e6.
+        Example:
+            vschvdc = "VSC ProjectA"
+            line = (1, 2, "1#")
+            location = 0.2
+            fault_r = 10
+            simulator.set_vsc_hvdc_dc_line_fault(vschvdc, line, location, shunt)
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        ibus, jbus, ickt = self.__extract_double_bus_device_id(line)
+        ickt = self.__get_c_char_p_of_string(ickt)
+        STEPS_LIB.api_set_vsc_hvdc_dc_line_fault(vschvdc, ibus, jbus, ickt, fault_location, fault_r, self.toolkit_index)
+        return
+
+    def clear_vsc_hvdc_dc_line_fault(self, vschvdc, line):
+        """
+        Clear VSC HVDC DC line fault without tripping the line.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+            (2) line: DC line device id in format of (ibus, jbus, ickt).
+        Rets: N/A
+        Tips: N/A
+        Example:
+            vschvdc = "VSC ProjectA"
+            line = (1, 2, "1#")
+            location = 0.2
+            fault_r = 10
+            simulator.set_vsc_hvdc_dc_line_fault(vschvdc, line, location, shunt)
+            simulator.clear_vsc_hvdc_dc_line_fault(vschvdc, line)
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        ibus, jbus, ickt = self.__extract_double_bus_device_id(line)
+        ickt = self.__get_c_char_p_of_string(ickt)
+        STEPS_LIB.api_clear_vsc_hvdc_dc_line_fault(vschvdc, ibus, jbus, ickt, self.toolkit_index)
+        return
+
+    def trip_vsc_hvdc_dc_line(self, vschvdc, line):
+        """
+        Trip VSC HVDC DC line. Breakers at the two sides of the line are both tripped.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+            (2) line: DC line device id in format of (ibus, jbus, ickt).
+        Rets: N/A
+        Example:
+            vschvdc = "VSC ProjectA"
+            line = (1, 2, "1#")
+            simulator.trip_vsc_hvdc_dc_line(vschvdc, line)
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        ibus, jbus, ickt = self.__extract_double_bus_device_id(line)
+        ickt = self.__get_c_char_p_of_string(ickt)
+        STEPS_LIB.api_trip_vsc_hvdc_dc_line(vschvdc, ibus, jbus, ickt, self.toolkit_index)
+        return
+
+    def close_vsc_hvdc_dc_line(self, vschvdc, line):
+        """
+        Close VSC HVDC DC line. Breakers at the two sides of the line are both closed.
+        Args:
+            (1) vschvdc: VSC HVDC project name in format of string.
+            (2) line: DC line device id in format of (ibus, jbus, ickt).
+        Rets: N/A
+        Example:
+            vschvdc = "VSC ProjectA"
+            line = (1, 2, "1#")
+            simulator.close_vsc_hvdc_dc_line(vschvdc, line)
+        """
+        global STEPS_LIB
+        vschvdc = self.__get_c_char_p_of_string(vschvdc)
+        ibus, jbus, ickt = self.__extract_double_bus_device_id(line)
+        ickt = self.__get_c_char_p_of_string(ickt)
+        STEPS_LIB.api_close_vsc_hvdc_dc_line(vschvdc, ibus, jbus, ickt, self.toolkit_index)
+        return
\ No newline at end of file