clear;

findmodes;
E=getresult("FDE::data::mode1","E");

# Get field and position vectors
Ex=pinch(E.Ex);
Ez=pinch(E.Ez);

x=E.x;
z=E.z;

# Specify the resolution of the plot  
res_x=30;
res_z=20;
 
# Define uniform vectors
x2=linspace(x(1),x(length(x)),res_x);
z2=linspace(z(1),z(length(z)),res_z);
  
xmesh = meshgridx(x2,z2);
zmesh = meshgridy(x2,z2);
 
# Interpolate fields on new uniform vectors
Ex2=interp(Ex, x, z, x2, z2);
Ez2=interp(Ez, x, z, x2, z2);
 
matlabsave("monitor_data",Ex2,Ez2,xmesh,zmesh);