/***************************************************************************
* nearname.c - search objects by its name with given radius,  ^D to exit
*
* History 
*
* 24-Aug-2004: include <stlib.h> to define atof() instead of local declaration
* 23-Jul-2001:
*   -Added double atof(char []); to insure main() knows that atof() 
*    returns a double. Joe Mazzarella
*   -Changed gets() to fgets() for compiler warnings about security. J.M.
* 
* 1996:
*   Original. Xiuqin Wu
*
***************************************************************************/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "ned_client.h"

extern int   ned_errno;

void
main(argc, argv)
int  argc;
char *argv[];
{
   
   int    st;
   int    no_obj;
   int    i, j;
   char   objname[101];
   char   sradius[101];
   double radius;
   CrossID *cp, *tmpcp;
   ObjInfo *op, *tmpop;
   MoreData     *mdatap;

   st = ned_connect();
   if (st < 0) {
      fprintf(stderr, "connection failed \n");
      exit(1);
      }
   fprintf(stdout, "input the objname:");
   fgets(objname, 101, stdin);
   fprintf(stdout, "input the radius(arcmin):");
   while(fgets(sradius, 101, stdin) != (char*)NULL) {

      radius = atof(sradius);
      fprintf(stdout, "DEBUG: You input sradius=%s, which atof translates into radius=%f", sradius, radius);
      st = ned_obj_nearname(objname, radius, &no_obj, &op, &cp);
      if (st < 0) {

	 /* for simple error message */
	 fprintf(stderr, "%s\n", ned_get_errmsg());

	 switch (ned_errno) {
	    case NE_NAME:
	       fprintf(stderr, 
	       "name %s can't be recognized by NED name interpreter\n",
		objname);
	       break;
	    case NE_RADIUS:
	       fprintf(stderr, "radius is out of range\n");
	       break;
	    case NE_AMBN:
	       fprintf(stderr, "%d ambiguous name: \n", no_obj);
	       for (i=0, tmpcp = cp; i<no_obj; i++, tmpcp++)
		  fprintf(stderr, "%s \n", tmpcp->objname);
	       break;
	    case NE_NOBJ:
	       fprintf(stderr, "object %s is not in NED database\n",
		  cp->objname);
	       break;
	    case NE_NOSPACE:
	       fprintf(stderr, "memory allocation error happened \n");
	       break;
	    case NE_QUERY:
	       fprintf(stderr, "Can't send query to the NED server\n");
	       break;
	    case NE_BROKENC:
	       fprintf(stderr, "The connection to server is broken\n");
	       break;
	    }
	 } /* -1 return code */
      else {
	 fprintf(stdout, "%d object(s) found in NED: \n", no_obj);
	 for (i=0, tmpop = op; i<no_obj; i++, tmpop++) {
	    fprintf(stdout, "\n\n%d crossid for object No. %d: \n\n", 
	       tmpop->no_crossid, i+1);
	    for (j=0, tmpcp = tmpop->cp; j<tmpop->no_crossid; j++, tmpcp++)
	       fprintf(stdout, "%s, %s \n", tmpcp->objname, tmpcp->objtype);
	    fprintf(stdout, "distance :    %f\n", tmpop->dist);
	    fprintf(stdout, "no_ref   :    %d\n", tmpop->no_ref);
	    fprintf(stdout, "no_note  :    %d\n", tmpop->no_note);
	    fprintf(stdout, "no_photom:    %d\n", tmpop->no_photom);
	    fprintf(stdout, "obj type :    %s\n", tmpop->objtype);
	    fprintf(stdout, "ra       :    %f\n", tmpop->ra);
	    fprintf(stdout, "dec      :    %f\n", tmpop->dec);
	    fprintf(stdout, "unc_maj  :    %f\n", tmpop->unc_maj);
	    fprintf(stdout, "unc_min  :    %f\n", tmpop->unc_min);
	    fprintf(stdout, "unc_ang  :    %f\n", tmpop->unc_ang);
	    fprintf(stdout, "refcode  :    %s\n", tmpop->refcode);
	    mdatap = tmpop->mdp;
	    while (mdatap) {
	       fprintf(stdout, "%-10s: %s\n", mdatap->data_typec,
				 mdatap->data);
	       mdatap = mdatap->next;
	       }
	    }
	 }

      if (cp)
	 ned_free_cp(cp);
      if (op)
	 ned_free_objp(op, no_obj);
      fprintf(stdout, "input objname:");
      fgets(objname, 101, stdin);
      fprintf(stdout, "input the radius(arcmin):");
      }
   ned_disconnect();
}
