anchor.diff:        omega = domega + 2.*PI*((i-1)/(PLFLT)NSHEAR)
         sin_omega = sin(omega)
         cos_omega = cos(omega)
  -      x_shear = 0.5*yrange*sin_omega
  +      x_shear = 0.5*xrange*sin_omega
         z_shear = 0.5*zrange*cos_omega
         call plptex3( &
            xmid, ymax, zmax -(zmax-0.2)*((PLFLT)i/(PLFLT)(NSHEAR-1)), &
  
positive.code: !   Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  !
  
     program x29f95
     use plplot, PI => PL_PI
  
     implicit none
  
negative.code:       sin_omega = sin(omega)
        cos_omega = cos(omega)
        y_shear = 0.5*yrange*cos_omega
        x_shear = 0.5*xrange*sin_omega
        call plptex3( &
           xmid, ymax -(ymax-0.2)*((PLFLT)i/(PLFLT)(NSHEAR-1)), zmin, &
           x_inclination, y_inclination, z_inclination, &
  