anchor.diff:        return
         end
   
  +
  
positive.code: c	The first angle is angle between cell and the satellite
  c	The second angle is the angle between the LOS and the cell
  c	The 3rd is for the fact towards high delta we have smaller surface
  		   proj_surf = surf*dcos(angle)*abs(prod/(R_Earth*d))*dcos(delta_grid)
  
                 if(include_pst) then
                   energy = energy + 
  
negative.code:             endif
           enddo
        enddo
  
        straylight_flux = energy * pixel !photons se px-1
  
        return
  