make_NEXRAD.c

Go to the documentation of this file.

 
#include "prototypes.h"
char fbuf[100];
 
void main(int argc, char **argv)
{
 
 struct nexrad nexrad;
 int i,j,k,ii,jj,kk,l,mm;
 struct HRAP hrap;
 struct HRAP irap;
 float cx,cy,lat,lon;
 FILE *fp,*fq,*fw;
 char *p,kbuf[120],dbuf[200],mbuf[200];
 float xlat,nlat,nlon,xlon;
 int nrapx[131][131],nrapy[131][131];
 char name[100];
 time_t ltim;
 float hmaxx,hmaxy,hminx,hminy;
 struct nexstruct nexstruct[25];
 float pw,pnexrad[255];
 int imaxx,iminx,imaxy,iminy;
 char fbuf[100];
 int maxi,maxj;
 time_t tim,stim;
 struct tm *gm;
 int irain;
 char tarbuf[1000];
 int nexarray[20][24];
 char *nexname[20];
 int ier;
 int m, maxdays=1;
 float frain;
 int  iflag; 
 int h;
 int ilen,slen;
 char ibuf[100];
 float val[500][500];
 int month,day,year,hour;
 
 month=atoi(argv[1]);
 day=atoi(argv[2]);
 year=atoi(argv[3]);
 hour=atoi(argv[4]);
 
 
 /* initialize val */
 
 for(i=0;i<500;i++) {
 for(j=0;j<500;j++) {
 
           val[i][j]=-9999;
           
           }
           
           }
 
 /* initialize name */
 
 for(i=0;i<20;i++)
        nexname[i]=calloc(20,1);
 
/* initialize pnexrad */
 
 for(i=0;i<255;i++) {
 
      pw=(-6) +( (float)(i-1))*.125;
      pw=pw/10;
      pw=.0397 * pow(10,pw);
 
      if(i==0)
             pnexrad[i]=0;
 
      else
             pnexrad[i]=pw;
 
    }
 
 /* pars argv for lat,lon,date and hour */
        
 sprintf(dbuf,"/usr/mapper/NEXRAD_radars");
 
 fp=fopen(dbuf,"r");
 
 if(fp==NULL) {
 
              printf("could not open /usr/mapper/NEXRAD_radars\n");
              exit(1);
 
              }
 
 p=fgets(kbuf,80,fp);
 if(p==NULL) {
 
              printf("error in format /usr/mapper/NEXRAD_radars\n");
              exit(1);
 
           }
 
 p=pars_line(kbuf,"x_lat=",dbuf);
 
 if(p==NULL) {
 
              printf("no x_lat\n");
              exit(1);
              
              }
              
 xlat=atof(dbuf);
   
 p=pars_line(kbuf,"n_lat=",dbuf);
 
 if(p==NULL) {
 
              printf("no n_lat\n");
              exit(1);
              
              }
              
 nlat=atof(dbuf);
 
 p=pars_line(kbuf,"x_lon=",dbuf);
 
 if(p==NULL) {
 
              printf("no x_lon\n");
              exit(1);
              
              }
              
 xlon=atof(dbuf);
   
 p=pars_line(kbuf,"n_lon=",dbuf);
 
 if(p==NULL) {
 
              printf("no n_lon\n");
              exit(1);
              
              }
              
 nlon=atof(dbuf);
  
 hmaxx=-1;hmaxy=-1;
 hminx=9999;hminy=9999;
  
 hrap=LatLongToHrap(xlat,xlon);
 
 if(hrap.x > hmaxx)
             hmaxx=hrap.x;
 
 
 if(hrap.x < hminx)
             hminx=hrap.x;
 
 
 if(hrap.y > hmaxy)
             hmaxy=hrap.y;
 
 
 if(hrap.y < hminy)
             hminy=hrap.y; 
 
 hrap=LatLongToHrap(xlat,nlon);
 
 if(hrap.x > hmaxx)
             hmaxx=hrap.x;
 
 
 if(hrap.x < hminx)
             hminx=hrap.x;
 
 
 if(hrap.y > hmaxy)
             hmaxy=hrap.y;
 
 
 if(hrap.y < hminy)
             hminy=hrap.y; 
 
 hrap=LatLongToHrap(nlat,xlon);
 
 if(hrap.x > hmaxx)
             hmaxx=hrap.x;
 
 
 if(hrap.x < hminx)
             hminx=hrap.x;
 
 
 if(hrap.y > hmaxy)
             hmaxy=hrap.y;
 
 
 if(hrap.y < hminy)
             hminy=hrap.y; 
 
 hrap=LatLongToHrap(nlat,nlon);
 
 if(hrap.x > hmaxx)
             hmaxx=hrap.x;
 
 
 if(hrap.x < hminx)
             hminx=hrap.x;
 
 
 if(hrap.y > hmaxy)
             hmaxy=hrap.y;
 
 
 if(hrap.y < hminy)
             hminy=hrap.y; 
 
 imaxx=hmaxx+.5;
 iminx=hminx+.5;
 
 imaxy=hmaxy+.5;
 iminy=hminy+.5;
 
 maxi=imaxx-iminx;
 maxj=imaxy-iminy;
   
 /* need times to get for each file */
 
 for(l=0;l<20;l++) {
 
     p=fgets(kbuf,80,fp);
     if(p==NULL) {
 
                 fclose(fp);
                 break;
 
           }
 
     p=pars_line(kbuf,"radar=",dbuf);
 
     if(p==NULL)
               continue;
                   
     strcpy(name,dbuf); 
 
     for(i=0;i<3;i++)
              name[i]=tolower(name[i]);
 
     name[3]=0;
 
     strcpy(nexname[l],name);
 
     p=pars_line(kbuf,"latitude=",dbuf);
 
     if(p==NULL)
               continue;
 
     lat=atof(p);
 
     p=pars_line(kbuf,"longitude=",dbuf);
 
     if(p==NULL)
               continue;
 
     lon=atof(p);
 
     if(lon > 0)
              lon=-lon; 
    
     /* get hrap coordinates of radar */
 
     hrap=LatLongToHrap(lat,-lon);
 
     cx=(int)hrap.x+1;
     cy=(int)hrap.y;
 
     for (i=0; i < 131; i++){     
     for (k=0; k < 131 ;k++){
 
 
        irap.x=cx + (float)(k - 65);
        irap.y=cy + (float)(i-65);
 
        nrapx[i][k]=irap.x;
        nrapy[i][k]=irap.y;        
 
         }
 
        }
 
     sprintf(fbuf,"/usr/mapper/nexrad/scratch/%s.%02d-%02d-%02d-%02d",
             name,year,month,day,hour);
 
     fq=fopen(fbuf,"rb");
 
             /* what if missing ?? */
 
     if(fq==NULL) {
 
                        printf("could not find %s\n",fbuf);
                        continue;
 
              }
 
     fread(&nexrad,sizeof(char),131*131,fq);
 
     fclose(fq);
             
     for (i=0; i < 131; i++) {
     for (j=0; j < 131 ;j++) {
 
                     ii=nrapx[i][j]-iminx;
                     jj=nrapy[i][j]-iminy;
 
                     if(ii < 0 || ii > maxi ||
                        jj < 0 || jj > maxj)
                                       continue;
                     
                     if(nexrad.value[i][j] == 255)
                                             continue;
 
                     if(pnexrad[nexrad.value[i][j]] > 
                        val[ii][jj])
                        val[ii][jj]=pnexrad[nexrad.value[i][j]];
 
           }
 
       }
 
       }
 
   }
 
/* now add up all numbers and divide by total contributions */
 
fclose(fp);
 
sprintf(fbuf,"/usr/mapper/nexrad/ngrid.%02d-%02d-%02d-%02d",year,month,day,hour);
 
printf("fbuf is %s\n",fbuf);
 
fp=fopen(fbuf,"w");
 
fprintf(fp,"%d %d %d %d 1\n",iminx,iminy,maxi,maxj);
 
for(i=0;i<maxi;i++) {
 
for(j=0;j<maxj;j++) { 
                       
 
                    if(val[i][j] >=0)
                    irain=val[i][j]* 100;
                   
                    else
                      irain=-9999;
                   
                    fprintf(fp," %5d",irain);
 
          }
 
       }
 
fclose(fp);
 
sprintf(tarbuf,"gzip -f %s",fbuf);
 
system(tarbuf);
 
}
 
}