;file: slopedy.dat rest slope0.sav def convert c_p = sqrt((b_mod + (4.0 * sh_mod / 3.0)) / m_dens) c_s = sqrt(sh_mod / m_dens) norm_str = -2.0 * m_dens * c_p shear_str = 2.0*(-2.0 * m_dens * c_s) end set m_dens=0.0025 b_mod=16667 sh_mod=10000 convert def blvel whilestepping if time > env_time wave = 0.0 else wave = ampl * sin(2.0*pi*freq*time) endif end set freq =10 ampl=1.0 env_time=0.1 blvel boundary ffield boundary xvisc ff_bulk=16667.0 ff_shear=10000.0 ff_density=0.0025 & range -0.1,22.1 -5.1,-4.9 boundary yvisc ff_bulk=16667.0 ff_shear=10000.0 ff_density=0.0025 & range -0.1,22.1 -5.1,-4.9 boundary stress 0.0,0.0,norm_str history=wave range -0.1,22.1 -5.1,-4.9 boundary stress 0.0,shear_str,0.0 history=wave range -0.1,22.1 -5.1,-4.9 reset hist time disp history wave history xvelocity 10.0,-5.0 history xvelocity 8.0,6.0 history xvelocity 18.0,-2.0 history yvelocity 10.0,-5.0 history yvelocity 22.0,10.0 damping 0.0010 10.0 save slopedy.sav cycle time 0.2 save slope02.sav cycle time 0.2 save slope03.sav cycle time 0.3 save slope04.sav