# -*- coding: utf-8 -*-
filename = 'device-pristine.hdf5'
filename_dynmat = 'ff-device-si.hdf5'
filesave = 'device-std-300K.hdf5'
device_configuration = nlread(filename, DeviceConfiguration)[-1]
dynamical_matrix = nlread(filename_dynmat,DynamicalMatrix)[-1]
T = 300

device_configuration_STD = SpecialThermalDisplacement(device_configuration, dynamical_matrix=dynamical_matrix, temperature=T*Kelvin)
device_configuration_STD.update()
nlsave(filesave,device_configuration_STD)