private  subroutine get_g_integral_contribution_source(g, source)  
        
    Layouts
Parameters
Grids
Calculations
assume there is only a single flux surface being simulated
integrate  over velocity space and sum over species within each processor
as v-space and species possibly spread over processors, wlil need to
gather sums from each proceessor and sum them all together below
gather sub-sums from each processor and add them together
store result in phi, which will be further modified below to account for polarization term
no longer need , so deallocate
    Arguments
        
    
      
        | Type | Intent | Optional | Attributes |  | Name |  | 
    
        
            | complex, | intent(in), |  | dimension(:, :, -nzgrid:, :, vmu_lo%llim_proc:) | :: | g |  | 
        
            | complex, | intent(inout), |  | dimension(:, :, -nzgrid:) | :: | source |  |