set dt 5e-8 

buoyancy -g 0 -1 -a 'Pr*Ra'

tol 1e-9
init
loop 200 
  loop 50
    step 20
    flux
  endl
  save
  tecp -f tecp_RBC_Ra1e8_Pr0.71.plt -vars vel p vort psi -s
endl
stop

