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

#define MAXBUF 81

main (argc,argv)
 int argc;
 char *argv[];
{  
   char filename[16];
   char id[4],sp[4];
   char t2[5],title[6],mon[3],day[3],hr[3],t3[6];
   char temp[MAXBUF],temp2[MAXBUF],t4[5];
   char loni[7],lonf[6],dlon[5],lati[6],latf[6],dlat[5];
   char  nlon[3],nlat[3],surp[4],surt[4];
   float longit[82],latit[50],fdata;
   float holdit[3970];
   char test[13];
   int i,j,l,kk, m,rec;
   FILE *f1;

/* the latitude and longitude range is from 35W to 5W, 20N to 45N at surface pressure 
   SP represents Surface pressure in Pascals and ST is the Sea Surface Temperature 
   in degrees Kelvin */

   banner(argc);
   if (argc == 1) {
     printf("\nPlease enter in the input file name\n\n"); 
     scanf("%s",filename);
   }
   else
     strcpy(filename,argv[1]);
   
   
   if((f1 = fopen(filename,"r")) == NULL)
   {  printf("Can't open the input file\n");
      exit(1);
   }

   
   fscanf(f1,"%s%s%s%s%s%s%s%s%s%s",title,t2,mon,day,hr,t3,t3,id,t3,sp);
   
   printf("Month is %s day %s hour %s ",mon,day,hr);
   printf("and the pressure level is %s Pascal.\n",sp);

   fscanf(f1,"%7s%7s%7s%7s%7s%7s%7s%7s",t4,t4,t4,t4,t4,t4,t4,t4);
   fscanf(f1,"%s%s%s%s%s%s%s%s", loni,lonf,dlon,nlon,lati,latf,dlat,nlat);
   if (atoi(lonf) == 355.00)
     strcpy(lonf,"5.00");
   printf("Initial longitude = %s and the final longitude = %s.\n\n\n", loni,lonf);
   printf("The longitude delta is %s and the number of longitude positions are %s.\n\n",dlon,nlon);
   printf("Initial latitude = %s and the final latitude = %s with a delta of %s\n",lati,latf,dlat);
   printf("providing %s latitude postions.\n",nlat);

/* The remaining data is a sweep of latitude numbers per longitude position */
   

   longit[0] = atof(loni);
   latit[0]  = atof(lati);
   fgets(&temp[0],MAXBUF,f1);
   i=0;
   j=0;
   while (feof(f1) == 0) {
     fgets(&temp[0],MAXBUF,f1);
     kk=1;
     for (l=0;l<6;l++) {
       strncpy(test,&temp[kk],13);
       holdit[m] = atof(test);
       kk = kk+13;
       m++;
     }
    } /* while */
   m = 0;
   for (i=0;i<atoi(nlon);i++) {
   printf("\n\n     Longitude         Latitude");
   
   switch (atoi(id)) {
      case 142:
        printf(" Large scale Precipitation (m/6hr)\n");
        break;
      case 143:
        printf(" Convective Precipitation (m/6hr)\n");
        break;
      case 145:
        printf(" Boundary layer dissipation (W/m^2)\n");
        break;
      case 146:
        printf(" Surface Sensible heat flux (W/m^2)\n");
        break;
      case 147:
        printf(" Surface latent heat flux (W/m^2)\n");
        break;
      case 164: 
        printf(" Total cloud cover \n");
        break;
      case 165: 
        printf(" 10 meter u (m/s) \n");
        break;
      case 166:
        printf(" 10 meter v (m/s) \n");
        break;
      case 167:
        printf(" 2 meter Temperature (K) \n");
        break;
      case 168:
        printf("2 meter Dewpoint Temperature (K) \n");
        break;
      case 176:
        printf("Surface Solar radiation (W/m^2) \n");
        break;
      case 177:
        printf("Surface thermal radiation (W/m^2) \n");
        break;
      case 178: 
        printf("Top Solar Readiation (W/m^2) \n");
        break;
      case 179:
        printf("Top Thermal Radiation (W/m^2) \n");
        break;
      case 180:
        printf("U Stress (N/m^2) \n");
        break;
      case 181:
        printf("V Stress (N/m^2) \n");
        break;
      case 182:
        printf("Evaporation (m H2O) \n");
        break;
      case 185:
        printf("Convective Cloud Cover \n");
        break;
      case 186:
        printf("Low Cloud Cover \n");
        break;
      case 187:
        printf("Medium Cloud Cover \n");
        break;
      case 188:
        printf("High Cloud Cover \n");
        break;
      case 208:
        printf("Top Solar Radiation Upward (W/m^2) \n");
        break;
      case 209:
        printf("Top Thermal Radiation Upward (W/m^2) \n");
        break;
      case 210:
        printf("Top Solar Radiation Upward Clear Sky \n");
        break;
      default:
        printf("\n\n Error: unknown field identifier \n");
        exit (0);
   } /* case */
   printf("---------------------------------------------------------------\n");
    for (j=0; j<atoi(nlat);j++) {
     printf( "     %f         %f         %13.6e\n", longit[i],latit[j], holdit[m]);
     m++;
     latit[j+1] = latit[j]-atof(dlat); 
    }     
   longit[i+1]=longit[i] + atof(dlon);
   }

   printf(" \n\nThe Read Program for ASTEX Data Set from ECMWF has completed.\n");
   rec = atoi(nlat)*atoi(nlon);
   printf(" There were %d records displayed\n\n", rec);

}

banner(argc)
{
  char cont;
 
  printf("**********************************************************\n\n");
  printf("  FIRE ASTEX Surface Diagnostic from ECMWF READ PROGRAM \n\n");
  printf("  Version 1.0                                  July 1995 \n\n");
  printf("  This program reads the Sfdiag Data generated in June 1992. The\n");
  printf("  input file name is built according to the naming convention formed\n");
  printf("  by the Langley DACC.  If you have changed the name of the\n");
  printf("  input file or need to specify a path, you will need to   \n");
  printf("  modify this program.  All data is expected to be in the  \n");
  printf("  current working directory.                               \n");
  printf("\n\n Platforms supported: SUN \n");
  printf("                      HP \n");
  printf("                      SGI \n");
  printf("                      Dec Alpha \n");
  printf(" \n\n  Contact: Langley DAAC User and Data Services Office \n");
  printf("           NASA Langley Research Center \n");
  printf("           Mail Stop 157D \n");
  printf("           Hampton, VA 23681-0001 \n");
  printf("           Phone (804)864-8656   Fax (804)864-8807 \n");
  printf("**********************************************************\n\n\n");
/*    printf("\n\n Hit a carriage control to continue. \n");
    scanf("%c",&cont);*/
return;
}