123 |
END FOR; |
END FOR; |
124 |
RUN values; |
RUN values; |
125 |
RUN solvernotes; |
RUN solvernotes; |
126 |
RUN solve; (* wish I knew why it needs three goes *) |
RUN solve; |
127 |
RUN ode_init; |
RUN ode_init; |
128 |
qdot_s := 10 {W/m}; |
qdot_s := 10 {W/m}; |
129 |
RUN solve; (* after fixing the states and freeing the derivatives *) |
RUN solve; (* after fixing the states and freeing the derivatives *) |
134 |
t.ode_type := -1; |
t.ode_type := -1; |
135 |
t := 0 {s}; |
t := 0 {s}; |
136 |
FOR i IN [2..n] DO |
FOR i IN [2..n] DO |
137 |
FREE drho_dt[i]; FIX node[i].rho; |
FREE drho_dt[i]; FREE node[i].rho; |
138 |
FREE dmdot_dt[i]; FIX node[i].mdot; |
FREE dmdot_dt[i]; FREE node[i].mdot; |
139 |
FREE drhou_dt[i]; FIX rhou[i]; |
FREE drhou_dt[i]; FREE rhou[i]; |
140 |
FREE dTw_dt[i]; FIX T_w[i]; |
FREE dTw_dt[i]; FREE T_w[i]; |
141 |
END FOR; |
END FOR; |
142 |
|
|
143 |
FOR i IN [2..n] DO |
FOR i IN [2..n] DO |
178 |
|
|
179 |
METHOD solvernotes; |
METHOD solvernotes; |
180 |
EXTERNAL solvernotes(SELF); |
EXTERNAL solvernotes(SELF); |
181 |
END setup_solver; |
END solvernotes; |
182 |
|
|
183 |
METHOD solve; |
METHOD solve; |
184 |
EXTERNAL solve(SELF); |
EXTERNAL solve(SELF); |