diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc index ad93f6f274dd00a857cb9279ca345a372a07228d..6c050198a78d0e5783faa117cfe707a24473e346 100644 --- a/cosserat-continuum.cc +++ b/cosserat-continuum.cc @@ -61,6 +61,9 @@ void dirichletValues(const FieldVector<double,dim>& in, FieldVector<double,3>& o double homotopy ) { + if (not vIt->geometry().corner(0)[0] > upper[0]-1e-3) + return; + double angle = 8*M_PI; angle *= homotopy;