anchor.diff:                     iWantDisp, iWantStrain, iWantStress, iWantTraction, &
                      nstations, stations, targets, &
                      lame, shearmod
  +use elast, only: strain2stress, max_shear_stress, calc_tractions
   use io, only: stderr, verbosity
   implicit none
   
  
positive.code:         endif
          ! Maximum (effective) shear stress: estress
          if (estress_file.ne.'') then
              call calc_estress(estress,stress)
              write(121,*) stations(iSta,:),estress
          endif
      endif
  
negative.code: 
  return
  end subroutine calc_strain
  
  !--------------------------------------------------------------------------------------------------!
  
  