anchor.diff:  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)
  +		   proj_surf = surf*dcos(angle)*abs(prod/(R_Earth*d))
  +     &		 *dcos(delta_grid)
   
                  if(include_pst) then
                    energy = energy + 
  
positive.code:       return
        end
  
  
negative.code:             endif
           enddo
        enddo
  
        straylight_flux = energy * pixel !photons se px-1
  
        return
  