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

Data files needed for look up table (hurricane version) 

cd_hur.dat:  contains drag coefficients
ch_hur.dat:  contains heat transfer coefficients
ce_hur.dat:  contains moisture transfer coefficients
table_data_hur.dat: contains data used for computing transfer coefficients
                     (wind speed, SST-Tair, SST, reference height)



Programs needed for flux table 

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

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

flux_hur.inc: used to share data via a COMMON block



TEST PROGRAM
-------------

(1) test_fluxtable_hur.f: (main program) Used to test 
                          the flux table. 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.f: (subroutine) Reads example input data (your_data.dat).                            Data 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_hur.dat: (data file) Used to compare output from flux
                             table. Compare to output_hur.dat

                            - output_hur.dat produced by test_fluxtable_hur.f 


Compile by typing: f77 read_data_hur.f read_your_data.f flux_table_hur.f                                test_fluxtable_hur.f 



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/2/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.