diff --git a/base/internals/Makefile b/base/internals/Makefile
index 66276c8a1..de55e4b10 100644
--- a/base/internals/Makefile
+++ b/base/internals/Makefile
@@ -20,7 +20,7 @@ lib: mpfobjs $(FOBJS) $(FOBJS2) $(COBJS) $(MPFOBJS2) $(MPFOBJS)
 
 $(FOBJS) $(FBOJS2): $(MODDIR)/psi_mod.o
 mpfobjs: 
-	(make $(MPFOBJS) F90="$(MPF90)" FC="$(MPF90)" FCOPT="$(F90COPT)")	
+	(make $(MPFOBJS) FC="$(MPFC)")	
 
 clean:
 	/bin/rm -f $(MPFOBJS) $(FOBJS) $(COBJS) $(FOBJS2) $(MPFOBJS2) *$(.mod)
diff --git a/util/psb_d_genpde_impl.f90 b/util/psb_d_genpde_impl.f90
index baaf92de2..abe9e9cb9 100644
--- a/util/psb_d_genpde_impl.f90
+++ b/util/psb_d_genpde_impl.f90
@@ -206,9 +206,9 @@ subroutine psb_d_gen_pde3d(ictxt,idim,a,bv,xv,desc_a,afmt,&
       endif
       iz = glob_row-(ix-1)*idim*idim-(iy-1)*idim
       ! x, y, x coordinates
-      x = ix*deltah
-      y = iy*deltah
-      z = iz*deltah
+      x = (ix-1)*deltah
+      y = (iy-1)*deltah
+      z = (iz-1)*deltah
       zt(k) = f_(x,y,z)
       ! internal point: build discretization
       !   
@@ -523,8 +523,8 @@ subroutine psb_d_gen_pde2d(ictxt,idim,a,bv,xv,desc_a,afmt,&
       endif
       iy = (glob_row-(ix-1)*idim)
       ! x, y
-      x = ix*deltah
-      y = iy*deltah
+      x = (ix-1)*deltah
+      y = (iy-1)*deltah
 
       zt(k) = f_(x,y)
       ! internal point: build discretization
diff --git a/util/psb_s_genpde_impl.f90 b/util/psb_s_genpde_impl.f90
index e517c3b5b..ba2f55895 100644
--- a/util/psb_s_genpde_impl.f90
+++ b/util/psb_s_genpde_impl.f90
@@ -191,9 +191,9 @@ subroutine psb_s_gen_pde3d(ictxt,idim,a,bv,xv,desc_a,afmt,&
       endif
       iz = glob_row-(ix-1)*idim*idim-(iy-1)*idim
       ! x, y, x coordinates
-      x = ix*deltah
-      y = iy*deltah
-      z = iz*deltah
+      x = (ix-1)*deltah
+      y = (iy-1)*deltah
+      z = (iz-1)*deltah
       zt(k) = f_(x,y,z)
       ! internal point: build discretization
       !   
@@ -493,8 +493,8 @@ subroutine psb_s_gen_pde2d(ictxt,idim,a,bv,xv,desc_a,afmt,&
       endif
       iy = (glob_row-(ix-1)*idim)
       ! x, y
-      x = ix*deltah
-      y = iy*deltah
+      x = (ix-1)*deltah
+      y = (iy-1)*deltah
 
       zt(k) = f_(x,y)
       ! internal point: build discretization