rc_sea_spray
============

.. py:module:: rc_sea_spray


Attributes
----------

.. autoapisummary::

   rc_sea_spray.A_w
   rc_sea_spray.h
   rc_sea_spray.U10
   rc_sea_spray.C
   rc_sea_spray.d_p
   rc_sea_spray.rho_p
   rc_sea_spray.k_a


Functions
---------

.. autoapisummary::

   rc_sea_spray.emission_factor
   rc_sea_spray.settling_velocity
   rc_sea_spray.aerosolization_flux
   rc_sea_spray.aerosolization_rate_constant


Module Contents
---------------

.. py:function:: emission_factor(dp_um)

   Computes the emission factor (E_f) based on particle diameter (dp) in micrometers. After fitting a power law distribution to empirical data.

   :param E_f_size_dict: dict
                         A dictionary containing the particle diameter (dp) in micrometers (dp_um) as the key and the corresponding emission factor (E_f) as the value.
   :param dp_um: float or array-like
                 Particle diameter in micrometers (µm)

   :returns:

             float or array-like
                 Emission factor in m³ m⁻² s⁻¹
   :rtype: E_f


.. py:function:: settling_velocity(d_p, rho_p, rho_w=1000, mu=0.001002, g=9.81)

   Compute the Stokes settling velocity for a particle in water.

   Parameters:
   d_p   : Particle diameter (m)
   rho_p : Particle density (kg/m³)
   rho_w : Water density (kg/m³) (default = 1000)
   mu    : Dynamic viscosity of water (Pa·s) (default = 1.002e-3)
   g     : Gravitational acceleration (m/s²) (default = 9.81)

   Returns:
   v_s   : Settling velocity (m/s)


.. py:function:: aerosolization_flux(U10, C, d_p, rho_p, E_f_base)

   Compute the aerosolization flux incorporating particle density.

   Parameters:
   U10       : Wind speed at 10 m height (m/s)
   C         : MNP concentration in water (m⁻³)
   d_p       : Particle diameter (m)
   rho_p     : Particle density (kg/m³)
   E_f_base  : Baseline emission factor (m³ m⁻² s⁻¹) from the paper

   Returns:
   F         : Aerosolization flux (m⁻² s⁻¹)


.. py:function:: aerosolization_rate_constant(U10, C, d_p, rho_p, A_w, h)

   Computes the first-order aerosolization rate constant k_a.

   Parameters:
   U10       : Wind speed at 10 m height (m/s)
   C         : MNP concentration in water (m⁻³)
   d_p       : Particle diameter (m)
   rho_p     : Particle density (kg/m³)
   E_f_base  : Baseline emission factor (m³ m⁻² s⁻¹)
   A_w       : Surface area of the water compartment (m²)
   h         : Depth of the water compartment (m)

   Returns:
   k_a       : Aerosolization rate constant (s⁻¹)


.. py:data:: A_w
   :value: 1000000.0


.. py:data:: h
   :value: 10


.. py:data:: U10
   :value: 5


.. py:data:: C
   :value: 1000000.0


.. py:data:: d_p
   :value: 1e-06


.. py:data:: rho_p
   :value: 900


.. py:data:: k_a