lambda1 1.064[um] "Fundamental wavelength" f1 c_const/lambda1 "Fundamental frequency" f2 2*f1 "Second harmonic frequency" lambda2 c_const/f2 "Second harmonic wavelength" sim_l lambda1*25 "Simulation length" sim_h lambda2/16 "Simulation height" d 1e-18[C/V^2] "Nonlinear coefficient" L "sim_l - 3*lambda1" "Length of nonlinear region" I1 "30[MW/m^2]" "Incident fundamental intensity" E1 sqrt(2*I1/c_const/epsilon0_const) "Incident fundamental electric field strength" offset 1.5*lambda1 "Start of nonlinear region"