diff --git a/PP/Makefile b/PP/Makefile
index 7161e8b7d..df01c77d6 100644
--- a/PP/Makefile
+++ b/PP/Makefile
@@ -46,6 +46,7 @@ MODULES = \
 ../Modules/mp_global.o \
 ../Modules/mp_wave.o \
 ../Modules/mp.o \
+../Modules/parallel_include.o \
 ../Modules/parameters.o \
 ../Modules/parser.o \
 ../Modules/pseudo_types.o \
diff --git a/PP/bands.f90 b/PP/bands.f90
index 7154cf445..c7e102604 100644
--- a/PP/bands.f90
+++ b/PP/bands.f90
@@ -23,7 +23,7 @@ program bands
   namelist / inputpp / outdir, prefix, filband
                                                                                 
   CHARACTER (LEN=80)  :: input_file
-  INTEGER             :: nargs, iiarg, ierr
+  INTEGER             :: nargs, iiarg, ierr, ilen
   INTEGER, EXTERNAL   :: iargc
 
   !
diff --git a/PP/chdens.f90 b/PP/chdens.f90
index 17ddc116b..aa0def5d1 100644
--- a/PP/chdens.f90
+++ b/PP/chdens.f90
@@ -58,7 +58,7 @@ program chdens
        makov, plot_out, output_format, fileout, epsilon, filepol
 
   CHARACTER (LEN=80)  :: input_file
-  INTEGER             :: nargs, iiarg, ierr
+  INTEGER             :: nargs, iiarg, ierr, ilen
   INTEGER, EXTERNAL   :: iargc
 
   !
@@ -1291,7 +1291,7 @@ subroutine makov_payne (ibrav, alat, at, nat, tau, ityp, zv, ntyp, x0, &
   !  differs from the "traditional" one found in the literature. See
   !  Lento, Mozos, Nieminen, J. Phys.: Condens. Matter 14 (2002), 2637-2645
 
-  real(kind=DP), parameter:: Madelung(3) = [2.8373, 2.8883, 2.885]
+  real(kind=DP), parameter:: Madelung(3) = (/ 2.8373, 2.8883, 2.885/)
 
   if (ibrav < 1 .or. ibrav > 3) then
      WRITE(stdout,'("Makov-Payne correction only for cubic lattices")')
@@ -1315,6 +1315,7 @@ subroutine makov_payne (ibrav, alat, at, nat, tau, ityp, zv, ntyp, x0, &
   dr = omega * deltax * deltay * deltaz
 
   rhomin = MAX ( MINVAL (rhor), 1.d-10 )
+  rhotot = 0.d0
   !
   !   compute electronic total charge, dipole and quadrupole moments
   !
diff --git a/PP/dos.f90 b/PP/dos.f90
index c71ec08bd..d513dd675 100644
--- a/PP/dos.f90
+++ b/PP/dos.f90
@@ -52,7 +52,7 @@ program dos
   logical :: minus_q
                                                                                 
   CHARACTER (LEN=80)  :: input_file
-  INTEGER             :: nargs, iiarg, ierr
+  INTEGER             :: nargs, iiarg, ierr, ilen
   INTEGER, EXTERNAL   :: iargc
                                                                                 
   !
diff --git a/PP/postproc.f90 b/PP/postproc.f90
index 2ca57a4e4..44b1bf6a6 100644
--- a/PP/postproc.f90
+++ b/PP/postproc.f90
@@ -44,7 +44,7 @@ program postproc
        filplot, lsign
 
   CHARACTER (LEN=80)  :: input_file
-  INTEGER             :: nargs, iiarg, ierr
+  INTEGER             :: nargs, iiarg, ierr, ilen
   INTEGER, EXTERNAL   :: iargc
 
   !
diff --git a/PP/projwfc.f90 b/PP/projwfc.f90
index c9bd58ddc..89683399e 100644
--- a/PP/projwfc.f90
+++ b/PP/projwfc.f90
@@ -42,7 +42,7 @@ program projwfc
              Emin, Emax, DeltaE, smoothing 
                                                                                 
   CHARACTER (LEN=80)  :: input_file
-  INTEGER             :: nargs, iiarg, ierr
+  INTEGER             :: nargs, iiarg, ierr, ilen
   INTEGER, EXTERNAL   :: iargc