obj = PCIC(hCarrierDensityEvolutionModel)
 
  'ODE Solution Method' = 'Implicit Euler' (default) /
  'Explicit Euler' / 'Const Rho'
 
  'Drude Model Diffraction' = 'off' (default) / 'on'
  in Drude model to calculate inverse Bremsstrahlung for the
  UPPE the divisiion by k_z is performed.
  For high values of k_x, k_y, k_r this can lead to division by
  k_z ~= 0 and loosing model stability. (maybe this can be
  avoided in non paraxial simulation with vectorial UPPE?)
  therefore setting 'Drude Model Diffraction' to 'off'
  (default) makes the code use division by wavevectro length
  (k) instead of k_z (a common assumption).