models/johnpye/moody.a4c

REQUIRE "atoms.a4l";

IMPORT "dsg";
IMPORT "johnpye/extpy/extpy";
IMPORT "johnpye/moody";
IMPORT "johnpye/solvernotes";

ADD NOTES IN moody;
	'QRSlv' iterationlimit {50}
END NOTES;
MODEL moody;
	f IS_A positive_factor;
	Re IS_A positive_factor;
	eps_on_D IS_A fraction;

	friction: dsg_friction_factor(
		Re, eps_on_D : INPUT;
		f : OUTPUT
	);
	
	eps IS_A distance;
	D IS_A distance;
	eps_on_D = eps / D;
METHODS
METHOD on_load;
	EXTERNAL solvernotes(SELF);
	EXTERNAL setsomevalue(SELF);
	RUN reset; 
	EXTERNAL checksomevalue(SELF);
	RUN values;
END on_load;
METHOD specify;
	FIX Re, eps, D;
END specify;
METHOD values;
	Re := 20000;
	eps := 0.007 {mm};
	D := 25.4 {mm};
END values;
METHOD plot;
	EXTERNAL moodyplot(SELF);
END plot;
END moody;
1	REQUIRE "atoms.a4l";
2
3	IMPORT "dsg";
4	IMPORT "johnpye/extpy/extpy";
5	IMPORT "johnpye/moody";
6	IMPORT "johnpye/solvernotes";
7
8	ADD NOTES IN moody;
9	'QRSlv' iterationlimit {50}
10	END NOTES;
11	MODEL moody;
12	f IS_A positive_factor;
13	Re IS_A positive_factor;
14	eps_on_D IS_A fraction;
15
16	friction: dsg_friction_factor(
17	Re, eps_on_D : INPUT;
18	f : OUTPUT
19	);
20
21	eps IS_A distance;
22	D IS_A distance;
23	eps_on_D = eps / D;
24	METHODS
25	METHOD on_load;
26	EXTERNAL solvernotes(SELF);
27	EXTERNAL setsomevalue(SELF);
28	RUN reset;
29	EXTERNAL checksomevalue(SELF);
30	RUN values;
31	END on_load;
32	METHOD specify;
33	FIX Re, eps, D;
34	END specify;
35	METHOD values;
36	Re := 20000;
37	eps := 0.007 {mm};
38	D := 25.4 {mm};
39	END values;
40	METHOD plot;
41	EXTERNAL moodyplot(SELF);
42	END plot;
43	END moody;