xdim = 1200;
ydim = 1200;
buffer = zeros(xdim,ydim);
a = 2;
f = 19+(a-1)*32;
IC_name = "IC_[-3,3][-3,3]" + xdim + "X" + ydim + "f" + f + "O80";
AT_name = "dat_files/AT[-3,3][-3,3]_" + xdim + "X" + ydim + "_f" + f + "_O80";
filename_IC = IC_name + ".dat";
IC = importdata(filename_IC);