anchor.diff:          res = res + 10000
         elseif(code == 'D') then
           res = res + 20000
  +      elseif(code == 'i') then
  +        res = res + 5000
         end if
   
       end function ind
  
positive.code:         call open_files('D')
          call write_devops()
          call close_files('D')
  
      end subroutine do_lindblad_fit_work
  
  
negative.code:       end do
      end subroutine read_evops
  
      subroutine write_evops(code)
        integer (i4b)          :: i, file_ios
        integer(i4b)           :: Uelement, Uelement2,Utnemele,Utnemele2
        character, intent(in)  :: code
  
        do i=1,size(Evops,5)
        do Uelement=1,Nl1
  