plot_pcpn_rsel.c

Go to the documentation of this file.

 
#include "prototypes_new.h"
 
void plot_pcpn(Pixmap pixm,int h, int display_flag,char *file,int num,int mnum)
 
{
extern long grey;
extern int dmvalue;
extern int pcpn_time_step;
extern int rsmode;
extern int wfo_all;
extern int wfo_in_use[20];
int m;
extern struct pcp *pcp;
extern struct pcp *spf;
extern int pcp_in_use[500];
extern Font font[10];
extern XFontStruct *info_font[10];
extern Display *display;
extern int npoint[2],zoom;
extern int kscale;
extern GC gc;
extern unsigned long xmap[16];
extern Widget drawing_area;
extern float delim[][16];
extern struct display_set display_set[4];
extern struct hrap_grid *hrap_grid;
int i,j,k;
int color,MY=10;
float xrat,yrat,mult,dmult,cmult;
int nx[4],ny[4];
Dimension width,height;
Arg args[10];
signed long XSIZE,YSIZE;
XPoint point[5];
long xmin,xmax,ymin,ymax;
signed long xcen,ycen,lint1,lint2;
int di[4],dj[4],dy,dx,yheight;
float value;
int i1;
 
XSetFont(display,gc,font[4]);  
yheight=info_font[4]->ascent+2;     
 
if(pcpn_time_step==0 && rsmode != 1) {
         
         i1=1;
     
     if(num==0)
             i1=0;          
 
                 if(pcp_in_use[num+mnum]!=-1 && pcp_in_use[num+mnum-i1]!=-1) {
         
                         read_file(file,num+mnum,spf);
             read_file(file,num+mnum-i1,pcp);
             
                         for(i=0;i<hrap_grid->maxi-1;i++) {
 
                         for(j=0;j<hrap_grid->maxj-1;j++) {
 
                         spf->value[i][j]=(spf->value[i][j] +
                               pcp->value[i][j])/2;
 
             
             
             }
             
             }
             
             
             }
             
             else if(pcp_in_use[num+mnum]==1) {
         
                         read_file(file,num+mnum,spf);
            
             
             }
             
             else if(pcp_in_use[num+mnum-i1]==1) {
         
                         read_file(file,num+mnum-i1,spf);
            
             
             }
             
             
                 
                   }
 
if(pcp_in_use[num]==-1)
                           return;
 
 
read_file(file,num,pcp);
 
dx=1;dy=1;
 
di[0]=0;dj[0]=0;
di[1]=1;dj[1]=0; 
di[2]=1;dj[2]=1;
di[3]=0;dj[3]=1;
 
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];
 
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;
 
cmult=dmult*mult;
lint1=(long)npoint[0]*(long)MY;
lint2=(long)npoint[1]*(long)MY;
 
XSetFillRule(display,gc,WindingRule);
 
for(i=0;i<hrap_grid->maxi-1;i++) {
 
for(j=0;j<hrap_grid->maxj-1;j++) {
 
             if(hrap_grid->owner[i][j]==-1)
                                 continue;
 
             if(wfo_all != 1) {
 
             for(m=0;m<20;m++) {
 
                if(wfo_in_use[m] == -1)
                                     break;
 
                if(hrap_grid->owner[i][j]==wfo_in_use[m])
                                     break;          
 
              }
 
            /* if(wfo_in_use[m]==-1)
                            continue; */
 
             }
         
            if(wfo_all ==1 || wfo_in_use[m]!=-1) {
      
            value=(float)pcp->value[i][j]/100;
            
            if(value<=0.01 && value >=-0.01)
                          continue;
 
            if(rsmode==1) {
 
 
             for(k=0;k<15;k++) {
 
                     if(value >= delim[kscale][k] &&
                        value <  delim[kscale][k+1]) {
 
                                          color=xmap[k];
                                          break;
 
                    }
                     
                  } 
 
              if(k==15)
                     color=xmap[15];
 
              if(value < -99.98)
                     color=xmap[7]; 
              
              if(color==0)
                        continue;
                
          }
          
          
          else {
          
             for(k=0;k<4;k++) {
 
                     if(value >= delim[kscale][k] &&
                        value <  delim[kscale][k+1]) {
 
                                          color=xmap[k];
                                          break;
 
                    }
                     
                  } 
 
              if(k==4)
                     color=xmap[4];
 
              if(value < -99.98)
                     color=xmap[7]; 
              
              if(color==0)
                        continue;
                    
          
          }
        
              if(rsmode != 1 && 
             (pcp_in_use[mnum+num] == 1 ||
         pcp_in_use[mnum+num-i1] == 1) ) {
         
                                  
                          if((spf->value[i][j] * 10 - dmvalue < 
                 hrap_grid->elev[i][j]) && spf->value[i][j] >= 0)
                                                      color=xmap[k+5];
                                                      
                           
                                                  
                           }                           
 
              XSetForeground(display,gc,color);
 
              }
 
              if(zoom <= 1) {
 
              for(k=0;k<4;k++) {               
 
                  nx[k]=xmin+(float)hrap_grid->coord[i+di[k]][j+dj[k]].x*cmult;                
                  ny[k]=ymin+(float)hrap_grid->coord[i+di[k]][j+dj[k]].y*cmult;
 
        }
        }
 
              else {
 
                 for(k=0;k<4;k++) {               
 
                     nx[k]=(xmin+((float)hrap_grid->coord[i+di[k]][j+dj[k]].x*
                           cmult-lint1))*zoom + xcen;
                     ny[k]=(ymin+((float)hrap_grid->coord[i+di[k]][j+dj[k]].y*
                           cmult-lint2))*zoom + ycen;
                   
 
           }
 
           }
 
           
              for(k=0;k<4;k++) {
 
                        point[k].x=(nx[k]+5)/MY;
                        point[k].y=(ny[k]+5)/MY;
 
 
              }
           
              point[4].x=point[0].x;
              point[4].y=point[0].y;
 
          if(wfo_all !=1 && wfo_in_use[m]==-1)         
                 XSetForeground(display,gc,grey);
 
              XFillPolygon(display,pixm,gc,point,5,Convex,CoordModeOrigin);
 
                                          
                 }
 
}
 
return;
 
}