from scipy import *
import shutil

alpha = 0.3

data = io.read_array('imp.0/Ac.inp')

aver = (data[:,1]+data[:,2]+data[:,3])/3.

data[:,1] = data[:,1]*alpha + (1-alpha)*aver
data[:,2] = data[:,2]*alpha + (1-alpha)*aver
data[:,3] = data[:,3]*alpha + (1-alpha)*aver
     
io.write_array("imp.0/Ac.inp_n", data)
shutil.copy2("imp.0/Ac.inp_n", "imp.0/Ac.inp")