|
From: Paulo C. <co...@as...> - 2004-10-15 01:25:13
|
Hi
thanks for the reply I wasn't that explicit because I wasn't sure if the
list was alive :). Here is little piece of my code.
PMap(i,j).getDensity() is an angle (polarization angle). I get the shape
that I should; except for a line like this:
| *
| *
| *
|*
_ _ _ _
in front of the vector field. Thanks!!!
Paulo
ps. Is it politically correct to include an image as an attachment?
so I can include my problem?
nx=33;ny=33;nz=33;
for (j = 0; j<ny; j++) {
for (i = 0; i<nx; i++) {
x = (i-nx/2)*dx;
y = (j-ny/2)*dy;
cgrid2.xg[i][j] = (i-nx/2);
cgrid2.yg[i][j] = (j-ny/2);
printf("y=%3.1f z=%3.1f Qmap=%3.8f Umap=%3.8f\n",y,z,
QMap(i,j).getDensity(),UMap(i,j).getDensity());
fpol =
sqrt(QMap(i,j).getDensity()*QMap(i,j).getDensity()+
UMap(i,j).getDensity()*UMap(i,j).getDensity())/
IMap(i,j).getDensity();
u[i][j]= fpol*cos(PMap(i,j).getDensity());
v[i][j]= fpol*sin(PMap(i,j).getDensity());
pp = sqrt(u[i][j]*u[i][j] + v[i][j]*v[i][j]);
}
}
pls->env(xmin, xmax, ymin, ymax, 0, 0);
pls->lab("(x)", "(y)", "#frPLplot B field");
pls->col0(2);
pls->vect(u,v,nx,ny,1.0,plstream::tr2,(void *)&cgrid2);
pls->col0(1);
Andrew Ross wrote:
> On Thu, Oct 14, 2004 at 04:20:53PM -0500, Paulo Cortes wrote:
>
>>Hi
>>
>>I am using plplot to plot a vector field, but
>>am having problems. I get the vector field, but I also
>>get a line that goes from the lower left corner to the
>>upper right corner of the plot, I can't get rid of it.
>>any ideas?
>
>
> Could it possibly be your array size is not quite right? Otherwise it
> is hard to tell from your description. Can so provide us with a simple
> example code?
>
> Andrew
>
>
> -------------------------------------------------------
> This SF.net email is sponsored by: IT Product Guide on ITManagersJournal
> Use IT products in your business? Tell us what you think of them. Give us
> Your Opinions, Get Free ThinkGeek Gift Certificates! Click to find out more
> http://productguide.itmanagersjournal.com/guidepromo.tmpl
> _______________________________________________
> Plplot-general mailing list
> Plp...@li...
> https://lists.sourceforge.net/lists/listinfo/plplot-general
>
|