FLUX TABLE 
----------

Data files needed for look up table (Jiping Liu hurricane version)

cd_hur_jiping.dat:  contains drag coefficients
ch_hur_jiping.dat:  contains heat transfer coefficients
ce_hur_jiping.dat:  contains moisture transfer coefficients
input_jiping.dat :  contains data used for computing transfer coefficients
                    (reference height, Air temp, Wind Speed, & SST - Tair)


Programs needed for flux table 

flux_table_hur.f90: (subroutine) returns 'u' stress, 'v' stress, 
                     sensible heat flux,and latent heat flux

read_data_hur.f90: (subroutine) reads data used to construct 
                     flux table and flux table transfer coefficients

flux_hur.f90     : (MODULE) used to share data between programs 



TEST PROGRAMS
--------------

(1) test_fluxtable_hur.f90: (main program) Used to test flux_table and shows 
                            the structure of main 
                            program and calling sequence of include file and 
                            subroutines. The correct calling sequence greatly 
                            improves the time efficiency of the flux table.
                            This is important when incorporting the flux table 
                            into a larger models.  

(2) read_your_data.f90: (subroutine) Reads example input data (your_data.dat).                                                                used in flux table to compute stress, sensible heat flux,
                         and latent heat flux

(3) your_data.dat: (data file) Contains example input data 


(4) example_output_hur90.dat: (data file) Used to compare output from flux
                              table.  Compare to output_hur90.dat

                            - output_hur90.dat produced by test_fluxtable_hur.f90


Compile by typing: f90 flux_hur.f90 read_data_hur.f90 read_your_data.f90 
                   flux_table_hur.f90 test_fluxtable_hur.f90 
 


DESCRIPTION OF FLUX TABLE SUBROUTINE
------------------------------------

SUBROUTINE flux_table_hur

      DEVELOPED BY: Dr. Mark A. Bourassa
                    Center for Ocean-Atmospheric Prediction Studies
                    Florida State University
                    Tallahassee, FL 32306-2840

                    Email: bourassa@coaps.fsu.edu
                    Phone: (850) 644-6923

                     
                    Paul Hughes
                    Center for Ocean-Atmospheric Prediction Studies
                    Florida State University
                    Tallahassee, FL 32306-2840 

                    Email: hughes@coaps.fsu.edu


      VERSION: 1.0

      MODIFICATIONS:

      Date       Programmer      Description of Change 
      ====       ==========      =====================
      12/19/05    P. Hughes       Original Code


      PURPOSE: Subroutine computes wind stress (u and v components), 
               sensible heat flux, and latent heat flux via the bulk flux 
               formulas.  Transfer coefficients calculated using a 
               surface flux model (Bourassa 2005)


      INPUT VARIABLES: 'u' component of wind speed  [m/s]
                       'v' component of wind speed  [m/s]
                       sea surface temperature  [Celsius]
                       air temperature  [Celsius]
                       reference height  [m]
                       significant wave height  [m]
                       'u' orbital velocity  [m/s]
                       'v' orbital velocity  [m/s]
                       specific humidity (surface)  [kg/kg]
                       specific humidity (reference height)  [kg/kg]
                       surface pressure  [Pa]

      OUTPUT VARIABLES: 'u' component of stress  [N/m^2]
                        'v' component of stress  [N/m^2]
                        sensible heat flux  [W/m^2]
                        latent heat flux  [W/m^2] 


     RANGE OF FLUX TABLE: 

     wind speed (0.0 < wind speed < 90.0 m/s)
     SST-Tair (-20.0 < SST-Tair < 20.0 degress C)
     reference height: (10.0 <= adjusted referecne height < 300.0 m)

     *** If additional ranges needed please contact ***


      ADDITIONAL COMMENTS:

     (1) Wind speed adjusted using orbital velocity
         Uorb=pi*Hsig/Peak Period
         adjusted wind speed = windspd - (0.8*Uorb)

     (1b) For more accurate fluxes the currents (if applicable) should
          also be subtracted from the wind speed
          adjusted wind speed = windspd -(0.8*Uorb) - current

     (2) Reference height adjusted using significant wave height
         adjusted reference height = ref_ht - (0.8*Hsig)

     (3) If no wave information then set significant wave height,
         'u' and 'v' orbital velocity equal to zero

     (4) If using total wind speed (orbital velocity) instead of 
         vector components then set 'v' ('v_orb') equal to 0.0 and 
         'u' ('u_orb') equal to the total wind speed (orbital velocity).
         The stress will be output in the 'tau_u' variable

     (5) The surface flux model used to determine transfer coefficients
         is a combination of Bourassa et al.(1999), Clayson et al.(1996),
         and Bourassa (2005) 

     (6) Postive sensible and latent heat flux => ocean to atmoshpere
         Postive Stress => atmoshpere to ocean

     (7) Validation of flux table in a paper in preparation


      REFERENCES:

      Bourassa, M.A., D.G. Vincent, and W.L. Wood, 1999: A flux 
         parameterization including the effects of capillary waves
         and sea state. J. Atmos. Sci., 56, 1123-1139.


      Bourassa, M.A., 2005, Satellite-based observations of surface
         turbulent stress during severe weather, Atmosphere-Ocean 
         Interactions, Vol. 2., ed., W. Perrie, Wessex Institute of 
         Technology, in press.

      Clayson, C.A., C.W. Fairall, and J.A. Curry, 1996: Evaluation of 
         turbulent fluxes at the ocean surface using surface renewal
         theory. J. Geophys. Res., 101, 28503-28513.