render_pcp.c

Go to the documentation of this file.

#include "prototypes.h"
 
void render_pcp(int pcpn_day,int pcpn_time,int pcpn_time_step)
 
{
 
extern Widget top_level;
extern int pcp_in_use[50];
extern int pcp_flag;
extern Cursor watch_cursor;
extern int method;
int i,j,h,hh,time_pos,htotal;
double distance,dist1,dist2,dist,value;
double temp;
extern struct hrap_grid *hrap_grid;
extern struct pcp *pcp;
extern Window window;
extern Display *display;
extern int max_stations;
 
XDefineCursor(display,window,watch_cursor); 
XmUpdateDisplay(top_level);
 
if(pcpn_time_step==0)
              time_pos=pcpn_time;
 
else
              time_pos=4;
 
for(i=0;i<hrap_grid->maxi;i++) {
 
for(j=0;j<hrap_grid->maxj;j++) {
 
           value=0.0;
           distance=0.0;
           htotal=0;
 
           for(h=0;h<20;h++) {
 
               hh=hrap_grid->gage[i][j].index[h];
 
               htotal++;
/*
               dist1=(double)(hrap_grid->coord[i][j].x-
                      station[hh].x);    
               dist2=(double)(hrap_grid->coord[i][j].y-
                  station[hh].y);  
*/
               dist=pow(dist1,2)+pow(dist2,2);
 
               if(dist < 1.0)
                        dist=1.0;
               
               dist=1/dist;
/*
               temp=pdata[pcpn_day].stn[hh].frain[time_pos].data*dist;
 
               if(method==2 && station[hh].isoh > 0)
                     temp=temp*(float)hrap_grid->isoh[i][j]/station[hh].isoh;
*/
               value=value+temp;
 
               distance=distance+dist;                    
 
         }
 
           if(htotal==0) {
 
           for(h=0;h<max_stations;h++) {
 
               htotal++;
/*
               dist1=(double)(hrap_grid->coord[i][j].x-
                      station[h].x);    
               dist2=(double)(hrap_grid->coord[i][j].y-
                      station[h].y);  
*/
               dist=pow(dist1,2)+pow(dist2,2);
 
               if(dist < 1.0)
                        dist=1.0;
               
               dist=1/dist;
/*
               temp=pdata[pcpn_day].stn[h].frain[time_pos].data*dist;
 
               if(method==2 && station[h].isoh > 0)
                     temp=temp*(float)hrap_grid->isoh[i][j]/station[h].isoh;
*/
               value=value+temp;
 
               distance=distance+dist;                    
 
         }
 
     }
 
           if(htotal == 0)
               pcp->value[i][j]=0;
  
           else
               pcp->value[i][j]=(int)(value/distance*100);
  
           if(pcp->value[i][j] < .01)
                    pcp->value[i][j]=0.0;
     
 
     }
 
}
 
 
if(pcpn_time_step==0)
              time_pos=pcp_flag;
 
else
              time_pos=pcp_flag+pcpn_day;
 
pcp_in_use[time_pos]=1;
 
write_file("pcp",time_pos,pcp);
 
}