#include <stdio.h>
#include <float.h>
/* #include <math.h> */
#include <ctype.h>
#include <stdlib.h>


main() 
{

/*  Reads data from the Scope2.new file.  Converts to MKS units  */

int   day, line_num, oph, test, tilt;   
float day_decimal, day_decimal_alt,  delta_q, delta_temp, freq_ave, freq_max, 
      freq_min, hour, hrms_ave, hrms_min, hrms_max, hsig, inv_ht, 
      lhf_ave, lhf_bulk, lhf_cov, lhf_iner, minutes, neg_z_over_L,  
      phase_spd, qa, qs, shf_ave, shf_bulk, shf_cov, shf_iner, ta, tau_ave, 
      tau_bulk, tau_cov, tau_iner, tw, u, ustar_cf, wave_age_ustar, wave_age_u, 
      wind_dir, wstar;  
char  month[4];  
FILE  *fp_data; 

/*  open data file  */
  if( ( fp_data = fopen( "Scope2.new", "r" ) ) == NULL )
     {
     printf( "Can not open the file Scope2.new\n" );
     exit( 1 );
     }

/*  skip first line  */
  fscanf(fp_data, "%*[^\n]\n" ); 

  line_num = 2; 
/*  read in data  */
   while ( ( test = fscanf( fp_data, "%i %s %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %i %i %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f\n",
   &day, month, &hour, &u, &tw, &ta, &delta_temp, &day_decimal, &qs, &qa,
   &delta_q, &shf_cov, &shf_iner, &shf_bulk, &shf_ave, &lhf_cov, &lhf_iner, 
   &lhf_bulk, &lhf_ave, &tau_cov, &tau_iner, &tau_bulk, &tau_ave, &ustar_cf,
   &oph, &tilt, &inv_ht, &neg_z_over_L, &wstar, &day_decimal_alt, &hrms_min,
   &hrms_max, &freq_min, &freq_max, &hrms_ave, &freq_ave, &phase_spd,
   &wave_age_ustar, &wave_age_u, &hsig, &wind_dir ) ) == 41 ) {         

    line_num++;   
 
/*  convert to MKS units  */
    hrms_min = 0.01 * hrms_min;
    hrms_max = 0.01 * hrms_max;   
    hrms_ave = 0.01 * hrms_ave; 
    qa = 0.001 * qa;
    qs = 0.001 * qs;

  }  
  if ( test != -1 ) printf( "Read error on line %i of Scope2.new; test = %i\n",
   line_num, test );   

  printf( "%i lines of data\n", line_num );


}