% area2rlc.m
% function area2RLC(ar, sl)
% used by syn_fire.m in synthesis
clear i;


	w = 2*pi * 4000.0;
	rho = 1.14e-3;
	mu = 1.86e-4;
	c = 3.53e4;
	eta = 1.4;
	lambda = 5.5e-5;
	epsilon = 0.24;
	mass = 1.5;
	viscous = 1500;
	compilance = 0;
	rdem = sqrt(rho*mu*w/2.0);
	rc2 = rho*c*c;
	Gpc = (eta-1.0)*sqrt(lambda*w/(2.0*epsilon*rho))/rc2;

	for j=1:NING,
 
%   Perimeter = 2.0 * sqrt(pi * A) <----- circle
%            = 4.0 * sqrt(pi * A) <----- ellipse

		peri = 4.0 * sqrt(pi*areatemp(j));

		Rs(j) = peri * rdem/((areatemp(j))^2);
		Ls(j) = rho / areatemp(j);
		Cp(j) = areatemp(j) / rc2;
		Gp(j) = peri * Gpc;

		wall = sqrt(pi)*areatemp(j)*lengtemp(j);        % Maeda */
%               wall = peri*sl->dl[i];  Flanagan */
		Rw(j) = viscous / wall;
		Lw(j) = mass / wall;
%               Cw[i] = wall / compilance;      */
	end
	Rr = (128.0*rho*c)/(9.0*pi*pi*areatemp(NING));
	Lr = 1.5 * ( 8.0*rho / (3.0*pi*sqrt(pi*areatemp(NING))) );