modified DL_calc to be callable and take the paths as attributes

applications/performDL/DL_calculation.py

@@ -1,39 +1,44 @@
 import sys
-
-sys.path.append('C:/Adam/pelicun')
-
 import numpy as np
 
 import pelicun
 from pelicun.control import FEMA_P58_Assessment
 from pelicun.file_io import write_SimCenter_DL_output
 
-# location of DL input file
-DL_input_path = "/Users/fmckenna/Documents/tmp.SimCenter/templatedir/dakota.json"
-EDP_input_path = "/Users/fmckenna/Documents/tmp.SimCenter/dakotaTab.out"
-CMP_data_path = '/Users/fmckenna/NHERI/resources/fragilities/'
-POP_data_path = '/Users/fmckenna/NHERI/resources/population.json'
+def run_pelicun(DL_input_path, EDP_input_path, CMP_data_path, POP_data_path):
+
+	A = FEMA_P58_Assessment()
+	A.read_inputs(DL_input_path, EDP_input_path, 
+		CMP_data_path, POP_data_path, verbose=False)
+
+	A.define_random_variables()
+
+	A.define_loss_model()
+
+	A.calculate_damage()
+
+	A.calculate_losses()
+
+	A.aggregate_results()
+
+	S = A._SUMMARY
 
-A = FEMA_P58_Assessment()
-A.read_inputs(DL_input_path, EDP_input_path, 
-	CMP_data_path, POP_data_path, verbose=False)
+	S.index.name = '#Num'
 
-A.define_random_variables()
+	S.columns = [('{}/{}'.format(s0, s1)).replace(' ', '_') 
+	             for s0, s1 in zip(S.columns.get_level_values(0),
+	                               S.columns.get_level_values(1))]
 
-A.define_loss_model()
+	write_SimCenter_DL_output('DL_summary.csv', S)
 
-A.calculate_damage()
+	S_stat = S.describe(np.arange(1, 100)/100.)
 
-A.calculate_losses()
+	S_stat.index.name = 'attribute'
 
-A.aggregate_results()
+	write_SimCenter_DL_output('DL_summary_stats.csv', S_stat)
 
-S = A._SUMMARY
+	return 0
 
-S.columns = [('{}/{}'.format(s0, s1)).replace(' ', '_') 
-             for s0, s1 in zip(S.columns.get_level_values(0),
-                               S.columns.get_level_values(1))]
+if __name__ == '__main__':
 
-write_SimCenter_DL_output('DL_summary.csv', S)
-write_SimCenter_DL_output('DL_summary_stats.csv', 
-                          S.describe(np.arange(1, 100)/100.))
+	sys.exit(run_pelicun(sys.argv[1], sys.argv[2], sys.argv[3], sys.argv[4]))