janus_fc_dat.py 文件源码

def calc_eff_area( self, v ) :

        # Note. #nvn is a vector similar to inflow particle bulk 
        # velocity and ndir is the look direction. 


        # Normalize the particle velocity.

        vn  = calc_arr_norm( v )
        nvn = tuple( [ -c for c in vn ] )

        # Calculate the particle inflow angle (in degrees) relative to
        # the cup normal (i.e., the cup pointing direction).

        psi = acos( calc_arr_dot( self['dir'], nvn ) )*pi/180.
        if ( psi > 90. ) :
            return 0. 

        # Return the effective collecting area corresponding to "psi".

        return interp( psi, self._spec._eff_deg, self._spec._eff_area )

    #-----------------------------------------------------------------------
    # DEFINE THE FUNCTION TO CALCULATE EXPECTED MAXWELLIAN CURRENT.
    #-----------------------------------------------------------------------