Para calcular la temperatura modificada del pinch, añadimos al final del programa estas líneas:
tPinch = tabla.getTemperaturaModificadaPinch() print 'T pinch = '+str(tPinch)+' ºC'
En cuanto al código para dibujar las curvas compuestas:
# Curva Compuesta Caliente caliente_file = open('curva_caliente','w') t = tabla.getCurvaCompuestaCalienteOriginal().getListaDeTemperaturas() h = tabla.getCurvaCompuestaCalienteOriginal().getListaEntalpias() for i in range(len(t)): caliente_file.write(str(t[i])+'\t'+str(h[i])+'\n') caliente_file.close() # Curva compuesta fría fria_file = open('curva_fria','w') t = tabla.getCurvaCompuestaFriaOriginal().getListaDeTemperaturas() h = tabla.getCurvaCompuestaFriaOriginal().getListaEntalpias() for i in range(len(t)): h[i] += tabla.getQC() fria_file.write(str(t[i])+'\t'+str(h[i])+'\n') fria_file.close()
Ahora, para dibujar la gráfica tenemos que hacer dos modificaciones. La primera, añadir la línea import os justo despúes de la línea para importar el módulo PinchPython. Después añadimos las siguientes líneas al final del programa:
# Comando para dibujar la gráfica os.system('echo "set terminal png small; set size .65; \ set grid; set output \'grafica.png\';\ set ylabel \'T (C)\'; set xlabel \'H (MW)\';\ plot \'curva_caliente\' using 2:1 with lines title \'Caliente\',\ \'curva_fria\' using 2:1 with lines title \'Fria\'"\ | gnuplot')
Si ejecutamos de nuevo el programa, los resultados son:
Qh = 11.5 MW Qc = 14.0 MW T pinch = 150.0 ºC
En este caso además deben haberse creado tres ficheros. El fichero grafica.png es una imagen de la gráfica de las curvas compuestas. En la figura se muestra la gráfica que se ha obtenido con el programa.
2004-05-30