grid_data.c

Go to the documentation of this file.

#include "prototypes.h"
 
void grid_data(int type)
 
{
extern struct dval dval;
extern struct grid *grid;
extern struct saved *fsaved;
extern Widget top_level;
extern Cursor watch_cursor;
extern int method;
int i,j,htotal;
double distance,dist1,dist2,dist,value;
double temp;
extern struct hrap_grid *hrap_grid;
extern Window window;
extern Display *display;
struct saved *saved;
int dir,speed,gust;
float fvalue;
 
XDefineCursor(display,window,watch_cursor); 
XmUpdateDisplay(top_level);
 
for(i=0;i<hrap_grid->maxi;i++) {
 
for(j=0;j<hrap_grid->maxj;j++) {
 
           value=0.0;
           distance=0.0;
    
           htotal=0;
           saved=fsaved;
 
           while( saved!=NULL) {
 
               if(saved->data < -98)
                              goto nextseg;
 
               if(saved->lat < 0 || saved->lon < 0)
                              goto nextseg;
 
               dist1=(double)(hrap_grid->coord[i][j].x-saved->r);
               dist2=(double)(hrap_grid->coord[i][j].y-saved->s);
 
               dist=pow(dist1,2)+pow(dist2,2);
 
               if(dist < 1.0)
                        dist=1.0;
               
               dist=1/dist;
 
               if(type==11) {
                           
                  dir=saved->data/1000000;
                  speed=(saved->data - dir*1000000)/1000;
                  gust=(saved->data - dir*1000000 - speed*1000);
                 
                  if(gust != 0)
                             fvalue=gust;
 
                  else
                             fvalue=speed;
 
                          } 
 
                else
                             fvalue=saved->data;
 
 
               temp=fvalue*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;    
               htotal++;
               
               nextseg:
 
               saved=saved->next;
 
               if(saved==NULL)
                          break;                
 
             }
 
           if(htotal == 0)
               grid->value[i][j]=0;
  
           else
               grid->value[i][j]=value/distance;
 
     }
 
 
}
 
}