read_and_sum_nexrad.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
void read_and_sum_nexrad(int type,int index)
 
{
 
extern char nexrad_file[100];
extern int nexrad_flag;
extern Display *display;
extern Window window;
extern Cursor watch_cursor;
extern Widget drawing_area;
extern struct hrap_grid *hrap_grid;
extern int ditime[];
extern time_t itime[500],dtime[500],xtime[500],htime[500];
extern struct pcp *pcp;
extern struct pcp *spf;
int gmini,gminj,gmaxi,gmaxj,ghrapi,ghrapj;
int minhrapi,minhrapj,maxhrapi,maxhrapj;
int hours;
time_t btim,tim;
int i,j,k;
struct tm *gm;
FILE *fq;
char fbuf[100];
char nbuf[100];
int ii,jj;
int ier;
char ibuf[100];
Cardinal argcount;
Arg args[10];
XmString msg;
Widget cerr;
char kbuf[10000];
int iflag;
int h;
char *p;
 
if(nexrad_flag==1)
              return;
 
XDefineCursor(display,window,watch_cursor);
 
msg=XmStringCreateLocalized("Reading and summing NEXRAD precipitation");
 
argcount=0;
XtSetArg(args[argcount],XmNmessageString,msg);argcount++;
XtSetArg(args[argcount],XmNdialogStyle,
           XmDIALOG_FULL_APPLICATION_MODAL);argcount++; 
cerr=XmCreateErrorDialog(drawing_area,"info",args,argcount);
 
XtUnmanageChild(XmMessageBoxGetChild(cerr,XmDIALOG_CANCEL_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(cerr,XmDIALOG_OK_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(cerr,XmDIALOG_HELP_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(cerr,XmDIALOG_DEFAULT_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(cerr,XmDIALOG_SEPARATOR));
      
XtManageChild(cerr);
XmStringFree(msg);
 
XmtWaitUntilMapped(cerr);
 
for(i=0;i<hrap_grid->maxi;i++) {
 
for(j=0;j<hrap_grid->maxj;j++) {
   
      spf->value[i][j]=-9999;
      
        }
        
        }     
 
if(type==13) {
 
/* instantaneous */
 
               hours=ditime[index]/60;               
               btim=htime[0];   
 
                         }
 
else if(type==14) {
 
/* hourly */
 
               hours=1;
               btim=htime[index];              
 
                    }
 
else if(type==15)  {
 
/* six hourly */
 
               if(index==0) {              
 
                  hours=6;                         
                  btim=xtime[0];
                  hours=(xtime[0]-xtime[1])/3600;
                  
                 }
                 
               else  {
               
                  hours=6;                         
                  btim=xtime[index];
                  
               
               
                   }
               
               
               }
 
else if(type==16)  {
 
/* daily */
printf("daly\n");
               if(index==0) {
                                  
                  btim=dtime[0];
                  hours=(dtime[0] - dtime[1])/3600;
                  
                  }
 
               else  {
 
                  hours=24;
                  btim=dtime[index];
      
                  }
                  
                  printf("%d %d %d %d\n",dtime[1],dtime[0],hours,index);
            
                  }
                  
for(k=0;k<hours;k++) {
      
    for(i=0;i<hrap_grid->maxi;i++) {
 
    for(j=0;j<hrap_grid->maxj;j++) {
   
        pcp->value[i][j]=-9999;
      
        }
        
        }     
   
    
    tim=btim - 3600*k;
    gm=gmtime(&tim);
 
    sprintf(fbuf,"/awips/hydroapps/rfc/local/apps/mapper/nexrad/%02d/%02d/ngrid.%02d-%02d-%02d-%02d",
            gm->tm_year,gm->tm_mon+1,gm->tm_year,
            gm->tm_mon+1,gm->tm_mday,gm->tm_hour);
            
    sprintf(nbuf,"%s/ngrid.%02d-%02d-%02d-%02d.%d",
            nexrad_file,
            gm->tm_year,gm->tm_mon+1,gm->tm_mday,gm->tm_hour,
            getpid());
            
    sprintf(ibuf,"cp %s.gz %s.gz",fbuf,nbuf);
            
    ier=system(ibuf);
            
    if(ier!=0)
                    continue;
 
    sprintf(ibuf,"gzip -t %s.gz",nbuf);
            
    ier=system(ibuf);
      
    if(ier != 0)
                continue;
    
    sprintf(ibuf,"gzip -d -f %s.gz",nbuf);
            
    ier=system(ibuf);
            
    if(ier != 0)
              continue;
            
    fq=fopen(nbuf,"r");
 
   /* what if missing ?? */
 
    if(fq==NULL) {
 
                        printf("could not find %s\n",nbuf);
                        continue;
 
              }
 
    minhrapi=hrap_grid->hrap_minx;
    minhrapj=hrap_grid->hrap_miny;
    maxhrapi=hrap_grid->hrap_minx + hrap_grid->maxi;
    maxhrapj=hrap_grid->hrap_miny + hrap_grid->maxj;
 
    p=fgets(kbuf,100,fq);
  
    ier=sscanf(kbuf,"%d %d %d %d %d\n",&gmini,&gminj,&gmaxi,&gmaxj,&iflag);
 
    if(ier==4)
        iflag=0;
 
    for(i=0;i<gmaxi;i++) {
 
     if(iflag==0)
     fread(kbuf,sizeof(char),gmaxj*4,fq);
 
     else
     fread(kbuf,sizeof(char),gmaxj*6,fq);
 
     /* get hrap coord of gridded data */
 
     ghrapi=gmini+i;
 
     ii=ghrapi-minhrapi;
 
     if(ghrapi >= minhrapi && ghrapi < maxhrapi) {  
 
     h=0;
 
     for(j=0;j<gmaxj;j++) {
 
            ghrapj=gminj+j;
 
            jj=ghrapj-minhrapj;
 
            if(ghrapj >= minhrapj && ghrapj < maxhrapj)  {
 
                        pcp->value[ii][jj]=atoi(&kbuf[h]);
                        
 
                }
 
            if(iflag==0)
                 h=h+4;
 
            else
                 h=h+6;
 
      }
 
    }
 
    }
 
    fclose(fq);
 
    for(i=0;i<hrap_grid->maxi;i++) {
 
    for(j=0;j<hrap_grid->maxj;j++) {
 
      if(spf->value[i][j] < 0 && pcp->value[i][j] >=0)
         spf->value[i][j]=pcp->value[i][j];
   
      else if(pcp->value[i][j] >=0)
         spf->value[i][j]=spf->value[i][j]+pcp->value[i][j];   
      
        }
        
        }     
    
}   
 
sprintf(ibuf,"rm %s/*%d*",nexrad_file,getpid());
 
system(ibuf);
 
XtDestroyWidget(cerr);
 
nexrad_flag=1;
 
}