plot_data.c

Go to the documentation of this file.

#include "prototypes.h"
void plot_data(Drawable pix,int plot_view,int type,int h,int display_flag)
{
extern int zoom_mode;
extern int rfs_flag;
extern float filter_value;
extern int elevation_filter_value;
extern time_t ditime[],xtime[500],htime[500],dtime[500];
extern char *legend[100];
extern int indexx;
extern time_t qtime;
extern struct saved *fsaved;
extern char find_station_buf[10];
extern struct dval dval;
extern int npoint[2],zoom;
extern GC gc;
extern unsigned long qmap[16];
extern Widget drawing_area;
extern Display *display;
extern Font font[10];
extern XFontStruct *info_font[10];
extern struct station station[1500];
extern int max_stations;
extern int dflag[10];
extern long amap[16];
int i,xc,yc,text_width,text_width2,length;
int x1,y1,xadd,yadd,MY=10;
float xrat,yrat,mult,dmult,cmult;
long xmin,xmax,ymin,ymax;
signed long xcen,ycen,lint1,lint2;
Dimension width,height;
Arg args[10];
signed long XSIZE,YSIZE;
extern struct display_set display_set[4];
long r,s;
float conv=.0174,fval;
float lat,lon;
char tbuf[100];
int yheight,ival;
struct tm *gm;
time_t ntime;
struct saved *saved;
int hour;
int year,yearp;
int mon;
int day;
char mbuf[100];
char *month[]=
{
"December","January","February","March","April",
"May","June","July","August","September",
"October","November","December"
};
 
 
if(type==-1)
  return;
 
XSetFont(display,gc,font[4]);
 
yheight=info_font[4]->ascent;
XSetForeground(display,gc,qmap[15]);
 
XtSetArg(args[0],XtNwidth,&width);
XtSetArg(args[1],XtNheight,&height);
XtGetValues(drawing_area,args,2);
 
XSIZE=(long)width*10L;
YSIZE=(long)height*10L;
 
xrat=(float)XSIZE/(float)12800;
yrat=(float)YSIZE/(float)9600;
 
if(xrat < yrat)
  mult=xrat;
else
  mult=yrat;
 
dmult=display_set[display_flag].dmult[h];
xmin= display_set[display_flag].xmin[h];
xmax= display_set[display_flag].xmax[h];
ymin= display_set[display_flag].ymin[h];
ymax= display_set[display_flag].ymax[h];
xcen= display_set[display_flag].xcen[h];
ycen= display_set[display_flag].ycen[h];
 
lint1=(long)npoint[0]*(long)MY;
lint2=(long)npoint[1]*(long)MY;
 
cmult=mult*dmult;
 
saved=fsaved;
 
while( saved!=NULL) 
{
  lat=saved->lat;
  lon=saved->lon;
 
  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;       
 
  if(zoom==1) 
  {
    x1=(xmin+(float)r * cmult)/MY;
    y1=(ymin+(float)s * cmult)/MY;
  }
  else      
  {
    x1=(((xmin+(float)r*cmult)-lint1)*zoom+xcen)/MY;
    y1=(((ymin+(float)s*cmult)-lint2)*zoom+ycen)/MY;
  }
    
  if(saved->pparm==1)
    XSetForeground(display,gc,qmap[11]);
  else if(saved->pparm==2 )
    XSetForeground(display,gc,qmap[13]);     
  else if(saved->data < -98)
    XSetForeground(display,gc,qmap[1]);
  else if(saved->parm==1)
    XSetForeground(display,gc,qmap[11]);
  else if(saved->parm==2 )
    XSetForeground(display,gc,qmap[13]);
  else 
    XSetForeground(display,gc,qmap[15]);
 
  if(strncmp(&saved->pc[4],"P",1)==0 && dflag[1] != 1)
    goto nextseg;
 
  if(strncmp(&saved->pc[4],"M",1)==0 && dflag[2] != 1)
    goto nextseg;
 
  if(strncmp(&saved->pc[4],"R",1)==0 && dflag[3] != 1)
    goto nextseg;
 
  if(strncmp(&saved->pc[4],"G",1)==0 && dflag[4] != 1)
    goto nextseg;
 
  if(strncmp(&saved->pc[4],"Z",1)==0 && dflag[5] != 1)
    goto nextseg;
                      
  if((type==3 || type==4 || type==5 || type==6) && rfs_flag==1 &&
      saved->rfs!=1)
    goto nextseg;
 
  if(elevation_filter_value > 0 && 
     elevation_filter_value > saved->elev)
    goto nextseg;
 
/* locate station in data stream */
 
  if(plot_view==6 && saved->data < -98 && 
     saved->fparm < 1 && saved->pparm < 1)
   goto nextseg;
 
  XDrawLine(display,pix,gc,x1+1,y1+1,x1+1,y1-1);
  XDrawLine(display,pix,gc,x1+1,y1-1,x1-1,y1-1);
  XDrawLine(display,pix,gc,x1-1,y1-1,x1-1,y1+1);
  XDrawLine(display,pix,gc,x1-1,y1+1,x1+1,y1-1);
 
 if(plot_view==1)
   strcpy(tbuf,saved->hb5);
 else if(plot_view==2) 
   strcpy(tbuf,saved->snum);
 else if(plot_view==3)
   strcpy(tbuf,saved->pc);
 else if(plot_view==4)
   strcpy(tbuf,saved->name);
 else if(plot_view==5) 
 {
   if(saved->elev < 0)
     goto nextseg;
   sprintf(tbuf,"%d",saved->elev);
 }
 else if(plot_view==6) 
 {
   if(filter_value > -9998 && saved->data < filter_value)
     goto nextseg;
         
   if(saved->data < -98) 
   {
     sprintf(tbuf,"M");                         
   }
 
   else if(type==0 || type==1 || type==2 || type==6 || type==7 ||
           type==8 || type==9 || type==10 || type==30) 
   {
     if(saved->data >= 0)
       ival=saved->data+.5; 
     else
       ival=saved->data -.5;
                     
     sprintf(tbuf,"%d",ival);
   }
 
   else if(type==3 || type==5) 
   {
     if(saved->data >= 0)
       fval=saved->data+.0005;
     else
       fval=saved->data-.0005; 
     sprintf(tbuf,"%5.2f",fval);
   }
 
   else if(type==4) 
   {
     if(saved->pc[1]=='G' || saved->pc[1]=='P' || saved->pc[0]=='Q')
       fval=(saved->data+0.0005)/1000.;
     else
       fval=saved->data;
     sprintf(tbuf,"%6.2f",fval);
   }
 
   else if(type==24 || type==26 || type==28 || type==29) 
   {
     fval=saved->data+.5; 
     sprintf(tbuf,"%4d",(int)fval);
   }
 
   else if(type==21 || type==22 || type==23 || type==25 || 
           type==13 || type==14 || type==15 || type==16) 
   {
     fval=saved->data+.0005; 
     sprintf(tbuf,"%5.2f",fval);
   }
 
   else if(type==17 || type==18 || type==27) 
   {
     if(saved->data >= 0)
       fval=saved->data+.05;
     else
       fval=saved->data-.05; 
     sprintf(tbuf,"%5.1f",fval);
   }
 
   else
     strcpy(tbuf,"dummy");
 
   for(i=0;i<strlen(tbuf);i++) 
     if(tbuf[i]!=' ')
       break;
               
   strcpy(tbuf,&tbuf[i]);
 }
 
 if(plot_view==6 && type==11) 
 {
   draw_wind(pix,x1,y1,saved->data);
   goto nextseg;
 } 
 
 length=strlen(tbuf);
 text_width=XTextWidth(info_font[4],tbuf,length); 
 text_width2= 0;
     
 if(saved->fparm==1 || saved->fparm==2) 
 {
   text_width2=XTextWidth(info_font[4],"^",1); 
 }
     
 xadd=saved->xadd;
 yadd=saved->yadd;
 
 if(xadd < 0) 
   xc=x1 - text_width;
 else
   xc=x1+3+text_width2;
 
 if(yadd < 0)
   yc=y1;
 else
   yc=y1+yheight;
 
 XDrawString(display,pix,gc,xc,yc,tbuf,length);
 
 if(saved->fparm==1 || saved->fparm==2) 
 {
   if(saved->fparm==1)
     XSetForeground(display,gc,qmap[11]);
 
   else if(saved->fparm==2)
     XSetForeground(display,gc,qmap[13]);
 
   if(xadd  < 0)
     XDrawString(display,pix,gc,x1-text_width-text_width2,yc,"^",1);          
 
   else
     XDrawString(display,pix,gc,x1+3,yc,"^",1);  
 }
 
 if(strcmp(find_station_buf,saved->hb5)==0) 
 {
   find_station_buf[0]=0;
   XDrawLine(display,pix,gc,xc,yc,xc+text_width,yc);
 }
 
 nextseg:
 saved=saved->next;
 if(saved==NULL)
   break;
}          
 
/* build legend &/
 
/* hourly instantaneous - 0,1,3,4,5,6,8,9,10,11 */
 
  if(type==0 || type==1 || type==3 || type==4  || type==5 ||
     type==6 || type==8 || type==9 || type==10 || type==11 || type==30) 
  {
    ntime=qtime - indexx*3600;
    gm=gmtime(&ntime);
    sprintf(tbuf,"%02d-%02d-%04d %02dZ",gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
  }
 
  /* daily instantaneous 2,7 */
  else if(type==2 || type==7) 
  {
    gm=gmtime(&qtime);
    ntime=qtime - indexx*3600*24;
    if(gm->tm_hour <= 8)
      ntime=ntime-3600*24;
                               
    gm=gmtime(&ntime);
    sprintf(tbuf,"%02d-%02d-%04d",gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
  }
         
  else if(type==13) 
  {
    if(indexx==0)
      sprintf(tbuf,"Last 30 minutes");
    else      
    {
      hour=ditime[indexx]/60;
      sprintf(tbuf,"Last %d hours",hour);
    }
  }
 
  else if(type==14) 
  {
    if(indexx > 0) 
    {
/*
      gm=gmtime(&htime[indexx+1]);
      sprintf(tbuf,"%02d-%02d-%04d %02dZ to ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
*/
      sprintf(tbuf,"Ending ");
      gm=gmtime(&htime[indexx]);
      sprintf(mbuf,"%02d-%02d-%04d %02dZ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
      strcat(tbuf,mbuf);
    }
    else         
    {
      gm=gmtime(&htime[indexx]);
      sprintf(tbuf,"Since %02d-%02d-%04d %02dZ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    }
  }
 
  else if(type==15)  
  {
    gm=gmtime(&xtime[indexx+1]);
    if(!indexx)
      sprintf(tbuf,"%02d-%02d-%04d %02dZ to ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    else
      sprintf(tbuf,"Ending ");
    gm=gmtime(&xtime[indexx]);
    sprintf(mbuf,"%02d-%02d-%04d %02dZ",
            gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    strcat(tbuf,mbuf);
  }
 
  else if(type==16)  
  {
    gm=gmtime(&dtime[indexx+1]);
    if(!indexx)
      sprintf(tbuf,"%02d-%02d-%04d 12Z to ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
    else
      sprintf(tbuf,"Ending ");
 
    gm=gmtime(&dtime[indexx]);
    if(!indexx)
      sprintf(mbuf,"%02d-%02d-%04d %02dZ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    else
      sprintf(mbuf,"%02d-%02d-%04d 12Z",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
    strcat(tbuf,mbuf);
  }
 
  else if(type==21)  
  {
    gm=gmtime(&xtime[indexx+1]);
    if(!indexx)
      sprintf(tbuf,"%02d-%02d-%04d %02dZ to ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    else
      sprintf(tbuf,"Ending ");
    gm=gmtime(&xtime[indexx]);
    sprintf(mbuf,"%02d-%02d-%04d %02dZ",
            gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    strcat(tbuf,mbuf);
  }
 
  else if(type==22)  
  {
    gm=gmtime(&dtime[indexx+1]);
    if(!indexx)
      sprintf(tbuf,"%02d-%02d-%04d 12Z to ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
    else
      sprintf(tbuf,"Ending ");
 
    gm=gmtime(&dtime[indexx]);
    if(!indexx)
      sprintf(mbuf,"%02d-%02d-%04d %02dZ",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour);
    else
      sprintf(mbuf,"%02d-%02d-%04d 12Z",
              gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
    strcat(tbuf,mbuf);
  }
 
  else if(type==23 || type==24)
  {
    gm=gmtime(&dtime[0]);
    year=gm->tm_year+1900;
    mon=gm->tm_mon+1;
    if(!indexx)
      sprintf(tbuf,"%s %4d (to date)",month[mon],year);
    else
    {
      mon=mon-indexx;
      while(mon<=0)
      {
        mon=mon+12;
        year=year-1;
      }
      sprintf(tbuf,"%s %4d",month[mon],year);
    }
  }
                        
  else if(type==25 || type==26)
  {
    gm=gmtime(&dtime[0]);
    year=gm->tm_year+1900;
    mon=gm->tm_mon+1;
    day=gm->tm_mday;
    if(indexx==0)
    {
      if(mon>=10) 
        yearp=year;
      else
        yearp=year-1;
      sprintf(tbuf,"October %04d thru %s %02d %04d",yearp,month[mon],day,year);
    }
 
    if(indexx==1)
    {
      mon=mon-1;
      if(mon<=0)
        year--;
      if(mon>=10) 
        yearp=year;
      else
        yearp=year-1;
      sprintf(tbuf,"October %04d thru %s %04d",yearp,month[mon],year);
    }
 
    if(indexx==2)
    {
      if(mon<=10)
        year--;
      sprintf(tbuf,"October %04d thru September %04d",year-1,year);
    }
  }
 
  else if(type==27 || type==28 || type==29)
  {
    gm=gmtime(&dtime[0]);
    year=gm->tm_year+1900;
    mon=gm->tm_mon+1;
    mon=mon-indexx;
    while(mon<=0)
    {
      mon=mon+12;
      year=year-1;
    }
    if(mon==12)
    {
      indexx=indexx+6;
      mon=6;
    }
    if(mon==7)
    {
      indexx=indexx-6;
      if(indexx < 0) 
        indexx=0;
      mon=1;
      year=year+1;
    }
    sprintf(tbuf,"%s 1 %4d",month[mon],year);
  }
                        
  else if(type==17)  
  {
    gm=gmtime(&dtime[indexx+1]);
    sprintf(tbuf,"%02d-%02d-%04d 12Z",gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900);
  }
   
  else if(type==18)  
  {
    sprintf(tbuf,"%2d day change",indexx+1);
  }
 
  XSetForeground(display,gc,qmap[15]);
  XSetBackground(display,gc,0);
  XSetFont(display,gc,font[2]);  
  yheight=info_font[2]->ascent+2;     
 
  XDrawString(display,pix,gc,30,ymax/10-yheight*2,tbuf,strlen(tbuf));
  XDrawString(display,pix,gc,30,ymax/10-yheight*3,legend[type],strlen(legend[type]));
 
  if(zoom_mode==1)
    strcpy(tbuf,"Mode=zoom");
  else
    strcpy(tbuf,"Mode=edit");
                       
  XDrawString(display,pix,gc,30,ymax/10-yheight,tbuf,strlen(tbuf));
  gm=gmtime(&qtime);
              
  sprintf(tbuf,"Last Update %02d-%02d-%04d %02d:%02dZ",
               gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,gm->tm_hour,gm->tm_min);
 
  XSetForeground(display,gc,qmap[15]);
  XSetBackground(display,gc,0);
  XSetFont(display,gc,font[4]);  
  yheight=info_font[4]->ascent+2;     
  XDrawString(display,pix,gc,30,yheight*2,tbuf,strlen(tbuf));
 
  return;
}