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

Annotation of /MITgcm/pkg/exf/exf_bulkcdn.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, 3 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_bulkcdn.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_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 ==================================================================
30     c FUNCTION exf_BulkCdn
31     c ==================================================================
32    
33     implicit none
34    
35 dimitri 1.3 #if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP))
36 heimbach 1.1
37     c == global variables ==
38    
39     #include "SIZE.h"
40     #include "exf_constants.h"
41    
42     c == function arguments ==
43     _RL umps
44    
45     c == end of interface ==
46    
47     exf_BulkCdn = cdrag_1/umps + cdrag_2 + cdrag_3*umps
48    
49     #else
50    
51     c == function arguments ==
52     _RL umps
53    
54     exf_BulkCdn = 0.0
55     #endif
56    
57     end

  ViewVC Help
Powered by ViewVC 1.1.22