def action(index, folder):
dataCurFolder = []
dataCurFolder.append([index, folder])
dp.getTrajactory()
dataCurFolder.append(pps.calculateRMSD())
dataCurFolder.append(dp.getEnergyTemperaturePressure())
return dataCurFolder
def action(index,folder): dataCurFolder = [] dataCurFolder.append([index, folder]) dp.getTrajactory() dataCurFolder.append(pps.calculateRMSD()) dataCurFolder.append(dp.getEnergyTemperaturePressure()) return dataCurFolder
def plotRMSD(ax, label=''):
dp.getTrajactory()
import pyramids.process.struct as pps
time, distance = pps.calculateRMSD()
ax.plot(time, distance, lw=2, label=label)
kargs=ma.getPropertyFromPosition(xlabel='Time (fs)', ylabel=r'$\langle u \rangle^\frac{1}{2}$ ($\AA$)',
title='RMSD')
ma.setProperty(ax,**kargs)
def plotRMSD(ax, label=''):
dp.getTrajactory()
import pyramids.process.struct as pps
time, distance = pps.calculateRMSD()
ax.plot(time, distance, lw=2, label=label)
kargs = ma.getPropertyFromPosition(
xlabel='Time (fs)',
ylabel=r'$\langle u \rangle^\frac{1}{2}$ ($\AA$)',
title='RMSD')
ma.setProperty(ax, **kargs)
#!/usr/bin/python
import numpy as np
from matplotlib import pyplot as plt
import pyramids.io.result as dp
import pyramids.plot.setting as ma
import pyramids.process.struct as pps
fig, axs = plt.subplots(2,1,sharex=True,sharey=False,figsize=(8,6))#
#--------------------------------------------------------------------------------------------
ax = axs[0]
dp.getTrajactory()
time, distance = pps.calculateRMSD()
ax.plot(time, distance, lw=3)
kargs=ma.getPropertyFromPosition(xlabel='Time (fs)', ylabel=r'$\langle u \rangle^\frac{1}{2}$ ($\AA$)',
title='RMSD')
ma.setProperty(ax,**kargs)
#ax.ticklabel_format(style='sci',axis='y',scilimits=[0,0])
#--------------------------------------------------------------------------------------------
ax = axs[1]
time, T, E_ks, E_tot, Vol, P = dp.getEnergyTemperaturePressure()
ax.plot(time, T, lw=3)
kargs=ma.getPropertyFromPosition(xlabel='Time (fs)', ylabel='T (K)',
title='Temperature')
ma.setProperty(ax,**kargs)
#ax.ticklabel_format(style='sci',axis='y',scilimits=[0,0])
#--------------------------------------------------------------------------------------------
#plt.style.use('ggplot')
plt.tight_layout()
SaveName = __file__.split('/')[-1].split('.')[0]