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