junk.c

Go to the documentation of this file.

 
void area_select(Widget w,XtPointer data)
 
{
 
Cardinal argcount;
Arg args[20];
Widget dialog,rowcol,pbutton;
Widget jwidget;
XmString text;
struct MenuItem items[50];
Widget main_bar;
int slider_size,scrollm,ivalue,increment,page_increment,ier;
FILE *fr,*fp;
char *q,*p, ibuf[100],dbuf[100],ebuf[100];
struct tm *gmtim;
int m,mer;
char preca[100],precb[100],precc[100],precd[100],prece[100],buf[100],databuf[50];
char zpointa[100],zpointb[100];
time_t tget;
struct tm *gm;
int k,kk,mk,num;
time_t sget;
int i;
int l;
 
XDefineCursor(display,window,watch_cursor);
XmUpdateDisplay(top_level);
 
argcount=0;
XtSetArg(args[argcount],XmNwidth,774);argcount++;
XtSetArg(args[argcount],XmNheight,580);argcount++;
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
menu_bar=XmCreateMenuBar(main_window,"Menubar",args,argcount);
text=XmStringCreateLocalized("Reading Map and Precipitation Data");
 
argcount=0;
XtSetArg(args[argcount],XmNwidth, 300);argcount++;
XtSetArg(args[argcount],XmNheight,300);argcount++;
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
XtSetArg(args[argcount],XmNmessageString,text);argcount++;
XtSetArg(args[argcount],XmNdialogStyle,
         XmDIALOG_FULL_APPLICATION_MODAL);argcount++;
dialog=XmCreateMessageDialog(drawing_area,"Map Backgrounds",args,argcount);
XtUnmanageChild(XmMessageBoxGetChild(dialog,XmDIALOG_CANCEL_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(dialog,XmDIALOG_HELP_BUTTON));
XtUnmanageChild(XmMessageBoxGetChild(dialog,XmDIALOG_SEPARATOR));
XtAddCallback(dialog,XmNokCallback,kill_dialog,NULL);
 
XtSetSensitive(XmMessageBoxGetChild(dialog,XmDIALOG_OK_BUTTON),False);
 
argcount=0; 
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
rowcol=XtCreateWidget("rowcolo",xmRowColumnWidgetClass,dialog,args,argcount);
 
argcount=0;
XtSetArg(args[argcount],XmNeditMode,XmMULTI_LINE_EDIT);argcount++;
XtSetArg(args[argcount],XmNscrollBarDisplayPolicy,XmSTATIC);argcount++;
XtSetArg(args[argcount],XmNeditable,False);argcount++;
XtSetArg(args[argcount],XmNcolumns,40);argcount++;
XtSetArg(args[argcount],XmNrows,10);argcount++;
XtSetArg(args[argcount],XmNscrollVertical,True);argcount++;
XtSetArg(args[argcount],XmNscrollHorizontal,True);argcount++;
XtSetArg(args[argcount],XmNcursorPositionVisible,False);argcount++;
XtSetArg(args[argcount],XmNscrollingPolicy,XmAPPLICATION_DEFINED);argcount++;
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
s_text=XmCreateScrolledText(rowcol,"text",args,argcount);
 
XtManageChild(s_text);
 
XtVaGetValues(XtParent(s_text),XmNverticalScrollBar,&scrollbar,NULL);
 
XtManageChild(rowcol);
XtManageChild(dialog);
XmStringFree(text);
 
XmtWaitUntilMapped(dialog);
 
posit=0;
 
clear_drawable(window);
 
legend_save=-1;
tlegend_save=-1;
add_topo_flag=-1;
zoom=1;
 
btim=ouptime(begin_year,begin_month,begin_day,12,0,0);
 
gmtim=gmtime(&btim);
 
emonth=gmtim->tm_mon;
 
tget=btim - 24L*3600L*(long)9;
 
gm=gmtime(&tget);
      
smonth=gm->tm_mon;
 
isom=emonth;
old_isom=emonth;
 
read_maps((char *)data,smonth,emonth);
 
basin_file[0]=0;hrap_file[0]=0;
grid_file[0]=0;
station_list_file[0]=0;rfc[0]=0;
zstation_list_file[0]=0;
map_out[0]=0;freezing_level_directory[100];
hrap_gage_file[0]=0;hrap_zgage_file[0]=0; 
topo_buf[0]=0;
ngrid_file[0]=0;
hrap_grid_mask_file[0]=0;
type[0]=0;
gif_file[0]=0;
gifbuf[0]=0;
station_climo_file[0]=0;
dbase_out[0]=0;zdbase_out[0]=0;
map_out[0]=0;zmap_out[0]=0;
 
fr=fopen((char *)data,"r");
 
for(;;) {
 
        p=fgets(ibuf,80,fr);
 
        if(p==NULL)
                break;
     
        p=strchr(ibuf,'\n');
        if(p!=NULL)
                 *p=0;
        p=pars_line(ibuf,"station_climo_list=",station_climo_file);
        p=pars_line(ibuf,"rfc=",rfc);
        p=pars_line(ibuf,"obs_archive_directory=",obs_archive_file);
        p=pars_line(ibuf,"scratch_directory=",scratch_file);
       
        p=pars_line(ibuf,"station_list=",station_list_file); 
        p=pars_line(ibuf,"stationz_list=",zstation_list_file);
        p=pars_line(ibuf,"hrap_gage_grid=",hrap_gage_file);
        p=pars_line(ibuf,"hrap_zgage_grid=",hrap_zgage_file);
        p=pars_line(ibuf,"basins=",basin_file);
        p=pars_line(ibuf,"hrap=",hrap_file);
        p=pars_line(ibuf,"type=",type);
        p=pars_line(ibuf,"OPOINT_out=",dbase_out);
        p=pars_line(ibuf,"OMAP_out=",map_out);
        p=pars_line(ibuf,"ZPOINT_out=",zdbase_out);
        p=pars_line(ibuf,"OMAZ_out=",zmap_out); 
    if((p=pars_line(ibuf,"pxtemp=",dbuf))!=NULL) {
    
                           pxtemp=atof(dbuf);
                   dmvalue=pxtemp*100*3.28/.55;
                   
    }
                               
        p=pars_line(ibuf,"freezing_level_directory=",freezing_level_directory);
        p=pars_line(ibuf,"topo_scale=",topo_buf);
        p=pars_line(ibuf,"hrap_grid_mask=",hrap_grid_mask_file);
        p=pars_line(ibuf,"gif_directory=",gif_file);
        
        if((p=pars_line(ibuf,"gif_location=",gifbuf))!=NULL)
                                 strcpy(gifbuf,p);
                                 
        if((p=pars_line(ibuf,"type_source=",dbuf))!=NULL) {
        
                           q=strchr(p,',');
                           
                           if(q==NULL)
                                     continue;
                                     
                           *q=0;
                           
                           strcpy(ts[tsmax].abr,p);
                           
                           strcpy(ts[tsmax].name,q+1);
                           
                           tsmax++;
                           
                           }
         
        
      }
      
if(gifbuf[0]!=0) {
 
                 m=0;
                 while(gifbuf[m]!=0) {
                 
                            if(gifbuf[m]==',')
                                       gifbuf[m]=' ';
                             
                            m++;    
                            
                            }      
                                       
                 ier=sscanf(gifbuf,"%f %f %f %f %f %f",
                            &giflat[1],&giflon[1],       
                            &giflat[2],&giflon[2], 
                            &giflat[3],&giflon[3]);
                            
                 maxgif=ier/2;
                            
                  }
                  
if(gif_file[0]==0)
          strcpy(gif_file,scratch_file);                  
 
if(type[0]==0)
            strcpy(type,"2");
 
if(topo_buf[0] != 0)
           topo_scale=atof(topo_buf);
           
fclose(fr);
 
sprintf(pcpn_file,"%s/opoint.%s.",scratch_file,rfc);
sprintf(proc_pcpn_file,"%s/opoint2.%s.",scratch_file,rfc);
sprintf(snow_file,"%s/snopoint.%s.",scratch_file,rfc);
 
sprintf(map_file,"%s/omap.%s.",scratch_file,rfc);
sprintf(grid_file,"%s/ogrid.%s.",scratch_file,rfc);
sprintf(ngrid_file,"%s/ngrid.",scratch_file);
 
sprintf(zpoint1_file,"%s/ozpoint.%s.",scratch_file,rfc);
sprintf(zpoint2_file,"%s/ozpoint2.%s.",scratch_file,rfc);
sprintf(zgrid_file,"%s/ozgrid.%s.",scratch_file,rfc);
sprintf(rsel_file,"%s/omaz.%s.",scratch_file,rfc);
 
if(zstation_list_file[0] != 0) {
 
                    strcpy(dbuf,"reading zstation list\n");
                    posit=posit+strlen(dbuf);
                    XmTextInsert(s_text,posit,dbuf);
 
                    XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
                    XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                                    &increment,&page_increment);
                    XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                                    slider_size,increment,page_increment,True);
 
                    XmUpdateDisplay(top_level);
                    get_zstation_list(zstation_list_file);
 
          }
 
if(station_list_file[0] != 0) {
 
                    strcpy(dbuf,"reading station list\n");
                    posit=posit+strlen(dbuf);
                    XmTextInsert(s_text,posit,dbuf);
 
                    XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
                    XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                                    &increment,&page_increment);
                    XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                                    slider_size,increment,page_increment,True);
 
                    XmUpdateDisplay(top_level);
                    get_station_list(station_list_file,smonth,emonth);
 
          }
 
for(m=0;m<max_stations;m++) {
 
         for(l=0;l<10;l++) {
     
        pdata[l].stn[m].scons=calloc(5,sizeof(short int));
        pdata[l].stn[m].rain=calloc(5,sizeof(struct rain));
        pdata[l].stn[m].frain=calloc(5,sizeof(struct rain));
        pdata[l].stn[m].frzlvl=calloc(5,sizeof(int));
        pdata[l].stn[m].snoflag=calloc(5,sizeof(short int));
        pdata[l].stn[m].sflag=calloc(5,sizeof(short int));
        pdata[l].stn[m].srain=calloc(5,sizeof(struct rain));
     
     }
     
}
 
get_hrap_coord(max_lat,min_lat,cen_lon,smonth,emonth);
 
get_zgage_grid(max_lat,min_lat,cen_lon); 
 
get_hrap_mask(hrap_grid_mask_file);
 
if(isohyets_used==0)
            method=1;
 
items[0].label=NULL;
items[0].callback=NULL;
items[0].callback_data=NULL;
 
main_bar=BuildPulldownPushMenu(menu_bar,"File",'M',items);
jwidget=AttachToCascade(menu_bar,"File",main_bar);
 
argcount=0;
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
 
pbutton=XmCreatePushButton(main_bar,"Send to Database",args,argcount);  
XtAddCallback(pbutton,XmNactivateCallback,send_dbase,NULL);
XtManageChild(pbutton);
 
argcount=0;
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
 
pbutton=XmCreatePushButton(main_bar,"Make GIF",args,argcount);  
XtAddCallback(pbutton,XmNactivateCallback,make_gif,NULL);
XtManageChild(pbutton);
 
argcount=0;
 
XtSetArg(args[argcount],XmNbackground,mmap[1]);argcount++;
XtSetArg(args[argcount],XmNforeground,mmap[0]);argcount++;
 
pbutton=XmCreatePushButton(main_bar,"Quit",args,argcount);  
XtAddCallback(pbutton,XmNactivateCallback,quit_all,NULL);
XtManageChild(pbutton);
 
strcpy(dbuf,"reading basin data\n");
posit=posit+strlen(dbuf);
XmTextInsert(s_text,posit,dbuf);
 
XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                    &increment,&page_increment);
XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                    slider_size,increment,page_increment,True);
        
XmUpdateDisplay(top_level);
 
get_basin_data(basin_file,NULL,hrap_file);
 
strcpy(dbuf,"reading Precipitation data\n");
posit=posit+strlen(dbuf);
XmTextInsert(s_text,posit,dbuf);
 
XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                    &increment,&page_increment);
XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                    slider_size,increment,page_increment,True);
        
XmUpdateDisplay(top_level);
 
ier=chdir(scratch_file);
 
if(ier==0) {
 
                             
           }
 
else {
 
     printf("could not change directory to %s\n",scratch_file);
     
     exit(1);
     
     }
 
for(m=0;m<maxdays;m++) {
 
     tget=btim - 24L*3600L*(long)m;
 
     read_archived_data(scratch_file,"obs",rfc,tget,obs_archive_file);
 
     read_archived_data(scratch_file,"nexrad","",tget,obs_archive_file);
 
     for(k=0;k<5;k++) {
 
          if(k < 2)
                 sget=tget-84600L;
 
          else
                 sget=tget;
          
          gmtim=gmtime(&sget);
 
          kk=3-k;
 
          sprintf(dbuf,"%s%02d-%02d-%02d.%s",ngrid_file,
                  gm->tm_mon+1,gm->tm_mday,gm->tm_year,timefile[k]);
               
          if(k < 4)
              num=50+m*4+kk;
 
          else
              num=90+m;
 
          mer=read_qpf_grids(num,dbuf);
 
           
  
    }
 
     pdata[m].stddev=3.0;
 
     tget=btim - 24L*3600L*(long)m;
 
     gm=gmtime(&tget);
      
     sprintf(preca,"%s%02d-%02d-%02d",
             pcpn_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(precb,"%s%02d-%02d-%02d",
             proc_pcpn_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(precc,"%sdev.%02d-%02d-%02d",
             proc_pcpn_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(precd,"%sbad.%02d-%02d-%02d",
             proc_pcpn_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(prece,"%s%02d-%02d-%02d",
             snow_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(zpointa,"%s%02d-%02d-%02d",
             zpoint1_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
     sprintf(zpointb,"%s%02d-%02d-%02d",
             zpoint2_file,gm->tm_mon+1,gm->tm_mday,gm->tm_year);
 
     sprintf(databuf,"%02d-%02d-%02d",gm->tm_mon+1,gm->tm_mday,gm->tm_year);
 
     sprintf(buf,"cp %s/opoint.%s.%02d-%02d-%02d %s",obs_archive_file,rfc,
             gm->tm_mon+1,gm->tm_mday,gm->tm_year,scratch_file);
 
     system(buf);    
     
     sprintf(buf,"cp %s/snopoint.%s.%02d-%02d-%02d %s",obs_archive_file,rfc,
             gm->tm_mon+1,gm->tm_mday,gm->tm_year,scratch_file);
 
     system(buf);
     
     sprintf(buf,"cp %s/ozpoint.%s.%02d-%02d-%02d %s",obs_archive_file,rfc,
             gm->tm_mon+1,gm->tm_mday,gm->tm_year,scratch_file);
 
     system(buf);
     
     read_snow(prece,m); 
          
     mer=read_precip_a(preca,tget,m); 
  
     read_bad_values(precd,m);
 
     ier=read_precip_b(precb,tget,m); 
 
     read_zlevel_a(zpointa,tget,m); 
 
     read_zlevel_b(zpointb,tget,m); 
  
     fp=fopen(precc,"r");
     if(fp!=NULL) {
 
                  p=fgets(ibuf,80,fp);
 
                  pdata[m].stddev=atof(ibuf);
 
                  if(pdata[m].stddev != 1.0 && 
                     pdata[m].stddev != 3.0 && 
                     pdata[m].stddev != 5.0)
                              pdata[m].stddev=3.0;
 
                fclose(fp);
 
        }
                   
     if(ier==1)
       sprintf(dbuf,"%s - level 2 data\n",databuf);
  
     else if(ier==-1 && mer==1)
       sprintf(dbuf,"%s - level 1 data\n",databuf);
 
     else if(ier==-2 && mer==1)
       sprintf(dbuf,"%s - level 1 data overwrite\n",databuf);
 
     else
     sprintf(dbuf,"%s - no data\n",databuf);
     posit=posit+strlen(dbuf); 
     XmTextInsert(s_text,posit,dbuf);
 
     XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
     XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                    &increment,&page_increment);
     XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                    slider_size,increment,page_increment,True);
        
     XmUpdateDisplay(top_level);
 
     estimate_daily_stations(m);
     
     estimate_partial_stations(m);
    
     quality_control_stations(m);
 
     check_consistency(m);
 
     restore_bad_values(m);
 
     /* read in grids and make maps if level 2 found */
 
     if(ier==1) {
 
        for(k=0;k<5;k++) {
 
          if(k < 2)
                 sget=tget-84600L;
 
          else
                 sget=tget;
          
          gmtim=gmtime(&sget);
 
          kk=3-k;
 
          sprintf(dbuf,"%s%02d-%02d-%02d.%s",ngrid_file,
                  gm->tm_mon+1,gm->tm_mday,gm->tm_year,timefile[k]);
               
          if(k < 4)
              num=50+m*4+kk;
 
          else
              num=90+m;
 
          mer=read_qpf_grids(num,dbuf);
 
          sprintf(dbuf,"%s%02d-%02d-%02d.%s",grid_file,
                  gm->tm_mon+1,gm->tm_mday,gm->tm_year,timefile[k]);
               
 
          if(k < 4)
              num=m*4+kk;
 
          else
              num=40+m;
 
          mer=read_qpf_grids(num,dbuf);
  
          if(mer==-1) {
 
              if(k==4)
                      mk=1;
 
              else
                      mk=0;
 
              if(pdata[m].used[k]!=0) {
              
              render_pcp(m,k,mk);
 
              write_qpf_grids(dbuf);
 
                 }
 
        }
     
        make_map(num);
  
        
    }
 
      }
   
for(k=0;k<4;k++) {
 
          if(zdata[m].level[k] != 2)
                   continue;
            
          if(k < 2)
                 sget=tget-84600L;
 
          else
                 sget=tget;
          
          gmtim=gmtime(&sget);
 
          kk=3-k;
 
          sprintf(dbuf,"%s%02d-%02d-%02d.%s",zgrid_file,
                  gm->tm_mon+1,gm->tm_mday,gm->tm_year,ztimefile[k]);
               
          num=100+m*4+kk;
 
          mer=read_qpf_grids(num,dbuf);
 
          if(mer==-1)
               render_z(pcpn_day,k,0);
 
          make_rsel(num);            
        
    }
 
}
 
get_bad_snotel(); 
 
calculate_zlevel();
 
for(m=0;m<10;m++)
         get_zlevel(m);
 
/* edit menus */
 
items[0].label="Precipitation";
items[0].callback=change_edit_stations;
items[0].callback_data=NULL;
items[1].label="Freezing level";
items[1].callback=change_zedit_stations;
items[1].callback_data=NULL;
items[2].label=NULL;
 
main_bar=BuildPulldownPushMenu(menu_bar,"Edit",'M',items);
jwidget=AttachToCascade(menu_bar,"Edit",main_bar);
 
main_bar=BuildPulldownToggleMenu(menu_bar,"Backgrounds",'M',ktems);
jwidget=AttachToCascade(menu_bar,"Backgrounds",main_bar);
 
items[0].label="Find station";
items[0].callback=other_pcpn_options;
items[0].callback_data=(XtPointer)1;
items[1].label=NULL;
 
main_bar=BuildPulldownPushMenu(menu_bar,"Options",'M',items);
jwidget=AttachToCascade(menu_bar,"Options",main_bar);
 
/* products menu */
 
XtManageChild(menu_bar);
 
XmMainWindowSetAreas(main_window,menu_bar,NULL,NULL,NULL,drawing_area);
 
strcpy(dbuf,"Done!!\n");
posit=posit+strlen(dbuf);
XmTextInsert(s_text,posit,dbuf);
 
XtVaGetValues(scrollbar,XmNmaximum,&scrollm,NULL);
XmScrollBarGetValues(scrollbar,&ivalue,&slider_size,
                    &increment,&page_increment);
XmScrollBarSetValues(scrollbar,scrollm-slider_size,
                    slider_size,increment,page_increment,True);
 
get_legend();
      
XmUpdateDisplay(top_level);
 
XtSetSensitive(XmMessageBoxGetChild(dialog,XmDIALOG_OK_BUTTON),True);
 
}