for key in Model_dict.keys():
    
    #Declare scripting name
    ScriptName = ScriptPrefix + key + '.in'
    
    #Generate the script
    ct.select_script('Argon_Sulfur_LinesModel_Ages', ScriptName, ScriptFolder, key = key, data_dict=Model_dict)
      
    #Run the cloudy script
#     ct.run_script(ScriptName, ScriptFolder)
#       
    #Plot the .out cloudy file
#     ct.open_output(ScriptName, ScriptFolder)

    #Generate lines dictionary with the output data
    Line_dict = ct.load_predicted_lines(ScriptName, ScriptFolder)
    
    x_values, y_values, TSIII, TOIII = Ar_S_model(Line_dict)
    dz.data_plot(x_values, y_values, color=Colors_dict[key], label=Legends_dict[key], markerstyle='o')
#     dz.data_plot(x_values, TSIII,  color=Colors_dict[key], label=Legends_dict[key], markerstyle='o')

    list_xvalues    = concatenate([list_xvalues, x_values])
    list_yvalues    = concatenate([list_yvalues, y_values])
    list_TSIII      = concatenate([list_TSIII, TSIII])
    list_TOIII      = concatenate([list_TOIII, TOIII])

Not_infinite = isfinite(list_yvalues)

list_xvalues_clean  = list_xvalues[Not_infinite]
list_yvalues_clean  = list_yvalues[Not_infinite]
list_TSIII_clean    = list_TSIII[Not_infinite]
  
Objects                     = Abundances_Matrix[:,0]
ArIII_HII_array             = Abundances_Matrix[:,1]
ArIV_HII_array              = Abundances_Matrix[:,2]
Temps                       = Abundances_Matrix[:,3]    
SIII_HII_array              = Abundances_Matrix[:,4]   
    
logArII_ArIII               = umlog10(ArIII_HII_array/ArIV_HII_array)

for key in Model_dict.keys():
    
    #Declare scripting name
    ScriptName = ScriptPrefix + key + '.in'
    
    #Generate lines dictionary with the output data
    Line_dict = ct.load_predicted_lines(ScriptName, ScriptFolder)
    
    print key, Colors_dict[key], Legends_dict[key]

    x_values, y_values, TSIII, TOIII = Ar_S_model(Line_dict, diags)
    dz.data_plot(x_values, y_values, color=Colors_dict[key], label=Legends_dict[key], markerstyle='o')
#     dz.data_plot(TSIII, y_values,  color=Colors_dict[key], label=Legends_dict[key], markerstyle='o')

    list_xvalues    = concatenate([list_xvalues, x_values])
    list_yvalues    = concatenate([list_yvalues, y_values])
    list_TSIII      = concatenate([list_TSIII, TSIII])
    list_TOIII      = concatenate([list_TOIII, TOIII])

Not_infinite = isfinite(list_yvalues)

list_xvalues_clean  = list_xvalues[Not_infinite]