/[MITgcm]/MITgcm/pkg/exf/exf_bulkqsat.F
ViewVC logotype

Annotation of /MITgcm/pkg/exf/exf_bulkqsat.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (hide annotations) (download)
Tue Feb 18 05:33:54 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint48i_post, checkpoint48h_post, checkpoint49, checkpoint48g_post, checkpoint50, checkpoint50a_post
Changes since 1.2: +2 -7 lines
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F

1 dimitri 1.3 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_bulkqsat.F,v 1.1.4.2 2003/01/12 08:20:10 dimitri 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     c ==================================================================
29     c FUNCTION exf_BulkqSat
30     c ==================================================================
31    
32     implicit none
33    
34 dimitri 1.3 #if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP))
35 heimbach 1.1
36     c == global variables ==
37    
38     #include "SIZE.h"
39     #include "exf_constants.h"
40    
41     c == function arguments ==
42     _RL tkelvin
43    
44     c == end of interface ==
45    
46     exf_BulkqSat = cvapor_fac/exp(cvapor_exp/tkelvin)
47    
48     #else
49    
50     c == function arguments ==
51     _RL tkelvin
52    
53     exf_BulkqSat = 0.0
54     #endif
55    
56     end
57    
58    
59    

  ViewVC Help
Powered by ViewVC 1.1.22