/* cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

  Name : readump_windii_l3lp_data.c

  Usage : This program reads UARS MLS Level 3LP data files distributed from
          the GSFC DAAC.  To run this program, type the program name
          followed by the metadata file name (*META extension).
          Example: 
          readump_windii_l3lp_data WINDII_L3LP_SPARAM_D0105.V0009_C01_META

  Description : This program first opens the metadata file (*.*META) to
                get the data file name (*.*PROD), and its physical record
                size.  The data file is then opened, and its contents are
                read and dumped to the screen.

  Author : KE - JUN SUN   Hughes_STX

  Date : February 2, 1994

  Limitations : This program and the data it reads are designed for machines
              running a UNIX operating system, and having 32-bit architecture.
              IEEE floating point and/or integer values will not translate
              correctly on other machines!

  Notes : Since this program uses the metadata file (*.*META) to get the
          data file name (*.*PROD) and record size, both the data file and
          its corresponding metadata file must be in the same directory.

  Algorithm information : No calculation is performed in this program.

  Version : 0
 
  To compile this progam : 
         cc readump_windii_l3lp_data.c -o readump_windii_l3lp_data
         executable code will be named readump_windii_l3lp_data


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc */
#include <stdio.h>
#include <string.h>

main(argc,argv)
int argc;
char *argv[];
{
 char line[1000];   
 char temp[30],data_file[180],str[80];
 int index,i,n,ret,rec,select,rec_size,*int4;
 int Total,Actual,T1,T2,Numb,Par_numb,yr,dy,ms;
 FILE *meta_ifp,*data_ifp;
 float *lat,*lon,*real4;
 float Lat,Lon,LST,SZA;

/* ******* Open meta data file first ******** */

  meta_ifp = fopen(argv[1],"r");
  if (meta_ifp == NULL)
   {
    printf(" Can not open meta data file \n");
    exit();
   }
   for(i=0;i<10;i++)
    fgets(line,80,meta_ifp);

/* get data file name from line 10 or from meta data file name  */

/*   strncpy(str,&line[14],50);
   sscanf(str,"%s",data_file);    */

   strncpy(data_file,argv[1],(strlen(argv[1])-5));
   strcat(data_file,"_PROD");

   fgets(line,80,meta_ifp);
   fgets(line,80,meta_ifp);

/* get record size from meta file */

   strncpy(temp,&line[16],5); 
   rec_size = atoi(temp);

   fclose(meta_ifp);

/* **** Now open data file **** */

   data_ifp=fopen(data_file,"r");
   if (data_ifp == NULL)
    {
     printf(" Can not open data file \n");
     exit(0);
    }

   printf(" ******************************************\n");
   printf(" Reading %s\n",data_file);
   printf(" ******************************************\n");
   memset(temp,' ',25);
   fread(line,1,rec_size,data_ifp);
   printf("-----------------------------\n");
   printf("   Descriptor       \n"); 
   printf("-----------------------------\n");
   printf(" SFDU T[Z] & L[Z] \n");
   printf("-----------------------------\n");
   strncpy(temp,&line[0],20);
   printf(" Record Key : %s\n",temp);
   memset(temp,' ',20);
   strncpy(temp,&line[20],4);
   printf(" Control ID : %s\n",temp);
   memset(temp,' ',8);   
   strncpy(temp,&line[24],1);
   printf(" Version ID : %s\n",temp);
   strncpy(temp,&line[25],1);
   printf(" Class ID   : %s\n",temp);
   strncpy(temp,&line[28],4);
   printf(" Data descriptive record ID : %s\n",temp);
   strncpy(temp,&line[32],8);
   printf(" Length : %s\n",temp);  
   printf("-----------------------------\n");
   printf(" SFDU T[I] & L[I] \n");
   printf("-----------------------------\n");
   memset(temp,' ',8);
   strncpy(temp,&line[40],4);
   printf(" Control ID : %s\n",temp);
   memset(temp,' ',8);
   strncpy(temp,&line[44],1);
   printf(" Version ID : %s\n",temp);
   strncpy(temp,&line[45],1);
   printf(" Class ID   : %s\n",temp);
   strncpy(temp,&line[48],4);
   printf(" Data descriptive record ID : %s\n",temp);
   strncpy(temp,&line[52],8);
   printf(" Length : %s\n",temp);
   memset(temp,' ',25);
   fread(line,1,rec_size,data_ifp); 
   printf("-----------------------------\n");
   printf("   Label record       \n");
   printf("-----------------------------\n");
   strncpy(temp,&line[0],20);
   printf(" Record Key    : %s\n",temp);
   memset(temp,' ',20);
   strncpy(temp,&line[20],4);
   printf(" Satellite ID  : %s\n",temp);
   memset(temp,' ',4);
   strncpy(temp,&line[24],2); 
   printf(" Record type   : %s\n",temp);
   strncpy(temp,&line[26],12);
   printf(" Instrument ID : %s\n",temp);
   strncpy(temp,&line[38],12);
   printf(" Subtype or Species : %s\n",temp);
   memset(temp,' ',12);
   strncpy(temp,&line[50],4);
   printf(" Format version number : %s\n",temp);
   strncpy(temp,&line[54],8);
   printf(" Physical record count : %s\n",temp);
   strncpy(temp,&line[66],8);
   printf(" Number of physical records in file : %s\n",temp);
   strncpy(temp,&line[74],23);
   printf(" File Creation Time in ASCII format : %s\n",temp);

   memset(temp,' ',23);
   printf(" Year day Ms for earliest data record : ");
   strncpy(temp,&line[97],3);
   sscanf(temp,"%d",&yr);
   strncpy(temp,&line[100],3);
   sscanf(temp,"%d",&dy);
   strncpy(temp,&line[103],8);
   sscanf(temp,"%d",&ms);
   printf("%3d %3d %8d\n",yr,dy,ms);
   memset(temp,' ',8);
   printf(" Year day Ms for last data record     : ");
   strncpy(temp,&line[111],3);
   sscanf(temp,"%d",&yr);
   strncpy(temp,&line[114],3);
   sscanf(temp,"%d",&dy);
   strncpy(temp,&line[117],8);
   sscanf(temp,"%d",&ms);
   printf("%3d %3d %8d\n",yr,dy,ms);
   memset(temp,' ',14);

   strncpy(temp,&line[125],3);
   printf(" Data level : %s\n",temp);
   strncpy(temp,&line[128],4);
   printf(" UARS day number : %s\n",temp);
   strncpy(temp,&line[132],4);
   printf(" Max. number of 32-bit words per record : %s\n",temp);
   strncpy(temp,&line[140],5);
   printf(" Record length in bytes : %s\n",temp);
   strncpy(temp,&line[145],6);
   printf(" Min/Max Latitude : %s\n",temp);
   strncpy(temp,&line[151],9);
   printf(" CCB version number : %s\n",temp);
   memset(temp,' ',12);
   strncpy(temp,&line[160],5);
   printf(" File cycle number : %s\n",temp);
   memset(temp,' ',12);
   strncpy(temp,&line[165],1);
   printf(" Virtual file flag : %s\n",temp);
   strncpy(temp,&line[166],4);
   printf(" Total number of time/version in file : %s\n",temp);
   strncpy(temp,&line[170],4);
   printf(" Number of time/version in record : %s\n",temp);
   rec = 2; 
   do 
   {
   memset(temp,' ',25);
   ret = fread(line,1,rec_size,data_ifp); 
   if (ret > 0)  
    {
      printf("-----------------------------\n");
      printf("  Record :  %d   data record\n",rec);
      printf("-----------------------------\n");
      strncpy(temp,&line[0],20);
      printf(" Record Key    : %s \n",temp);
      memset(temp,' ',20);
      strncpy(temp,&line[20],4);
      printf(" Satellite ID  : %s\n",temp);
      memset(temp,' ',8);
      strncpy(temp,&line[24],2);
      printf(" Record type   : %s\n",temp);
      strncpy(temp,&line[26],12);
      printf(" Instrument ID : %s\n",temp);
      strncpy(temp,&line[38],10);
      printf(" Physical record count : %s\n",temp);
      memset(temp,' ',10);

      int4 = (int *)&line[48];
      Total = *int4;
      printf(" Total number of points in record : %d\n",Total); 

      int4 = (int *)&line[52];
      Actual = *int4;
/*      printf(" Number of actual points : %d\n",Actual); */

      int4 = (int *)&line[56];
      printf(" Starting index of first actual point : %d\n",*int4);

      int4 = (int *)&line[60];
      T1 = *int4;
      int4 = (int *)&line[64];
      T2 = *int4;
      printf(" UDTF time : %d    %d (year day, milliseconds)\n",T1,T2);

      real4 = (float *)&line[68];
      Lat = *real4;
      printf(" Latitude  : %f (degress)\n",Lat);

      real4 = (float *)&line[72];
      Lon = *real4;
      printf(" Longitude : %f (degress)\n",Lon);

      int4 = (int *)&line[84];
      Par_numb = *int4;
      printf(" Number of 32-bit Parameter words : %d\n",Par_numb);

      strncpy(temp,&line[88],8);
      printf(" JOB Version : %s\n",temp);

      strncpy(temp,&line[96],8);
      printf(" CDB Version : %s\n",temp);

      memset(temp,' ',10);

      printf(" Flag Quantity : %X\n",line[104]);
      printf(" Filter Number : %X\n",line[105]);
      printf(" Quality       : %X\n",line[106]); 
 
       rec++;

    }  
   /* if ret != EOF */
 
  } while (ret > 0); 


  fclose(data_ifp);
}