render_data.c

Go to the documentation of this file.

#include "prototypes.h"
 
void render_data(int type,int index)
{
extern float rainfall_level_1[20];
extern float rainfall_level_2[20];
extern int hpflag;
extern int hgflag;
extern int dflag[10];
extern struct nexrad nexrad_define[5];
extern struct dval dval;
extern int regrid_flag;
/*
extern struct custom custom[6000]; 
*/
extern struct custom custom[9000]; 
extern struct datalayer *sdatalayer;
extern struct saved *fsaved;
struct datalayer *datalayer;
struct saved *saved;
struct saved *ssaved;
int xadd,yadd;
int first;
float lat,lon;
int r,s;
float conv=.0174;
int h;
 
regrid_flag=1;
 
if(fsaved != NULL) 
{
  saved=fsaved;
  while(saved != NULL) 
  {
    ssaved=saved->next;
    free(saved);
    saved=ssaved;
  }
}
 
saved=NULL;
fsaved=NULL;
first=0;
datalayer=sdatalayer;
 
while( datalayer!=NULL) 
{
  if(datalayer->type != type) 
  {
    datalayer=datalayer->next;
    if(datalayer==NULL) 
    {
      if(first != 0)
        saved->next=NULL;
      break;
    }
    continue;
  }
                                    
  if(type==3 || type==4 || type==5 || type==6) 
  {
    if((hgflag==-1 && 
      ((datalayer->pc[0]=='H' && datalayer->pc[1]=='G') ||
      (datalayer->pc[0]=='Q' && datalayer->pc[1]=='T')))
      ||
      (hpflag==-1 &&
      ((datalayer->pc[0]=='H' && datalayer->pc[1]=='P') ||
      (datalayer->pc[0]=='Q' && datalayer->pc[1]=='I')))) 
    {
      datalayer=datalayer->next;
      if(datalayer==NULL) 
      {
        if(first != 0)
          saved->next=NULL;
        break;
      }
      continue;
    }
  }
 
  if(first==0) 
  {
    saved=(struct saved *) calloc(1,sizeof(struct saved));
    fsaved=saved;
    first=1;
  }
  else 
  {
    saved->next=(struct saved *) calloc(1,sizeof(struct saved));
    saved=saved->next;
  }
 
  lat=datalayer->lat;
  lon=datalayer->lon;
 
  if(lat <= 0 || lon <= 0) 
  {
    r=-99;
    s=-99;
  }
  else 
  {
    r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
             * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
    s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
            * sin((lon-dval.lo-90)*conv) + dval.yo + .5;       
  }
 
  saved->r=r;
  saved->s=s;
  saved->lat=lat;
  saved->lon=lon;
 
  strcpy(saved->name,datalayer->name);
  strcpy(saved->hb5,datalayer->hb5);
  strcpy(saved->pc,datalayer->pc);
  strcpy(saved->snum,datalayer->snum);
  saved->elev=(int)datalayer->elev;
  saved->rfs=-1;
     
  if(datalayer->lvl4 > 0)
    saved->rfs=1;
 
  if(index >= datalayer->num)
    saved->data=-9999;
  else           
    saved->data=datalayer->data[index];
     
  saved->parm=0;
  saved->fparm=0;
     
  if(type==14) 
  {
    if(saved->data >= rainfall_level_1[1])
      saved->parm=1;
    if(saved->data >= rainfall_level_2[1])
      saved->parm=2;
  }                    
               
  else if(type==15) 
  {
    if(saved->data >= rainfall_level_1[5])
      saved->parm=1;
    if(saved->data >= rainfall_level_2[5])
      saved->parm=2;
  }                    
                    
  else if(type==16) 
  {
    if(saved->data >= rainfall_level_1[8])
      saved->parm=1;
    if(saved->data >= rainfall_level_2[8])
      saved->parm=2;
  }                    
                                 
  else if(type==13) 
  {
    if(saved->data >= rainfall_level_1[index])
      saved->parm=1;
    if(saved->data >= rainfall_level_2[index])
      saved->parm=2;
  }                    
                                    
  else 
  {
    if(datalayer->lvl1 > -10 && 
       datalayer->data[index] >= datalayer->lvl1)
      saved->parm=1;
    if(datalayer->lvl2 > -10 && 
       datalayer->data[index] >= datalayer->lvl2)
      saved->parm=2;
  }
           
  if(datalayer->lvl3 ==1.0)
    saved->fparm=1;
  if(datalayer->lvl3 ==2.0)
    saved->fparm=2;
 
  saved->pparm = -1;
    
  if(index < 6 && (type == 3 || type == 4 || type == 5 || type==6) &&
     datalayer->data[index] < -98) 
  {
    for(h=index; h < 6; h++) 
    {
      if(datalayer->data[h] >= 0) 
      {
        if(datalayer->lvl1 > -10 && 
           datalayer->data[h] >= datalayer->lvl1)
          saved->pparm=1;
  
        if(datalayer->lvl2 > -10 && 
           datalayer->data[h] >= datalayer->lvl2)
          saved->pparm=2;
                                             
        break;
      }
    }
  }
         
  get_offset(datalayer->hb5,datalayer->pc,&xadd,&yadd);     
 
  saved->xadd=xadd;
  saved->yadd=yadd;
  datalayer=datalayer->next;
 
  if(datalayer==NULL) 
  {
    saved->next=NULL;
    break;
  }
}
 
 
/*
if((type==13 || type==14 || type==15 || type==16) &&
    dflag[0]==1)     
  read_and_sum_nexrad(type,index);
*/
 
return;
 
}
 
void get_offset(char *hb5,char *pc,int *xadd,int *yadd)
{
/*
extern struct custom custom[6000];
*/
extern struct custom custom[9000];
int i;
 
i=0;
 
while(custom[i].hb5[0]!=0) 
{
  if(strcmp(custom[i].hb5,hb5)==0 &&
     strcmp(custom[i].pc,pc)==0) 
  {
    *xadd=custom[i].xadd;
    *yadd=custom[i].yadd;
    return;
  }
  i++;
}
*xadd=0;
*yadd=0;
}