plot  = []

show_plot = True
font_family= 'serif'
font_type= 'Times New Roman'

wall = 'tracking_wall/mix_line_p00000_t10.000E+03.res'
moleflux = 'tracking_moleFlux/mix_line_p00000_t10.000E+03.res'

wall_time = 'tracking_wall/mix_point_p00000.res'
moleflux_time = 'tracking_moleFlux/mix_point_p00000.res'

plot.append(dict(kind='xy', data=[wall], 
                 row=[1,-1],col=[1,4], startplot = True,
                 label='Acetone: wall BC',ls = 'k-') ) 

plot.append(dict(kind='xy', data=[wall],
                 row=[1,-1],col=[1,5], label='Methanol: wall BC',ls = 'k--'))

plot.append(dict(kind='xy', data=[wall],
                 row=[1,-1],col=[1,6], label='Air: wall BC',ls = 'k-.'))

plot.append(dict(kind='xy', data=[moleflux],
                 row=[1,-1],col=[1,4], label='Methanol: wall BC',ls = 'r-'))

plot.append(dict(kind='xy', data=[moleflux],
                 row=[1,-1],col=[1,5], label='Methanol: wall BC',ls = 'r--'))

plot.append(dict(kind='xy', data=[moleflux],
                 row=[1,-1],col=[1,6], endplot = True,
                 label='Air: moleFlux BC',ls = 'r-.', xlabel='length (m)', 
                 ylabel ='mole fraction ($\chi_i =n_i/n_t$)',
                 format = 'eps', figname = 'stefantube', dpi=100))

plot.append(dict(kind='xy', data=[wall_time], 
                 row=[1,-1],col=[1,5], startplot = True,
                 label='Acetone: wall BC',ls = 'k-') ) 

plot.append(dict(kind='xy', data=[wall_time], 
                 row=[1,-1],col=[1,8], 
                 label='Methanol: wall BC',ls = 'k--') )

plot.append(dict(kind='xy', data=[moleflux_time], 
                 row=[1,-1],col=[1,5], 
                 label='Acetone: moleFlux BC',ls = 'r-') ) 

plot.append(dict(kind='xy', data=[moleflux_time], 
                 row=[1,-1],col=[1,8], endplot = True,
                 label='Methanol: moleFlux BC',ls = 'r--',
                 xlabel='time (s)', ylabel =' mole flux [$mol/m^2/s$]',
                 title = 'mole flux of acetone and methanol at top of the tube',
                 format = 'eps', figname = 'moleFluxOverTime', dpi=100))

