
// gravity
gx #eval { -9.81*sin($theta*pi()/180.0) };
gy #eval { -9.81*cos($theta*pi()/180.0) };

/// Tank Dimensions
x1 #eval  {-1.0*$R*$x1OverR};
x2 #eval { 1.0*$R*$x2OverR };
y1 #eval  { -1.0*$R*$y1OverR };
y2 #eval  { 1.0*$R*$y2OverR };
z1 #eval { -1.0*$R };
z2 $R;

// clast
xc1 #eval  { -1.0*$R*$LxOverR };
xc2 #eval { 1.0*$R*$LxOverR };
yc1 #eval  { -1.0*$R*$LyOverR };
yc2 #eval  { 1.0*$R*$LyOverR };

// shell
xs1 #eval  { $xc1-$dsOnLx };
xs2 #eval  { $xc2+$dsOnLx };
ys1 #eval  { $yc1-$dsOnLy };
ys2 #eval  { $yc2+$dsOnLy };

Nx #eval { $N*($x2OverR+$x1OverR) };
Ny #eval { $N*($y2OverR+$y1OverR) };

/// Refine boxes
r1x1 #eval { -1.0*$R*$r1x1OverR };
r1x2 #eval { 1.0*$R*$r1x2OverR };
r1y1 #eval { (-1.0*$R*$r1y1OverR) }; 
r1y2 #eval { (1.0*$R*$r1y2OverR) };
r1z1 $r1x1;
r1z2 $r1x2;

r2x1 #eval { -1.0*$R*$r2x1OverR };
r2x2 #eval { 1.0*$R*$r2x2OverR };
r2y1 #eval { (-1.0*$R*$r2y1OverR) }; 
r2y2 #eval { (1.0*$R*$r2y2OverR) };
r2z1 $r1x1;
r2z2 $r1x2;

r3x1 #eval { -1.0*$R*$r3x1OverR };
r3x2 #eval { 1.0*$R*$r3x2OverR };
r3y1 #eval { (-1.0*$R*$r3y1OverR) }; 
r3y2 #eval { (1.0*$R*$r3y2OverR) };
r3z1 $r1x1;
r3z2 $r1x2;

r4x1 #eval { -1.0*$R*$r4x1OverR };
r4x2 #eval { 1.0*$R*$r4x2OverR };
r4y1 #eval { (-1.0*$R*$r4y1OverR) }; 
r4y2 #eval { (1.0*$R*$r4y2OverR) };
r4z1 $r1x1;
r4z2 $r1x2;

