diff --git a/test/fileread/getp.f90 b/test/fileread/getp.f90
index db8207b6..2a32f7de 100644
--- a/test/fileread/getp.f90
+++ b/test/fileread/getp.f90
@@ -67,7 +67,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=psb_inp_unit
+        write(psb_err_unit,*) 'Usage: psb_d_fileread ctrl-file '
+        call psb_abort(ctxt)
+        stop        
       end if
       ! Read Input Parameters
       read(inp_unit,*) ip
@@ -187,7 +189,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=inp_unit
+        write(psb_err_unit,*) 'Usage: psb_s_fileread ctrl-file '
+        call psb_abort(ctxt)
+        stop        
       end if
       ! Read Input Parameters
       read(inp_unit,*) ip
diff --git a/test/pdegen/psb_d_pde2d.F90 b/test/pdegen/psb_d_pde2d.F90
index e04ec678..780a183d 100644
--- a/test/pdegen/psb_d_pde2d.F90
+++ b/test/pdegen/psb_d_pde2d.F90
@@ -877,7 +877,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=psb_inp_unit
+        call pr_usage(psb_err_unit)
+        call psb_abort(ctxt)
+        stop        
       end if
       read(inp_unit,*) ip
       if (ip >= 3) then
@@ -1012,12 +1014,10 @@ contains
   !
   subroutine pr_usage(iout)
     integer(psb_ipk_) :: iout
-    write(iout,*)'incorrect parameter(s) found'
-    write(iout,*)' usage:  pde2d90 methd prec dim &
-         &[ipart istop itmax itrace]'
-    write(iout,*)' where:'
-    write(iout,*)'     methd:    cgstab cgs rgmres bicgstabl'
-    write(iout,*)'     prec :    bjac diag none'
+    write(iout,*)' usage:  psb_d_pde2d cntrl-file '
+    write(iout,*)' where ctrl-file contains the following:'
+    write(iout,*)'     methd     BICGSTAB CGS  BICG BICGSTABL RGMRES FCG CGR RICHARDSON'
+    write(iout,*)'     prec      bjac diag none'
     write(iout,*)'     dim       number of points along each axis'
     write(iout,*)'               the size of the resulting linear '
     write(iout,*)'               system is dim**2'
diff --git a/test/pdegen/psb_d_pde3d.F90 b/test/pdegen/psb_d_pde3d.F90
index cac1c413..567ec029 100644
--- a/test/pdegen/psb_d_pde3d.F90
+++ b/test/pdegen/psb_d_pde3d.F90
@@ -933,7 +933,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=psb_inp_unit
+        call pr_usage(psb_err_unit)
+        call psb_abort(ctxt)
+        stop        
       end if
       read(inp_unit,*) ip
       if (ip >= 3) then
@@ -1070,12 +1072,10 @@ contains
   !
   subroutine pr_usage(iout)
     integer(psb_ipk_) :: iout
-    write(iout,*)'incorrect parameter(s) found'
-    write(iout,*)' usage:  pde3d90 methd prec dim &
-         &[istop itmax itrace]'
-    write(iout,*)' where:'
-    write(iout,*)'     methd:    cgstab cgs rgmres bicgstabl'
-    write(iout,*)'     prec :    bjac diag none'
+    write(iout,*)' usage:  psb_d_pde3d cntrl-file '
+    write(iout,*)' where ctrl-file contains the following:'
+    write(iout,*)'     methd     BICGSTAB CGS  BICG BICGSTABL RGMRES FCG CGR RICHARDSON'
+    write(iout,*)'     prec      bjac diag none'
     write(iout,*)'     dim       number of points along each axis'
     write(iout,*)'               the size of the resulting linear '
     write(iout,*)'               system is dim**3'
diff --git a/test/pdegen/psb_s_pde2d.F90 b/test/pdegen/psb_s_pde2d.F90
index 121980bf..179994e8 100644
--- a/test/pdegen/psb_s_pde2d.F90
+++ b/test/pdegen/psb_s_pde2d.F90
@@ -877,7 +877,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=psb_inp_unit
+        call pr_usage(psb_err_unit)
+        call psb_abort(ctxt)
+        stop        
       end if
       read(inp_unit,*) ip
       if (ip >= 3) then
@@ -1012,12 +1014,10 @@ contains
   !
   subroutine pr_usage(iout)
     integer(psb_ipk_) :: iout
-    write(iout,*)'incorrect parameter(s) found'
-    write(iout,*)' usage:  pde2d90 methd prec dim &
-         &[ipart istop itmax itrace]'
-    write(iout,*)' where:'
-    write(iout,*)'     methd:    cgstab cgs rgmres bicgstabl'
-    write(iout,*)'     prec :    bjac diag none'
+    write(iout,*)' usage:  psb_s_pde2d cntrl-file '
+    write(iout,*)' where ctrl-file contains the following:'
+    write(iout,*)'     methd     BICGSTAB CGS  BICG BICGSTABL RGMRES FCG CGR RICHARDSON'
+    write(iout,*)'     prec      bjac diag none'
     write(iout,*)'     dim       number of points along each axis'
     write(iout,*)'               the size of the resulting linear '
     write(iout,*)'               system is dim**2'
diff --git a/test/pdegen/psb_s_pde3d.F90 b/test/pdegen/psb_s_pde3d.F90
index 49b7e979..00edf046 100644
--- a/test/pdegen/psb_s_pde3d.F90
+++ b/test/pdegen/psb_s_pde3d.F90
@@ -933,7 +933,9 @@ contains
           write(psb_err_unit,*) 'Opened file ',trim(filename),' for input'
         end if
       else
-        inp_unit=psb_inp_unit
+        call pr_usage(psb_err_unit)
+        call psb_abort(ctxt)
+        stop        
       end if
       read(inp_unit,*) ip
       if (ip >= 3) then
@@ -1070,12 +1072,10 @@ contains
   !
   subroutine pr_usage(iout)
     integer(psb_ipk_) :: iout
-    write(iout,*)'incorrect parameter(s) found'
-    write(iout,*)' usage:  pde3d90 methd prec dim &
-         &[istop itmax itrace]'
-    write(iout,*)' where:'
-    write(iout,*)'     methd:    cgstab cgs rgmres bicgstabl'
-    write(iout,*)'     prec :    bjac diag none'
+    write(iout,*)' usage:  psb_s_pde3d cntrl-file '
+    write(iout,*)' where ctrl-file contains the following:'
+    write(iout,*)'     methd     BICGSTAB CGS  BICG BICGSTABL RGMRES FCG CGR RICHARDSON'
+    write(iout,*)'     prec      bjac diag none'
     write(iout,*)'     dim       number of points along each axis'
     write(iout,*)'               the size of the resulting linear '
     write(iout,*)'               system is dim**3'