calculate_NEXRAD.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
void calculate_NEXRAD(int inex,int time_step)
 
{
 
extern struct nexrad *nexrad[240];
extern struct nexrad nexrad_final;
int i,j,k,kk;
double pw,tpw;
 
if(time_step==1) {
 
                 memcpy(&nexrad_final,nexrad[inex],sizeof(struct nexrad));
                 return;
 
           }
 
strcpy(nexrad_final.date,nexrad[inex]->date);
strcat(nexrad_final.date," thru ");
strcat(nexrad_final.date,nexrad[inex+time_step]->date);
 
for(i=0;i<131;i++) {
for(j=0;j<131;j++) {
     
     nexrad_final.value[i][j]=0;
     tpw=0;
 
     for(k=inex;k<time_step+inex;k++) {                                
 
             if(nexrad[k]==NULL) {
 
                           printf("error in allocation\n");
                           continue;
 
             }
 
             if(nexrad[k]->flag < 0)
                               continue;
 
             pw=(-170+((float)nexrad[k]->value[i][j]*5))/100.;
             pw=.03937 * pow(10,pw);
 
             if(pw < .01)
                       continue;
 
             tpw=tpw+pw;       
            
 
                       }
 
     if(tpw <= .01)
              continue;
 
     tpw=tpw/.03937;
 
     pw=log10(tpw);
 
     kk=(int)(100*pw+170)/5;            
     nexrad_final.value[i][j]=kk;
 
     }
 
}
 
return;
 
}