#include <stdio.h> /*c-style io */
#include <math.h>

main(int argc, char *argv[])
{
 FILE *input;
 float p[53],t[53],z[53];
 char line[80];
 int i,nlevel;
 float rgas,grav,tzero;
 float selv,tavg;
 float pei,pei1,e;
                                   
 rgas = 287.;
 grav=9.80665;
 tzero=273.;
 nlevel=52;
 sscanf(argv[1],"%f",&e);
 i=0;
 input = fopen("/home/meteor/x/class/jhorel/3410/PT.341","r");

   while ((fgets(line,sizeof(line),input)) != NULL )
 {
   sscanf(line,"%f %f",&p[i],&t[i]);
   p[i]=p[i]*100.;
   t[i]=t[i]+tzero;
   i++;
 }

  selv = 1288.;
  z[0] = selv;                                                               

  for  (i=1; i <= nlevel; i++)
 {
   pei = p[i] + e * ( p[0] - p[i] ) ;
   pei1 = p[i-1] + e * ( p[0] - p[i-1] ) ;
   tavg = 0.5 * ( t[i] + t[i-1] );                                            
   z[i] = z[i-1] - ( rgas * tavg / grav ) * log ( pei / pei1 );
 }
  for ( i=0; i <= nlevel; i++ )
 {
   printf("%f %f %f\n",p[i]/100.,t[i],z[i]);
 }
}