time_series.c

Go to the documentation of this file.

#include "prototypes.h" 
 
void time_series(int win_x,int win_y,unsigned int win_button,char *disp)
 
{
extern Cursor watch_cursor;
extern Cursor dotbox_cursor;
extern Widget top_level;
extern Window window;
extern Widget drawing_area;
extern Display *display;
int isave;
extern struct pdata pdata[10];
extern int qflag[10];
extern struct station station[3000];
extern int max_stations;
extern int npoint[2],zoom;
extern struct display_set display_set[4];
extern struct dval dval;
extern int dflag[10];
extern int pcpn_day,pcpn_time;
extern int pcpn_time_step;
Widget rowcol,rowcol3,pbutton,rowcol1,dialog,sep;
float conv=.0174;
int MY=10;
Arg args[10];
Cardinal argcount;
int i;
double testdist,maxdist;
int x,y;
XmString t;
Dimension width,height;
signed long XSIZE,YSIZE;
float xrat,yrat,mult;
int xmin,xmax,ymin,ymax,display_flag,xcen,ycen,h;
long lint1,lint2;
float dmult,lat,lon;
char cmd[100],pc[10],buf[100],muf[10];
long r,s;
int time_pos,naflag,m;
char tbuf[100];
 
if(pcpn_time_step==0)
                  time_pos=pcpn_time;
 
else 
                  time_pos=4;
 
/* need to set flag to plot */
 
/* need to set resource in widget also */
 
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;
 
display_flag=0;
h=0;
 
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;
 
x=win_x;
y=win_y;
 
if(zoom==1) {
 
     win_x=x*MY/mult;
     win_y=y*MY/mult;
 
   }
 
else         {
 
 
     win_x=((x*MY-xcen)/zoom + lint1)/mult;
     win_y=((y*MY-ycen)/zoom + lint2)/mult;
 
   }
 
 
if(zoom==1) {
 
             x=(win_x*mult)/MY;
             y=(win_y*mult)/MY;
 
             }
 
else         {
 
             x=(((win_x*mult)-lint1)*zoom+xcen)/MY;
             y=(((win_y*mult)-lint2)*zoom+ycen)/MY;
 
             }
 
isave=-1;
maxdist=9999;
 
for(i=0;i<max_stations;i++){
                    
     lat=station[i].lat;
     lon=station[i].lon;
 
     if(strncmp(&station[i].parm[4],"P",1)==0 && dflag[1] != 1)
                                                  continue;
 
     if(strncmp(&station[i].parm[4],"M",1)==0 && dflag[2] != 1)
                                                  continue;
 
     if(strncmp(&station[i].parm[4],"R",1)==0 && dflag[3] != 1)
                                                  continue;
 
     if(strncmp(&station[i].parm[4],"G",1)==0 && dflag[4] != 1)
                                                  continue;
 
     if(strncmp(&station[i].parm[4],"Z",1)==0 && dflag[5] != 1)
                                                  continue; 
 
     for(m=0;m<8;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==8) 
              continue;
 
     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;
          
     testdist= pow((double)(win_x-r),2) + pow((double)(win_y-s),2);
     testdist= pow(testdist,.5);
 
     if(testdist < maxdist) {
 
                        isave=i;
                        maxdist=testdist;
 
               }
                            
                          } 
               
if(isave==-1)
            return;
 
strcpy(cmd,"rsh wasatch -l oper /usr/local/bin/subdisplay ");
strcat(cmd,station[isave].hb5);
strcat(cmd," ");
strcpy(pc,station[isave].parm);
 
pc[3]='R';
if(pc[4]!='Z') {
 
           pc[1]='C';
           pc[2]='I';
 
         }
 
else
          pc[2]='D';
       
strcat(cmd,pc);
strcat(cmd," ");
 
strcat(cmd,"15 ");
sprintf(tbuf,"-display %s",disp);
strcat(cmd,tbuf);
 
XDefineCursor(display,window,watch_cursor);
XmUpdateDisplay(top_level);
 
XUndefineCursor(display,window);
XDefineCursor(display,window,dotbox_cursor);
system(cmd);
 
}