gbe1 = 0.0400928287784106; %eV/ang^2;
ev2j = 1.60218e-19;
gbe2 = gbe1*ev2j*(1e10)^2 ; % J/m^2
gbe3 = gbe2*1e3 %mJ/m^2;