group_edit_stations.c

Go to the documentation of this file.

#include "prototypes_new.h" 
 
void group_edit_stations(int win_x,int win_y,unsigned int win_button)
 
{
extern int pcp_in_use[500];
extern Widget diswidget[3];
extern float filter_value;
extern struct ts ts[20];
extern int tsmax;
extern Widget edit_dialog;
extern int isom;
extern int gage_char[2];
extern int method;
extern int qflag[10];
extern struct pdata pdata[10];
extern unsigned long cmap[16];
extern Widget drawing_area;
extern struct station station[3000];
extern int max_stations;
extern Font font[10];
extern int npoint[2],zoom;
extern struct display_set display_set[4];
extern Display *display;
extern GC gc;
extern int mmap[16];
extern void kill_widget();
extern struct dval dval;
extern int dflag[10];
void change_station_location();
void change_custom_file();
void change_station_quality();
void reset_station_quality();
void cancel_edit();
void set_snotel();
extern int pcpn_day,pcpn_time;
extern int pcpn_time_step;
Widget rowcol,rowcol3,pbutton,rowcol1,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 *st[10],buf[100],muf[10];
long r,s;
int time_pos,naflag,m;
XmString help_string;
extern int group_qual;
int isave;
int k;
 
 
if(pcpn_time_step==0)
                  time_pos=pcpn_time;
 
else 
                  time_pos=4;
 
XSetFont(display,gc,font[3]);
 
/* need to set flag to plot */
 
/* need to set resource in widget also */
 
XSetFunction(display,gc,GXcopy);
 
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++){
 
 
     if(pdata[pcpn_day].stn[i].frain[time_pos].data < 0)
                                          continue;
                                          
     if(pdata[pcpn_day].stn[i].frain[time_pos].data < filter_value)
                                          continue;
                    
     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;
 
 
pdata[pcpn_day].stn[isave].frain[time_pos].qual=group_qual;
 
if(group_qual==1 &&
   time_pos==4 && pdata[pcpn_day].stn[isave].sflag[time_pos]==1) {
      
                  pdata[pcpn_day].stn[isave].frain[time_pos].data=
                  pdata[pcpn_day].stn[isave].rain[time_pos].data;
                  
                  pdata[pcpn_day].stn[isave].sflag[time_pos]=-1;
                  
                            }
 
if(time_pos==4 && ( group_qual==1 || group_qual==0 || group_qual==8)) {
 
            for(k=0;k<4;k++) 
                   pdata[pcpn_day].stn[isave].frain[k].qual=group_qual;
 
      }
 
    /* 6 hour data set bad set 24 hour bad too */
 
if(time_pos != 4 && group_qual==1 && 
       pdata[pcpn_day].stn[isave].frain[4].qual != 5 &&
       pdata[pcpn_day].stn[isave].frain[4].qual != 4)
                 pdata[pcpn_day].stn[isave].frain[4].qual=group_qual;
   
         
for(k=0;k<5;k++) {
      
             if(k < 4)
                             time_pos=pcpn_day*4 +k;
 
             else
                             time_pos=40+pcpn_day;
 
             pcp_in_use[time_pos]=-1;
      
 
             if(pdata[pcpn_day].used[k]!=0)
                           pdata[pcpn_day].used[k]=2;
   }
    
XtSetSensitive(diswidget[1],False);
XtSetSensitive(diswidget[2],False);
XtSetArg(args[0],XmNset,False);        
XtSetValues(diswidget[1],args,1);
XtSetArg(args[0],XmNset,False);        
XtSetValues(diswidget[2],args,1);
 
write_screen();
 
return;
 
}