anchor.diff:      ! If there is only 1 value, set the full constraint array to that value
       if (rake_constraint%nrows.eq.1) then
           dp1 = rake_constraint%array(1,1)
  -        dp2 = rake_constraint%array(1,2)
  +        if (rake_constraint%ncols.eq.2) then
  +            dp2 = rake_constraint%array(1,2)
  +        endif
           if (allocated(rake_constraint%array)) then
               deallocate(rake_constraint%array)
           endif
  
positive.code:         endif
          rake_constraint%nrows = fault%nrows
          rake_constraint%array(:,1) = dp1
          rake_constraint%array(:,2) = dp2
      endif
  else
      if (inversion_mode.eq.'anneal') then
  
negative.code: !         i = i + 1
  !         call get_command_argument(i,tag)
  !         read(tag,*) los_weight
  !     elseif (trim(tag).eq.'-los:misfit') then
  !         i = i + 1
  !         call get_command_argument(i,los_misfit_file)
  