plot_stations.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
 
void plot_stations(Drawable pix,int type,int h,int display_flag)
 
{
extern int find_station_flag;
extern struct pdata pdata[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 int qflag[10];
extern long amap[16];
extern int pcpn_time,pcpn_day;
extern int pcpn_time_step;
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];
int yheight;   
struct tm *gm;
int time_pos;
time_t ltime;
char ast[300];
 
for(i=0;i<255;i++)
           ast[i]=i;
           ast[255]=0;
 
if(pcpn_time_step==0)
              time_pos=pcpn_time;
 
else
              time_pos=4;
 
XSetFont(display,gc,font[4]);
 
yheight=info_font[4]->ascent;
 XDrawString(display,pix,gc,0,50,ast,255);
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_stations;i++) {
     
     lat=station[i].lat;
     lon=station[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(strncmp(&station[i].parm[3],"RP",2)==0 && dflag[1] != 1)
                        continue;
 
     if(strncmp(&station[i].parm[3],"RM",2)==0 && dflag[2] != 1)
                        continue;
 
     if(strncmp(&station[i].parm[3],"RR",2)==0 && dflag[3] != 1)
                        continue;
 
     if(strncmp(&station[i].parm[3],"RG",2)==0 && dflag[4] != 1)
                        continue;
 
     if(strncmp(&station[i].parm[3],"RZ",2)==0 && dflag[5] != 1)
                        continue;
 
/* locate station in data stream */
 
 
     if(type==4 && pdata[pcpn_day].stn[i].frain[time_pos].data == -1)
                                  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,station[i].hb5);
 
     else if(type==2) {
 
              strcpy(tbuf,&station[i].parm[3]);
              tbuf[2]=0;
 
        }
 
     else if(type==3)
              strcpy(tbuf,station[i].name);
 
     else if(type==4) {
 
              if(pdata[pcpn_day].used[time_pos]==0)
                                       continue;
 
              if(pdata[pcpn_day].stn[i].frain[time_pos].data == -2)
                                      continue;
 
              for(m=0;m<5;m++) {
 
                  if((m==pdata[pcpn_day].stn[i].frain[time_pos].qual &&
                      qflag[m]==1) || 
                     (m==(pdata[pcpn_day].stn[i].frain[time_pos].qual-100) &&
                      qflag[m]==1) )
                                               break;
 
               }
 
              if(m==5) 
                     continue;
 
              sprintf(tbuf,"%5.2f",pdata[pcpn_day].stn[i].frain[time_pos].data);
              XSetForeground(display,gc,amap[m]);   
 
        }
       
     length=strlen(tbuf);
     text_width=XTextWidth(info_font[4],tbuf,length); 
 
     xadd=station[i].xadd;
     yadd=station[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);
 
        }
 
       }          
 
if(type==4) {
 
         if(pcpn_time < 2)
                        ltime=pdata[pcpn_day].data_time - 24L*3600L;
 
         else
                        ltime=pdata[pcpn_day].data_time;
             
         gm=gmtime(&ltime);
 
         sprintf(tbuf,"%02d-%02d-%02d",gm->tm_mon+1,gm->tm_mday,gm->tm_year);
 
         if(pcpn_time_step==0) {
 
             if(pcpn_time==3)
                    strcat(tbuf," 06z-12z");
             else if(pcpn_time==2)
                    strcat(tbuf," 00z-06z");
             else if(pcpn_time==1)
                    strcat(tbuf," 18z-00z");
             else if(pcpn_time==0)
                    strcat(tbuf," 12z-18z");
       }
 
         else
             strcat(tbuf," ending at 12z");
 
         XSetForeground(display,gc,qmap[15]);
         XSetBackground(display,gc,0);
 
         if(pdata[pcpn_day].used[time_pos]==3) 
                                     strcat(tbuf," Modified - Saved");
 
         else if(pdata[pcpn_day].used[time_pos]==2) 
                                     strcat(tbuf," Modified - Not Saved");
 
         else if(pdata[pcpn_day].used[time_pos]==1) 
                                     strcat(tbuf," Not Modified");
 
         else 
                                     strcat(tbuf," No Data");
 
         XDrawString(display,pix,gc,30,ymax/10-30,tbuf,strlen(tbuf));
 
 
       }
 
return;
}