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