set dt 0.005

init

rand

loop 50
  loop 100
    step 10
    int -f integral_vel_minus_frestream_circ_cyl_2d_re100.dat -u '(u-1)^2+v^2' 
    forces 2 forces_circ_cyl_2d_re100.dat
  endl
  save -f save_circ_cyl_2d_re100.dat
  tecp -f tecp_circ_cyl_2d_re100.plt
endl

