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	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_bulkcdn.F,v 1.1 2001/05/14 22:08:40 heimbach Exp $
2
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	c Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
30	c
31	c - Added EXF_NO_BULK_COMPUTATIONS option.
32	c
33	c ==================================================================
34	c FUNCTION exf_BulkCdn
35	c ==================================================================
36
37	implicit none
38
39	#if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP) \
40	&& ~defined (EXF_NO_BULK_COMPUTATIONS) )
41
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