00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034 #include "plumeinput.h"
00035 #include "plumevars.h"
00036 #include <utils/utils.h>
00037
00038 int plumedist( Plume_grid *grid , Plume_options *opt )
00039 {
00040 void itoa(int,char[]);
00041 #ifdef DBG
00042 FILE *fiddbg;
00043 char fndbg[11], ss[3];
00044 #endif
00045 int ii, jj, kk, k1, k2, kin;
00046 double dst1, dst0, rdis;
00047
00048
00049 dst0 = sqrt( sq((grid->ymax-grid->ymin)) + sq((grid->xmax-grid->xmin)) );
00050
00051 if ( opt->fjrd || opt->strt )
00052 {
00053
00054 for ( ii=0 ; ii<grid->lx ; ii++ )
00055 {
00056 for( jj=0 ; jj<grid->ly ; jj++ )
00057 {
00058 grid->dist[ii][jj][0] = (double)ii;
00059 grid->dist[ii][jj][1] = (double)(grid->zy);
00060 grid->dist[ii][jj][2] = grid->xval[ii];
00061 grid->dist[ii][jj][3] = fabs(grid->yval[jj]);
00062 grid->dist[ii][jj][4] = grid->pcent[ii][3];
00063 }
00064 }
00065
00066 }
00067 else
00068 {
00069
00070 for( ii=0 ; ii<grid->lx ; ii++ )
00071 {
00072 jj = 0;
00073
00074 #ifdef DBG
00075 fprintf( stderr ,
00076 "\r PlumeDist ii = %d, lx = %d, ly = %d " ,
00077 ii ,
00078 grid->lx ,
00079 grid->ly );
00080 #endif
00081
00082
00083
00084
00085 dst1 = dst0;
00086 for ( kk=0 ; kk<grid->lc ; kk++ )
00087 {
00088 rdis = sqrt( sq(grid->pcent[kk][0]-grid->xval[ii])
00089 + sq(grid->pcent[kk][1]-grid->yval[jj]) );
00090 if ( rdis < dst1 )
00091 {
00092 dst1 = rdis;
00093 kin = kk;
00094 }
00095 }
00096
00097 grid->dist[ii][jj][0] = grid->pcent[kin][0]/grid->dx;
00098 grid->dist[ii][jj][1] = grid->pcent[kin][1]/grid->dy;
00099 grid->dist[ii][jj][2] = grid->pcent[kin][2];
00100 grid->dist[ii][jj][3] = dst1;
00101 grid->dist[ii][jj][4] = grid->pcent[kin][3];
00102
00103 for ( jj=1 ; jj<grid->ly ; jj++ )
00104 {
00105
00106 k1 = mx( kin-npts, 0 );
00107 k2 = mn( kin+npts+1, grid->lc);
00108
00109 dst1 = dst0;
00110 for ( kk=k1 ; kk<k2 ; kk++ )
00111 {
00112 rdis = sqrt( sq(grid->pcent[kk][0]-grid->xval[ii])
00113 + sq(grid->pcent[kk][1]-grid->yval[jj]) );
00114 if ( rdis < dst1 )
00115 {
00116 dst1 = rdis;
00117 kin = kk;
00118 }
00119 }
00120
00121
00122
00123 grid->dist[ii][jj][0] = grid->pcent[kin][0]/grid->dx;
00124 grid->dist[ii][jj][1] = grid->pcent[kin][1]/grid->dy;
00125 grid->dist[ii][jj][2] = grid->pcent[kin][2];
00126 grid->dist[ii][jj][3] = dst1;
00127 grid->dist[ii][jj][4] = grid->pcent[kin][3];
00128 }
00129 }
00130
00131 #ifdef DBG
00132 fprintf(stderr,"\n");
00133 #endif
00134 }
00135
00136 #ifdef DBG
00137 for ( kk=0 ; kk<5 ; kk++ )
00138 {
00139 strcpy( fndbg, "pdst" );
00140 itoa( kk, ss );
00141 strcat( fndbg, ss );
00142 strcat( fndbg, ".dat" );
00143 fprintf( stderr," Writing dist[ii][jj][%d] to %s \n", kk, fndbg );
00144
00145 fiddbg = fopen(fndbg,"w");
00146 for( ii=0; ii<grid->lx; ii++)
00147 {
00148 for( jj=0; jj<ly; jj++ )
00149 {
00150 fprintf(fiddbg,"%5e ", grid->dist[ii][jj][kk]);
00151 }
00152 fprintf(fiddbg,"\n");
00153 }
00154 fclose(fiddbg);
00155 }
00156 #endif
00157
00158 return 0;
00159
00160 }
00161