[ayoung94]

I am having trouble generating point geometry based around a curve equation. I wish to create several points to describe the curve, and want to do this via a for loop, as part of a project on automating geometry and mesh for given input parameters.

However, my current code returns several errors, including the For loop not being closed, and the 'integralExpr' function not being a valid command.

ic_geo_new_family GEOM
ic_empty_tetin
set d_throat 0.1288
set d_ef 0.1524
set n_profiles 6
set x_duct 0.4572
set l_curve 0.481613
for {set i 0} {$i<$n_profiles} {set i [expr $i+1]} {
set x_l [expr $i / ($n_profiles - 1)]
set x_curve [expr $l_curve * $x_l]
set x_global [integralExpr 0 $x_curve 10 sqrt(1 + (0.0225 * 3.14159 * 3.14159 * sin(3.14159 * $x_l) * sin(3.14159 * $x_l)))]
set y_global 0
set z_global [expr -0.15 * $x_duct * (1 - cos(3.14159 * $x_l))]
ic_point {} GEOM pnt.1 $x_global,$y_global,$z_global
}

Any suggestions and corrections are most welcome. Thanks in advance.