diff --git a/input/example_user_continue.cfg b/input/example_user_continue.cfg
new file mode 100644
index 0000000000000000000000000000000000000000..4fe7cceb32cb0514da3a88600c1969b8ddc4ccd7
--- /dev/null
+++ b/input/example_user_continue.cfg
@@ -0,0 +1,506 @@
+name = "default";
+param = ( 
+  {
+    type = "string";
+    name = "prog_id";
+    comment = "An identification name for this run";
+    value = "example";
+  }, 
+  {
+    type = "string";
+    name = "output_dir";
+    comment = "Output directory";
+    value = "output";
+  }, 
+  {
+    type = "string";
+    name = "kin_input";
+    comment = "Kinetics input file";
+    value = "input/kinetic_example.cfg";
+  }, 
+  {
+    type = "int";
+    name = "restart";
+    comment = "If restart is 1 (TRUE), the simulation will continue with data from a previous run";
+    value = 0;
+  }, 
+  {
+    type = "string";
+    name = "load_file";
+    comment = "If restart is TRUE, the name of the file with data from previous run, otherwise empty";
+    value = "";
+  }, 
+  {
+    type = "double";
+    name = "output_dt";
+    comment = "Time interval for output to be written to disk";
+    value = 0.1;
+  }, 
+  {
+    type = "int";
+    name = "pois_output";
+    comment = "Output of the Poisson grids, including the potential?";
+    value = 0;
+  }, 
+  {
+    type = "int";
+    name = "cdr_output_margin";
+    comment = "Margin outside the grids in the output of the cdr equation";
+    value = 0;
+  }, 
+  {
+    type = "int";
+    name = "pois_output_margin";
+    comment = "Margin outside the grids in the output of the poisson equation";
+    value = 1;
+  }, 
+  {
+    type = "double";
+    name = "warn_min_timestep";
+    comment = "If the time steps are smaller than this number, the program issues a warning";
+    value = 1e-06;
+  }, 
+  {
+    type = "int";
+    name = "max_disk_space_mb";
+    comment = "Maximum disk space, in Mb, to use";
+    value = 1048576;
+  }, 
+  {
+    type = "int";
+    name = "gridpoints_r";
+    comment = "Number of R gridpoints at level 0";
+    value = 600;
+  }, 
+  {
+    type = "int";
+    name = "gridpoints_z";
+    comment = "Number of Z gridpoints at level 0";
+    value = 600;
+  }, 
+  {
+    type = "int";
+    name = "max_ntheta";
+    comment = "Number of azimuthal gridcells and modes";
+    value = 1;
+  }, 
+  {
+    type = "double";
+    name = "start_t";
+    comment = "Initial time";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "end_t";
+    comment = "End time";
+    value = 0.6;
+  }, 
+  {
+    type = "double";
+    name = "attempt_dt";
+    comment = "Attempted timestep.  The actual timestep may be larger";
+    value = 50.0;
+  }, 
+  {
+    type = "int";
+    name = "extra_pois_levels";
+    comment = "Extra levels for the Poisson solver";
+    value = 2;
+  }, 
+  {
+    type = "int";
+    name = "max_levels";
+    comment = "Maximum level of refinement. Use a big number here";
+    value = 64;
+  }, 
+  {
+    type = "double";
+    name = "pois_max_error";
+    comment = "Error threshold that leads to refinement in the Poisson code.";
+    value = 0.001;
+  }, 
+  {
+    type = "int";
+    name = "pois_max_level";
+    comment = "Maximum level of refinement in the Poisson equation.";
+    value = 3;
+  }, 
+  {
+    type = "int";
+    name = "extra_photo_levels";
+    comment = "Extra levels for the photo-ionization solver";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_max_level";
+    comment = "Maximum level of refinement in the photo-ionization solver.";
+    value = 4;
+  }, 
+  {
+    type = "double";
+    name = "photo_max_error";
+    comment = "Error threshold that leads to refinement in the photo-ionization code.";
+    value = 0.01;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_right";
+    comment = "Photo-ionization boundary condition at r = L_r.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_bottom";
+    comment = "Photo-ionization boundary condition at z = 0.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_top";
+    comment = "Photo-ionization boundary condition at z = L_z.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "extra_photo_levels_2";
+    comment = "Extra levels for the photo-ionization solver";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_max_level_2";
+    comment = "Maximum level of refinement in the photo-ionization solver.";
+    value = 4;
+  }, 
+  {
+    type = "double";
+    name = "photo_max_error_2";
+    comment = "Error threshold that leads to refinement in the photo-ionization code.";
+    value = 0.01;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_right_2";
+    comment = "Photo-ionization boundary condition at r = L_r.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_bottom_2";
+    comment = "Photo-ionization boundary condition at z = 0.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "photo_bnd_top_2";
+    comment = "Photo-ionization boundary condition at z = L_z.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "cdr_bnd_bottom";
+    comment = "Particles boundary condition at z = 0.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = 1;
+  }, 
+  {
+    type = "int";
+    name = "cdr_bnd_top";
+    comment = "Particles boundary condition at z = L_z.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = 1;
+  }, 
+  {
+    type = "int";
+    name = "cdr_bnd_right";
+    comment = "Particles boundary condition at r = L_r.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = 1;
+  }, 
+  {
+    type = "int";
+    name = "pois_bnd_right";
+    comment = "Potential boundary condition at r = L_r.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "pois_bnd_bottom";
+    comment = "Potential boundary condition at z = 0.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "int";
+    name = "pois_bnd_top";
+    comment = "Potential boundary condition at z = L_z.  1 for Hom. Neumann, -1 for Hom. Dirichlet";
+    value = -1;
+  }, 
+  {
+    type = "double";
+    name = "nu_a";
+    comment = "Maximum advection Courant number";
+    value = 0.2;
+  }, 
+  {
+    type = "double";
+    name = "nu_d";
+    comment = "Maximum diffusion Courant number";
+    value = 0.2;
+  }, 
+  {
+    type = "double";
+    name = "nu_rt";
+    comment = "Maximum ratio of dt/relaxation time";
+    value = 0.2;
+  }, 
+  {
+    type = "double";
+    name = "nu_f";
+    comment = "Maximum ratio of change of the densities (set to a very large number to ignore)";
+    value = 1e+20;
+  }, 
+  {
+    type = "double";
+    name = "ref_threshold_eabs";
+    comment = "Refinement threshold for the electric field";
+    value = 0.2;
+  }, 
+  {
+    type = "int";
+    name = "ref_level_eabs";
+    comment = "Maximum refinement level reached through ref_threshold_eabs";
+    value = 4;
+  }, 
+  {
+    type = "double";
+    name = "ref_threshold_charge";
+    comment = "Refinement threshold for the curvature of the charge";
+    value = 0.004;
+  }, 
+  {
+    type = "double";
+    name = "ref_threshold_dens";
+    comment = "Refinement threshold for the curvature of the densities";
+    value = 0.004;
+  }, 
+  {
+    type = "double";
+    name = "ref_threshold_edge";
+    comment = "Refinement threshold for the densities in the leading edge";
+    value = 10000.0;
+  }, 
+  {
+    type = "int";
+    name = "cdr_brick_dr";
+    comment = "r-length of the minimal refinement area in the cdr equation";
+    value = 8;
+  }, 
+  {
+    type = "int";
+    name = "cdr_brick_dz";
+    comment = "z-length of the minimal refinement area in the cdr equation";
+    value = 8;
+  }, 
+  {
+    type = "int";
+    name = "cdr_max_level";
+    comment = "Maximum level of refinement in the Fluid equation.";
+    value = 3;
+  }, 
+  {
+    type = "int";
+    name = "cdr_interp_in";
+    comment = "Interpolation method for the grid interior (0=zero_masses, 1=quadratic_masses [default], 2=wackers_masses, 3=quadlog";
+    value = 1;
+  }, 
+  {
+    type = "int";
+    name = "cdr_interp_bnd";
+    comment = "Interpolation method for the grid boundaries (0=zero_masses, 1=quadratic_masses [default], 2=wackers_masses, 3=quadlog";
+    value = 1;
+  }, 
+  {
+    type = "double";
+    name = "L_r";
+    comment = "Length in r of the complete domain";
+    value = 13044.0;
+  }, 
+  {
+    type = "double";
+    name = "L_z";
+    comment = "Length in z of the complete domain";
+    value = 13044.0;
+  }, 
+  {
+    type = "double";
+    name = "diffusion_coeff";
+    comment = "Isotropic difussion coefficient";
+    value = 0.1;
+  }, 
+  {
+    type = "int";
+    name = "has_photoionization";
+    comment = "Whether the code includes photoionization or not";
+    value = 1;
+  }, 
+  {
+    type = "string";
+    name = "photoionization_file";
+    comment = "The name of a file from which we can read the photoionization parameters";
+    value = "input/air760torr.photo";
+  }, 
+  {
+    type = "double";
+    name = "attachment_rate";
+    comment = "Rate of dissociative attachment";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "attachment_E0";
+    comment = "E0 in the exp(-E0/E) factor in the attachment expression.";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "E0_x";
+    comment = "x component of the external electric field";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "E0_y";
+    comment = "y component of the external electric field";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "E0_z";
+    comment = "z component of the external electric field";
+    value = -0.06;
+  }, 
+  {
+    type = "double";
+    name = "rise_time";
+    comment = "Rise time of the electric field (0 for instantaneous rise)";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "off_time";
+    comment = "Time to switch off the electric field (0.0 means never)";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "seed_sigma_x";
+    comment = "x width of the initial seed";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "seed_sigma_y";
+    comment = "y width of the initial seed";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "seed_sigma_z";
+    comment = "z width of the initial seed";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "seed_N";
+    comment = "Number of electrons in the initial seed";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "background_ionization";
+    comment = "Initial at z=0 densities of electrons and ions";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "background_increase_length";
+    comment = "Length of exponential increase of the pre-ionization (for atmospherical models)";
+    value = 0.0;
+  }, 
+  {
+    type = "int";
+    name = "pois_inhom";
+    comment = "Use the point-plane geometry?";
+    value = 1;
+  }, 
+  {
+    type = "int";
+    name = "pois_inhom_reflections";
+    comment = "Number of mirror charges to use";
+    value = 4;
+  }, 
+  {
+    type = "double";
+    name = "needle_length";
+    comment = "Length of the needle";
+    value = 2500.0;
+  }, 
+  {
+    type = "double";
+    name = "needle_radius";
+    comment = "Radius of the needle";
+    value = 400.0;
+  }, 
+  {
+    type = "double";
+    name = "pois_inhom_fixed_q";
+    comment = "If nonzero, the charge is fixed, not floating (simulation of charged clouds close to the earth surface)";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "constant_source";
+    comment = "Constant ionization rate";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "perturb_epsilon";
+    comment = "Initial perturbation to the axisymmetric configuration";
+    value = 0.0;
+  }, 
+  {
+    type = "int";
+    name = "perturb_max_k";
+    comment = "Perturb only modes up to perturb_max_k (large number to perturb all)";
+    value = 1024;
+  }, 
+  {
+    type = "int";
+    name = "sprite_module";
+    comment = "1 if the sprite module is activated, 0 otherwise";
+    value = 0;
+  }, 
+  {
+    type = "double";
+    name = "dens_decay_len";
+    comment = "Lenght of exponential decay of the density w/r to altitude";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "sprite_dens_0";
+    comment = "Density at z = 0";
+    value = 0.0;
+  }, 
+  {
+    type = "double";
+    name = "sprite_dens_q";
+    comment = "Quenching density";
+    value = 0.0;
+  }, 
+  {
+    type = "int";
+    name = "sprite_sign";
+    comment = "Sign of the sprite head that we are following (the other will not be reliable";
+    value = -1;
+  } );