pkg/exf/exf_bulkcdn.F

c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_bulkcdn.F,v 1.1 2001/05/14 22:08:40 heimbach Exp $

#include "EXF_CPPOPTIONS.h"


      _RL function exf_BulkCdn(
     I                          umps
     &                        )

c     ==================================================================
c     FUNCTION exf_BulkCdn
c     ==================================================================
c
c     o Compute the neutral drag coefficient as a function of the wind
c       speed.
c
c     started: Christian Eckert eckert@mit.edu  27-Aug-1999
c
c     changed: Christian Eckert eckert@mit.edu  14-Jan-2000
c
c              - Restructured the code in order to create a package
c                for the MITgcmUV.
c
c              Patrick Heimbach heimbach@mit.edu  05-May-2000
c
c              - Included appropriate CPP options
c                ALLOW_BULKFORMULAE, ALLOW_ATM_TEMP
c
c              Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
c
c              - Added EXF_NO_BULK_COMPUTATIONS option.
c
c     ==================================================================
c     FUNCTION exf_BulkCdn
c     ==================================================================

      implicit none

#if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP) \
    && ~defined (EXF_NO_BULK_COMPUTATIONS) )

c     == global variables ==

#include "SIZE.h"
#include "exf_constants.h"

c     == function arguments ==
      _RL umps

c     == end of interface ==

      exf_BulkCdn = cdrag_1/umps + cdrag_2 + cdrag_3*umps

#else

c     == function arguments ==
      _RL umps

      exf_BulkCdn = 0.0
#endif

      end
1	dimitri	1.2	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_bulkcdn.F,v 1.1 2001/05/14 22:08:40 heimbach Exp $
2	heimbach	1.1
3			#include "EXF_CPPOPTIONS.h"
4
5
6			_RL function exf_BulkCdn(
7			I umps
8			& )
9
10			c ==================================================================
11			c FUNCTION exf_BulkCdn
12			c ==================================================================
13			c
14			c o Compute the neutral drag coefficient as a function of the wind
15			c speed.
16			c
17			c started: Christian Eckert eckert@mit.edu 27-Aug-1999
18			c
19			c changed: Christian Eckert eckert@mit.edu 14-Jan-2000
20			c
21			c - Restructured the code in order to create a package
22			c for the MITgcmUV.
23			c
24			c Patrick Heimbach heimbach@mit.edu 05-May-2000
25			c
26			c - Included appropriate CPP options
27			c ALLOW_BULKFORMULAE, ALLOW_ATM_TEMP
28			c
29	dimitri	1.2	c Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
30			c
31			c - Added EXF_NO_BULK_COMPUTATIONS option.
32			c
33	heimbach	1.1	c ==================================================================
34			c FUNCTION exf_BulkCdn
35			c ==================================================================
36
37			implicit none
38
39	dimitri	1.2	#if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP) \
40			&& ~defined (EXF_NO_BULK_COMPUTATIONS) )
41	heimbach	1.1
42			c == global variables ==
43
44			#include "SIZE.h"
45			#include "exf_constants.h"
46
47			c == function arguments ==
48			_RL umps
49
50			c == end of interface ==
51
52			exf_BulkCdn = cdrag_1/umps + cdrag_2 + cdrag_3*umps
53
54			#else
55
56			c == function arguments ==
57			_RL umps
58
59			exf_BulkCdn = 0.0
60			#endif
61
62			end