Add if __name__ == '__main__' statement
@@ -494,6 +494,9 @@ def calcStationaryAgeDstn(LivPrb,terminal_period):
AgeDstn = (x/np.sum(x))
return AgeDstn
+####################################################################################################
+
+if __name__ == '__main__':
# Set targets for K/Y and the Lorenz curve based on the data
if Params.do_liquid:
