/* param.c   Calcul des divers parametres et indices de la photometrie
             de Geneve, avec une sortie sur un fichier utilisable par
             mongo
 */
#include <stdio.h>

main(argc,argv)
  int argc;
  char *argv[];
 {
  short p1,p2,no;
  char buf[BUFSIZ];
  float mv,c[7],x,y;
  extern float parmf();
  
  p1=atoi(argv[1]);
  p2=atoi(argv[2]);
  c[0]=0.0;

  while (fgets(buf, sizeof buf, stdin) != NULL)
   {
    sscanf(buf,"%hd%*f%f%*f%f%f%f%f%f%f",&no,&mv,
           &c[1],&c[2],&c[3],&c[4],&c[5],&c[6]);

   if(p1 == 0 )
     {
      x=parmf(c,p2);
      y=mv;
      switch (p2)
        {
      case 31:
      if( x < 3.)
         printf(" %04d %7.3f %+6.3f\n",no,y,x);
      break;
      default:
         printf(" %04d %7.3f %+6.3f\n",no,y,x);
        }
     }
   else
     {
      x=parmf(c,p2);
      y=parmf(c,p1);
      if(x != 9. && y != 9.)
         printf(" %04d %+6.3f %+6.3f\n",no,y,x);
     }
    }
 }