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

int main(int argc,char **argv) {
  int i,j,k;
  int hr,mi,se;
  int ms;
  static float pwr1[900];  // RH
  static float pwr2[900];  // LH
  static float dum[900];
  FILE *fp;
  static char buf[256];
  int c0,c1,c2,c3;
  int year,mon,day;
  float f;
  int i1;
  int i2;

  year=(argv[1][0]-'0')*1000+(argv[1][1]-'0')*100+
    (argv[1][2]-'0')*10+(argv[1][3]-'0');
  mon=(argv[1][4]-'0')*10+(argv[1][5]-'0');
  day=(argv[1][6]-'0')*10+(argv[1][7]-'0');

  hr=(argv[2][0]-'0')*10+(argv[2][1]-'0');
  mi=(argv[2][2]-'0')*10+(argv[2][3]-'0');
  se=(argv[2][4]-'0')*10+(argv[2][5]-'0');
  i1=(hr*60+mi)/10;
  i2=((mi%10)*60+se)*2;

  for (i=0;i<900;i++) dum[i]=-999;

  sprintf(buf,"/home/jupiter/wdp5j/10min/%s/%s%02d%02d00.10min",
	    argv[1],argv[1],i1/6,(i1%6)*10);
  fp=fopen(buf,"r");
  for (j=0;j<1200;j++) {
    //hr=i/6;
    //mi=(i%6)*10+j/120;
    //se=(j%120)/2;
    //ms=(j%2)*50;
    for (k=0;k<900;k++) {
      c0=fgetc(fp);c1=fgetc(fp);
      c2=fgetc(fp);c3=fgetc(fp);
      pwr1[k]=(c0+256*c1)/100.0;  // Ch1 -> dBJy
      pwr2[k]=(c2+256*c3)/100.0;  // Ch2 -> dBJy
    }
    if (j==i2) break;
  }
  fclose(fp);
  for (k=0;k<900;k++) {
    f=(452+k)/4096.0*125;  // [MHz]
    printf("%.3lf %.2lf %.2lf\n",f,pwr1[k],pwr2[k]);
  }
  printf("\n");
  exit(0);
}