! ============================================================
! M09 - Punto de partida
! Barrido parametrico con seleccion automatica de diseno.
! ============================================================
/CLEAR,START
/FILNAME,m09_start,1
/UNITS,SI

! Contrato previsto:
! /INPUT,09_parametric_study,mac

beam_l=1.0
beam_b=0.05
young=210E9
nu=0.30
tip_force=-1000
uy_limit=0.00025
stress_limit=150E6
n_cases=3

! TODO 1: define arrays de alturas y mesh_h = h/2 para cada caso.
! TODO 2: implementa *DO que invoque 09_build, 09_solve y 09_extract.
! TODO 3: reconstruye el modelo en cada iteracion (sin estado residual).
! TODO 4: distingue validation_pass (M07-M08) de feasible (limites de diseno).
! TODO 5: comprueba tendencia: |uy| decrece al aumentar h.
! TODO 6: selecciona min(beam_h) con feasible=1.
! TODO 7: escribe m09_results.csv una sola vez desde el driver.
!
! Casos del barrido base:
! case_id=1 -> h=0.08, mesh_h=0.04, 156 nodos, 50 elementos
! case_id=2 -> h=0.10, mesh_h=0.05, 126 nodos, 40 elementos
! case_id=3 -> h=0.12, mesh_h=0.06, 108 nodos, 34 elementos
!
! Referencias analiticas:
! uy_ref_1=-7.44047619E-4 m ; sigma_ref_1=1.50000E7 Pa ; feasible=0
! uy_ref_2=-3.80952381E-4 m ; sigma_ref_2=9.60000E6 Pa ; feasible=0
! uy_ref_3=-2.20458598E-4 m ; sigma_ref_3=6.86275E6 Pa (x alineado a malla) ; feasible=1
!
! Seleccion esperada: case_id=3, beam_h=0.12 m

FINISH