anchor.diff:         end if
          ! Need to change ispot if sphere and falls inside spot
       case(4)
  -       call emit_from_map(src,p,icell)
  +       call emit_from_map(src,p)
       case(5)
          call emit_from_extern_sph(src,p)
       case(6)
  
positive.code:           case(2)
               call random_planck_frequency(p%nu, src%spot(ispot)%temperature)
            case(3)
               p%dust_id = select_dust_specific_energy_rho(icell)
               p%emiss_var_id = jnu_var_id(icell%ic, p%dust_id)
               p%emiss_var_frac = jnu_var_frac(icell%ic, p%dust_id)
               call dust_sample_j_nu(d(p%dust_id),p%emiss_var_id,p%emiss_var_frac,p%nu)
            end select
         end if
  
negative.code:   ! emit_from_map : emit a photon from a luminosity map
    !**********************************************************************!
  
    subroutine emit_from_map(src,p)
  
      implicit none
  
  