group_edit_tstations.c

Go to the documentation of this file.

#include "prototypes_new.h" 
 
void group_edit_tstations(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 tdata tdata[10];
extern unsigned long cmap[16];
extern Widget drawing_area;
extern struct station tstation[1000];
extern int max_tstations;
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_tstation_location();
void change_custom_file();
void change_tstation_quality();
void reset_tstation_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;
 
printf("start\n");
 
if(pcpn_time_step==0)
                  time_pos=pcpn_time;
 
else if(pcpn_time_step==1)
                  time_pos=4;
 
else if(pcpn_time_step==2)
                  time_pos=5;
 
XSetFont(display,gc,font[3]);
 
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_tstations;i++){
 
 
     if(tdata[pcpn_day].stn[i].tlevel2[time_pos].data ==-99 )
                                          continue;
    
     lat=tstation[i].lat;
     lon=tstation[i].lon;
 
                                    
     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) 
              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;
 
tdata[pcpn_day].stn[isave].tlevel2[time_pos].qual=group_qual;
 
if(pcpn_time_step==0)
                  time_pos=150+pcpn_day*4+pcpn_time;
 
else if(pcpn_time_step==1)
                  time_pos=190+pcpn_day;
 
else if(pcpn_time_step==2)
                  time_pos=200+pcpn_day;
          
pcp_in_use[time_pos]=-1;
 
tdata[pcpn_day].used[pcpn_time]=2;
 
if(pcpn_time_step==1 || pcpn_time_step==2) {
            
for(k=0;k<4;k++) {
                 
             time_pos=150+pcpn_day*4+k;
 
             pcp_in_use[time_pos]=-1;
      
             if(tdata[pcpn_day].used[k]!=0)
                           tdata[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;
 
}