#!python3
# read the evolution of layer thickness 1 for Bump config
import sys

prfx='/home/abozec/PYTHON/'
sys.path.append(prfx)
import myenv as my

# directory
io='/nexsan/people/abozec/NUMO/MICHAL/BUMP/expt_10.0/data/output/'
idm=41; jdm=41
baclin=2.0

# initial condition
t=8640
imin='    0'
it='{:05d}'.format(t)
file1='archv.000000'+it+'.a'
thk=my.hio.sub_var2(io+file1,idm,jdm,13)
oneta=my.hio.sub_var2(io+file1,idm,jdm,3)
thk=thk*oneta/9806.0
print(thk.shape)
fig, axs = my.plot.subplots(nrows=1,width='10cm')
axs[0].format(title='Interface depth (m) Bump Time:'+imin+'min',ylim=(22,18),xlim=(0,jdm),xlabel='y points')
m = axs[0].plot(thk[:,20],lw=2,color='blue7')
ps_dir=io+'Figures/'
file_ps='bump_intdepth_'+it+'.png'
fig.savefig(ps_dir+file_ps,dpi=100,\
    facecolor='w', edgecolor='w',transparent=False) ## .pdf,.eps
my.plot.close()


#for d in np.arange(30)+2:
t=43230
#for i in [0]:
#for i in my.np.arange(720): ## 12h
for i in my.np.arange(120): ## 2h
  print(t)
  imin='{:5d}'.format(i+1)
  it='{:05d}'.format(t)
  file1='archv.000000'+it+'.a'
  thk=my.hio.sub_var2(io+file1,idm,jdm,13)
  oneta=my.hio.sub_var2(io+file1,idm,jdm,3)
  thk=thk*oneta/9806.0
  print(thk.shape)
  fig, axs = my.plot.subplots(nrows=1,width='10cm')
  axs[0].format(title='Interface depth (m) Bump Time:'+imin+'min',ylim=(22,18),xlim=(0,jdm),xlabel='y points')
  m = axs[0].plot(thk[:,20],lw=2,color='blue7')
  ps_dir=io+'Figures/'
  file_ps='bump_intdepth_'+it+'.png'
  fig.savefig(ps_dir+file_ps,dpi=100,\
        facecolor='w', edgecolor='w',transparent=False) ## .pdf,.eps
  my.plot.close()
  #my.plot.show()
  t=t+30