plot_tstations.c

Go to the documentation of this file.

 
#include "prototypes_new.h"
 
 
void plot_tstations(Drawable pix,int type,int h,int display_flag)
 
{
extern struct station tstation[1000];
extern int elevation_filter_value;
extern struct ts ts[20];
extern int tsmax;
extern int find_station_flag;
extern struct tdata tdata[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 int max_tstations;
extern long amap[16];
extern int pcpn_time,pcpn_day;
extern int pcpn_time_step;
int typemax=9;
extern int dflag[10];
extern int qflag[10];
char *typename[10];
int i,xc,yc,text_width,length,m;
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;
float lat,lon;
char tbuf[100],mbuf[100];
int yheight;   
int time_pos;
extern int funct[];
int mm;
 
typename[0]="Verified";
typename[1]="Screened";
typename[2]="Time Distributed";
typename[3]="Manual";
typename[4]="Questionable";
typename[5]="Partial";
typename[6]="Calculated";
typename[7]="Bad";
typename[8]="Missing";
 
if(pcpn_time_step==1)
                  time_pos=4;
          
else if(pcpn_time_step==2)
                  time_pos=5;
          
else
                  
    time_pos=pcpn_time;
 
XSetFont(display,gc,font[4]);
 
yheight=info_font[4]->ascent;
 
if(type==0)
         return;
 
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;
 
for(i=0;i<max_tstations;i++) {
     
     lat=tstation[i].lat;
     lon=tstation[i].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;
 
     }
 
     XSetForeground(display,gc,qmap[15]);
 
     if(tstation[i].elev > 0 && tstation[i].elev < elevation_filter_value)
                  continue;
                                                     
     for(m=0;m<tsmax;m++) {
     
           if(strncmp(&tstation[i].parm[3],ts[m].abr,2)==0 && dflag[m+1] == 1)
                                   break;
 
                           }
                         
     if(m==tsmax)
                continue;
 
     for(m=0;m<9;m++) {
 
             if(m==tdata[pcpn_day].stn[i].tlevel2[time_pos].qual &&
                 qflag[m]==1)
                                               break;
 
               }
   
     if(m==9 && (type==4 || type==5))
              continue;
 
/* locate station in data stream */
    
     if((type==4 || type==5) 
         && tdata[pcpn_day].stn[i].tlevel2[time_pos].data == -99)
                                  continue;
 
     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(type==1)
              strcpy(tbuf,tstation[i].hb5);
 
     else if(type==2) {
 
              strcpy(tbuf,&tstation[i].parm[4]);
              tbuf[1]=0;
 
        }
 
     else if(type==3)
              strcpy(tbuf,tstation[i].name);
 
     else if(type==4) {
     
              if(tdata[pcpn_day].used[time_pos]==0)
                                       continue;
 
              sprintf(mbuf,"%d",tdata[pcpn_day].stn[i].tlevel2[time_pos].data);
 
              strcpy(tbuf,mbuf);
 
        }
        
     else if(type==5) {
     
              if(tstation[i].tip==0)
                            continue;
                
              strcpy(tbuf,tstation[i].hb5);
          
          }
 
     m=tdata[pcpn_day].stn[i].tlevel2[time_pos].qual;
                                             
     XSetForeground(display,gc,amap[m]);  
     length=strlen(tbuf);
     text_width=XTextWidth(info_font[4],tbuf,length); 
 
     xadd=tstation[i].xadd;
     yadd=tstation[i].yadd;
 
     if(xadd < 0) 
                    xc=x1 - text_width;
                        
     else
                    xc=x1+3;
 
     if(yadd < 0)
                    yc=y1;
 
     else
                    yc=y1+yheight;
 
     XDrawString(display,pix,gc,xc,yc,tbuf,length);
 
     if(i==find_station_flag) {
 
              find_station_flag=-1;
 
              XDrawLine(display,pix,gc,xc,yc,xc+text_width,yc);
 
        }
       
       }          
 
XSetFont(display,gc,font[4]);
 
yheight=info_font[4]->ascent;
   
i=0;    
for(m=0;m<typemax;m++) {
 
              if(m==5 || m==6)
                         continue; 
          
              i++;
    
              XSetForeground(display,gc,amap[funct[m]]);   
 
              XDrawLine(display,pix,gc,15,ymin/10+(i+1)*yheight-
                        yheight/2,35,ymin/10+(i+1)*yheight-yheight/2);
 
              XSetForeground(display,gc,qmap[15]);
              XDrawString(display,pix,gc,35,ymin/10+(i+1)*yheight,
                           typename[m],strlen(typename[m]));
 
        }
 
 
return;
 
}