time_series.c

Go to the documentation of this file.

#include "prototypes_new.h" 
 
void time_series(int win_x,int win_y,unsigned int win_button)
 
{
extern struct ts ts[20];
extern int tsmax;
extern char scratch_file[100];
extern int gage_char[2];
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;
float conv=.0174;
int MY=10;
Arg args[10];
int i;
double testdist,maxdist;
int x,y;
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];
long r,s;
int time_pos,m;
char tbuf[50];
time_t ctime;
struct tm *gm;
 
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(station[i].tip==0 && gage_char[0]==-1)
                                    continue;
 
     if(station[i].tip==1 && gage_char[1]==-1)
                                    continue;
                         
     for(m=0;m<tsmax;m++) {
     
           if(strncmp(&station[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==pdata[pcpn_day].stn[i].frain[time_pos].qual &&
                 qflag[m]==1) 
                                               break;
 
               }
       
     if(m==9) 
              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(pc,station[isave].parm);
 
pc[3]='R';
if(pc[4]!='Z') {
 
           pc[1]='C';
           pc[2]='I';
 
         }
 
else
          pc[2]='D';
 
ctime=pdata[0].data_time;  
 
gm=gmtime(&ctime);
 
/*display_data(station[isave].hb5,pc,gm->tm_mon+1,gm->tm_mday,gm->tm_year+1900,12,15,scratch_file);*/
 
}