#from simParam import *
import math
plot  = []

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

#plot.append(dict(kind='xy', data=['tracking_n40/mix_numDensity*_t10*.res'],
#                 row=[1,-1],col=[1,4], startplot = True, label='n40',ls='-' ))
#plot.append(dict(kind='xy', data=['tracking_n40/mix_numDensity*_t0*.res'],
#                 row=[1,-1],col=[1,4], label='n40 t0',ls='-' ))
time = ['0.000E+00', '100.000E+00', '200.000E+00','300.000E+00','400.000E+00','500.000E+00','600.000E+00','700.000E+00','800.000E+00','900.000E+00','1.000E+03']

for itime in range(len(time)):
  splot = False
  eplot = False
  if itime == 0:
    splot = True
    eplot = False
  if itime == len(time)-1:
    splot = False
    eplot = True
  plot.append(dict(kind='xy', data=['tracking/mix_numDensity*_t'+time[itime]+'.res'],
                   row=[1,-1],col=[1,4], startplot = splot, endplot = eplot,
                   label='H20, n20, t'+time[itime],ls='-', 
                   xlabel = 'length (m)', ylabel = 'molar densiy [$mol/m^3$]',
                   format = 'png', figname = 'H20_concenOnLengthOverTime', dpi=100))

for itime in range(len(time)):
  splot = False
  eplot = False
  if itime == 0:
    splot = True
    eplot = False
  if itime == len(time)-1:
    splot = False
    eplot = True
  plot.append(dict(kind='xy', data=['tracking/mix_numDensity*_t'+time[itime]+'.res'],
                   row=[1,-1],col=[1,5], startplot = splot, endplot = eplot,
                   label='Na, n20, t'+time[itime],ls='-', 
                   xlabel = 'length (m)', ylabel = 'molar densiy [$mol/m^3$]',
                   format = 'png', figname = 'Na_concenOnLengthOverTime', dpi=100))

plot.append(dict(kind='xy', data=['tracking/mix_point_p00000.res'],
                   row=[1,-1],col=[1,11], startplot = True, endplot = True,
                   label='H2O, n20 concentration' ,ls='-', 
                   xlabel = 'time (s)', ylabel = 'molar densiy [$mol/m^3$]',
                   format = 'png', figname = 'H2OOverTime', dpi=100))

plot.append(dict(kind='xy', data=['tracking/mix_point_p00000.res'],
                   row=[1,-1],col=[1,12], startplot = True, endplot = True,
                   label='Na, n20 concentration' , ls='-',
                   xlabel = 'time (s)', ylabel = 'molar densiy [$mol/m^3$]',
                   format = 'png', figname = 'NaOverTime', dpi=100))


