pkg/exf/exf_bulkqsat.F

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

#include "EXF_CPPOPTIONS.h"


      _RL function exf_BulkqSat(
     I                           tkelvin
     &                         )

c     ==================================================================
c     FUNCTION exf_BulkqSat
c     ==================================================================
c
c     o Compute the saturation specific humidity.
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_BulkqSat
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 tkelvin

c     == end of interface ==

      exf_BulkqSat = cvapor_fac/exp(cvapor_exp/tkelvin)

#else

c     == function arguments ==
      _RL tkelvin

      exf_BulkqSat = 0.0
#endif

      end


1	dimitri	1.2	c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_bulkqsat.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_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	dimitri	1.2	c Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002
29			c
30			c - Added EXF_NO_BULK_COMPUTATIONS option.
31			c
32	heimbach	1.1	c ==================================================================
33			c FUNCTION exf_BulkqSat
34			c ==================================================================
35
36			implicit none
37
38	dimitri	1.2	#if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP) \
39			&& ~defined (EXF_NO_BULK_COMPUTATIONS) )
40	heimbach	1.1
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