show_plot = True
font_size = 16
figsize=[8,6]
plot = []
#plot.append(dict(kind='function',data=['tracking/mix_line_p*_t0.*.res'],
#                 row=[1,-1],col=[1,4], startplot = True,
#                 fun = lambda *x:x[0]*x[1]))
files = ['tracking/mix_line_p*_t1.000E+03.res']
#files_2 = ['/home/UNISIEGEN/gk779/apes_default/musubi/testsuite/validation/liquidModel/stefanTube/tracking/mix_line_p*_t10.000E+03.res']
files_2 = ['tracking_dt1e-2/mix_line_p*_t1.000E+03.res']
#files_3 = ['tracking_eq_dt1e-1/mix_line_p*_t1.000E+03.res']
files_3 = ['tracking/mix_line_p*_t1.000E+03.res']

reference = ['stefantube_reference/maxwell_stefan_reference.res']

plot.append(dict(kind='xy', data=files, 
                 row=[1,-1],col=[1,4], startplot = True,
                 convert=True, facs=[1.0/0.238,1.0],
                 label='Species 1 (LBM)',ls = 'k-') ) 

plot.append(dict(kind='xy', data=files,
                 convert=True, facs=[1.0/0.238,1.0],
                 row=[1,-1],col=[1,5], label='Species 2 (LBM)',ls = 'k--'))

plot.append(dict(kind='xy', data=files,
                 convert=True, facs=[1.0/0.238,1.0],
                 row=[1,-1],col=[1,6], label='Species 3 (LBM)',ls = 'k-.'))

#plot.append(dict(kind='xy', data=files_2, 
#                 row=[1,-1],col=[1,4], 
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 label='Species 1 (LBM) wall dt = 1e-2',ls = 'r-') ) 
#
#plot.append(dict(kind='xy', data=files_2,
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 row=[1,-1],col=[1,5], label='Species 2 (LBM) wall dt = 1e-2',ls = 'r--'))
#
#plot.append(dict(kind='xy', data=files_2,
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 row=[1,-1],col=[1,6], label='Species 3 (LBM) wall dt = 1e-2',ls = 'r-.'))
# 
#plot.append(dict(kind='xy', data=files_3, 
#                 row=[1,-1],col=[1,4], 
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 label='Species 1 (LBM) eq dt = 1e-1',ls = 'b-') ) 
#
#plot.append(dict(kind='xy', data=files_3,
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 row=[1,-1],col=[1,5], label='Species 2 (LBM) eq dt = 1e-1',ls = 'b--'))
#
#plot.append(dict(kind='xy', data=files_3,
#                 convertPhy=True, facs=[1.0/0.238,1.0],
#                 row=[1,-1],col=[1,6], label='Species 3 (LBM) eq dt = 1e-1',ls = 'b-.'))

plot.append(dict(kind='xy', data=reference, 
                 row=[1,-1],col=[1,2], npoints = 20,
                 label='Species 1 (shooting)',ls = 'kx') ) 

plot.append(dict(kind='xy', data=reference,
                 row=[1,-1],col=[1,3], npoints = 20,
                 label='Species 2 (shooting)',ls = 'kd'))

plot.append(dict(kind='xy', data=reference,
                 row=[1,-1],col=[1,4], npoints = 20,
                 label='Species 3 (shooting)',ls = 'kv', 
                 legend=dict(loc = 2),
                 endplot=True, figname = 'StefanTube', format='png', dpi=100,
                 xlabel = 'x/Length', ylabel='Mole fraction ($n_i/n$)'))
               
