anchor.diff:        f1=spa+spd
         f2=f0**2+four*det                                                  !hr06
         if(f2 .lt. zero) then 
  -        write(lout,*) 'ERROR in betalf() - f2 < 0: ',  f2, ' f0: ', f0, ' det: ', det
  +        write(lout,'(a,F12.5, a, F12.5, a, F12.5)') 'ERROR in betalf() - f2 < 0: ',  f2, ' f0: ', f0, ' det: ', det
           goto 160
         end if
         f2=sqrt(f2)
  
positive.code:       rclam2=(egwg2+rca2)*half
        yclam2=yca2*half
        if(egwg1**2 .ge. four) then 
          write(lout,*) 'ERROR in betalf() - egwg1**2 > 4: ',  egwg1**2, ' f0: ', spa-spd, &
          ' f1: ', spa+spd, ' f2: ', f0**2+four*det, ' spa: ', spa, ' spd: ', spd
          write(lout,*) 'ERROR in betalf() - am: ',  am
          goto 160                                    !hr06
        end if
        if(egwg2**2 .ge. four) then
          write(lout,*) 'ERROR in betalf() - egwg2**2 > 4: ',  egwg2**2
          goto 160                                    !hr06
        end if
     50 continue
  
negative.code:      &+ta(3,1)*ta(4,2))-ta(4,1)*ta(3,2)                                  !hr06
        if(rn1.lt.zero) goto 70                                             !hr06
        if(rn1.eq.zero) then 
          write(lout,'(a,F12.5)') 'ERROR in betalf() - rn1 = 0: ', rn1
          goto 160                                            !hr06
        end if
        if(rn1.gt.zero) goto 90                                             !hr06
  