anchor.diff:      write(stdout,'(X,"(DEBUG)",X,A,F12.6)')   'read_halfspace: poissons_ratio= ',poisson
       write(stdout,'(X,"(DEBUG)",X,A,1PE12.4)') 'read_halfspace: lame_parameter= ',lame
   endif
  -if (verbosity.ge.3) then
  +if (verbosity.ge.2) then
       write(stdout,'(X,A,1PE12.4)') 'read_halfspace: shear_modulus=  ',shearmod
       write(stdout,'(X,A,F12.6)')   'read_halfspace: poissons_ratio= ',poisson
       write(stdout,'(X,A,1PE12.4)') 'read_halfspace: lame_parameter= ',lame
  
positive.code: endif
  
  
  if (verbosity.ge.3) then
      if (coord_type.eq.'geographic') then
          write(stdout,*) '   #         Lon         Lat   Dep(km)     Str     Dip     Rak   '//&
  
negative.code: endif
  
  
  if (verbosity.ge.2) then
      write(stdout,*) 'read_stations: read',nstations,' stations'
  endif
  if (verbosity.ge.3) then
      if (auto_mode.eq.'') then
          if (coord_type.eq.'geographic') then
  