#include <stdio.h> /*c-style io */
#include <math.h>
#define infile "/home/meteor/x/class/jhorel/3410/hat.990909"

main(int argc, char *argv[])
{
  FILE *input;
  char line[80];
  char stn[3],date[12];
  float tmpf,tsoi;
  float tm,ts,tlm,tls,cor,ct,slope,inter,px,py;


/**************  OPEN FILE AND GET DATA ************************/

   input = fopen(infile,"r"); 
   printf(" Stn Time Temp TLake\n"); 
   tm=0;
   ts=0;
   tlm=0;
   tls=0;
   cor=0;
   ct=0;


   while ((fgets(line,sizeof(line),input)) != NULL )
{
   sscanf(line,"     %s    %s %f %f " ,stn,date,&tmpf,&tsoi);
 
  if ( strncmp("HAT",stn,3) == 0 )
  {
        ct=ct+1;
        tm=tm+tmpf;
        ts=ts+tmpf*tmpf;
        tlm=tlm+tsoi;
        tls=tls+tsoi*tsoi;
        cor=cor+tmpf*tsoi;

     printf("%f %f \n", tmpf, tsoi);
  }
}
     tm=tm/ct;
     ts=ts/ct-tm*tm;
     tlm=tlm/ct;
     tls=tls/ct-tlm*tlm;
     cor=cor/ct-tm*tlm;
     cor=cor/sqrt(ts*tls);
     slope=cor*tls/ts;
     inter=tlm-slope*tm;
     px=74;
     py=inter+slope*px;
     printf("%f %f %f %f %f %f %f %f %f \n", tm,ts,tlm,tls,cor,slope,inter,px,py);
     return(0);
}