diff --git a/NN/volume.c b/NN/volume.c
index 79a606b30b2d5b4eda411fee671b7318d35a0639..3fbfa3983f4d574525711485bdeb6aafad305f34 100644
--- a/NN/volume.c
+++ b/NN/volume.c
@@ -216,6 +216,7 @@ return(v);
 
 ----------------------------------------------------------------*/
 
+/*
 volume_(a,b,m,n,mmax,nmax,result)
      int   *n, *m, *mmax, *nmax;
      float *a, *b;
@@ -223,6 +224,7 @@ volume_(a,b,m,n,mmax,nmax,result)
 {
   volume (a,b,m,n,mmax,nmax,result);
 }
+*/
 
 void volume (a,b,m,n,mmax,nmax,result)
      int   *n, *m, *mmax, *nmax;
@@ -433,6 +435,7 @@ return(v);
 
 ----------------------------------------------------------------*/
 
+/*
 volumeb_ (a,b,m,n,mmax,nmax,opt,volume,dvdb) 
      int   *n, *m, *mmax, *nmax, *opt;
      float *a, *b;
@@ -441,6 +444,7 @@ volumeb_ (a,b,m,n,mmax,nmax,opt,volume,dvdb)
 {
   volumeb(a,b,m,n,mmax,nmax,opt,volume,dvdb);
 }
+*/
 
 void volumeb (a,b,m,n,mmax,nmax,opt,volume,dvdb)
 
@@ -668,6 +672,7 @@ return(v);
    A dummy routine used to call cvolumebj from a fortran routine
 
 ----------------------------------------------------------------*/
+/*
 volumebj_ (a,b,m,n,mmax,nmax,con,volume,dvdb) 
      int   *n, *m, *mmax, *nmax,*con;
      float *a, *b;
@@ -676,6 +681,7 @@ volumebj_ (a,b,m,n,mmax,nmax,con,volume,dvdb)
 {
   volumebj(a,b,m,n,mmax,nmax,con,volume,dvdb);
 }
+*/
 
 void volumebj (a,b,m,n,mmax,nmax,con,volume,dvdb)
 
@@ -1249,7 +1255,7 @@ return(v);
    A dummy routine used to call cvolume from a fortran routine 
 
 ----------------------------------------------------------------*/
-
+/*
 volumef_ (a,b,m,n,mmax,nmax,volume) 
      int   *n, *m, *mmax, *nmax;
      float *a, *b;
@@ -1257,6 +1263,7 @@ volumef_ (a,b,m,n,mmax,nmax,volume)
 {
   volumef(a,b,m,n,mmax,nmax,volume);
 }
+*/
 
 void volumef (a,b,m,n,mmax,nmax,volume)
 
@@ -1628,7 +1635,7 @@ return(v);
    A dummy routine used to call cdvda from a fortran routine 
 
 ----------------------------------------------------------------*/
-
+/*
 dvdaf_ (a,b,m,n,mmax,nmax,idim,dvda,code) 
      int   *n, *m, *mmax, *nmax, *idim, *code;
      float *a, *b;
@@ -1636,6 +1643,7 @@ dvdaf_ (a,b,m,n,mmax,nmax,idim,dvda,code)
 {
   dvdaf(a,b,m,n,mmax,nmax,idim,dvda,code);
 }
+*/
 
 void dvdaf (a,b,m,n,mmax,nmax,idim,dvda,code)