Code covered by the BSD License
- BipolarCoordinates(Xi,Eta...a is radii that defies height of the domain
- EllipticCylinderCoordinat...a is radius of ellips
row is aspect ratio of the major to minor axes of ellipse
R is aspect ratio of the major to minor axes of ellipse
- PolarCoordinates(Xi,Eta,t...theta0, theta1 are two angles.
- plotgrid(X,Y)plotgrid: To plot structured grid.
View all files
04 Mar 2013
To demonstrate grid generation using analytic coordinates systems.
function [x y] = PolarCoordinates(Xi,Eta,theta0,theta1,r0,r1)
% theta0, theta1 are two angles.
% r0, r1 are two radii
% Check for the input arguments and set default values
if nargin == 2
theta0 = 0. ;
theta1 = -pi/4 ;
r0 = 1. ;
r1 = 2. ;
theta = (theta0-theta1)*Xi ;
r = r0+(r1-r0)*Eta ;
x = r*cos(theta) ;
y = r*sin(theta) ;