! NAME
! neqDMFT
! DESCRIPTION
!   Run the non-equilibrium DMFT in presence of an external electric field E. 
!   The field is controlled by few flags in the nml/cmd options. It can be 
!   constant, pulse or switched off smoothly. Many more fields can be added by  
!   simply coding them in ELECTRIC_FIELD.f90. The executable read the file
!   *inputFILE.ipt, if not found dump default values to a defualt file.
!
!   The output consist of very few data files that contain all the information,
!   these are eventually read by a second program *get_data_neqDMFT to extract 
!   all the relevant information.
!   In this version the impurity solver is: IPT
! OPTIONS (important)
!  dt=[0.1]            -- Time step for solution of KB equations
!  beta=[100.0]             -- Inverse temperature 
!  U=[6]                    -- Hubbard local interaction value
!  Efield=[0]               -- Strenght of the electric field
!  Vbath=[0]                -- Strenght of the coupling to bath (Lambda=Vbath^2/Wbath)
!  Wbath=[10.0]             -- Bandwidth of the fermionic thermostat
!  ts=[1]                   -- Hopping parameter
!  nstep=[50]               -- Number of time steps: T_max = dt*nstep
!  nloop=[30]               -- Maximum number of DMFT loops allowed (then exit)
!  eps_error=[1.d-4]        -- Tolerance on convergence
!  weight=[0.9]             -- Mixing parameter
!  Nsuccess =[2]            -- Number of consecutive success for convergence to be true
!  Ex=[1]                   -- X-component of the Electric field vector
!  Ey=[1]                   -- Y-component of the Electric field vector
!  t0=[0]                   -- Switching on time parameter for the Electric field
!  t1=[10^6]                -- Switching off time parameter for the Electric field
!  Ncycles=[3]              -- Number of cycles in the  gaussian packect envelope for the impulsive field. fix width.
!  omega0=[pi]            -- Frequency of the of the Oscillating Electric field        
!  E1=[0]                   -- Strenght of the electric field for the AC+DC case, to be tuned to resonate
!  field_type=[dc]       -- Type of electric field profile (dc,ac,ac+dc,etc..)
!  bath_type=[flat]     -- Fermionic thermostat type (constant,gaussian,bethe,etc..)
!  int_method=["trapz"]     -- 
!  data_dir=[DATAneq]       -- Name of the directory containing data files
!  plot_dir=[PLOT]          -- Name of the directory containing plot files
!  fchi=[F]                 -- Flag for the calculation of the optical response
!  L=[1024]                 -- A large number for whatever reason
!  Ltau=[200]               -- A large number for whatever reason
!  P=[5]                    -- Uniform Power mesh power-mesh parameter
!  Q=[5]                    -- Uniform Power mesh uniform-mesh parameter
!  eps=[0.05d0]             -- Broadening on the real-axis
!  Nx=[50]                  -- Size of the K-points grid
!  solve_wfftw =[F]         -- 
!  plot3D=[F]       -- 
!  Lkreduced=[200]  -- 
!  eps=[0.05d0]         -- 
!  irdSFILE=[restartSigma]-- 
!  irdNkFILE=[restartNk]-- 
!#####################################################################