jnk.c

Go to the documentation of this file.

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;
 
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,"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/zpoint.%s.",scratch_file,rfc);
sprintf(zpoint2_file,"%s/zpoint2.%s.",scratch_file,rfc);
sprintf(zgrid_file,"%s/zgrid.%s.",scratch_file,rfc);
 
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);
 
          }
 
/*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);
 
          }
 
*/
max_zstations=0;
 
get_hrap_coord(max_lat,min_lat,cen_lon,smonth,emonth);
printf("pcp %d\n",pcp_in_use[103]);
 
/*get_zgage_grid(max_lat,min_lat,cen_lon); */
 
get_hrap_mask(hrap_grid_mask_file);
printf("2pcp %d\n",pcp_in_use[103]);
 
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/zpoint.%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); 
 
/*
     zier=-1;
 
     zier=read_zlevel_a(zpointa,tget,m); 
 
     zier=read_zlevel_b(zpointb,tget,m); 
printf("pcp %d\n",pcp_in_use[103]);
  
  */
     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);
 
/*     get_frzlvl_data(m); */
     
     estimate_daily_stations(m);
     
     estimate_partial_stations(m);
    
     quality_control_stations(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);
  
        
    }
 
      }
     printf("pcp %d\n",pcp_in_use[103]); 
 if(zier==1) {
 
        for(k=0;k<4;k++) {
 
          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);            
        
    }
 
      }
 
   
   }
 
printf("pcp %d\n",pcp_in_use[103]);
/*get_bad_snotel(); */
 
if(menu_flag==1) 
       XtDestroyWidget(save_widget);
 
/* 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);
 
menu_flag=1;
 
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);
 
}
void kill_dialog(Widget w,XtPointer data,XmDrawingAreaCallbackStruct *call_data)
{
 
first_through=1;
XtDestroyWidget(w);
send_expose();