.. +
..
.. TODO
.. ====
..
.. anchors
..
.. better ReST
..
.. pb with math directives on aedon : probably some latex installation pb
..
.. ref to tropflux.sh covering the 3 last steps
..
..
.. EVOLUTIONS
.. ==========
..
.. $Id$
..
.. $URL$
..
.. - fplod 20110826T133731Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * fix bad looking in html of 2m indice : m is too hight
..
.. - fplod 20110811T162125Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * revision of steps with pk
..   * revison of ISCCP paragraph
..
.. - fplod 20101222T173813Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * add ref to interpolation tools
..
.. - fplod 20101217T105140Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * add anchors
..
.. - fplod 20101216T160822Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * add ref to OAFLUX tools
..
.. - fplod 20101216T120541Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * add ref to compute_erai_daily_region_2d.sh
..
.. - fplod 20101125T084625Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * use :math: to improve notation of geophysical parameter thanks to
..     http://sphinx.pocoo.org/latest/ext/math.html#module-sphinx.ext.mathbase
..     and ftp://ftp.ams.org/pub/tex/doc/amsmath/amsldoc.pdf
..
.. - fplod 20101123T145240Z aedon.locean-ipsl.upmc.fr (Darwin)
..
..   * creation from mail send by Jérôme Vialard
..
.. -

Generate Tropflux heat fluxes, 4 steps
======================================

.. _gather_data:

Gather data
-----------

OAFLUX data

   for sea land mask

   tools:

   - :ref:`get_oaflux.sh`
   - :func:`oaflux_mask_30n30s`

ERA-I data
    (wind components at 10 m, air temperature at 2 m, dewpoint temperature at 2 m, net SW, net LW)

   tools:

   - :ref:`compute_erai_daily_region_2d.sh`

ISCCP shortwave

   tools:
    - :ref:`get_swr.sh`


OLR
   get NOAA interpolated and uninterpolated OLR data

   tools:
    - :ref:`get_olr.sh`
    - :ref:`get_olrnint.sh`

.. _interpolate_data:

Interpolate datasets onto the OAFLUX grid
-----------------------------------------

tools :
 - :func:`interp_erai_dewt`
 - :func:`interp_erai_lwr`
 - :func:`interp_erai_sst`
 - :func:`interp_erai_msl.pro`
 - :func:`interp_erai_t2m`
 - :func:`interp_erai_ws`


The OLR interpolation also involves combining the delayed time and realtime OLR
datasets, which is  not done extremely cleanly now.

 - :func:`interp_olr_30n30s`

There's also one additional step to convert ERA-I dewpoint T :math:`d_{2m}` into specific humidity :math:`q_{2m}` :
 
 - :func:`d2m_to_q2m_erai` (:func:`rh_to_spechum`, :func:`qsat`)

Result after that, 7 netcdf files for: NOAA OLR, ERA-I :math:`T_{2m}`, ERA-I :math:`q_{2m}`, ERA-I wind speed, ERA-I SST, ISCCP net SW, ERA-I net LW.

.. _mooring_corrections:

Apply corrections which have been previously estimated from mooring data (see :ref:`generate_statistics`)
---------------------------------------------------------------------------------------------------------

There are two steps for the wind correction: the statistical correction from
observations, and the gustiness correction.

We apply statistical correction at this stage.

tools :
 - :func:`tropflux_swr_dt`, :func:`tropflux_swr_nrt`, :func:`tropflux_swr_blnd`
 - :func:`lwr_correction_ncdf`
 - :func:`q2m_correction_ncdf`
 - :func:`sst_correction_ncdf`
 - :func:`t2m_correction_ncdf`
 - :func:`ws_correction_ncdf`

Result after that, 7 netcdf files for: NOAA OLR, ERA-I :math:`T_a`, ERA-I :math:`q_a`, ERA-I wind speed, ERA-I SST, ISCCP net SW, ERA-I net LW.

.. _apply_coare:

Apply the COARE v3 algorithm to compute turbulent fluxes
--------------------------------------------------------

Wind gustiness correction is applied here.

tools :
 - :func:`cronin_gustiness_ncdf`
 - :func:`tropflux` (:func:`cor30a`), :func:`tropflux_nrt_ncdf` 

After that step, there are 5 files: Tropflux :math:`T_{2m}`, :math:`q_{2m}`, SST, wind, and one file with the 4 components of the net flux (sensible, latent, SWR, LWR)

.. ++ cf pk paper 20110811 to be add in th document with better form
.. math:`T_{2m}`, :math:`q_{2m}`, SST, wind, swr and lwr came from correction mooring
.. final wind product is not corrected for wind gustiness.
.. sensible, latent shf and lhf came from coare correction 

.. _momentum_fluxes:

Momentum Fluxes
===============

tools :

 - :ref:`TropFlux_wind_stress_19890101_20101231.pro`

After generation of heat fluxes we process Momentum Fluxes using ERA-I wind components, COARE gustiness and wind stress magnitude returned from COARE algorithm.