add the _ versions

NN/stack.c

@@ -4,7 +4,7 @@
 	{ int key; struct node *next; };
 	static struct node *head, *z, *t;
 
-/*        stackinit_() {stackinit();} */
+        void stackinit_() {stackinit();}
         void stackinit() 
 	   {
 	     head = (struct node *) malloc(sizeof *head);
@@ -14,7 +14,7 @@
 	     z->key = 0;
 	   }
 
-/*        push_(p) int *p; {push(p);} */
+        void push_(p) int *p; {push(p);}
 	void push(p)
            int *p;
 	   {
@@ -25,7 +25,7 @@
 	     head->next =t;	
 	   }
 
-/*        pop_(x) int *x; {pop(x);} */
+        void pop_(x) int *x; {pop(x);}
 	void pop(x)
            int *x;
 	   {
@@ -34,7 +34,7 @@
 	     free(t);
 	   }
 
-/*        stackempty_(i) int *i; {stackempty(i);} */
+        void stackempty_(i) int *i; {stackempty(i);}
 	void stackempty(i)
           int *i;
 	  { 
@@ -42,7 +42,7 @@
             if(head->next == z) *i = 1;
           }
 
-/*        stackflush_() {stackflush();} */
+        void stackflush_() {stackflush();}
         void stackflush()
            {
              free(head);

NN/stackpair.c

@@ -5,7 +5,7 @@
 	static struct node *heada, *za, *ta;
 	static struct node *headb, *zb, *tb;
 
-/*        stackpairinit_() {stackpairinit();} */
+        void stackpairinit_() {stackpairinit();}
 	void stackpairinit () 
 	   {
 	     heada = (struct node *) malloc(sizeof *heada);
@@ -21,7 +21,7 @@
 	     zb->key = 0;
 	   }
 
-/*        stackpairflush_() {stackpairflush();} */
+        void stackpairflush_() {stackpairflush();}
 	void stackpairflush () 
 	   {
              free(heada);
@@ -30,7 +30,7 @@
              free(zb);
 	   }
 
-/*        pushpair_(pa, pb)  int *pa; int *pb; {pushpair(pa, pb);} */
+        void pushpair_(pa, pb)  int *pa; int *pb; {pushpair(pa, pb);}
 	void pushpair (pa, pb)
              int *pa;
              int *pb;
@@ -48,7 +48,7 @@
 	     headb->next =tb;	
 	   }
 
-/*        poppair_(xa,xb)  int *xa; int *xb; {poppair(xa,xb);} */
+        void poppair_(xa,xb)  int *xa; int *xb; {poppair(xa,xb);}
 	void poppair (xa,xb)
              int *xa;
              int *xb;
@@ -61,7 +61,7 @@
 	     free(tb);
 	   }
 
-/*        stackpairempty_(i) int *i; {stackpairempty(i);} */
+        void stackpairempty_(i) int *i; {stackpairempty(i);}
 	void stackpairempty (i)
             int *i;
 	  { 

NN/volume.c

@@ -435,8 +435,8 @@ return(v);
 
 ----------------------------------------------------------------*/
 
-/*
-volumeb_ (a,b,m,n,mmax,nmax,opt,volume,dvdb) 
+
+void volumeb_ (a,b,m,n,mmax,nmax,opt,volume,dvdb) 
      int   *n, *m, *mmax, *nmax, *opt;
      float *a, *b;
      float *volume; 
@@ -444,7 +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)
 
@@ -672,8 +672,8 @@ return(v);
    A dummy routine used to call cvolumebj from a fortran routine
 
 ----------------------------------------------------------------*/
-/*
-volumebj_ (a,b,m,n,mmax,nmax,con,volume,dvdb) 
+
+void volumebj_ (a,b,m,n,mmax,nmax,con,volume,dvdb) 
      int   *n, *m, *mmax, *nmax,*con;
      float *a, *b;
      float *volume; 
@@ -681,7 +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)
 
@@ -1051,16 +1051,16 @@ return(v);
    A dummy routine used to call cdvda from a fortran routine 
 
 ----------------------------------------------------------------*/
-/*  Not used...won't compile...function name and variable name identical */
-/*
-dvda_ (a,b,m,n,mmax,nmax,idim,result,code) 
+/*  Not used...won't compile...function name and variable name (were) identical */
+
+void dvda_ (a,b,m,n,mmax,nmax,idim,result,code) 
      int   *n, *m, *mmax, *nmax, *idim, *code;
      float *a, *b;
      float *result; 
 {
   dvda(a,b,m,n,mmax,nmax,idim,result,code);
 }
-*/
+
 void dvda (a,b,m,n,mmax,nmax,idim,result,code)
 
 int   *n, *m, *mmax, *nmax, *idim, *code;
@@ -1255,15 +1255,15 @@ return(v);
    A dummy routine used to call cvolume from a fortran routine 
 
 ----------------------------------------------------------------*/
-/*
-volumef_ (a,b,m,n,mmax,nmax,volume) 
+
+void volumef_ (a,b,m,n,mmax,nmax,volume) 
      int   *n, *m, *mmax, *nmax;
      float *a, *b;
      float *volume; 
 {
   volumef(a,b,m,n,mmax,nmax,volume);
 }
-*/
+
 
 void volumef (a,b,m,n,mmax,nmax,volume)
 
@@ -1635,15 +1635,15 @@ return(v);
    A dummy routine used to call cdvda from a fortran routine 
 
 ----------------------------------------------------------------*/
-/*
-dvdaf_ (a,b,m,n,mmax,nmax,idim,dvda,code) 
+
+void dvdaf_ (a,b,m,n,mmax,nmax,idim,dvda,code) 
      int   *n, *m, *mmax, *nmax, *idim, *code;
      float *a, *b;
      float *dvda; 
 {
   dvdaf(a,b,m,n,mmax,nmax,idim,dvda,code);
 }
-*/
+
 
 void dvdaf (a,b,m,n,mmax,nmax,idim,dvda,code)