How to plot a data in spherical coordinates?

Let me explain, what I'm actually solving.

I'm solving Diffusion-Reaction equation(2D) in spherical coordinates(axis-symmetric case).

I'm solving this equation numerically. So at the end I'm getting , r, θ.

So I want to plot this function in r and θ plane instead of palne

I'm attaching data.mat file to this question

Can anyone help in this regard?

Thanks in advance.

Answer 1

[David Hill]

[x,y,z]=sph2cart(th,conc,r);
surf(x,y,z);

[Jagadeesh Korukonda]

No it is giving different result, conc is not elavation angle. conc is a function of r and theta.

ψ(r, θ) is not ϕ.

Actually ψ is a function of ψ(r, θ, ϕ). I solved my solution for axis-symmetry case.

Let me explain, what I'm actually solving.

I'm solving Diffusion-Reaction equation(2D) in spherical coordinates(axis-symmetric case).

I'm solving this equation numerically. So at the end I'm getting , r, θ.

So I want to plot this function in r and θ plane instead of palne

[Torsten]

If you solve the equation in r and theta, you also get the solution C in r and theta. So plotting in the r-theta plane would simply mean

[R,THETA] = ndgrid(r,theta);
surf(R,THETA,C)

where C(i,j) = C(r(i),theta(j)).

So you plot over a rectangle with x-axis r and y-axis theta (0<=r<=R, 0<=theta<=2*pi).

If you want to plot in the x-y-plane (thus over the circle with radius R), you must convert to x and y coordinates via x = r*cos(theta), y= r*sin(theta).

[Jagadeesh Korukonda]

Thanks @Torsten