- FPMD/CP more merging

- new module for merged FPMD/CP restart subroutines - starting work on XML restart/punch (cp_restart.f90) carlo git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@1443 c92efa57-630b-4861-b058-cf58834340f0
2004-11-15 12:13:22 +00:00 · 2004-11-15 12:13:22 +00:00 · 3b22ff0108
parent 4aee60f71e
commit 3b22ff0108
28 changed files with 764 additions and 509 deletions
--- a/CPV/cpr.f90
+++ b/CPV/cpr.f90
@ -80,7 +80,7 @@
      use reciprocal_vectors, only: gstart
      use ions_base, only: na, nat, pmass, nas => nax, nsp, rcmax
      use ions_base, only: ind_srt, ions_cofmass, ions_kinene, ions_temp, ions_thermal_stress
-      use ions_base, only: ions_vrescal, fricp, greasp
+      use ions_base, only: ions_vrescal, fricp, greasp, iforce
      use grid_dimensions, only: nnr => nnrx, nr1, nr2, nr3
      use cell_base, only: ainv, a1, a2, a3, frich, greash
      use cell_base, only: omega, alat, ibrav, celldm
@ -125,7 +125,7 @@
      use cp_electronic_mass, only: emass_precond
      use cpr_subroutines
      use ions_positions, only: tau0, taum, taup, taus, tausm, tausp, vels, velsm, velsp
-      use ions_positions, only: ions_hmove, ions_move, iforce
+      use ions_positions, only: ions_hmove, ions_move
      use ions_nose, only: gkbt, qnp, vnhp, xnhp0, xnhpm, xnhpp, ions_nosevel, &
                           ions_noseupd, tempw
      use electrons_nose, only: qne, ekincw
--- a/CPV/input.f90
+++ b/CPV/input.f90
@ -88,14 +88,11 @@ CONTAINS
           ekin_conv_thr, etot_conv_thr, max_seconds, na_inp, rd_pos, atom_label, rd_vel, &
           smd_polm, smd_kwnp, smd_linr, smd_stcd, smd_stcd1, smd_stcd2, smd_stcd3, smd_codf, &
           smd_forf, smd_smwf, smd_lmfreq, smd_tol, smd_maxlm, smd_smcp, smd_smopt, smd_smlm, &
-           num_of_images, smd_ene_ini, smd_ene_fin
-
-      use read_namelists_module, only: read_namelists
-      use read_cards_module, only: read_cards
+           num_of_images, smd_ene_ini, smd_ene_fin, title

      use constants, only: pi, scmass, factem, eps8, uma_au, terahertz 

-      use parameters, only: nsx, natx, nbndxx
+      use parameters, only: natx

      use io_global, only: ionode, stdout

@ -140,7 +137,6 @@ CONTAINS
           fricp_ => fricp, &
           greasp_ => greasp
      USE ions_positions, ONLY: &
-           iforce_ => iforce, &
           tau0_ => tau0
      ! 
      USE cell_base, ONLY: cell_base_init, a1, a2, a3, &
@ -210,6 +206,8 @@ CONTAINS
           maxlm_ => smd_maxlm, &
           ene_ini_ => smd_ene_ini, &
           ene_fin_ => smd_ene_fin
+      USE printout_base, ONLY: &
+           title_ => title

      !
      implicit none
@ -218,7 +216,7 @@ CONTAINS
      ! local variables
      !

-      real(kind=8) :: taus( 3, natx ), ocp, fsum
+      real(kind=8) :: ocp, fsum
      integer :: i, ia, is, iss, in, isa
      real(kind=8) :: alat_

@ -226,6 +224,8 @@ CONTAINS
      ! Subroutine body
      !

+      title_ = title   ! simulation title
+
      IF( TRIM( calculation ) == 'nscf' ) trhor_ = .true.
     
      ! 
@ -284,10 +284,6 @@ CONTAINS
        CALL ions_base_init( ntyp , nat , na_inp , sp_pos , rd_pos , rd_vel, atom_mass, &
             atom_label, if_pos, atomic_positions , alat_ , a1, a2, a3  )
      end if
-      iforce_ = 0
-      DO isa = 1, nat
-        iforce_ ( :, isa ) = if_pos( :, ind_srt( isa ) ) 
-      END DO

      ! ...   Set Values for bands and spin

--- a/CPV/make.depend
+++ b/CPV/make.depend
@ -5,12 +5,13 @@ compute_scf.o : ../Modules/input_parameters.o
 compute_scf.o : ../Modules/io_files.o
 compute_scf.o : ../Modules/io_global.o
 compute_scf.o : ../Modules/kind.o
-compute_scf.o : ../Modules/mp.o
 compute_scf.o : ../Modules/mp_global.o
+compute_scf.o : ../Modules/mp.o
 compute_scf.o : ../Modules/parser.o
 compute_scf.o : ../Modules/path_formats.o
 compute_scf.o : ../Modules/path_variables.o
 compute_scf.o : restart.o
+cplib.o : fft_cp.o
 cplib.o : ../Modules/atom.o
 cplib.o : ../Modules/cell_base.o
 cplib.o : ../Modules/constants.o
@ -26,14 +27,19 @@ cplib.o : ../Modules/io_global.o
 cplib.o : ../Modules/ions_base.o
 cplib.o : ../Modules/kind.o
 cplib.o : ../Modules/mp.o
+cplib.o : modules.o
 cplib.o : ../Modules/parameters.o
 cplib.o : ../Modules/recvec.o
 cplib.o : ../Modules/smallbox.o
 cplib.o : ../Modules/stick_base.o
 cplib.o : ../Modules/uspp.o
-cplib.o : fft_cp.o
-cplib.o : modules.o
 cplib.o : para.o
+cpr_mod.o : ../Modules/io_global.o
+cpr_mod.o : ../Modules/mp.o
+cpr_mod.o : ../Modules/recvec.o
+cpr_mod.o : ../Modules/wave_base.o
+cpr.o : cpr_mod.o
+cpr.o : input.o
 cpr.o : ../Modules/cell_base.o
 cpr.o : ../Modules/check_stop.o
 cpr.o : ../Modules/constants.o
@ -47,8 +53,9 @@ cpr.o : ../Modules/input_parameters.o
 cpr.o : ../Modules/io_files.o
 cpr.o : ../Modules/io_global.o
 cpr.o : ../Modules/ions_base.o
-cpr.o : ../Modules/mp.o
 cpr.o : ../Modules/mp_global.o
+cpr.o : ../Modules/mp.o
+cpr.o : modules.o
 cpr.o : ../Modules/parameters.o
 cpr.o : ../Modules/printout_base.o
 cpr.o : ../Modules/recvec.o
@ -56,31 +63,25 @@ cpr.o : ../Modules/timestep.o
 cpr.o : ../Modules/uspp.o
 cpr.o : ../Modules/wave_base.o
 cpr.o : ../Modules/wavefunctions.o
-cpr.o : cpr_mod.o
-cpr.o : input.o
-cpr.o : modules.o
 cpr.o : para.o
 cpr.o : restart.o
 cpr.o : restart_sub.o
 cpr.o : wannier.o
-cpr_mod.o : ../Modules/io_global.o
-cpr_mod.o : ../Modules/mp.o
-cpr_mod.o : ../Modules/recvec.o
-cpr_mod.o : ../Modules/wave_base.o
+cprstart.o : input.o
 cprstart.o : ../Modules/control_flags.o
 cprstart.o : ../Modules/input_parameters.o
 cprstart.o : ../Modules/io_files.o
 cprstart.o : ../Modules/io_global.o
 cprstart.o : ../Modules/ions_base.o
 cprstart.o : ../Modules/kind.o
-cprstart.o : ../Modules/mp.o
 cprstart.o : ../Modules/mp_global.o
+cprstart.o : ../Modules/mp.o
 cprstart.o : ../Modules/path_base.o
 cprstart.o : ../Modules/path_io_routines.o
 cprstart.o : ../Modules/path_variables.o
-cprstart.o : input.o
 cprstart.o : para.o
 cprstart.o : path_routines.o
+cprsub.o : cpr_mod.o
 cprsub.o : ../Modules/atom.o
 cprsub.o : ../Modules/cell_base.o
 cprsub.o : ../Modules/constants.o
@ -91,12 +92,11 @@ cprsub.o : ../Modules/griddim.o
 cprsub.o : ../Modules/io_global.o
 cprsub.o : ../Modules/ions_base.o
 cprsub.o : ../Modules/kind.o
+cprsub.o : modules.o
 cprsub.o : ../Modules/parameters.o
 cprsub.o : ../Modules/recvec.o
 cprsub.o : ../Modules/smallbox.o
 cprsub.o : ../Modules/uspp.o
-cprsub.o : cpr_mod.o
-cprsub.o : modules.o
 cprsub.o : restart.o
 errore.o : ../Modules/io_global.o
 errore.o : ../Modules/parallel_include.o
@ -119,6 +119,7 @@ input.o : ../Modules/ions_base.o
 input.o : ../Modules/mp.o
 input.o : ../Modules/parameters.o
 input.o : ../Modules/path_variables.o
+input.o : ../Modules/printout_base.o
 input.o : ../Modules/read_cards.o
 input.o : ../Modules/read_namelists.o
 input.o : ../Modules/recvec.o
@ -147,8 +148,8 @@ path_routines.o : ../Modules/input_parameters.o
 path_routines.o : ../Modules/io_global.o
 path_routines.o : ../Modules/ions_base.o
 path_routines.o : ../Modules/kind.o
-path_routines.o : ../Modules/mp.o
 path_routines.o : ../Modules/mp_global.o
+path_routines.o : ../Modules/mp.o
 path_routines.o : ../Modules/parser.o
 path_routines.o : ../Modules/path_base.o
 path_routines.o : ../Modules/path_variables.o
@ -170,25 +171,27 @@ restart.o : ../Modules/input_parameters.o
 restart.o : ../Modules/io_base.o
 restart.o : ../Modules/io_global.o
 restart.o : ../Modules/ions_base.o
-restart.o : ../Modules/mp.o
 restart.o : ../Modules/mp_global.o
+restart.o : ../Modules/mp.o
+restart.o : modules.o
 restart.o : ../Modules/parameters.o
 restart.o : ../Modules/parser.o
 restart.o : ../Modules/recvec.o
-restart.o : modules.o
+restart_sub.o : cpr_mod.o
 restart_sub.o : ../Modules/cell_base.o
 restart_sub.o : ../Modules/control_flags.o
 restart_sub.o : ../Modules/cp_emass.o
 restart_sub.o : ../Modules/cp_fpmd.o
 restart_sub.o : ../Modules/io_global.o
 restart_sub.o : ../Modules/ions_base.o
+restart_sub.o : modules.o
 restart_sub.o : ../Modules/recvec.o
 restart_sub.o : ../Modules/uspp.o
 restart_sub.o : ../Modules/wave_base.o
 restart_sub.o : ../Modules/wavefunctions.o
-restart_sub.o : cpr_mod.o
-restart_sub.o : modules.o
 restart_sub.o : wannier.o
+smcp.o : cpr_mod.o
+smcp.o : input.o
 smcp.o : ../Modules/atom.o
 smcp.o : ../Modules/cell_base.o
 smcp.o : ../Modules/check_stop.o
@ -203,34 +206,32 @@ smcp.o : ../Modules/input_parameters.o
 smcp.o : ../Modules/io_files.o
 smcp.o : ../Modules/io_global.o
 smcp.o : ../Modules/ions_base.o
-smcp.o : ../Modules/mp.o
 smcp.o : ../Modules/mp_global.o
+smcp.o : ../Modules/mp.o
+smcp.o : modules.o
 smcp.o : ../Modules/parameters.o
 smcp.o : ../Modules/path_variables.o
 smcp.o : ../Modules/recvec.o
 smcp.o : ../Modules/timestep.o
 smcp.o : ../Modules/uspp.o
 smcp.o : ../Modules/wave_base.o
-smcp.o : cpr_mod.o
-smcp.o : input.o
-smcp.o : modules.o
 smcp.o : para.o
 smcp.o : restart.o
 smcp.o : smd_modules.o
+smd_modules.o : ../Modules/parameters.o
+smd.o : cpr_mod.o
 smd.o : ../Modules/cell_base.o
 smd.o : ../Modules/constants.o
 smd.o : ../Modules/control_flags.o
 smd.o : ../Modules/input_parameters.o
 smd.o : ../Modules/io_global.o
 smd.o : ../Modules/ions_base.o
+smd.o : modules.o
 smd.o : ../Modules/parameters.o
 smd.o : ../Modules/path_variables.o
 smd.o : ../Modules/recvec.o
-smd.o : cpr_mod.o
-smd.o : modules.o
 smd.o : restart.o
 smd.o : smd_modules.o
-smd_modules.o : ../Modules/parameters.o
 smlam.o : ../Modules/ions_base.o
 smlam.o : ../Modules/parameters.o
 smlam.o : ../Modules/path_variables.o
@ -241,9 +242,9 @@ wannier.o : ../Modules/griddim.o
 wannier.o : ../Modules/io_global.o
 wannier.o : ../Modules/ions_base.o
 wannier.o : ../Modules/mp.o
+wannier.o : modules.o
 wannier.o : ../Modules/recvec.o
 wannier.o : ../Modules/wave_base.o
-wannier.o : modules.o
 wannier.o : restart.o
 wf.o : ../Modules/cell_base.o
 wf.o : ../Modules/constants.o
@ -253,12 +254,12 @@ wf.o : ../Modules/fft_scalar.o
 wf.o : ../Modules/griddim.o
 wf.o : ../Modules/io_global.o
 wf.o : ../Modules/ions_base.o
+wf.o : modules.o
 wf.o : ../Modules/parallel_include.o
 wf.o : ../Modules/parameters.o
 wf.o : ../Modules/recvec.o
 wf.o : ../Modules/smallbox.o
 wf.o : ../Modules/uspp.o
-wf.o : modules.o
 wf.o : para.o
 wf.o : wannier.o
 wrapper.o : ../Modules/io_global.o
--- a/CPV/smcp.f90
+++ b/CPV/smcp.f90
@ -78,7 +78,7 @@ subroutine smdmain( tau, fion_out, etot_out, nat_out )
  use reciprocal_vectors, only: ng0 => gstart
  use ions_base, only: na, nat, pmass, nas => nax, nsp, rcmax
  use ions_base, only: ind_srt, ions_vel, ions_cofmass, ions_kinene, ions_temp
-  use ions_base, only: ions_thermal_stress, ions_vrescal, fricp, greasp
+  use ions_base, only: ions_thermal_stress, ions_vrescal, fricp, greasp, iforce
  use grid_dimensions, only: nnr => nnrx, nr1, nr2, nr3
  use cell_base, only: ainv, a1, a2, a3, r_to_s, celldm, ibrav
  use cell_base, only: omega, alat, frich, greash, press
@ -117,7 +117,7 @@ subroutine smdmain( tau, fion_out, etot_out, nat_out )
  USE control_flags, ONLY : conv_elec, tconvthrs
  USE check_stop, ONLY : check_stop_now
  USE cpr_subroutines
-  use ions_positions, only: tau0, velsp, iforce
+  use ions_positions, only: tau0, velsp
  use ions_positions, only: ions_hmove, ions_move
  use ions_nose, only: gkbt, qnp, ions_nosevel, ions_noseupd, tempw
  USE cell_base, ONLY: cell_kinene, cell_move, cell_gamma, cell_hmove
--- a/D3/Makefile
+++ b/D3/Makefile
@ -190,6 +190,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/readpseudo.o \
 ../Modules/recvec.o \
--- a/D3/make.depend
+++ b/D3/make.depend
@ -1,22 +1,23 @@
+allocate_d3.o : d3com.o
 allocate_d3.o : ../Modules/ions_base.o
 allocate_d3.o : ../PH/phcom.o
 allocate_d3.o : ../PW/pwcom.o
-allocate_d3.o : d3com.o
+bcast_d3_input.o : d3com.o
 bcast_d3_input.o : ../Modules/control_flags.o
 bcast_d3_input.o : ../Modules/io_files.o
 bcast_d3_input.o : ../Modules/mp.o
 bcast_d3_input.o : ../PH/phcom.o
 bcast_d3_input.o : ../PW/pwcom.o
-bcast_d3_input.o : d3com.o
 ch_psi_all2.o : ../Modules/kind.o
 ch_psi_all2.o : ../PH/phcom.o
 ch_psi_all2.o : ../PW/becmod.o
 ch_psi_all2.o : ../PW/pwcom.o
+close_open.o : d3com.o
 close_open.o : ../Modules/io_files.o
 close_open.o : ../Modules/io_global.o
 close_open.o : ../PH/phcom.o
 close_open.o : ../PW/pwcom.o
-close_open.o : d3com.o
+d0rhod2v.o : d3com.o
 d0rhod2v.o : ../Modules/io_files.o
 d0rhod2v.o : ../Modules/io_global.o
 d0rhod2v.o : ../Modules/ions_base.o
@ -26,35 +27,50 @@ d0rhod2v.o : ../Modules/uspp.o
 d0rhod2v.o : ../Modules/wavefunctions.o
 d0rhod2v.o : ../PH/phcom.o
 d0rhod2v.o : ../PW/pwcom.o
-d0rhod2v.o : d3com.o
 d2mxc.o : ../Modules/kind.o
+d3com.o : ../Modules/kind.o
+d3dyn_cc.o : d3com.o
+d3dyn_cc.o : ../Modules/ions_base.o
+d3dyn_cc.o : ../Modules/kind.o
+d3dyn_cc.o : ../PH/phcom.o
+d3dyn_cc.o : ../PW/pwcom.o
+d3_exc.o : d3com.o
 d3_exc.o : ../Modules/io_global.o
 d3_exc.o : ../Modules/ions_base.o
 d3_exc.o : ../Modules/kind.o
-d3_exc.o : ../Modules/mp.o
 d3_exc.o : ../Modules/mp_global.o
+d3_exc.o : ../Modules/mp.o
 d3_exc.o : ../PH/phcom.o
 d3_exc.o : ../PW/pwcom.o
-d3_exc.o : d3com.o
+d3_init.o : d3com.o
 d3_init.o : ../Modules/atom.o
 d3_init.o : ../Modules/ions_base.o
 d3_init.o : ../Modules/mp.o
 d3_init.o : ../Modules/uspp.o
 d3_init.o : ../PH/phcom.o
 d3_init.o : ../PW/pwcom.o
-d3_init.o : d3com.o
+d3ionq.o : ../Modules/constants.o
+d3ionq.o : ../Modules/io_global.o
+d3ionq.o : ../Modules/kind.o
+d3matrix.o : d3com.o
+d3matrix.o : ../Modules/control_flags.o
+d3matrix.o : ../Modules/ions_base.o
+d3matrix.o : ../Modules/kind.o
+d3matrix.o : ../PH/phcom.o
+d3matrix.o : ../PW/pwcom.o
+d3_readin.o : d3com.o
 d3_readin.o : ../Modules/control_flags.o
 d3_readin.o : ../Modules/io_files.o
 d3_readin.o : ../Modules/io_global.o
 d3_readin.o : ../Modules/ions_base.o
 d3_readin.o : ../PH/phcom.o
 d3_readin.o : ../PW/pwcom.o
-d3_readin.o : d3com.o
+d3_recover.o : d3com.o
 d3_recover.o : ../Modules/io_global.o
 d3_recover.o : ../Modules/mp.o
 d3_recover.o : ../PH/phcom.o
 d3_recover.o : ../PW/pwcom.o
-d3_recover.o : d3com.o
+d3_setup.o : d3com.o
 d3_setup.o : ../Modules/atom.o
 d3_setup.o : ../Modules/constants.o
 d3_setup.o : ../Modules/control_flags.o
@ -64,7 +80,7 @@ d3_setup.o : ../Modules/kind.o
 d3_setup.o : ../Modules/mp_global.o
 d3_setup.o : ../PH/phcom.o
 d3_setup.o : ../PW/pwcom.o
-d3_setup.o : d3com.o
+d3_summary.o : d3com.o
 d3_summary.o : ../Modules/atom.o
 d3_summary.o : ../Modules/control_flags.o
 d3_summary.o : ../Modules/io_global.o
@ -73,30 +89,10 @@ d3_summary.o : ../Modules/kind.o
 d3_summary.o : ../Modules/uspp.o
 d3_summary.o : ../PH/phcom.o
 d3_summary.o : ../PW/pwcom.o
-d3_summary.o : d3com.o
 d3_symdyn.o : ../Modules/kind.o
 d3_symdynph.o : ../Modules/constants.o
 d3_symdynph.o : ../Modules/kind.o
-d3_valence.o : ../Modules/ions_base.o
-d3_valence.o : ../Modules/kind.o
-d3_valence.o : ../PH/phcom.o
-d3_valence.o : ../PW/pwcom.o
-d3_valence.o : d3com.o
-d3com.o : ../Modules/kind.o
-d3dyn_cc.o : ../Modules/ions_base.o
-d3dyn_cc.o : ../Modules/kind.o
-d3dyn_cc.o : ../PH/phcom.o
-d3dyn_cc.o : ../PW/pwcom.o
-d3dyn_cc.o : d3com.o
-d3ionq.o : ../Modules/constants.o
-d3ionq.o : ../Modules/io_global.o
-d3ionq.o : ../Modules/kind.o
-d3matrix.o : ../Modules/control_flags.o
-d3matrix.o : ../Modules/ions_base.o
-d3matrix.o : ../Modules/kind.o
-d3matrix.o : ../PH/phcom.o
-d3matrix.o : ../PW/pwcom.o
-d3matrix.o : d3com.o
+d3toten.o : d3com.o
 d3toten.o : ../Modules/io_files.o
 d3toten.o : ../Modules/io_global.o
 d3toten.o : ../Modules/ions_base.o
@ -104,7 +100,12 @@ d3toten.o : ../Modules/kind.o
 d3toten.o : ../Modules/version.o
 d3toten.o : ../PH/phcom.o
 d3toten.o : ../PW/pwcom.o
-d3toten.o : d3com.o
+d3_valence.o : d3com.o
+d3_valence.o : ../Modules/ions_base.o
+d3_valence.o : ../Modules/kind.o
+d3_valence.o : ../PH/phcom.o
+d3_valence.o : ../PW/pwcom.o
+d3vrho.o : d3com.o
 d3vrho.o : ../Modules/io_files.o
 d3vrho.o : ../Modules/ions_base.o
 d3vrho.o : ../Modules/kind.o
@ -112,29 +113,29 @@ d3vrho.o : ../Modules/uspp.o
 d3vrho.o : ../Modules/wavefunctions.o
 d3vrho.o : ../PH/phcom.o
 d3vrho.o : ../PW/pwcom.o
-d3vrho.o : d3com.o
 davcio_drho2.o : ../Modules/io_global.o
 davcio_drho2.o : ../Modules/kind.o
-davcio_drho2.o : ../Modules/mp.o
 davcio_drho2.o : ../Modules/mp_global.o
+davcio_drho2.o : ../Modules/mp.o
 davcio_drho2.o : ../PH/phcom.o
 davcio_drho2.o : ../PW/para.o
 davcio_drho2.o : ../PW/pwcom.o
+dpsi_corr.o : d3com.o
 dpsi_corr.o : ../Modules/kind.o
 dpsi_corr.o : ../PH/phcom.o
 dpsi_corr.o : ../PW/pwcom.o
-dpsi_corr.o : d3com.o
+dpsidpsidv.o : d3com.o
 dpsidpsidv.o : ../Modules/ions_base.o
 dpsidpsidv.o : ../Modules/kind.o
 dpsidpsidv.o : ../PH/phcom.o
 dpsidpsidv.o : ../PW/pwcom.o
-dpsidpsidv.o : d3com.o
+dpsidvdpsi.o : d3com.o
 dpsidvdpsi.o : ../Modules/io_files.o
 dpsidvdpsi.o : ../Modules/ions_base.o
 dpsidvdpsi.o : ../Modules/kind.o
 dpsidvdpsi.o : ../PH/phcom.o
 dpsidvdpsi.o : ../PW/pwcom.o
-dpsidvdpsi.o : d3com.o
+dqrhod2v.o : d3com.o
 dqrhod2v.o : ../Modules/io_files.o
 dqrhod2v.o : ../Modules/ions_base.o
 dqrhod2v.o : ../Modules/kind.o
@ -143,51 +144,51 @@ dqrhod2v.o : ../Modules/uspp.o
 dqrhod2v.o : ../Modules/wavefunctions.o
 dqrhod2v.o : ../PH/phcom.o
 dqrhod2v.o : ../PW/pwcom.o
-dqrhod2v.o : d3com.o
+drho_cc.o : d3com.o
 drho_cc.o : ../Modules/kind.o
 drho_cc.o : ../PH/phcom.o
 drho_cc.o : ../PW/pwcom.o
-drho_cc.o : d3com.o
+drhod2v.o : d3com.o
+drhod2v.o : ../Modules/ions_base.o
+drhod2v.o : ../Modules/kind.o
+drhod2v.o : ../PH/phcom.o
+drhod2v.o : ../PW/pwcom.o
+drho_drc.o : d3com.o
 drho_drc.o : ../Modules/atom.o
 drho_drc.o : ../Modules/ions_base.o
 drho_drc.o : ../Modules/kind.o
 drho_drc.o : ../Modules/mp.o
 drho_drc.o : ../PH/phcom.o
 drho_drc.o : ../PW/pwcom.o
-drho_drc.o : d3com.o
-drhod2v.o : ../Modules/ions_base.o
-drhod2v.o : ../Modules/kind.o
-drhod2v.o : ../PH/phcom.o
-drhod2v.o : ../PW/pwcom.o
-drhod2v.o : d3com.o
+dvdpsi.o : d3com.o
 dvdpsi.o : ../Modules/ions_base.o
 dvdpsi.o : ../Modules/uspp.o
 dvdpsi.o : ../PH/phcom.o
 dvdpsi.o : ../PW/pwcom.o
-dvdpsi.o : d3com.o
+dvscf.o : d3com.o
 dvscf.o : ../Modules/atom.o
 dvscf.o : ../Modules/ions_base.o
 dvscf.o : ../Modules/kind.o
 dvscf.o : ../PH/phcom.o
 dvscf.o : ../PW/pwcom.o
-dvscf.o : d3com.o
+gen_dpdvp.o : d3com.o
 gen_dpdvp.o : ../Modules/io_files.o
 gen_dpdvp.o : ../Modules/ions_base.o
 gen_dpdvp.o : ../Modules/kind.o
 gen_dpdvp.o : ../Modules/wavefunctions.o
 gen_dpdvp.o : ../PH/phcom.o
 gen_dpdvp.o : ../PW/pwcom.o
-gen_dpdvp.o : d3com.o
+gen_dwfc.o : d3com.o
 gen_dwfc.o : ../Modules/io_global.o
 gen_dwfc.o : ../PH/phcom.o
 gen_dwfc.o : ../PW/pwcom.o
-gen_dwfc.o : d3com.o
 incdrhoscf2.o : ../Modules/ions_base.o
 incdrhoscf2.o : ../Modules/kind.o
 incdrhoscf2.o : ../Modules/uspp.o
 incdrhoscf2.o : ../Modules/wavefunctions.o
 incdrhoscf2.o : ../PH/phcom.o
 incdrhoscf2.o : ../PW/pwcom.o
+openfild3.o : d3com.o
 openfild3.o : ../Modules/control_flags.o
 openfild3.o : ../Modules/io_files.o
 openfild3.o : ../Modules/io_global.o
@ -195,37 +196,37 @@ openfild3.o : ../Modules/mp_global.o
 openfild3.o : ../PH/phcom.o
 openfild3.o : ../PW/pwcom.o
 openfild3.o : ../PW/restart.o
-openfild3.o : d3com.o
-print_clock_d3.o : ../Modules/io_global.o
 print_clock_d3.o : d3com.o
+print_clock_d3.o : ../Modules/io_global.o
+psymd0rho.o : d3com.o
 psymd0rho.o : ../Modules/ions_base.o
 psymd0rho.o : ../Modules/kind.o
 psymd0rho.o : ../Modules/mp_global.o
 psymd0rho.o : ../PH/phcom.o
 psymd0rho.o : ../PW/para.o
 psymd0rho.o : ../PW/pwcom.o
-psymd0rho.o : d3com.o
 qstar_d3.o : ../Modules/kind.o
+read_ef.o : d3com.o
 read_ef.o : ../Modules/io_global.o
 read_ef.o : ../Modules/mp.o
 read_ef.o : ../PW/pwcom.o
-read_ef.o : d3com.o
 rotate_and_add_d3.o : ../Modules/kind.o
+set_d3irr.o : d3com.o
 set_d3irr.o : ../Modules/control_flags.o
 set_d3irr.o : ../Modules/ions_base.o
 set_d3irr.o : ../Modules/kind.o
 set_d3irr.o : ../PH/phcom.o
 set_d3irr.o : ../PW/pwcom.o
-set_d3irr.o : d3com.o
+set_efsh.o : d3com.o
 set_efsh.o : ../Modules/io_global.o
 set_efsh.o : ../Modules/kind.o
 set_efsh.o : ../PH/phcom.o
 set_efsh.o : ../PW/pwcom.o
-set_efsh.o : d3com.o
 set_sym_irr.o : ../Modules/kind.o
-set_sym_irr.o : ../Modules/mp.o
 set_sym_irr.o : ../Modules/mp_global.o
+set_sym_irr.o : ../Modules/mp.o
 sgama_d3.o : ../Modules/kind.o
+solve_linter_d3.o : d3com.o
 solve_linter_d3.o : ../Modules/io_files.o
 solve_linter_d3.o : ../Modules/io_global.o
 solve_linter_d3.o : ../Modules/ions_base.o
@ -233,32 +234,31 @@ solve_linter_d3.o : ../Modules/kind.o
 solve_linter_d3.o : ../Modules/wavefunctions.o
 solve_linter_d3.o : ../PH/phcom.o
 solve_linter_d3.o : ../PW/pwcom.o
-solve_linter_d3.o : d3com.o
+stop_d3.o : d3com.o
 stop_d3.o : ../Modules/control_flags.o
 stop_d3.o : ../Modules/io_files.o
-stop_d3.o : ../Modules/mp.o
 stop_d3.o : ../Modules/mp_global.o
+stop_d3.o : ../Modules/mp.o
 stop_d3.o : ../PH/phcom.o
 stop_d3.o : ../PW/pwcom.o
-stop_d3.o : d3com.o
+symd0rho.o : ../Modules/kind.o
+sym_def1.o : d3com.o
 sym_def1.o : ../Modules/kind.o
 sym_def1.o : ../PH/phcom.o
 sym_def1.o : ../PW/pwcom.o
-sym_def1.o : d3com.o
-symd0rho.o : ../Modules/kind.o
 trntnsc_3.o : ../Modules/kind.o
 w_1gauss.o : ../Modules/kind.o
+write_aux.o : d3com.o
 write_aux.o : ../Modules/ions_base.o
 write_aux.o : ../PH/phcom.o
 write_aux.o : ../PW/pwcom.o
-write_aux.o : d3com.o
-write_d3dyn.o : ../Modules/kind.o
-write_igk.o : ../Modules/io_files.o
-write_igk.o : ../PH/phcom.o
-write_igk.o : ../PW/pwcom.o
+writed3dyn_5.o : d3com.o
 writed3dyn_5.o : ../Modules/io_global.o
 writed3dyn_5.o : ../Modules/ions_base.o
 writed3dyn_5.o : ../Modules/kind.o
 writed3dyn_5.o : ../PH/phcom.o
 writed3dyn_5.o : ../PW/pwcom.o
-writed3dyn_5.o : d3com.o
+write_d3dyn.o : ../Modules/kind.o
+write_igk.o : ../Modules/io_files.o
+write_igk.o : ../PH/phcom.o
+write_igk.o : ../PW/pwcom.o
--- a/Gamma/Makefile
+++ b/Gamma/Makefile
@ -251,6 +251,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/readpseudo.o \
 ../Modules/recvec.o \
--- a/Gamma/make.depend
+++ b/Gamma/make.depend
@ -1,3 +1,4 @@
+a_h.o : cgcom.o
 a_h.o : ../Modules/cell_base.o
 a_h.o : ../Modules/constants.o
 a_h.o : ../Modules/functionals.o
@ -6,14 +7,19 @@ a_h.o : ../Modules/uspp.o
 a_h.o : ../Modules/wavefunctions.o
 a_h.o : ../PW/becmod.o
 a_h.o : ../PW/pwcom.o
-a_h.o : cgcom.o
+cgcom.o : ../Modules/kind.o
+cg_readin.o : cgcom.o
 cg_readin.o : ../Modules/io_files.o
 cg_readin.o : ../Modules/io_global.o
 cg_readin.o : ../Modules/ions_base.o
 cg_readin.o : ../Modules/kind.o
 cg_readin.o : ../Modules/mp.o
 cg_readin.o : ../PW/pwcom.o
-cg_readin.o : cgcom.o
+cg_setupdgc.o : cgcom.o
+cg_setupdgc.o : ../Modules/functionals.o
+cg_setupdgc.o : ../Modules/kind.o
+cg_setupdgc.o : ../PW/pwcom.o
+cg_setup.o : cgcom.o
 cg_setup.o : ../Modules/atom.o
 cg_setup.o : ../Modules/functionals.o
 cg_setup.o : ../Modules/io_files.o
@ -23,64 +29,59 @@ cg_setup.o : ../Modules/mp_global.o
 cg_setup.o : ../Modules/wavefunctions.o
 cg_setup.o : ../PW/pwcom.o
 cg_setup.o : ../PW/restart.o
-cg_setup.o : cgcom.o
-cg_setupdgc.o : ../Modules/functionals.o
-cg_setupdgc.o : ../Modules/kind.o
-cg_setupdgc.o : ../PW/pwcom.o
-cg_setupdgc.o : cgcom.o
+cgsolve.o : ../Modules/io_global.o
+cgsolve.o : ../Modules/kind.o
+cg_summary.o : cgcom.o
 cg_summary.o : ../Modules/io_global.o
 cg_summary.o : ../Modules/ions_base.o
 cg_summary.o : ../Modules/uspp.o
 cg_summary.o : ../PW/pwcom.o
-cg_summary.o : cgcom.o
-cgcom.o : ../Modules/kind.o
-cgsolve.o : ../Modules/io_global.o
-cgsolve.o : ../Modules/kind.o
 d2ion.o : ../Modules/io_global.o
 d2ion.o : ../Modules/kind.o
 dgcxc.o : ../Modules/kind.o
 dgcxc_spin.o : ../Modules/kind.o
 dgradcorr.o : ../Modules/kind.o
+dielec.o : cgcom.o
 dielec.o : ../Modules/ions_base.o
 dielec.o : ../PW/pwcom.o
-dielec.o : cgcom.o
 dmxc.o : ../Modules/functionals.o
 dmxc.o : ../Modules/kind.o
-drhodv.o : ../PW/pwcom.o
 drhodv.o : cgcom.o
+drhodv.o : ../PW/pwcom.o
+dvpsi_e.o : cgcom.o
 dvpsi_e.o : ../Modules/ions_base.o
 dvpsi_e.o : ../Modules/kind.o
 dvpsi_e.o : ../Modules/uspp.o
 dvpsi_e.o : ../Modules/wavefunctions.o
 dvpsi_e.o : ../PW/becmod.o
 dvpsi_e.o : ../PW/pwcom.o
-dvpsi_e.o : cgcom.o
+dvpsi_kb.o : cgcom.o
 dvpsi_kb.o : ../Modules/atom.o
 dvpsi_kb.o : ../Modules/ions_base.o
 dvpsi_kb.o : ../Modules/kind.o
 dvpsi_kb.o : ../Modules/uspp.o
 dvpsi_kb.o : ../Modules/wavefunctions.o
 dvpsi_kb.o : ../PW/pwcom.o
-dvpsi_kb.o : cgcom.o
 dyndiar.o : ../Modules/io_global.o
 dyndiar.o : ../Modules/kind.o
-dynmat_init.o : ../Modules/ions_base.o
-dynmat_init.o : ../PW/pwcom.o
-dynmat_init.o : cgcom.o
+dynmatcc.o : cgcom.o
 dynmatcc.o : ../Modules/atom.o
 dynmatcc.o : ../Modules/ions_base.o
 dynmatcc.o : ../Modules/wavefunctions.o
 dynmatcc.o : ../PW/pwcom.o
-dynmatcc.o : cgcom.o
+dynmat_init.o : cgcom.o
+dynmat_init.o : ../Modules/ions_base.o
+dynmat_init.o : ../PW/pwcom.o
 generate_dynamical_matrix.o : ../Modules/kind.o
 generate_effective_charges.o : ../Modules/kind.o
+h_h.o : cgcom.o
 h_h.o : ../Modules/kind.o
 h_h.o : ../Modules/uspp.o
 h_h.o : ../PW/becmod.o
 h_h.o : ../PW/pwcom.o
-h_h.o : cgcom.o
-macro.o : ../PW/pwcom.o
 macro.o : cgcom.o
+macro.o : ../PW/pwcom.o
+phcg.o : cgcom.o
 phcg.o : ../Modules/check_stop.o
 phcg.o : ../Modules/control_flags.o
 phcg.o : ../Modules/functionals.o
@ -90,25 +91,24 @@ phcg.o : ../Modules/ions_base.o
 phcg.o : ../Modules/mp.o
 phcg.o : ../Modules/version.o
 phcg.o : ../PW/pwcom.o
-phcg.o : cgcom.o
 pw_dot.o : ../Modules/kind.o
 pw_dot.o : ../PW/pwcom.o
+rhod2vkb.o : cgcom.o
 rhod2vkb.o : ../Modules/ions_base.o
 rhod2vkb.o : ../Modules/uspp.o
 rhod2vkb.o : ../Modules/wavefunctions.o
 rhod2vkb.o : ../PW/pwcom.o
-rhod2vkb.o : cgcom.o
+solve_e.o : cgcom.o
 solve_e.o : ../Modules/io_global.o
 solve_e.o : ../Modules/wavefunctions.o
 solve_e.o : ../PW/becmod.o
 solve_e.o : ../PW/pwcom.o
-solve_e.o : cgcom.o
+solve_ph.o : cgcom.o
 solve_ph.o : ../Modules/io_files.o
 solve_ph.o : ../Modules/io_global.o
 solve_ph.o : ../Modules/wavefunctions.o
 solve_ph.o : ../PW/becmod.o
 solve_ph.o : ../PW/pwcom.o
-solve_ph.o : cgcom.o
+writedyn.o : cgcom.o
 writedyn.o : ../Modules/ions_base.o
 writedyn.o : ../PW/pwcom.o
-writedyn.o : cgcom.o
--- a/Modules/Makefile
+++ b/Modules/Makefile
@ -13,6 +13,7 @@ clocks.o \
 constants.o \
 control_flags.o \
 cp_emass.o \
+cp_restart.o \
 cp_fpmd.o \
 descriptors.o \
 electrons_base.o \
--- a/Modules/cp_restart.f90
+++ b/Modules/cp_restart.f90
@ -0,0 +1,230 @@
+!
+! Copyright (C) 2002 FPMD group
+! This file is distributed under the terms of the
+! GNU General Public License. See the file `License'
+! in the root directory of the present distribution,
+! or http://www.gnu.org/copyleft/gpl.txt .
+!
+
+!------------------------------------------------------------------------------!
+  MODULE cp_restart
+!------------------------------------------------------------------------------!
+
+    !  This module contains subroutines to write and read data required to
+    !  restart a calculation from the disk  
+!
+      IMPLICIT NONE
+      SAVE
+
+!------------------------------------------------------------------------------!
+  CONTAINS
+!------------------------------------------------------------------------------!
+
+    SUBROUTINE create_directories( dirname )
+      USE io_global, ONLY: ionode, ionode_id, stdout
+      USE mp, ONLY: mp_bcast
+      !
+      IMPLICIT NONE
+      !
+      CHARACTER(LEN=256) :: dirname
+      INTEGER :: ios
+      !
+      INTEGER, EXTERNAL :: C_MKDIR
+
+      IF ( ionode ) THEN
+          WRITE( stdout, * ) 'Creating dir : ',TRIM( dirname )
+          !
+          ios = C_MKDIR( TRIM( dirname ), LEN_TRIM( dirname ) )
+      END IF
+      CALL mp_bcast( ios, ionode_id )
+      IF( ios /= 0 ) THEN
+        CALL errore(' cp_writefile ', ' unable to create directory '//dirname , ios )
+      END IF  
+
+      IF ( ionode ) THEN
+        OPEN( UNIT = 4, FILE = TRIM( dirname ) // '/' // 'tst' , &
+            STATUS = 'UNKNOWN', FORM = 'UNFORMATTED', IOSTAT = ios )
+        CLOSE( UNIT = 4, STATUS = 'DELETE' )
+      END IF
+      !
+      CALL mp_bcast( ios, ionode_id )
+      !
+      IF ( ios /= 0 ) &
+           CALL errore( ' cp_writefile : ', TRIM( dirname ) // &
+                      & ' non existent or non writable', 1 )
+
+
+      RETURN
+    END SUBROUTINE
+
+!------------------------------------------------------------------------------!
+
+    SUBROUTINE cp_writefile( scradir, ascii, nfi, simtime, acc, nk, xk, wk )
+      USE iotk_module
+      USE io_global, ONLY: ionode, ionode_id, stdout
+      USE parser, ONLY: int_to_char
+      USE control_flags, ONLY: ndw, gamma_only
+      USE io_files, ONLY: iunpun, xmlpun
+      USE printout_base, ONLY: title
+      USE kinds, ONLY: dbl
+      USE grid_dimensions, ONLY: nr1, nr2, nr3
+      USE smooth_grid_dimensions, ONLY: nr1s, nr2s, nr3s
+      USE smallbox_grid_dimensions, ONLY: nr1b, nr2b, nr3b
+      USE gvecp, ONLY: ngm, ecutp
+      USE gvecs, ONLY: ngs
+      USE gvecw, ONLY: ngs, ecutw
+      USE electrons_base, ONLY: nspin, nbnd, nelt, nel, nupdwn, iupdwn
+      USE cell_base, ONLY: ibrav, alat, celldm
+      USE ions_base, ONLY: nsp, nat, na, atm, zv, pmass
+      
+      IMPLICIT NONE
+      CHARACTER(LEN=*), INTENT(IN) :: scradir
+      LOGICAL, INTENT(IN) :: ascii
+      INTEGER, INTENT(IN) :: nfi    !  index of the current step
+      REAL(dbl), INTENT(IN) :: simtime   !  simulated time
+      REAL(dbl), INTENT(IN) :: acc(:)   !  
+      INTEGER, INTENT(IN) :: nk    !  number of kpoints
+      REAL(dbl), INTENT(IN) :: xk(:,:) ! k points coordinates 
+      REAL(dbl), INTENT(IN) :: wk(:)   ! k points weights
+      !
+      CHARACTER(LEN=256) :: dirname
+      INTEGER :: strlen
+      CHARACTER(iotk_attlenx) :: attr
+      INTEGER :: kunit
+      INTEGER :: k1, k2, k3
+      INTEGER :: nk1, nk2, nk3
+      INTEGER :: i
+
+      dirname = 'RESTART' // int_to_char( ndw )
+      IF ( LEN( scradir ) > 1 ) THEN
+         strlen  = index(scradir,' ') - 1
+         dirname = scradir(1:strlen) // '/' // dirname
+      END IF 
+
+      CALL create_directories( dirname )
+
+      !  Some default values
+      !
+      kunit = 1
+      k1 = 0
+      k2 = 0
+      k3 = 0
+      nk1 = 1
+      nk2 = 1
+      nk3 = 1
+
+
+      !  Open XML descriptor
+
+      IF ( ionode ) THEN
+        write(stdout,*) "Opening file "//trim(xmlpun)
+        IF( ascii ) THEN
+          call iotk_open_write(iunpun,file=TRIM(dirname)//'/'//TRIM(xmlpun),binary=.false.)
+        ELSE
+          call iotk_open_write(iunpun,file=TRIM(dirname)//'/'//TRIM(xmlpun),binary=.true.)
+        ENDIF
+      END IF
+
+      IF( ionode ) THEN
+        !
+        call iotk_write_begin(iunpun,"STATUS")
+        call iotk_write_attr (attr,"nfi",nfi,first=.true.)
+        call iotk_write_empty(iunpun,"STEP",attr)
+        call iotk_write_dat (iunpun, "TIME", simtime)
+        call iotk_write_dat (iunpun, "TITLE", TRIM(title) )
+        call iotk_write_end(iunpun,"STATUS")
+        !
+        call iotk_write_begin(iunpun,"DIMENSIONS")
+        call iotk_write_attr (attr,"nx",nr1,first=.true.)
+        call iotk_write_attr (attr,"ny",nr2)
+        call iotk_write_attr (attr,"nz",nr3)
+        call iotk_write_empty(iunpun,"FFT_GRID",attr)
+        call iotk_write_attr (attr,"nx",nr1s,first=.true.)
+        call iotk_write_attr (attr,"ny",nr2s)
+        call iotk_write_attr (attr,"nz",nr3s)
+        call iotk_write_empty(iunpun,"SMOOTH_FFT_GRID",attr)
+        call iotk_write_attr (attr,"nx",nr1b,first=.true.)
+        call iotk_write_attr (attr,"ny",nr2b)
+        call iotk_write_attr (attr,"nz",nr3b)
+        call iotk_write_empty(iunpun,"SMALLBOX_FFT_GRID",attr)
+        call iotk_write_attr (attr,"ng",ngm,first=.true.)
+        call iotk_write_empty(iunpun,"GVECTORS",attr)
+        call iotk_write_attr (attr,"ng",ngs,first=.true.)
+        call iotk_write_empty(iunpun,"SMOOTH_GVECTORS",attr)
+        call iotk_write_end(iunpun,"DIMENSIONS")
+        !
+        call iotk_write_begin(iunpun,"ELECTONS")
+        call iotk_write_attr (attr,"nspin",nspin,first=.true.)
+        call iotk_write_attr (attr,"nbnd",nbnd)
+        call iotk_write_attr (attr,"nel",nelt)
+        call iotk_write_empty(iunpun,"BANDS",attr)
+        call iotk_write_attr (attr,"nel",nel(1),first=.true.)
+        call iotk_write_attr (attr,"nbnd",nupdwn(1))
+        call iotk_write_empty(iunpun,"SPINUP",attr)
+        call iotk_write_attr (attr,"nel",nel(2),first=.true.)
+        call iotk_write_attr (attr,"nbnd",nupdwn(2))
+        call iotk_write_empty(iunpun,"SPINDW",attr)
+        call iotk_write_end(iunpun,"ELECTONS")
+        ! 
+        call iotk_write_begin(iunpun,"AVERAGES")
+        call iotk_write_dat (iunpun, "ACCUMULATORS", acc)
+        call iotk_write_end(iunpun,"AVERAGES")
+        !
+        call iotk_write_attr (attr,"unit","Hartree",first=.true.)
+        call iotk_write_begin(iunpun,"CUTOFFS",attr)
+        call iotk_write_dat (iunpun, "ECUTRHO", ecutp/2.0d0)
+        call iotk_write_dat (iunpun, "ECUTWFC", ecutw/2.0d0)
+        call iotk_write_end(iunpun,"CUTOFFS")
+        !
+        call iotk_write_begin(iunpun,"K_POINTS")
+        call iotk_write_attr (attr, "nk",nk,first=.true.)
+        call iotk_write_attr (attr, "gamma_only", gamma_only )
+        call iotk_write_attr (attr, "kunit", kunit )
+        call iotk_write_attr (attr, "k1", k1 )
+        call iotk_write_attr (attr, "k2", k2 )
+        call iotk_write_attr (attr, "k3", k3 )
+        call iotk_write_attr (attr, "nk1", nk1 )
+        call iotk_write_attr (attr, "nk2", nk2 )
+        call iotk_write_attr (attr, "nk3", nk3 )
+        call iotk_write_empty(iunpun,"GRID",attr)
+        DO i = 1, nk
+          call iotk_write_attr (attr, "ik",i,first=.true.)
+          call iotk_write_attr (attr,"xyz",xk(:,i))
+          call iotk_write_attr (attr,"w",wk(i))
+          call iotk_write_empty(iunpun,"K",attr)
+        END DO
+        call iotk_write_end(iunpun,"K_POINTS")
+        !
+        call iotk_write_begin(iunpun,"LATTICE")
+        call iotk_write_attr (attr, "ibrav",ibrav,first=.true.)
+        call iotk_write_attr (attr, "alat", alat )
+        call iotk_write_dat (iunpun, "CELLDM", celldm, attr)
+        call iotk_write_end(iunpun,"LATTICE")
+        ! 
+        call iotk_write_attr (attr, "nsp", nsp, first=.true.)
+        call iotk_write_attr (attr, "nat", nat )
+        call iotk_write_begin(iunpun,"IONS", attr)
+        DO i = 1, nsp
+          call iotk_write_attr (attr, "label",atm(i),first=.true.)
+          call iotk_write_attr (attr, "na",na(i) )
+          call iotk_write_attr (attr, "zv", zv(i) )
+          call iotk_write_attr (attr, "mass", pmass(i) )
+          call iotk_write_empty(iunpun,"SPECIE",attr)
+        END DO
+        call iotk_write_end(iunpun,"IONS")
+        
+
+      
+      END IF
+
+      if( ionode ) then
+        call iotk_close_write( iunpun )
+      end if 
+
+      RETURN
+    END SUBROUTINE
+
+
+!------------------------------------------------------------------------------!
+  END MODULE cp_restart
+!------------------------------------------------------------------------------!
--- a/Modules/io_files.f90
+++ b/Modules/io_files.f90
@ -54,6 +54,7 @@ MODULE io_files
    axsf_file = 'os.axsf'    ! file containing coordinates of all images in axsf format
  CHARACTER (LEN=261) :: &
    exit_file = "os.EXIT"    ! file required for a soft exit  
+  CHARACTER (LEN=11), PARAMETER :: xmlpun = 'restart.xml'
  !
  ! ... The units where various variables are saved
  !
--- a/Modules/ions_base.f90
+++ b/Modules/ions_base.f90
@ -48,10 +48,10 @@
      CHARACTER(LEN=3  ) :: atm( ntypx ) 
      CHARACTER(LEN=80 ) :: tau_units

-      ! if if_pos( x, i ) = 0 then 
-      !    x coordinate of the i-th atom will be kept fixed

-      INTEGER, ALLOCATABLE :: if_pos(:,:)  
+      INTEGER, ALLOCATABLE :: if_pos(:,:)     ! if if_pos( x, i ) = 0 then  x coordinate of 
+                                              ! the i-th atom will be kept fixed
+      INTEGER, ALLOCATABLE :: iforce(:,:)     ! if_pos sorted by specie 
      INTEGER :: fixatom                      !!! to be removed

      INTEGER :: ind_localisation(natx) = 0   ! true if we want to know the localization arount the atom
@ -217,6 +217,7 @@
      ALLOCATE( vel_srt( 3, nat ) )
      ALLOCATE( ind_srt( nat ) )
      ALLOCATE( if_pos( 3, nat ) )
+      ALLOCATE( iforce( 3, nat ) )

      ityp( 1:nat )      = ityp_ ( 1:nat )
      vel( : , 1:nat )   = vel_ ( : , 1:nat )
@ -309,6 +310,13 @@
      if_pos = 1
      if_pos(:,:) = if_pos_ (:,1:nat)

+      iforce = 0
+      DO ia = 1, nat
+        iforce ( :, ia ) = if_pos ( :, ind_srt( ia ) )
+      END DO
+
+
+
      IF( PRESENT( sic_ ) ) THEN
        select case ( TRIM( sic_ ) )
        case ( 'sic_pz' ) 
@ -373,6 +381,7 @@
      IF( ALLOCATED( vel_srt ) ) DEALLOCATE( vel_srt )
      IF( ALLOCATED( ind_srt ) ) DEALLOCATE( ind_srt )
      IF( ALLOCATED( if_pos ) ) DEALLOCATE( if_pos )
+      IF( ALLOCATED( iforce ) ) DEALLOCATE( iforce )
      IF( ALLOCATED( pos_localisation ) ) DEALLOCATE( pos_localisation )
      tions_base_init = .FALSE.
      RETURN
@ -705,8 +714,6 @@
      REAL(dbl) :: taus(3,natx), tausm(3,natx), tausp(3,natx)
      REAL(dbl) :: vels(3,natx), velsm(3,natx), velsp(3,natx)

-      INTEGER :: iforce(3,natx)
-
 !------------------------------------------------------------------------------!
  CONTAINS 
 !------------------------------------------------------------------------------!
--- a/Modules/make.depend
+++ b/Modules/make.depend
@ -2,8 +2,8 @@ atom.o : kind.o
 atom.o : parameters.o
 basic_algebra_routines.o : kind.o
 berry_phase.o : io_global.o
-berry_phase.o : mp.o
 berry_phase.o : mp_global.o
+berry_phase.o : mp.o
 berry_phase.o : recvec.o
 berry_phase.o : stick_base.o
 bfgs_module.o : basic_algebra_routines.o
@ -16,12 +16,12 @@ cell_base.o : kind.o
 check_stop.o : io_files.o
 check_stop.o : io_global.o
 check_stop.o : kind.o
-check_stop.o : mp.o
 check_stop.o : mp_global.o
+check_stop.o : mp.o
 clocks.o : io_global.o
 clocks.o : kind.o
-clocks.o : mp.o
 clocks.o : mp_global.o
+clocks.o : mp.o
 constants.o : kind.o
 control_flags.o : kind.o
 control_flags.o : parameters.o
@ -40,6 +40,14 @@ cp_fpmd.o : parameters.o
 cp_fpmd.o : recvec.o
 cp_fpmd.o : smallbox.o
 cp_fpmd.o : uspp.o
+cp_restart.o : control_flags.o
+cp_restart.o : io_files.o
+cp_restart.o : io_global.o
+cp_restart.o : iotk_module.o
+cp_restart.o : kind.o
+cp_restart.o : mp.o
+cp_restart.o : parser.o
+cp_restart.o : printout_base.o
 descriptors.o : io_global.o
 descriptors.o : parallel_types.o
 electrons_base.o : kind.o
@ -47,9 +55,9 @@ energies.o : io_global.o
 energies.o : kind.o
 fft_base.o : fft_types.o
 fft_base.o : kind.o
-fft_base.o : mp.o
 fft_base.o : mp_buffers.o
 fft_base.o : mp_global.o
+fft_base.o : mp.o
 fft_base.o : parallel_include.o
 fft_scalar.o : kind.o
 fft_types.o : io_global.o
@ -57,14 +65,14 @@ griddim.o : fft_scalar.o
 griddim.o : fft_types.o
 griddim.o : io_global.o
 griddim.o : kind.o
-griddim.o : mp.o
 griddim.o : mp_global.o
+griddim.o : mp.o
 input_parameters.o : kind.o
 input_parameters.o : parameters.o
 io_base.o : io_global.o
 io_base.o : kind.o
-io_base.o : mp.o
 io_base.o : mp_global.o
+io_base.o : mp.o
 io_base.o : mp_wave.o
 io_base.o : parameters.o
 io_files.o : parameters.o
@ -75,23 +83,23 @@ ions_base.o : io_global.o
 ions_base.o : kind.o
 ions_base.o : parameters.o
 kind.o : io_global.o
-mp.o : io_global.o
-mp.o : kind.o
-mp.o : parallel_include.o
 mp_buffers.o : io_global.o
 mp_buffers.o : kind.o
 mp_buffers.o : parallel_include.o
 mp_buffers.o : shmem_include.o
 mp_global.o : parallel_include.o
 mp_global.o : shmem_include.o
+mp.o : io_global.o
+mp.o : kind.o
+mp.o : parallel_include.o
 mp_wave.o : kind.o
 mp_wave.o : parallel_include.o
 parallel_include.o : kind.o
 parallel_types.o : kind.o
 parser.o : io_global.o
 parser.o : kind.o
-parser.o : mp.o
 parser.o : mp_global.o
+parser.o : mp.o
 path_base.o : basic_algebra_routines.o
 path_base.o : cell_base.o
 path_base.o : check_stop.o
@ -118,8 +126,8 @@ path_io_routines.o : io_files.o
 path_io_routines.o : io_global.o
 path_io_routines.o : ions_base.o
 path_io_routines.o : kind.o
-path_io_routines.o : mp.o
 path_io_routines.o : mp_global.o
+path_io_routines.o : mp.o
 path_io_routines.o : parser.o
 path_io_routines.o : path_formats.o
 path_io_routines.o : path_variables.o
@ -131,12 +139,12 @@ path_opt_routines.o : path_variables.o
 path_variables.o : kind.o
 printout_base.o : io_global.o
 printout_base.o : kind.o
-printout_base.o : mp.o
 printout_base.o : mp_global.o
-pseudo_types.o : kind.o
-pseudo_types.o : parameters.o
+printout_base.o : mp.o
 pseudodata.o : kind.o
 pseudodata.o : parameters.o
+pseudo_types.o : kind.o
+pseudo_types.o : parameters.o
 ptoolkit.o : io_global.o
 ptoolkit.o : kind.o
 ptoolkit.o : parallel_include.o
@ -144,20 +152,20 @@ read_cards.o : constants.o
 read_cards.o : input_parameters.o
 read_cards.o : io_global.o
 read_cards.o : kind.o
-read_cards.o : mp.o
 read_cards.o : mp_global.o
+read_cards.o : mp.o
 read_cards.o : parser.o
 read_namelists.o : constants.o
 read_namelists.o : input_parameters.o
 read_namelists.o : io_global.o
 read_namelists.o : kind.o
-read_namelists.o : mp.o
 read_namelists.o : mp_global.o
+read_namelists.o : mp.o
 readpseudo.o : kind.o
 readpseudo.o : pseudo_types.o
 recvec.o : kind.o
-recvec.o : mp.o
 recvec.o : mp_global.o
+recvec.o : mp.o
 shmem_include.o : kind.o
 smallbox.o : constants.o
 smallbox.o : io_global.o
@ -167,8 +175,8 @@ stick_base.o : control_flags.o
 stick_base.o : fft_types.o
 stick_base.o : io_global.o
 stick_base.o : kind.o
-stick_base.o : mp.o
 stick_base.o : mp_global.o
+stick_base.o : mp.o
 supercell.o : cell_base.o
 supercell.o : kind.o
 timestep.o : kind.o
@ -176,6 +184,6 @@ uspp.o : constants.o
 uspp.o : kind.o
 uspp.o : parameters.o
 wave_base.o : kind.o
-wave_base.o : mp.o
 wave_base.o : mp_global.o
+wave_base.o : mp.o
 wavefunctions.o : kind.o
--- a/Modules/printout_base.f90
+++ b/Modules/printout_base.f90
@ -14,13 +14,15 @@ MODULE printout_base
  IMPLICIT NONE
  SAVE

+  CHARACTER(LEN=75) :: title
+  ! ...  title of the simulation
+
  CHARACTER(LEN=256) :: fort_unit(30:40)
  ! ...  fort_unit = fortran units for saving physical quantity

  CHARACTER(LEN=256) :: pprefix
  ! ...  prefix combined with the outpur path

-
 CONTAINS


--- a/Nmr/make.depend
+++ b/Nmr/make.depend
@ -7,6 +7,9 @@ magnetic_kkterm.o : ../Modules/kind.o
 magnetic_kkterm.o : ../Modules/wavefunctions.o
 magnetic_kkterm.o : ../PH/phcom.o
 magnetic_kkterm.o : ../PW/pwcom.o
+nmr_mod.o : ../Modules/kind.o
+nmr_mod.o : ../Modules/uspp.o
+nmr_mod.o : ../PW/pwcom.o
 nmr.o : ../Modules/cell_base.o
 nmr.o : ../Modules/io_files.o
 nmr.o : ../Modules/ions_base.o
@ -15,14 +18,11 @@ nmr.o : ../Modules/parameters.o
 nmr.o : ../Modules/uspp.o
 nmr.o : ../Modules/version.o
 nmr.o : ../Modules/wavefunctions.o
+nmr.o : nmr_mod.o
 nmr.o : ../PH/phcom.o
 nmr.o : ../PW/becmod.o
 nmr.o : ../PW/paw.o
 nmr.o : ../PW/pwcom.o
-nmr.o : nmr_mod.o
-nmr_mod.o : ../Modules/kind.o
-nmr_mod.o : ../Modules/uspp.o
-nmr_mod.o : ../PW/pwcom.o
 solve_cg.o : ../Modules/kind.o
 solve_cg.o : ../Modules/uspp.o
 solve_cg.o : ../Modules/wavefunctions.o
@ -33,7 +33,7 @@ take_nloc_k_kq.o : ../Modules/ions_base.o
 take_nloc_k_kq.o : ../Modules/kind.o
 take_nloc_k_kq.o : ../Modules/uspp.o
 take_nloc_k_kq.o : ../Modules/wavefunctions.o
+take_nloc_k_kq.o : nmr_mod.o
 take_nloc_k_kq.o : ../PH/phcom.o
 take_nloc_k_kq.o : ../PW/becmod.o
 take_nloc_k_kq.o : ../PW/pwcom.o
-take_nloc_k_kq.o : nmr_mod.o
--- a/PH/Makefile
+++ b/PH/Makefile
@ -141,6 +141,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/read_cards.o \
 ../Modules/read_namelists.o \
--- a/PH/make.depend
+++ b/PH/make.depend
@ -1,136 +1,136 @@
-add_dkmds.o : ../Modules/ions_base.o
-add_dkmds.o : ../Modules/kind.o
-add_dkmds.o : ../Modules/uspp.o
-add_dkmds.o : ../Modules/wavefunctions.o
-add_dkmds.o : ../PW/pwcom.o
-add_dkmds.o : phcom.o
-add_for_charges.o : ../Modules/ions_base.o
-add_for_charges.o : ../Modules/kind.o
-add_for_charges.o : ../Modules/uspp.o
-add_for_charges.o : ../PW/pwcom.o
-add_for_charges.o : phcom.o
-add_zstar_ue.o : ../Modules/io_files.o
-add_zstar_ue.o : ../Modules/kind.o
-add_zstar_ue.o : ../Modules/wavefunctions.o
-add_zstar_ue.o : ../PW/pwcom.o
-add_zstar_ue.o : phcom.o
-add_zstar_ue_us.o : ../Modules/io_files.o
-add_zstar_ue_us.o : ../Modules/kind.o
-add_zstar_ue_us.o : ../Modules/wavefunctions.o
-add_zstar_ue_us.o : ../PW/pwcom.o
-add_zstar_ue_us.o : phcom.o
 addcore.o : ../Modules/atom.o
 addcore.o : ../Modules/cell_base.o
 addcore.o : ../Modules/ions_base.o
 addcore.o : ../Modules/kind.o
-addcore.o : ../PW/pwcom.o
 addcore.o : phcom.o
+addcore.o : ../PW/pwcom.o
+add_dkmds.o : ../Modules/ions_base.o
+add_dkmds.o : ../Modules/kind.o
+add_dkmds.o : ../Modules/uspp.o
+add_dkmds.o : ../Modules/wavefunctions.o
+add_dkmds.o : phcom.o
+add_dkmds.o : ../PW/pwcom.o
 adddvepsi_us.o : ../Modules/ions_base.o
 adddvepsi_us.o : ../Modules/kind.o
 adddvepsi_us.o : ../Modules/uspp.o
-adddvepsi_us.o : ../PW/pwcom.o
 adddvepsi_us.o : phcom.o
+adddvepsi_us.o : ../PW/pwcom.o
 adddvscf.o : ../Modules/ions_base.o
 adddvscf.o : ../Modules/kind.o
 adddvscf.o : ../Modules/uspp.o
-adddvscf.o : ../PW/pwcom.o
 adddvscf.o : phcom.o
+adddvscf.o : ../PW/pwcom.o
+add_for_charges.o : ../Modules/ions_base.o
+add_for_charges.o : ../Modules/kind.o
+add_for_charges.o : ../Modules/uspp.o
+add_for_charges.o : phcom.o
+add_for_charges.o : ../PW/pwcom.o
 addnlcc.o : ../Modules/functionals.o
 addnlcc.o : ../Modules/ions_base.o
 addnlcc.o : ../Modules/kind.o
-addnlcc.o : ../PW/pwcom.o
 addnlcc.o : phcom.o
+addnlcc.o : ../PW/pwcom.o
 addnlcc_zstar_eu_us.o : ../Modules/functionals.o
 addnlcc_zstar_eu_us.o : ../Modules/kind.o
 addnlcc_zstar_eu_us.o : ../Modules/mp_global.o
-addnlcc_zstar_eu_us.o : ../PW/pwcom.o
 addnlcc_zstar_eu_us.o : phcom.o
+addnlcc_zstar_eu_us.o : ../PW/pwcom.o
 addusdbec.o : ../Modules/ions_base.o
 addusdbec.o : ../Modules/kind.o
 addusdbec.o : ../Modules/uspp.o
-addusdbec.o : ../PW/pwcom.o
 addusdbec.o : phcom.o
+addusdbec.o : ../PW/pwcom.o
+addusddense.o : ../Modules/ions_base.o
+addusddense.o : ../Modules/kind.o
+addusddense.o : ../Modules/uspp.o
+addusddense.o : phcom.o
+addusddense.o : ../PW/pwcom.o
 addusddens.o : ../Modules/ions_base.o
 addusddens.o : ../Modules/kind.o
 addusddens.o : ../Modules/uspp.o
 addusddens.o : ../Modules/wavefunctions.o
-addusddens.o : ../PW/pwcom.o
 addusddens.o : phcom.o
-addusddense.o : ../Modules/ions_base.o
-addusddense.o : ../Modules/kind.o
-addusddense.o : ../Modules/uspp.o
-addusddense.o : ../PW/pwcom.o
-addusddense.o : phcom.o
+addusddens.o : ../PW/pwcom.o
 addusdynmat.o : ../Modules/ions_base.o
 addusdynmat.o : ../Modules/kind.o
 addusdynmat.o : ../Modules/uspp.o
-addusdynmat.o : ../PW/pwcom.o
 addusdynmat.o : phcom.o
+addusdynmat.o : ../PW/pwcom.o
 addusldos.o : ../Modules/ions_base.o
 addusldos.o : ../Modules/uspp.o
 addusldos.o : ../Modules/wavefunctions.o
 addusldos.o : ../PW/pwcom.o
+add_zstar_ue.o : ../Modules/io_files.o
+add_zstar_ue.o : ../Modules/kind.o
+add_zstar_ue.o : ../Modules/wavefunctions.o
+add_zstar_ue.o : phcom.o
+add_zstar_ue.o : ../PW/pwcom.o
+add_zstar_ue_us.o : ../Modules/io_files.o
+add_zstar_ue_us.o : ../Modules/kind.o
+add_zstar_ue_us.o : ../Modules/wavefunctions.o
+add_zstar_ue_us.o : phcom.o
+add_zstar_ue_us.o : ../PW/pwcom.o
 allocate_part.o : ../Modules/ions_base.o
 allocate_part.o : ../Modules/kind.o
-allocate_part.o : ../PW/pwcom.o
 allocate_part.o : phcom.o
+allocate_part.o : ../PW/pwcom.o
+allocate_phq.o : elph.o
 allocate_phq.o : ../Modules/ions_base.o
 allocate_phq.o : ../Modules/kind.o
 allocate_phq.o : ../Modules/uspp.o
 allocate_phq.o : ../Modules/wavefunctions.o
+allocate_phq.o : phcom.o
 allocate_phq.o : ../PW/becmod.o
 allocate_phq.o : ../PW/pwcom.o
-allocate_phq.o : elph.o
-allocate_phq.o : phcom.o
+bcast_ph_input1.o : ../Modules/mp.o
+bcast_ph_input1.o : phcom.o
+bcast_ph_input1.o : ../PW/pwcom.o
 bcast_ph_input.o : ../Modules/control_flags.o
 bcast_ph_input.o : ../Modules/io_files.o
 bcast_ph_input.o : ../Modules/mp.o
-bcast_ph_input.o : ../PW/pwcom.o
 bcast_ph_input.o : phcom.o
-bcast_ph_input1.o : ../Modules/mp.o
-bcast_ph_input1.o : ../PW/pwcom.o
-bcast_ph_input1.o : phcom.o
+bcast_ph_input.o : ../PW/pwcom.o
 cg_psi.o : ../Modules/kind.o
 cgsolve_all.o : ../Modules/kind.o
 ch_psi_all.o : ../Modules/kind.o
+ch_psi_all.o : phcom.o
 ch_psi_all.o : ../PW/becmod.o
 ch_psi_all.o : ../PW/pwcom.o
-ch_psi_all.o : phcom.o
 close_phq.o : ../Modules/control_flags.o
 close_phq.o : ../Modules/io_files.o
 close_phq.o : ../Modules/mp_global.o
-close_phq.o : ../PW/pwcom.o
 close_phq.o : phcom.o
+close_phq.o : ../PW/pwcom.o
 compute_alphasum.o : ../Modules/ions_base.o
 compute_alphasum.o : ../Modules/kind.o
 compute_alphasum.o : ../Modules/uspp.o
-compute_alphasum.o : ../PW/pwcom.o
 compute_alphasum.o : phcom.o
+compute_alphasum.o : ../PW/pwcom.o
 compute_becalp.o : ../Modules/io_files.o
 compute_becalp.o : ../Modules/kind.o
-compute_becalp.o : ../PW/pwcom.o
 compute_becalp.o : phcom.o
+compute_becalp.o : ../PW/pwcom.o
 compute_becsum.o : ../Modules/ions_base.o
 compute_becsum.o : ../Modules/kind.o
 compute_becsum.o : ../Modules/uspp.o
-compute_becsum.o : ../PW/pwcom.o
 compute_becsum.o : phcom.o
+compute_becsum.o : ../PW/pwcom.o
 compute_drhous.o : ../Modules/io_files.o
 compute_drhous.o : ../Modules/ions_base.o
 compute_drhous.o : ../Modules/kind.o
 compute_drhous.o : ../Modules/uspp.o
 compute_drhous.o : ../Modules/wavefunctions.o
-compute_drhous.o : ../PW/pwcom.o
 compute_drhous.o : phcom.o
+compute_drhous.o : ../PW/pwcom.o
 compute_dvloc.o : ../Modules/ions_base.o
 compute_dvloc.o : ../Modules/kind.o
-compute_dvloc.o : ../PW/pwcom.o
 compute_dvloc.o : phcom.o
+compute_dvloc.o : ../PW/pwcom.o
 compute_nldyn.o : ../Modules/ions_base.o
 compute_nldyn.o : ../Modules/kind.o
 compute_nldyn.o : ../Modules/uspp.o
-compute_nldyn.o : ../PW/pwcom.o
 compute_nldyn.o : phcom.o
+compute_nldyn.o : ../PW/pwcom.o
 compute_qdipol.o : ../Modules/atom.o
 compute_qdipol.o : ../Modules/constants.o
 compute_qdipol.o : ../Modules/ions_base.o
@ -138,23 +138,23 @@ compute_qdipol.o : ../Modules/kind.o
 compute_qdipol.o : ../Modules/uspp.o
 compute_qdipol.o : phcom.o
 compute_weight.o : ../Modules/kind.o
-compute_weight.o : ../PW/pwcom.o
 compute_weight.o : phcom.o
+compute_weight.o : ../PW/pwcom.o
 d2ionq.o : ../Modules/constants.o
 d2ionq.o : ../Modules/io_global.o
 d2ionq.o : ../Modules/kind.o
 davcio_drho.o : ../Modules/kind.o
-davcio_drho.o : ../Modules/mp.o
 davcio_drho.o : ../Modules/mp_global.o
+davcio_drho.o : ../Modules/mp.o
 davcio_drho.o : ../Modules/parallel_include.o
+davcio_drho.o : phcom.o
 davcio_drho.o : ../PW/para.o
 davcio_drho.o : ../PW/pwcom.o
-davcio_drho.o : phcom.o
 deallocate_part.o : phcom.o
-deallocate_phq.o : ../Modules/wavefunctions.o
-deallocate_phq.o : ../PW/becmod.o
 deallocate_phq.o : elph.o
+deallocate_phq.o : ../Modules/wavefunctions.o
 deallocate_phq.o : phcom.o
+deallocate_phq.o : ../PW/becmod.o
 dgcxc.o : ../Modules/kind.o
 dgcxc_spin.o : ../Modules/kind.o
 dgradcorr.o : ../Modules/constants.o
@ -162,125 +162,125 @@ dgradcorr.o : ../Modules/kind.o
 dielec.o : ../Modules/io_files.o
 dielec.o : ../Modules/io_global.o
 dielec.o : ../Modules/kind.o
-dielec.o : ../PW/pwcom.o
 dielec.o : phcom.o
+dielec.o : ../PW/pwcom.o
 dmxc.o : ../Modules/functionals.o
 dmxc.o : ../Modules/kind.o
 dmxc_spin.o : ../Modules/functionals.o
 dmxc_spin.o : ../Modules/kind.o
-drho.o : ../Modules/ions_base.o
-drho.o : ../Modules/kind.o
-drho.o : ../Modules/uspp.o
-drho.o : ../PW/pwcom.o
-drho.o : phcom.o
-drhodv.o : ../Modules/io_files.o
-drhodv.o : ../Modules/ions_base.o
-drhodv.o : ../Modules/kind.o
-drhodv.o : ../PW/pwcom.o
-drhodv.o : phcom.o
 drhodvloc.o : ../Modules/ions_base.o
 drhodvloc.o : ../Modules/kind.o
-drhodvloc.o : ../PW/pwcom.o
 drhodvloc.o : phcom.o
+drhodvloc.o : ../PW/pwcom.o
 drhodvnl.o : ../Modules/ions_base.o
 drhodvnl.o : ../Modules/kind.o
 drhodvnl.o : ../Modules/uspp.o
-drhodvnl.o : ../PW/pwcom.o
 drhodvnl.o : phcom.o
+drhodvnl.o : ../PW/pwcom.o
+drhodv.o : ../Modules/io_files.o
+drhodv.o : ../Modules/ions_base.o
+drhodv.o : ../Modules/kind.o
+drhodv.o : phcom.o
+drhodv.o : ../PW/pwcom.o
 drhodvus.o : ../Modules/ions_base.o
 drhodvus.o : ../Modules/kind.o
-drhodvus.o : ../PW/pwcom.o
 drhodvus.o : phcom.o
-dv_of_drho.o : ../Modules/functionals.o
-dv_of_drho.o : ../Modules/kind.o
-dv_of_drho.o : ../PW/pwcom.o
-dv_of_drho.o : phcom.o
+drhodvus.o : ../PW/pwcom.o
+drho.o : ../Modules/ions_base.o
+drho.o : ../Modules/kind.o
+drho.o : ../Modules/uspp.o
+drho.o : phcom.o
+drho.o : ../PW/pwcom.o
 dvanqq.o : ../Modules/ions_base.o
 dvanqq.o : ../Modules/kind.o
 dvanqq.o : ../Modules/uspp.o
-dvanqq.o : ../PW/pwcom.o
 dvanqq.o : phcom.o
+dvanqq.o : ../PW/pwcom.o
 dvkb3.o : ../Modules/ions_base.o
 dvkb3.o : ../Modules/kind.o
 dvkb3.o : ../Modules/uspp.o
 dvkb3.o : ../Modules/wavefunctions.o
-dvkb3.o : ../PW/pwcom.o
 dvkb3.o : phcom.o
+dvkb3.o : ../PW/pwcom.o
+dv_of_drho.o : ../Modules/functionals.o
+dv_of_drho.o : ../Modules/kind.o
+dv_of_drho.o : phcom.o
+dv_of_drho.o : ../PW/pwcom.o
 dvpsi_e.o : ../Modules/io_global.o
 dvpsi_e.o : ../Modules/ions_base.o
 dvpsi_e.o : ../Modules/kind.o
 dvpsi_e.o : ../Modules/uspp.o
 dvpsi_e.o : ../Modules/wavefunctions.o
+dvpsi_e.o : phcom.o
 dvpsi_e.o : ../PW/becmod.o
 dvpsi_e.o : ../PW/pwcom.o
-dvpsi_e.o : phcom.o
 dvqpsi_us.o : ../Modules/atom.o
 dvqpsi_us.o : ../Modules/ions_base.o
 dvqpsi_us.o : ../Modules/kind.o
 dvqpsi_us.o : ../Modules/wavefunctions.o
-dvqpsi_us.o : ../PW/pwcom.o
-dvqpsi_us.o : phcom.o
 dvqpsi_us_only.o : ../Modules/ions_base.o
 dvqpsi_us_only.o : ../Modules/kind.o
 dvqpsi_us_only.o : ../Modules/uspp.o
-dvqpsi_us_only.o : ../PW/pwcom.o
 dvqpsi_us_only.o : phcom.o
+dvqpsi_us_only.o : ../PW/pwcom.o
+dvqpsi_us.o : phcom.o
+dvqpsi_us.o : ../PW/pwcom.o
 dyndia.o : ../Modules/io_global.o
 dyndia.o : ../Modules/kind.o
 dynmat0.o : ../Modules/control_flags.o
 dynmat0.o : ../Modules/ions_base.o
 dynmat0.o : ../Modules/kind.o
-dynmat0.o : ../PW/pwcom.o
 dynmat0.o : phcom.o
+dynmat0.o : ../PW/pwcom.o
+dynmatcc.o : ../Modules/ions_base.o
+dynmatcc.o : ../Modules/kind.o
+dynmatcc.o : phcom.o
+dynmatcc.o : ../PW/pwcom.o
+dynmatrix.o : ../Modules/control_flags.o
+dynmatrix.o : ../Modules/io_global.o
+dynmatrix.o : ../Modules/ions_base.o
+dynmatrix.o : ../Modules/kind.o
+dynmatrix.o : phcom.o
+dynmatrix.o : ../PW/pwcom.o
 dynmat_us.o : ../Modules/io_files.o
 dynmat_us.o : ../Modules/ions_base.o
 dynmat_us.o : ../Modules/kind.o
 dynmat_us.o : ../Modules/mp_global.o
 dynmat_us.o : ../Modules/uspp.o
 dynmat_us.o : ../Modules/wavefunctions.o
-dynmat_us.o : ../PW/pwcom.o
 dynmat_us.o : phcom.o
-dynmatcc.o : ../Modules/ions_base.o
-dynmatcc.o : ../Modules/kind.o
-dynmatcc.o : ../PW/pwcom.o
-dynmatcc.o : phcom.o
-dynmatrix.o : ../Modules/control_flags.o
-dynmatrix.o : ../Modules/io_global.o
-dynmatrix.o : ../Modules/ions_base.o
-dynmatrix.o : ../Modules/kind.o
-dynmatrix.o : ../PW/pwcom.o
-dynmatrix.o : phcom.o
+dynmat_us.o : ../PW/pwcom.o
 ef_shift.o : ../Modules/cell_base.o
 ef_shift.o : ../Modules/io_global.o
 ef_shift.o : ../Modules/kind.o
 ef_shift.o : ../Modules/wavefunctions.o
-ef_shift.o : ../PW/pwcom.o
 ef_shift.o : phcom.o
+ef_shift.o : ../PW/pwcom.o
 elph.o : ../Modules/kind.o
+elphon.o : elph.o
 elphon.o : ../Modules/io_files.o
 elphon.o : ../Modules/ions_base.o
 elphon.o : ../Modules/kind.o
 elphon.o : ../Modules/mp_global.o
 elphon.o : ../Modules/wavefunctions.o
-elphon.o : ../PW/pwcom.o
-elphon.o : elph.o
 elphon.o : phcom.o
+elphon.o : ../PW/pwcom.o
 h_psiq.o : ../Modules/kind.o
 h_psiq.o : ../Modules/wavefunctions.o
+h_psiq.o : phcom.o
 h_psiq.o : ../PW/becmod.o
 h_psiq.o : ../PW/pwcom.o
-h_psiq.o : phcom.o
 incdrhoscf.o : ../Modules/ions_base.o
 incdrhoscf.o : ../Modules/kind.o
 incdrhoscf.o : ../Modules/uspp.o
 incdrhoscf.o : ../Modules/wavefunctions.o
-incdrhoscf.o : ../PW/pwcom.o
 incdrhoscf.o : phcom.o
+incdrhoscf.o : ../PW/pwcom.o
 incdrhous.o : ../Modules/ions_base.o
 incdrhous.o : ../Modules/kind.o
 incdrhous.o : ../Modules/uspp.o
-incdrhous.o : ../PW/pwcom.o
 incdrhous.o : phcom.o
+incdrhous.o : ../PW/pwcom.o
 io_pattern.o : ../Modules/io_global.o
 io_pattern.o : ../Modules/ions_base.o
 io_pattern.o : ../Modules/kind.o
@ -290,22 +290,22 @@ localdos.o : ../Modules/ions_base.o
 localdos.o : ../Modules/kind.o
 localdos.o : ../Modules/uspp.o
 localdos.o : ../Modules/wavefunctions.o
-localdos.o : ../PW/pwcom.o
 localdos.o : phcom.o
+localdos.o : ../PW/pwcom.o
 newdq.o : ../Modules/ions_base.o
 newdq.o : ../Modules/kind.o
 newdq.o : ../Modules/uspp.o
-newdq.o : ../PW/pwcom.o
 newdq.o : phcom.o
+newdq.o : ../PW/pwcom.o
 openfilq.o : ../Modules/control_flags.o
 openfilq.o : ../Modules/io_files.o
 openfilq.o : ../Modules/io_global.o
 openfilq.o : ../Modules/kind.o
 openfilq.o : ../Modules/mp_global.o
 openfilq.o : ../Modules/uspp.o
+openfilq.o : phcom.o
 openfilq.o : ../PW/pwcom.o
 openfilq.o : ../PW/restart.o
-openfilq.o : phcom.o
 phcom.o : ../Modules/kind.o
 phcom.o : ../Modules/parameters.o
 phonon.o : ../Modules/control_flags.o
@ -313,12 +313,12 @@ phonon.o : ../Modules/io_files.o
 phonon.o : ../Modules/io_global.o
 phonon.o : ../Modules/ions_base.o
 phonon.o : ../Modules/kind.o
-phonon.o : ../Modules/mp.o
 phonon.o : ../Modules/mp_global.o
+phonon.o : ../Modules/mp.o
 phonon.o : ../Modules/parser.o
 phonon.o : ../Modules/version.o
-phonon.o : ../PW/pwcom.o
 phonon.o : phcom.o
+phonon.o : ../PW/pwcom.o
 phq_init.o : ../Modules/atom.o
 phq_init.o : ../Modules/constants.o
 phq_init.o : ../Modules/io_files.o
@ -327,8 +327,8 @@ phq_init.o : ../Modules/ions_base.o
 phq_init.o : ../Modules/kind.o
 phq_init.o : ../Modules/uspp.o
 phq_init.o : ../Modules/wavefunctions.o
-phq_init.o : ../PW/pwcom.o
 phq_init.o : phcom.o
+phq_init.o : ../PW/pwcom.o
 phq_readin.o : ../Modules/check_stop.o
 phq_readin.o : ../Modules/constants.o
 phq_readin.o : ../Modules/control_flags.o
@ -337,91 +337,91 @@ phq_readin.o : ../Modules/io_global.o
 phq_readin.o : ../Modules/ions_base.o
 phq_readin.o : ../Modules/kind.o
 phq_readin.o : ../Modules/mp.o
-phq_readin.o : ../PW/pwcom.o
 phq_readin.o : phcom.o
+phq_readin.o : ../PW/pwcom.o
 phq_recover.o : ../Modules/control_flags.o
 phq_recover.o : ../Modules/io_global.o
 phq_recover.o : ../Modules/ions_base.o
 phq_recover.o : ../Modules/kind.o
-phq_recover.o : ../PW/pwcom.o
 phq_recover.o : phcom.o
+phq_recover.o : ../PW/pwcom.o
+phqscf.o : ../Modules/check_stop.o
+phqscf.o : ../Modules/io_global.o
+phqscf.o : ../Modules/kind.o
+phqscf.o : phcom.o
+phqscf.o : ../PW/pwcom.o
 phq_setup.o : ../Modules/atom.o
 phq_setup.o : ../Modules/constants.o
 phq_setup.o : ../Modules/control_flags.o
 phq_setup.o : ../Modules/io_global.o
 phq_setup.o : ../Modules/ions_base.o
 phq_setup.o : ../Modules/kind.o
-phq_setup.o : ../PW/pwcom.o
 phq_setup.o : phcom.o
+phq_setup.o : ../PW/pwcom.o
 phq_summary.o : ../Modules/atom.o
 phq_summary.o : ../Modules/control_flags.o
 phq_summary.o : ../Modules/io_global.o
 phq_summary.o : ../Modules/ions_base.o
 phq_summary.o : ../Modules/uspp.o
-phq_summary.o : ../PW/pwcom.o
 phq_summary.o : phcom.o
-phqscf.o : ../Modules/check_stop.o
-phqscf.o : ../Modules/io_global.o
-phqscf.o : ../Modules/kind.o
-phqscf.o : ../PW/pwcom.o
-phqscf.o : phcom.o
+phq_summary.o : ../PW/pwcom.o
 print_clock_ph.o : ../Modules/io_global.o
-print_clock_ph.o : ../PW/pwcom.o
 print_clock_ph.o : phcom.o
+print_clock_ph.o : ../PW/pwcom.o
 psidspsi.o : ../Modules/ions_base.o
 psidspsi.o : ../Modules/kind.o
 psidspsi.o : ../Modules/uspp.o
 psidspsi.o : ../Modules/wavefunctions.o
-psidspsi.o : ../PW/pwcom.o
 psidspsi.o : phcom.o
+psidspsi.o : ../PW/pwcom.o
 psymdvscf.o : ../Modules/kind.o
 psymdvscf.o : ../Modules/mp_global.o
+psymdvscf.o : phcom.o
 psymdvscf.o : ../PW/para.o
 psymdvscf.o : ../PW/pwcom.o
-psymdvscf.o : phcom.o
 psyme.o : ../Modules/kind.o
 psyme.o : ../Modules/mp_global.o
+psyme.o : phcom.o
 psyme.o : ../PW/para.o
 psyme.o : ../PW/pwcom.o
-psyme.o : phcom.o
 punch_plot_e.o : ../Modules/io_global.o
 punch_plot_e.o : ../Modules/ions_base.o
 punch_plot_e.o : ../Modules/kind.o
-punch_plot_e.o : ../PW/pwcom.o
 punch_plot_e.o : phcom.o
+punch_plot_e.o : ../PW/pwcom.o
 punch_plot_ph.o : ../Modules/io_global.o
 punch_plot_ph.o : ../Modules/ions_base.o
 punch_plot_ph.o : ../Modules/kind.o
-punch_plot_ph.o : ../PW/pwcom.o
 punch_plot_ph.o : phcom.o
+punch_plot_ph.o : ../PW/pwcom.o
 q2qstar_ph.o : ../Modules/kind.o
 q_points.o : ../Modules/cell_base.o
 q_points.o : ../Modules/io_global.o
 q_points.o : ../Modules/kind.o
-q_points.o : ../PW/pwcom.o
 q_points.o : phcom.o
+q_points.o : ../PW/pwcom.o
 random_matrix.o : ../Modules/kind.o
 rotate_and_add_dyn.o : ../Modules/kind.o
 set_drhoc.o : ../Modules/atom.o
 set_drhoc.o : ../Modules/ions_base.o
 set_drhoc.o : ../Modules/kind.o
 set_drhoc.o : ../Modules/parameters.o
-set_drhoc.o : ../PW/pwcom.o
 set_drhoc.o : phcom.o
+set_drhoc.o : ../PW/pwcom.o
+set_irr_mode.o : ../Modules/constants.o
+set_irr_mode.o : ../Modules/kind.o
+set_irr_nosym.o : ../Modules/kind.o
 set_irr.o : ../Modules/constants.o
 set_irr.o : ../Modules/io_global.o
 set_irr.o : ../Modules/kind.o
 set_irr.o : ../Modules/mp.o
-set_irr_mode.o : ../Modules/constants.o
-set_irr_mode.o : ../Modules/kind.o
-set_irr_nosym.o : ../Modules/kind.o
 setlocq.o : ../Modules/constants.o
 setlocq.o : ../Modules/kind.o
 setqmod.o : ../Modules/kind.o
 setup_dgc.o : ../Modules/functionals.o
 setup_dgc.o : ../Modules/kind.o
-setup_dgc.o : ../PW/pwcom.o
 setup_dgc.o : phcom.o
+setup_dgc.o : ../PW/pwcom.o
 smallgq.o : ../Modules/kind.o
 solve_e.o : ../Modules/check_stop.o
 solve_e.o : ../Modules/control_flags.o
@ -431,9 +431,9 @@ solve_e.o : ../Modules/ions_base.o
 solve_e.o : ../Modules/kind.o
 solve_e.o : ../Modules/uspp.o
 solve_e.o : ../Modules/wavefunctions.o
+solve_e.o : phcom.o
 solve_e.o : ../PW/becmod.o
 solve_e.o : ../PW/pwcom.o
-solve_e.o : phcom.o
 solve_linter.o : ../Modules/check_stop.o
 solve_linter.o : ../Modules/constants.o
 solve_linter.o : ../Modules/control_flags.o
@ -443,9 +443,9 @@ solve_linter.o : ../Modules/ions_base.o
 solve_linter.o : ../Modules/kind.o
 solve_linter.o : ../Modules/uspp.o
 solve_linter.o : ../Modules/wavefunctions.o
+solve_linter.o : phcom.o
 solve_linter.o : ../PW/becmod.o
 solve_linter.o : ../PW/pwcom.o
-solve_linter.o : phcom.o
 star_q.o : ../Modules/io_global.o
 star_q.o : ../Modules/kind.o
 stop_ph.o : ../Modules/kind.o
@ -454,20 +454,20 @@ stop_ph.o : ../Modules/parallel_include.o
 sym_and_write_zue.o : ../Modules/io_global.o
 sym_and_write_zue.o : ../Modules/ions_base.o
 sym_and_write_zue.o : ../Modules/kind.o
-sym_and_write_zue.o : ../PW/pwcom.o
 sym_and_write_zue.o : phcom.o
+sym_and_write_zue.o : ../PW/pwcom.o
 sym_def.o : ../Modules/kind.o
-sym_def.o : ../PW/pwcom.o
 sym_def.o : phcom.o
+sym_def.o : ../PW/pwcom.o
 symdvscf.o : ../Modules/kind.o
-symdvscf.o : ../PW/pwcom.o
 symdvscf.o : phcom.o
+symdvscf.o : ../PW/pwcom.o
 symdyn_munu.o : ../Modules/kind.o
 symdynph_gq.o : ../Modules/constants.o
 symdynph_gq.o : ../Modules/kind.o
 syme.o : ../Modules/kind.o
-syme.o : ../PW/pwcom.o
 syme.o : phcom.o
+syme.o : ../PW/pwcom.o
 tra_write_matrix.o : ../Modules/io_global.o
 tra_write_matrix.o : ../Modules/kind.o
 trntnsc.o : ../Modules/kind.o
@ -481,12 +481,12 @@ zstar_eu.o : ../Modules/io_global.o
 zstar_eu.o : ../Modules/ions_base.o
 zstar_eu.o : ../Modules/kind.o
 zstar_eu.o : ../Modules/wavefunctions.o
-zstar_eu.o : ../PW/pwcom.o
 zstar_eu.o : phcom.o
+zstar_eu.o : ../PW/pwcom.o
 zstar_eu_us.o : ../Modules/io_files.o
 zstar_eu_us.o : ../Modules/ions_base.o
 zstar_eu_us.o : ../Modules/kind.o
 zstar_eu_us.o : ../Modules/uspp.o
 zstar_eu_us.o : ../Modules/wavefunctions.o
-zstar_eu_us.o : ../PW/pwcom.o
 zstar_eu_us.o : phcom.o
+zstar_eu_us.o : ../PW/pwcom.o
--- a/PP/Makefile
+++ b/PP/Makefile
@ -53,6 +53,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/readpseudo.o \
 ../Modules/recvec.o \
--- a/PP/make.depend
+++ b/PP/make.depend
@ -34,8 +34,8 @@ bands.o : ../Modules/constants.o
 bands.o : ../Modules/io_files.o
 bands.o : ../Modules/io_global.o
 bands.o : ../Modules/ions_base.o
-bands.o : ../Modules/mp.o
 bands.o : ../Modules/mp_global.o
+bands.o : ../Modules/mp.o
 bands.o : ../Modules/uspp.o
 bands.o : ../Modules/wavefunctions.o
 bands.o : ../PW/pwcom.o
@ -58,6 +58,7 @@ do_initial_state.o : ../Modules/kind.o
 do_initial_state.o : ../Modules/parameters.o
 do_initial_state.o : ../Modules/uspp.o
 do_initial_state.o : ../PW/pwcom.o
+dosg.o : ../Modules/kind.o
 do_shift_ew.o : ../Modules/kind.o
 dos.o : ../Modules/constants.o
 dos.o : ../Modules/io_files.o
@ -65,7 +66,6 @@ dos.o : ../Modules/io_global.o
 dos.o : ../Modules/kind.o
 dos.o : ../Modules/mp.o
 dos.o : ../PW/pwcom.o
-dosg.o : ../Modules/kind.o
 efg.o : ../Modules/atom.o
 efg.o : ../Modules/cell_base.o
 efg.o : ../Modules/constants.o
@ -97,6 +97,12 @@ initial_state.o : ../Modules/parameters.o
 initial_state.o : ../Modules/uspp.o
 initial_state.o : ../Modules/wavefunctions.o
 initial_state.o : ../PW/pwcom.o
+local_dos1d.o : ../Modules/cell_base.o
+local_dos1d.o : ../Modules/ions_base.o
+local_dos1d.o : ../Modules/uspp.o
+local_dos1d.o : ../Modules/wavefunctions.o
+local_dos1d.o : ../PW/becmod.o
+local_dos1d.o : ../PW/pwcom.o
 local_dos.o : ../Modules/cell_base.o
 local_dos.o : ../Modules/io_files.o
 local_dos.o : ../Modules/ions_base.o
@ -105,12 +111,6 @@ local_dos.o : ../Modules/mp.o
 local_dos.o : ../Modules/uspp.o
 local_dos.o : ../Modules/wavefunctions.o
 local_dos.o : ../PW/pwcom.o
-local_dos1d.o : ../Modules/cell_base.o
-local_dos1d.o : ../Modules/ions_base.o
-local_dos1d.o : ../Modules/uspp.o
-local_dos1d.o : ../Modules/wavefunctions.o
-local_dos1d.o : ../PW/becmod.o
-local_dos1d.o : ../PW/pwcom.o
 openfil_pp.o : ../Modules/control_flags.o
 openfil_pp.o : ../Modules/io_files.o
 openfil_pp.o : ../Modules/io_global.o
@ -126,12 +126,12 @@ plan_avg.o : ../Modules/uspp.o
 plan_avg.o : ../Modules/wavefunctions.o
 plan_avg.o : ../PW/becmod.o
 plan_avg.o : ../PW/pwcom.o
+plotband.o : ../Modules/io_global.o
 plot_io.o : ../Modules/io_global.o
 plot_io.o : ../Modules/kind.o
-plot_whole_cell.o : ../Modules/kind.o
-plotband.o : ../Modules/io_global.o
 plotrho.o : ../Modules/io_global.o
 plotrho.o : ../Modules/kind.o
+plot_whole_cell.o : ../Modules/kind.o
 poormanwannier.o : ../Modules/atom.o
 poormanwannier.o : ../Modules/cell_base.o
 poormanwannier.o : ../Modules/constants.o
@ -187,8 +187,8 @@ pw2wan.o : ../Modules/io_base.o
 pw2wan.o : ../Modules/io_files.o
 pw2wan.o : ../Modules/io_global.o
 pw2wan.o : ../Modules/ions_base.o
-pw2wan.o : ../Modules/mp.o
 pw2wan.o : ../Modules/mp_global.o
+pw2wan.o : ../Modules/mp.o
 pw2wan.o : ../Modules/wavefunctions.o
 pw2wan.o : ../PW/pwcom.o
 pw_export.o : ../Modules/io_files.o
@ -196,8 +196,8 @@ pw_export.o : ../Modules/io_global.o
 pw_export.o : ../Modules/ions_base.o
 pw_export.o : ../Modules/iotk_module.o
 pw_export.o : ../Modules/kind.o
-pw_export.o : ../Modules/mp.o
 pw_export.o : ../Modules/mp_global.o
+pw_export.o : ../Modules/mp.o
 pw_export.o : ../Modules/mp_wave.o
 pw_export.o : ../Modules/parameters.o
 pw_export.o : ../Modules/wavefunctions.o
--- a/PW/Makefile
+++ b/PW/Makefile
@ -278,6 +278,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/read_cards.o \
 ../Modules/read_namelists.o \
--- a/PW/make.depend
+++ b/PW/make.depend
@ -4,15 +4,10 @@ add_efield.o : ../Modules/control_flags.o
 add_efield.o : ../Modules/io_global.o
 add_efield.o : ../Modules/ions_base.o
 add_efield.o : ../Modules/kind.o
-add_efield.o : ../Modules/mp.o
 add_efield.o : ../Modules/mp_global.o
+add_efield.o : ../Modules/mp.o
 add_efield.o : para.o
 add_efield.o : pwcom.o
-add_vuspsi.o : ../Modules/ions_base.o
-add_vuspsi.o : ../Modules/kind.o
-add_vuspsi.o : ../Modules/uspp.o
-add_vuspsi.o : becmod.o
-add_vuspsi.o : pwcom.o
 addusdens.o : ../Modules/ions_base.o
 addusdens.o : ../Modules/kind.o
 addusdens.o : ../Modules/uspp.o
@ -28,6 +23,11 @@ addusstress.o : ../Modules/ions_base.o
 addusstress.o : ../Modules/kind.o
 addusstress.o : ../Modules/uspp.o
 addusstress.o : pwcom.o
+add_vuspsi.o : becmod.o
+add_vuspsi.o : ../Modules/ions_base.o
+add_vuspsi.o : ../Modules/kind.o
+add_vuspsi.o : ../Modules/uspp.o
+add_vuspsi.o : pwcom.o
 allocate_fft.o : ../Modules/io_global.o
 allocate_fft.o : ../Modules/ions_base.o
 allocate_fft.o : ../Modules/wavefunctions.o
@ -70,6 +70,13 @@ bfgs.o : ../Modules/io_global.o
 bfgs.o : ../Modules/ions_base.o
 bfgs.o : ../Modules/kind.o
 bfgs.o : pwcom.o
+bp_calc_btq.o : ../Modules/atom.o
+bp_calc_btq.o : ../Modules/cell_base.o
+bp_calc_btq.o : ../Modules/constants.o
+bp_calc_btq.o : ../Modules/ions_base.o
+bp_calc_btq.o : ../Modules/kind.o
+bp_calc_btq.o : ../Modules/parameters.o
+bp_calc_btq.o : ../Modules/uspp.o
 bp_c_phase.o : ../Modules/cell_base.o
 bp_c_phase.o : ../Modules/constants.o
 bp_c_phase.o : ../Modules/io_files.o
@ -80,18 +87,14 @@ bp_c_phase.o : ../Modules/parameters.o
 bp_c_phase.o : ../Modules/uspp.o
 bp_c_phase.o : ../Modules/wavefunctions.o
 bp_c_phase.o : pwcom.o
-bp_calc_btq.o : ../Modules/atom.o
-bp_calc_btq.o : ../Modules/cell_base.o
-bp_calc_btq.o : ../Modules/constants.o
-bp_calc_btq.o : ../Modules/ions_base.o
-bp_calc_btq.o : ../Modules/kind.o
-bp_calc_btq.o : ../Modules/parameters.o
-bp_calc_btq.o : ../Modules/uspp.o
 bp_qvan3.o : ../Modules/ions_base.o
 bp_qvan3.o : ../Modules/kind.o
 bp_qvan3.o : ../Modules/uspp.o
 bp_qvan3.o : pwcom.o
 bp_strings.o : ../Modules/kind.o
+c_bands.o : becmod.o
+c_bands.o : complex_diis_module.o
+c_bands.o : g_psi_mod.o
 c_bands.o : ../Modules/cell_base.o
 c_bands.o : ../Modules/constants.o
 c_bands.o : ../Modules/control_flags.o
@ -100,25 +103,24 @@ c_bands.o : ../Modules/io_global.o
 c_bands.o : ../Modules/kind.o
 c_bands.o : ../Modules/uspp.o
 c_bands.o : ../Modules/wavefunctions.o
-c_bands.o : becmod.o
-c_bands.o : complex_diis_module.o
-c_bands.o : g_psi_mod.o
 c_bands.o : pwcom.o
 c_bands.o : real_diis_module.o
 ccalbec.o : ../Modules/kind.o
 ccalbec.o : pwcom.o
 ccgdiagg.o : ../Modules/io_global.o
 ccgdiagg.o : ../Modules/kind.o
-cdiagh.o : ../Modules/kind.o
-cdiagh.o : ../Modules/mp.o
-cdiagh.o : ../Modules/mp_global.o
 cdiaghg.o : ../Modules/kind.o
-cdiaghg.o : ../Modules/mp.o
 cdiaghg.o : ../Modules/mp_global.o
+cdiaghg.o : ../Modules/mp.o
+cdiagh.o : ../Modules/kind.o
+cdiagh.o : ../Modules/mp_global.o
+cdiagh.o : ../Modules/mp.o
 cegterg.o : ../Modules/io_global.o
 cegterg.o : ../Modules/kind.o
 cft3.o : ../Modules/fft_base.o
+cft_3.o : ../Modules/fft_scalar.o
 cft3.o : ../Modules/fft_scalar.o
+cft_3.o : ../Modules/kind.o
 cft3.o : ../Modules/kind.o
 cft3.o : ../Modules/mp_global.o
 cft3.o : para.o
@ -129,13 +131,11 @@ cft3s.o : ../Modules/kind.o
 cft3s.o : ../Modules/mp_global.o
 cft3s.o : para.o
 cft3s.o : pwcom.o
-cft_3.o : ../Modules/fft_scalar.o
-cft_3.o : ../Modules/kind.o
-cft_fftw.o : ../Modules/kind.o
 cft_fftw.o : fftw.o
+cft_fftw.o : ../Modules/kind.o
+cfts_3.o : ../Modules/kind.o
 cft_sgi.o : ../Modules/kind.o
 cft_sun.o : ../Modules/kind.o
-cfts_3.o : ../Modules/kind.o
 cgramg1.o : ../Modules/constants.o
 cgramg1.o : ../Modules/io_global.o
 cgramg1.o : ../Modules/kind.o
@ -143,6 +143,8 @@ cgramg1.o : pwcom.o
 checkallsym.o : ../Modules/kind.o
 checksym.o : ../Modules/kind.o
 cinitcgg.o : ../Modules/kind.o
+clean_pw.o : cft_3.o
+clean_pw.o : constraints_module.o
 clean_pw.o : ../Modules/berry_phase.o
 clean_pw.o : ../Modules/fft_types.o
 clean_pw.o : ../Modules/ions_base.o
@ -150,19 +152,17 @@ clean_pw.o : ../Modules/recvec.o
 clean_pw.o : ../Modules/stick_base.o
 clean_pw.o : ../Modules/uspp.o
 clean_pw.o : ../Modules/wavefunctions.o
-clean_pw.o : cft_3.o
-clean_pw.o : constraints_module.o
 clean_pw.o : noncol.o
 clean_pw.o : pwcom.o
 close_files.o : ../Modules/control_flags.o
 close_files.o : ../Modules/io_files.o
-close_files.o : ../Modules/mp.o
 close_files.o : ../Modules/mp_global.o
+close_files.o : ../Modules/mp.o
 close_files.o : pwcom.o
-complex_diis_module.o : ../Modules/constants.o
-complex_diis_module.o : ../Modules/control_flags.o
 complex_diis_module.o : diis_base.o
 complex_diis_module.o : g_psi_mod.o
+complex_diis_module.o : ../Modules/constants.o
+complex_diis_module.o : ../Modules/control_flags.o
 compute_dip.o : ../Modules/cell_base.o
 compute_dip.o : ../Modules/constants.o
 compute_dip.o : ../Modules/io_global.o
@ -179,8 +179,8 @@ compute_scf.o : ../Modules/io_files.o
 compute_scf.o : ../Modules/io_global.o
 compute_scf.o : ../Modules/ions_base.o
 compute_scf.o : ../Modules/kind.o
-compute_scf.o : ../Modules/mp.o
 compute_scf.o : ../Modules/mp_global.o
+compute_scf.o : ../Modules/mp.o
 compute_scf.o : ../Modules/parser.o
 compute_scf.o : ../Modules/path_formats.o
 compute_scf.o : ../Modules/path_variables.o
@ -194,16 +194,13 @@ constraints_module.o : pwcom.o
 coset.o : ../Modules/kind.o
 cryst_to_car.o : ../Modules/kind.o
 cubicsym.o : ../Modules/kind.o
-d_matrix.o : ../Modules/cell_base.o
-d_matrix.o : ../Modules/kind.o
-d_matrix.o : pwcom.o
 data_structure.o : ../Modules/cell_base.o
 data_structure.o : ../Modules/fft_scalar.o
 data_structure.o : ../Modules/fft_types.o
 data_structure.o : ../Modules/io_global.o
 data_structure.o : ../Modules/kind.o
-data_structure.o : ../Modules/mp.o
 data_structure.o : ../Modules/mp_global.o
+data_structure.o : ../Modules/mp.o
 data_structure.o : ../Modules/stick_base.o
 data_structure.o : para.o
 data_structure.o : pwcom.o
@ -217,44 +214,47 @@ diis_base.o : ../Modules/kind.o
 diropn.o : ../Modules/io_files.o
 diropn.o : ../Modules/kind.o
 diropn.o : ../Modules/mp_global.o
-divide.o : ../Modules/mp_global.o
 divide_et_impera.o : ../Modules/io_global.o
 divide_et_impera.o : ../Modules/kind.o
 divide_et_impera.o : ../Modules/mp_global.o
+divide.o : ../Modules/mp_global.o
+d_matrix.o : ../Modules/cell_base.o
+d_matrix.o : ../Modules/kind.o
+d_matrix.o : pwcom.o
+dndepsilon.o : becmod.o
 dndepsilon.o : ../Modules/atom.o
 dndepsilon.o : ../Modules/io_files.o
 dndepsilon.o : ../Modules/ions_base.o
 dndepsilon.o : ../Modules/kind.o
 dndepsilon.o : ../Modules/uspp.o
 dndepsilon.o : ../Modules/wavefunctions.o
-dndepsilon.o : becmod.o
 dndepsilon.o : pwcom.o
+dndtau.o : becmod.o
+dndtau_of_k.o : ../Modules/ions_base.o
+dndtau_of_k.o : ../Modules/kind.o
+dndtau_of_k.o : pwcom.o
 dndtau.o : ../Modules/atom.o
 dndtau.o : ../Modules/io_files.o
 dndtau.o : ../Modules/ions_base.o
 dndtau.o : ../Modules/kind.o
 dndtau.o : ../Modules/uspp.o
 dndtau.o : ../Modules/wavefunctions.o
-dndtau.o : becmod.o
 dndtau.o : pwcom.o
-dndtau_of_k.o : ../Modules/ions_base.o
-dndtau_of_k.o : ../Modules/kind.o
-dndtau_of_k.o : pwcom.o
+dprojdepsilon.o : becmod.o
 dprojdepsilon.o : ../Modules/atom.o
 dprojdepsilon.o : ../Modules/cell_base.o
 dprojdepsilon.o : ../Modules/ions_base.o
 dprojdepsilon.o : ../Modules/kind.o
 dprojdepsilon.o : ../Modules/uspp.o
 dprojdepsilon.o : ../Modules/wavefunctions.o
-dprojdepsilon.o : becmod.o
 dprojdepsilon.o : pwcom.o
+dprojdtau.o : becmod.o
 dprojdtau.o : ../Modules/atom.o
 dprojdtau.o : ../Modules/cell_base.o
 dprojdtau.o : ../Modules/ions_base.o
 dprojdtau.o : ../Modules/kind.o
 dprojdtau.o : ../Modules/uspp.o
 dprojdtau.o : ../Modules/wavefunctions.o
-dprojdtau.o : becmod.o
 dprojdtau.o : pwcom.o
 dqvan2.o : ../Modules/kind.o
 dqvan2.o : ../Modules/parameters.o
@ -289,12 +289,12 @@ electrons.o : para.o
 electrons.o : pwcom.o
 eqvect.o : ../Modules/kind.o
 estimate.o : ../Modules/kind.o
-ewald.o : ../Modules/kind.o
 ewald_dipole.o : ../Modules/cell_base.o
 ewald_dipole.o : ../Modules/constants.o
 ewald_dipole.o : ../Modules/ions_base.o
 ewald_dipole.o : ../Modules/kind.o
 ewald_dipole.o : pwcom.o
+ewald.o : ../Modules/kind.o
 force_cc.o : ../Modules/atom.o
 force_cc.o : ../Modules/cell_base.o
 force_cc.o : ../Modules/constants.o
@ -310,6 +310,7 @@ force_corr.o : ../Modules/kind.o
 force_corr.o : ../Modules/wavefunctions.o
 force_corr.o : pwcom.o
 force_ew.o : ../Modules/kind.o
+force_hub.o : becmod.o
 force_hub.o : ../Modules/atom.o
 force_hub.o : ../Modules/cell_base.o
 force_hub.o : ../Modules/io_files.o
@ -318,9 +319,13 @@ force_hub.o : ../Modules/kind.o
 force_hub.o : ../Modules/mp_global.o
 force_hub.o : ../Modules/uspp.o
 force_hub.o : ../Modules/wavefunctions.o
-force_hub.o : becmod.o
 force_hub.o : pwcom.o
 force_lc.o : ../Modules/kind.o
+forces.o : ../Modules/cell_base.o
+forces.o : ../Modules/io_global.o
+forces.o : ../Modules/ions_base.o
+forces.o : ../Modules/kind.o
+forces.o : pwcom.o
 force_us.o : ../Modules/cell_base.o
 force_us.o : ../Modules/io_files.o
 force_us.o : ../Modules/ions_base.o
@ -328,17 +333,8 @@ force_us.o : ../Modules/kind.o
 force_us.o : ../Modules/uspp.o
 force_us.o : ../Modules/wavefunctions.o
 force_us.o : pwcom.o
-forces.o : ../Modules/cell_base.o
-forces.o : ../Modules/io_global.o
-forces.o : ../Modules/ions_base.o
-forces.o : ../Modules/kind.o
-forces.o : pwcom.o
 functionals.o : ../Modules/functionals.o
 functionals.o : ../Modules/kind.o
-g_psi.o : ../Modules/constants.o
-g_psi.o : ../Modules/kind.o
-g_psi.o : g_psi_mod.o
-g_psi_mod.o : ../Modules/kind.o
 gen_at_dj.o : ../Modules/atom.o
 gen_at_dj.o : ../Modules/cell_base.o
 gen_at_dj.o : ../Modules/constants.o
@ -378,16 +374,15 @@ ggen.o : pwcom.o
 gk_sort.o : ../Modules/constants.o
 gk_sort.o : ../Modules/kind.o
 gk_sort.o : pwcom.o
+g_psi_mod.o : ../Modules/kind.o
+g_psi.o : g_psi_mod.o
+g_psi.o : ../Modules/constants.o
+g_psi.o : ../Modules/kind.o
 gradcorr.o : ../Modules/functionals.o
 gradcorr.o : ../Modules/kind.o
 gradcorr.o : pwcom.o
 gweights.o : ../Modules/kind.o
 h_1psi.o : ../Modules/kind.o
-h_psi.o : ../Modules/kind.o
-h_psi.o : ../Modules/uspp.o
-h_psi.o : ../Modules/wavefunctions.o
-h_psi.o : becmod.o
-h_psi.o : pwcom.o
 hexsym.o : ../Modules/kind.o
 hinit0.o : ../Modules/cell_base.o
 hinit0.o : ../Modules/io_files.o
@ -397,6 +392,11 @@ hinit1.o : ../Modules/cell_base.o
 hinit1.o : ../Modules/control_flags.o
 hinit1.o : ../Modules/ions_base.o
 hinit1.o : pwcom.o
+h_psi.o : becmod.o
+h_psi.o : ../Modules/kind.o
+h_psi.o : ../Modules/uspp.o
+h_psi.o : ../Modules/wavefunctions.o
+h_psi.o : pwcom.o
 init_at_1.o : ../Modules/atom.o
 init_at_1.o : ../Modules/cell_base.o
 init_at_1.o : ../Modules/constants.o
@ -422,8 +422,8 @@ init_paw_2.o : ../Modules/ions_base.o
 init_paw_2.o : ../Modules/kind.o
 init_paw_2.o : paw.o
 init_paw_2.o : pwcom.o
-init_pool.o : ../Modules/mp.o
 init_pool.o : ../Modules/mp_global.o
+init_pool.o : ../Modules/mp.o
 init_pool.o : ../Modules/parallel_include.o
 init_run.o : ../Modules/ions_base.o
 init_run.o : ../Modules/mp_global.o
@ -447,6 +447,7 @@ init_vloc.o : ../Modules/cell_base.o
 init_vloc.o : ../Modules/ions_base.o
 init_vloc.o : ../Modules/uspp.o
 init_vloc.o : pwcom.o
+input.o : constraints_module.o
 input.o : ../Modules/basic_algebra_routines.o
 input.o : ../Modules/bfgs_module.o
 input.o : ../Modules/cell_base.o
@ -457,23 +458,22 @@ input.o : ../Modules/input_parameters.o
 input.o : ../Modules/io_files.o
 input.o : ../Modules/io_global.o
 input.o : ../Modules/ions_base.o
-input.o : ../Modules/mp.o
 input.o : ../Modules/mp_global.o
+input.o : ../Modules/mp.o
 input.o : ../Modules/parser.o
 input.o : ../Modules/path_variables.o
 input.o : ../Modules/read_cards.o
 input.o : ../Modules/read_namelists.o
-input.o : constraints_module.o
 input.o : noncol.o
 input.o : pwcom.o
 interpolate.o : ../Modules/kind.o
 interpolate.o : pwcom.o
-io_pot.o : ../Modules/kind.o
-io_pot.o : ../Modules/mp.o
-io_pot.o : ../Modules/mp_global.o
-io_pot.o : pwcom.o
 ions.o : ../Modules/control_flags.o
 ions.o : pwcom.o
+io_pot.o : ../Modules/kind.o
+io_pot.o : ../Modules/mp_global.o
+io_pot.o : ../Modules/mp.o
+io_pot.o : pwcom.o
 irrek.o : ../Modules/kind.o
 iweights.o : ../Modules/kind.o
 iweights.o : noncol.o
@ -498,8 +498,8 @@ memory.o : ../Modules/control_flags.o
 memory.o : ../Modules/io_files.o
 memory.o : ../Modules/ions_base.o
 memory.o : ../Modules/kind.o
-memory.o : ../Modules/mp.o
 memory.o : ../Modules/mp_global.o
+memory.o : ../Modules/mp.o
 memory.o : ../Modules/parameters.o
 memory.o : ../Modules/uspp.o
 memory.o : ../Modules/version.o
@ -514,6 +514,7 @@ mix_rho.o : ../Modules/kind.o
 mix_rho.o : ../Modules/wavefunctions.o
 mix_rho.o : pwcom.o
 mode_group.o : ../Modules/kind.o
+move_ions.o : constraints_module.o
 move_ions.o : ../Modules/basic_algebra_routines.o
 move_ions.o : ../Modules/bfgs_module.o
 move_ions.o : ../Modules/cell_base.o
@ -523,19 +524,10 @@ move_ions.o : ../Modules/io_files.o
 move_ions.o : ../Modules/io_global.o
 move_ions.o : ../Modules/ions_base.o
 move_ions.o : ../Modules/kind.o
-move_ions.o : ../Modules/mp.o
 move_ions.o : ../Modules/mp_global.o
-move_ions.o : constraints_module.o
+move_ions.o : ../Modules/mp.o
 move_ions.o : pwcom.o
 multable.o : ../Modules/kind.o
-n_plane_waves.o : ../Modules/kind.o
-new_ns.o : ../Modules/atom.o
-new_ns.o : ../Modules/io_files.o
-new_ns.o : ../Modules/io_global.o
-new_ns.o : ../Modules/ions_base.o
-new_ns.o : ../Modules/kind.o
-new_ns.o : ../Modules/wavefunctions.o
-new_ns.o : pwcom.o
 newd.o : ../Modules/cell_base.o
 newd.o : ../Modules/ions_base.o
 newd.o : ../Modules/kind.o
@ -543,8 +535,16 @@ newd.o : ../Modules/uspp.o
 newd.o : ../Modules/wavefunctions.o
 newd.o : noncol.o
 newd.o : pwcom.o
+new_ns.o : ../Modules/atom.o
+new_ns.o : ../Modules/io_files.o
+new_ns.o : ../Modules/io_global.o
+new_ns.o : ../Modules/ions_base.o
+new_ns.o : ../Modules/kind.o
+new_ns.o : ../Modules/wavefunctions.o
+new_ns.o : pwcom.o
 noncol.o : ../Modules/kind.o
 noncol.o : ../Modules/parameters.o
+n_plane_waves.o : ../Modules/kind.o
 ns_adj.o : ../Modules/io_global.o
 ns_adj.o : ../Modules/ions_base.o
 ns_adj.o : ../Modules/kind.o
@ -557,20 +557,20 @@ openfil.o : ../Modules/mp_global.o
 openfil.o : noncol.o
 openfil.o : pwcom.o
 openfil.o : restart.o
-ortho.o : ../Modules/io_files.o
-ortho.o : ../Modules/kind.o
-ortho.o : ../Modules/uspp.o
-ortho.o : ../Modules/wavefunctions.o
-ortho.o : becmod.o
-ortho.o : pwcom.o
+orthoatwfc.o : becmod.o
 orthoatwfc.o : ../Modules/io_files.o
 orthoatwfc.o : ../Modules/io_global.o
 orthoatwfc.o : ../Modules/ions_base.o
 orthoatwfc.o : ../Modules/kind.o
 orthoatwfc.o : ../Modules/parameters.o
 orthoatwfc.o : ../Modules/uspp.o
-orthoatwfc.o : becmod.o
 orthoatwfc.o : pwcom.o
+ortho.o : becmod.o
+ortho.o : ../Modules/io_files.o
+ortho.o : ../Modules/kind.o
+ortho.o : ../Modules/uspp.o
+ortho.o : ../Modules/wavefunctions.o
+ortho.o : pwcom.o
 output_tau.o : ../Modules/cell_base.o
 output_tau.o : ../Modules/constants.o
 output_tau.o : ../Modules/io_global.o
@ -578,8 +578,8 @@ output_tau.o : ../Modules/ions_base.o
 output_tau.o : ../Modules/kind.o
 output_tau.o : pwcom.o
 para.o : ../Modules/kind.o
-para.o : ../Modules/mp.o
 para.o : ../Modules/mp_global.o
+para.o : ../Modules/mp.o
 para.o : ../Modules/parallel_include.o
 paw.o : ../Modules/atom.o
 paw.o : ../Modules/io_global.o
@ -593,8 +593,8 @@ potinit.o : ../Modules/io_files.o
 potinit.o : ../Modules/io_global.o
 potinit.o : ../Modules/ions_base.o
 potinit.o : ../Modules/kind.o
-potinit.o : ../Modules/mp.o
 potinit.o : ../Modules/mp_global.o
+potinit.o : ../Modules/mp.o
 potinit.o : pwcom.o
 print_clock_pw.o : ../Modules/control_flags.o
 print_clock_pw.o : ../Modules/io_global.o
@ -611,16 +611,17 @@ punch.o : ../Modules/wavefunctions.o
 punch.o : noncol.o
 punch.o : pwcom.o
 punch.o : restart.o
-pw_gemm.o : ../Modules/kind.o
-pw_gemm.o : pwcom.o
 pwcom.o : ../Modules/cell_base.o
 pwcom.o : ../Modules/constants.o
 pwcom.o : ../Modules/fft_types.o
 pwcom.o : ../Modules/ions_base.o
 pwcom.o : ../Modules/kind.o
 pwcom.o : ../Modules/parameters.o
+pwcom.o : ../Modules/printout_base.o
 pwcom.o : ../Modules/recvec.o
 pwcom.o : ../Modules/uspp.o
+pw_gemm.o : ../Modules/kind.o
+pw_gemm.o : pwcom.o
 pwscf.o : ../Modules/control_flags.o
 pwscf.o : ../Modules/io_files.o
 pwscf.o : ../Modules/io_global.o
@ -635,8 +636,8 @@ qvan2.o : ../Modules/kind.o
 qvan2.o : ../Modules/uspp.o
 qvan2.o : pwcom.o
 rdiaghg.o : ../Modules/kind.o
-rdiaghg.o : ../Modules/mp.o
 rdiaghg.o : ../Modules/mp_global.o
+rdiaghg.o : ../Modules/mp.o
 read_conf_from_file.o : ../Modules/cell_base.o
 read_conf_from_file.o : ../Modules/io_files.o
 read_conf_from_file.o : ../Modules/io_global.o
@ -652,20 +653,17 @@ read_file.o : ../Modules/mp_global.o
 read_file.o : noncol.o
 read_file.o : pwcom.o
 read_file.o : restart.o
+readin.o : ../Modules/atom.o
+readin.o : ../Modules/functionals.o
+readin.o : ../Modules/io_files.o
+readin.o : ../Modules/ions_base.o
+readin.o : ../Modules/uspp.o
 read_ncpp.o : ../Modules/atom.o
 read_ncpp.o : ../Modules/functionals.o
 read_ncpp.o : ../Modules/kind.o
 read_ncpp.o : ../Modules/parameters.o
 read_ncpp.o : ../Modules/uspp.o
 read_ncpp.o : pwcom.o
-read_pseudo.o : ../Modules/pseudo_types.o
-read_pseudo.o : ../Modules/readpseudo.o
-read_pseudo.o : upf_to_internal.o
-readin.o : ../Modules/atom.o
-readin.o : ../Modules/functionals.o
-readin.o : ../Modules/io_files.o
-readin.o : ../Modules/ions_base.o
-readin.o : ../Modules/uspp.o
 readnewvan.o : ../Modules/atom.o
 readnewvan.o : ../Modules/constants.o
 readnewvan.o : ../Modules/functionals.o
@ -673,16 +671,19 @@ readnewvan.o : ../Modules/kind.o
 readnewvan.o : ../Modules/parameters.o
 readnewvan.o : ../Modules/uspp.o
 readnewvan.o : pwcom.o
+read_pseudo.o : ../Modules/pseudo_types.o
+read_pseudo.o : ../Modules/readpseudo.o
+read_pseudo.o : upf_to_internal.o
 readvan.o : ../Modules/atom.o
 readvan.o : ../Modules/functionals.o
 readvan.o : ../Modules/kind.o
 readvan.o : ../Modules/parameters.o
 readvan.o : ../Modules/uspp.o
 readvan.o : pwcom.o
-real_diis_module.o : ../Modules/constants.o
-real_diis_module.o : ../Modules/control_flags.o
 real_diis_module.o : diis_base.o
 real_diis_module.o : g_psi_mod.o
+real_diis_module.o : ../Modules/constants.o
+real_diis_module.o : ../Modules/control_flags.o
 real_diis_module.o : pwcom.o
 regterg.o : ../Modules/io_global.o
 regterg.o : ../Modules/kind.o
@ -692,26 +693,6 @@ remove_atomic_rho.o : ../Modules/kind.o
 remove_atomic_rho.o : pwcom.o
 reset_k_points.o : ../Modules/input_parameters.o
 reset_k_points.o : pwcom.o
-restart.o : ../Modules/atom.o
-restart.o : ../Modules/cell_base.o
-restart.o : ../Modules/constants.o
-restart.o : ../Modules/control_flags.o
-restart.o : ../Modules/functionals.o
-restart.o : ../Modules/io_base.o
-restart.o : ../Modules/io_files.o
-restart.o : ../Modules/io_global.o
-restart.o : ../Modules/ions_base.o
-restart.o : ../Modules/kind.o
-restart.o : ../Modules/mp.o
-restart.o : ../Modules/mp_global.o
-restart.o : ../Modules/parameters.o
-restart.o : ../Modules/pseudo_types.o
-restart.o : ../Modules/readpseudo.o
-restart.o : ../Modules/uspp.o
-restart.o : ../Modules/wavefunctions.o
-restart.o : noncol.o
-restart.o : pwcom.o
-restart.o : upf_to_internal.o
 restart_from_file.o : ../Modules/control_flags.o
 restart_from_file.o : ../Modules/io_files.o
 restart_from_file.o : ../Modules/io_global.o
@ -731,23 +712,37 @@ restart_in_ions.o : ../Modules/io_global.o
 restart_in_ions.o : ../Modules/kind.o
 restart_in_ions.o : ../Modules/wavefunctions.o
 restart_in_ions.o : pwcom.o
+restart.o : ../Modules/atom.o
+restart.o : ../Modules/cell_base.o
+restart.o : ../Modules/constants.o
+restart.o : ../Modules/control_flags.o
+restart.o : ../Modules/functionals.o
+restart.o : ../Modules/io_base.o
+restart.o : ../Modules/io_files.o
+restart.o : ../Modules/io_global.o
+restart.o : ../Modules/ions_base.o
+restart.o : ../Modules/kind.o
+restart.o : ../Modules/mp_global.o
+restart.o : ../Modules/mp.o
+restart.o : ../Modules/parameters.o
+restart.o : ../Modules/pseudo_types.o
+restart.o : ../Modules/readpseudo.o
+restart.o : ../Modules/uspp.o
+restart.o : ../Modules/wavefunctions.o
+restart.o : noncol.o
+restart.o : pwcom.o
+restart.o : upf_to_internal.o
 rgen.o : ../Modules/kind.o
 rho2zeta.o : ../Modules/io_global.o
 rho2zeta.o : ../Modules/kind.o
-rotate_wfc.o : ../Modules/kind.o
 rotate_wfc_gamma.o : ../Modules/kind.o
 rotate_wfc_gamma.o : pwcom.o
+rotate_wfc.o : ../Modules/kind.o
 ruotaijk.o : ../Modules/kind.o
+s_1psi.o : becmod.o
 s_1psi.o : ../Modules/kind.o
 s_1psi.o : ../Modules/uspp.o
-s_1psi.o : becmod.o
 s_1psi.o : pwcom.o
-s_axis_to_ca.o : ../Modules/kind.o
-s_psi.o : ../Modules/ions_base.o
-s_psi.o : ../Modules/kind.o
-s_psi.o : ../Modules/uspp.o
-s_psi.o : becmod.o
-s_psi.o : pwcom.o
 save_in_cbands.o : ../Modules/control_flags.o
 save_in_cbands.o : ../Modules/io_files.o
 save_in_cbands.o : ../Modules/kind.o
@ -760,10 +755,11 @@ save_in_ions.o : ../Modules/control_flags.o
 save_in_ions.o : ../Modules/io_files.o
 save_in_ions.o : ../Modules/kind.o
 save_in_ions.o : pwcom.o
-scala_cdiag.o : ../Modules/io_global.o
-scala_cdiag.o : ../Modules/kind.o
+s_axis_to_ca.o : ../Modules/kind.o
 scala_cdiaghg.o : ../Modules/io_global.o
 scala_cdiaghg.o : ../Modules/kind.o
+scala_cdiag.o : ../Modules/io_global.o
+scala_cdiag.o : ../Modules/kind.o
 scale_h.o : ../Modules/cell_base.o
 scale_h.o : ../Modules/io_global.o
 scale_h.o : ../Modules/ions_base.o
@ -779,17 +775,17 @@ set_fft_dim.o : pwcom.o
 set_kplusb.o : ../Modules/kind.o
 set_kplusq.o : ../Modules/kind.o
 set_kup_and_kdw.o : ../Modules/kind.o
+setlocal.o : ../Modules/ions_base.o
+setlocal.o : ../Modules/kind.o
+setlocal.o : pwcom.o
+setqf.o : ../Modules/kind.o
 set_rhoc.o : ../Modules/atom.o
 set_rhoc.o : ../Modules/cell_base.o
 set_rhoc.o : ../Modules/io_global.o
 set_rhoc.o : ../Modules/ions_base.o
 set_rhoc.o : ../Modules/kind.o
 set_rhoc.o : pwcom.o
-set_vrs.o : ../Modules/kind.o
-setlocal.o : ../Modules/ions_base.o
-setlocal.o : ../Modules/kind.o
-setlocal.o : pwcom.o
-setqf.o : ../Modules/kind.o
+setupkpt.o : ../Modules/kind.o
 setup.o : ../Modules/atom.o
 setup.o : ../Modules/cell_base.o
 setup.o : ../Modules/constants.o
@ -802,19 +798,24 @@ setup.o : ../Modules/parameters.o
 setup.o : ../Modules/uspp.o
 setup.o : noncol.o
 setup.o : pwcom.o
-setupkpt.o : ../Modules/kind.o
-sgam_at.o : ../Modules/io_global.o
-sgam_at.o : ../Modules/kind.o
+set_vrs.o : ../Modules/kind.o
+sgama.o : ../Modules/kind.o
 sgam_at_mag.o : ../Modules/io_global.o
 sgam_at_mag.o : ../Modules/kind.o
+sgam_at.o : ../Modules/io_global.o
+sgam_at.o : ../Modules/kind.o
 sgam_ph.o : ../Modules/kind.o
-sgama.o : ../Modules/kind.o
 smallg_q.o : ../Modules/kind.o
 sph_ind.o : ../Modules/kind.o
 spinor.o : ../Modules/kind.o
+s_psi.o : becmod.o
+s_psi.o : ../Modules/ions_base.o
+s_psi.o : ../Modules/kind.o
+s_psi.o : ../Modules/uspp.o
+s_psi.o : pwcom.o
 startup.o : ../Modules/io_global.o
-startup.o : ../Modules/mp.o
 startup.o : ../Modules/mp_global.o
+startup.o : ../Modules/mp.o
 startup.o : para.o
 stop_pw.o : ../Modules/control_flags.o
 stop_pw.o : ../Modules/input_parameters.o
@ -856,13 +857,6 @@ stres_loc.o : ../Modules/kind.o
 stres_loc.o : ../Modules/uspp.o
 stres_loc.o : ../Modules/wavefunctions.o
 stres_loc.o : pwcom.o
-stres_us.o : ../Modules/constants.o
-stres_us.o : ../Modules/ions_base.o
-stres_us.o : ../Modules/kind.o
-stres_us.o : ../Modules/mp_global.o
-stres_us.o : ../Modules/uspp.o
-stres_us.o : ../Modules/wavefunctions.o
-stres_us.o : pwcom.o
 stress.o : ../Modules/cell_base.o
 stress.o : ../Modules/constants.o
 stress.o : ../Modules/control_flags.o
@ -870,13 +864,20 @@ stress.o : ../Modules/io_global.o
 stress.o : ../Modules/ions_base.o
 stress.o : ../Modules/kind.o
 stress.o : pwcom.o
+stres_us.o : ../Modules/constants.o
+stres_us.o : ../Modules/ions_base.o
+stres_us.o : ../Modules/kind.o
+stres_us.o : ../Modules/mp_global.o
+stres_us.o : ../Modules/uspp.o
+stres_us.o : ../Modules/wavefunctions.o
+stres_us.o : pwcom.o
 struct_fact.o : ../Modules/kind.o
 sum_band.o : ../Modules/cell_base.o
 sum_band.o : ../Modules/io_files.o
 sum_band.o : ../Modules/ions_base.o
 sum_band.o : ../Modules/kind.o
-sum_band.o : ../Modules/mp.o
 sum_band.o : ../Modules/mp_global.o
+sum_band.o : ../Modules/mp.o
 sum_band.o : ../Modules/uspp.o
 sum_band.o : ../Modules/wavefunctions.o
 sum_band.o : pwcom.o
@ -923,9 +924,6 @@ usnldiag.o : ../Modules/ions_base.o
 usnldiag.o : ../Modules/kind.o
 usnldiag.o : ../Modules/uspp.o
 usnldiag.o : pwcom.o
-v_of_rho.o : ../Modules/io_global.o
-v_of_rho.o : ../Modules/kind.o
-v_of_rho.o : pwcom.o
 vcsmd.o : ../Modules/cell_base.o
 vcsmd.o : ../Modules/constants.o
 vcsmd.o : ../Modules/control_flags.o
@ -947,20 +945,23 @@ vloc_of_g.o : ../Modules/kind.o
 vloc_psi.o : ../Modules/kind.o
 vloc_psi.o : ../Modules/wavefunctions.o
 vloc_psi.o : pwcom.o
+v_of_rho.o : ../Modules/io_global.o
+v_of_rho.o : ../Modules/kind.o
+v_of_rho.o : pwcom.o
 vpack.o : ../Modules/kind.o
 w0gauss.o : ../Modules/kind.o
 w1gauss.o : ../Modules/kind.o
+wfcinit.o : becmod.o
 wfcinit.o : ../Modules/cell_base.o
 wfcinit.o : ../Modules/constants.o
 wfcinit.o : ../Modules/control_flags.o
 wfcinit.o : ../Modules/io_files.o
 wfcinit.o : ../Modules/io_global.o
 wfcinit.o : ../Modules/kind.o
-wfcinit.o : ../Modules/mp.o
 wfcinit.o : ../Modules/mp_global.o
+wfcinit.o : ../Modules/mp.o
 wfcinit.o : ../Modules/uspp.o
 wfcinit.o : ../Modules/wavefunctions.o
-wfcinit.o : becmod.o
 wfcinit.o : pwcom.o
 wgauss.o : ../Modules/kind.o
 write_config_to_file.o : ../Modules/control_flags.o
--- a/PW/pwcom.f90
+++ b/PW/pwcom.f90
@ -408,9 +408,10 @@ MODULE char
  !
  ! ... The names of the system and of the symmetries
  !  
+  USE printout_base, ONLY: title    ! title of the run
+  !
  SAVE
  !
-  CHARACTER(LEN=75) ::  title       ! title of the run
  CHARACTER(LEN=20) ::  crystal     ! type of the solid
  CHARACTER(LEN=45) ::  sname(48)   ! name of the symmetries
  !
--- a/PWCOND/Makefile
+++ b/PWCOND/Makefile
@ -70,6 +70,7 @@ MODULES = \
 ../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
+../Modules/printout_base.o \
 ../Modules/pseudo_types.o \
 ../Modules/readpseudo.o \
 ../Modules/recvec.o \
--- a/PWCOND/make.depend
+++ b/PWCOND/make.depend
@ -1,87 +1,87 @@
+allocate_cond_2.o : condcom.o
+allocate_cond_2.o : ../PW/noncol.o
+allocate_cond.o : condcom.o
 allocate_cond.o : ../Modules/atom.o
 allocate_cond.o : ../Modules/ions_base.o
 allocate_cond.o : ../Modules/uspp.o
 allocate_cond.o : ../PW/noncol.o
 allocate_cond.o : ../PW/pwcom.o
-allocate_cond.o : condcom.o
-allocate_cond_2.o : ../PW/noncol.o
-allocate_cond_2.o : condcom.o
 bessj.o : ../Modules/kind.o
+compbs_2.o : ../Modules/kind.o
+compbs.o : condcom.o
 compbs.o : ../Modules/uspp.o
 compbs.o : ../PW/noncol.o
 compbs.o : ../PW/pwcom.o
-compbs.o : condcom.o
-compbs_2.o : ../Modules/kind.o
 condcom.o : ../Modules/kind.o
+do_cond.o : condcom.o
 do_cond.o : ../Modules/io_files.o
 do_cond.o : ../Modules/io_global.o
 do_cond.o : ../Modules/ions_base.o
 do_cond.o : ../Modules/mp.o
 do_cond.o : ../PW/pwcom.o
-do_cond.o : condcom.o
 eigenchnl.o : ../Modules/kind.o
 form_zk.o : ../Modules/kind.o
+four.o : condcom.o
 four.o : ../Modules/atom.o
 four.o : ../Modules/kind.o
 four.o : ../Modules/parameters.o
 four.o : ../Modules/uspp.o
 four.o : ../PW/pwcom.o
-four.o : condcom.o
 free_mem.o : condcom.o
-gep_gep.o : ../Modules/kind.o
 gep_gep.o : condcom.o
+gep_gep.o : ../Modules/kind.o
 gep_x.o : ../Modules/kind.o
 gramsh.o : ../Modules/kind.o
 hev_ab.o : ../Modules/kind.o
+init_cond.o : condcom.o
 init_cond.o : ../Modules/io_global.o
 init_cond.o : ../Modules/ions_base.o
 init_cond.o : ../Modules/uspp.o
 init_cond.o : ../PW/noncol.o
 init_cond.o : ../PW/pwcom.o
-init_cond.o : condcom.o
+init_gper.o : condcom.o
 init_gper.o : ../Modules/io_global.o
 init_gper.o : ../PW/pwcom.o
-init_gper.o : condcom.o
 integrals.o : ../Modules/kind.o
 jbloch.o : ../Modules/kind.o
 jbloch.o : ../PW/noncol.o
 kbloch.o : ../Modules/constants.o
 kbloch.o : ../Modules/kind.o
+local.o : condcom.o
 local.o : ../Modules/io_files.o
 local.o : ../Modules/io_global.o
 local.o : ../Modules/mp_global.o
 local.o : ../Modules/parallel_include.o
 local.o : ../PW/noncol.o
 local.o : ../PW/pwcom.o
-local.o : condcom.o
+poten.o : condcom.o
 poten.o : ../Modules/mp_global.o
 poten.o : ../PW/noncol.o
 poten.o : ../PW/para.o
 poten.o : ../PW/pwcom.o
-poten.o : condcom.o
+print_clock_pwcond.o : condcom.o
 print_clock_pwcond.o : ../Modules/io_global.o
 print_clock_pwcond.o : ../Modules/mp_global.o
-print_clock_pwcond.o : condcom.o
 rotate.o : ../Modules/kind.o
 rotproc.o : ../Modules/kind.o
 rotproc.o : ../Modules/mp_global.o
 rotproc.o : ../Modules/parallel_include.o
 rotproc.o : ../PW/noncol.o
+scatter_back.o : condcom.o
 scatter_back.o : ../PW/noncol.o
 scatter_back.o : ../PW/pwcom.o
-scatter_back.o : condcom.o
+scatter_forw.o : condcom.o
 scatter_forw.o : ../PW/noncol.o
 scatter_forw.o : ../PW/pwcom.o
-scatter_forw.o : condcom.o
 slabcpu.o : ../Modules/kind.o
+summary_band.o : condcom.o
 summary_band.o : ../Modules/io_files.o
 summary_band.o : ../Modules/io_global.o
 summary_band.o : ../PW/noncol.o
 summary_band.o : ../PW/pwcom.o
-summary_band.o : condcom.o
 summary_tran.o : ../Modules/kind.o
+transmit.o : condcom.o
 transmit.o : ../Modules/io_global.o
 transmit.o : ../Modules/uspp.o
 transmit.o : ../PW/noncol.o
 transmit.o : ../PW/pwcom.o
-transmit.o : condcom.o
--- a/atomic/Makefile
+++ b/atomic/Makefile
@ -86,6 +86,7 @@ MODULES = \
 ../Modules/functionals.o \
 ../Modules/io_global.o \
 ../Modules/kind.o \
+../Modules/io_files.o \
 ../PW/bachel.o \
 ../PW/functionals.o \
 ../PW/lsda_functionals.o
--- a/atomic/make.depend
+++ b/atomic/make.depend
@ -9,53 +9,53 @@ descreening.o : ld1inc.o
 el_config.o : ../Modules/kind.o
 el_config.o : parameters.o
 elsd.o : ../Modules/functionals.o
-elsdps.o : ../Modules/functionals.o
 elsdps.o : ld1inc.o
+elsdps.o : ../Modules/functionals.o
 esic.o : ld1inc.o
 find_qi.o : ld1inc.o
 gener_pseudo.o : ld1inc.o
-ld1.o : ld1inc.o
-ld1_readin.o : ../Modules/functionals.o
-ld1_readin.o : ../Modules/kind.o
-ld1_readin.o : ld1inc.o
-ld1_readin.o : parameters.o
-ld1_setup.o : ../Modules/functionals.o
-ld1_setup.o : ld1inc.o
-ld1_writeout.o : ../Modules/functionals.o
-ld1_writeout.o : ld1inc.o
 ld1inc.o : ../Modules/kind.o
 ld1inc.o : parameters.o
+ld1.o : ld1inc.o
+ld1_readin.o : ld1inc.o
+ld1_readin.o : ../Modules/functionals.o
+ld1_readin.o : ../Modules/kind.o
+ld1_readin.o : parameters.o
+ld1_setup.o : ld1inc.o
+ld1_setup.o : ../Modules/functionals.o
+ld1_writeout.o : ld1inc.o
+ld1_writeout.o : ../Modules/functionals.o
 lderiv.o : ld1inc.o
 lderivps.o : ld1inc.o
-new_potential.o : ../Modules/functionals.o
 newd_at.o : ld1inc.o
+new_potential.o : ../Modules/functionals.o
 normalize.o : ld1inc.o
 openfil.o : ld1inc.o
 pseudo_types.o : ../Modules/kind.o
 pseudovloc.o : ld1inc.o
-read_newpseudo.o : ../Modules/functionals.o
 read_newpseudo.o : ld1inc.o
+read_newpseudo.o : ../Modules/functionals.o
 read_pseudo_module.o : ../Modules/kind.o
 read_pseudo_module.o : pseudo_types.o
-read_pseudoupf.o : ../Modules/functionals.o
 read_pseudoupf.o : ld1inc.o
+read_pseudoupf.o : ../Modules/functionals.o
 read_pseudoupf.o : pseudo_types.o
 read_pseudoupf.o : read_pseudo_module.o
 run_pseudo.o : ld1inc.o
 run_test.o : ld1inc.o
-scf.o : ../Modules/constants.o
+scf.o : constants.o
 scf.o : ld1inc.o
 set_rho_core.o : ld1inc.o
-sic_correction.o : ../Modules/constants.o
-sic_correction.o : ../Modules/functionals.o
+sic_correction.o : constants.o
 sic_correction.o : ld1inc.o
+sic_correction.o : ../Modules/functionals.o
 start_potps.o : ld1inc.o
 vpack.o : ../Modules/kind.o
-write_paw_recon.o : ../Modules/functionals.o
 write_paw_recon.o : ld1inc.o
-write_results.o : ../Modules/functionals.o
+write_paw_recon.o : ../Modules/functionals.o
 write_results.o : ld1inc.o
-write_resultsps.o : ../Modules/functionals.o
+write_results.o : ../Modules/functionals.o
 write_resultsps.o : ld1inc.o
-write_upf.o : ../Modules/functionals.o
+write_resultsps.o : ../Modules/functionals.o
 write_upf.o : ld1inc.o
+write_upf.o : ../Modules/functionals.o
--- a/flib/make.depend
+++ b/flib/make.depend
@ -10,16 +10,16 @@ linpack.o : ../Modules/kind.o
 ngnr_set.o : ../Modules/fft_scalar.o
 ngnr_set.o : ../Modules/io_global.o
 ngnr_set.o : ../Modules/kind.o
-ngnr_set.o : ../Modules/mp.o
 ngnr_set.o : ../Modules/mp_global.o
+ngnr_set.o : ../Modules/mp.o
 recips.o : ../Modules/kind.o
 rranf.o : ../Modules/kind.o
 scnds.o : ../Modules/kind.o
 simpsn.o : ../Modules/kind.o
-sort.o : ../Modules/constants.o
-sort.o : ../Modules/kind.o
 sort_gvec.o : ../Modules/constants.o
 sort_gvec.o : ../Modules/kind.o
+sort.o : ../Modules/constants.o
+sort.o : ../Modules/kind.o
 sph_bes.o : ../Modules/kind.o
 volume.o : ../Modules/kind.o
 ylmr2.o : ../Modules/kind.o
--- a/pwtools/make.depend
+++ b/pwtools/make.depend
@ -1,6 +1,6 @@
 matdyn.o : ../Modules/kind.o
-matdyn.o : ../Modules/mp.o
 matdyn.o : ../Modules/mp_global.o
+matdyn.o : ../Modules/mp.o
 path_int.o : ../Modules/basic_algebra_routines.o
 path_int.o : ../Modules/constants.o
 path_int.o : ../Modules/io_files.o
@ -10,5 +10,5 @@ path_int.o : ../Modules/splinelib.o
 path_int.o : ../Modules/supercell.o
 q2r.o : ../Modules/fft_scalar.o
 q2r.o : ../Modules/kind.o
-q2r.o : ../Modules/mp.o
 q2r.o : ../Modules/mp_global.o
+q2r.o : ../Modules/mp.o