1d_fd_simple_shear

continuum model for granular flows with pore-pressure dynamics
git clone git://src.adamsgaard.dk/1d_fd_simple_shear
Log | Files | Refs | README | LICENSE Back to index

shear_flux (176B)

      1 #!/usr/bin/awk -f
      2 BEGIN {
      3 	getline
      4 	integral = 0.0
      5 }
      6 NF{
      7 	if (NR > 1)
      8 	{
      9 		integral += ($2 + v_prev)/2.0*($1 - z_prev);
     10 	}
     11 	z_prev = $1;
     12 	v_prev = $2;
     13 }
     14 END{
     15 	print integral
     16 }