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

Annotation of /MITgcm/pkg/bulk_force/exf_bulkcdn.F

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


Revision 1.1 - (hide annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 5 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint47j_post, checkpoint48d_pre, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint47i_post, checkpoint48i_post, checkpoint47f_post, checkpoint48d_post, checkpoint47c_post, checkpoint50e_post, checkpoint50c_post, checkpoint47d_post, checkpoint47a_post, checkpoint51f_pre, checkpoint48e_post, checkpoint48h_post, checkpoint50c_pre, branchpoint-genmake2, checkpoint50d_pre, checkpoint47h_post, checkpoint48c_post, checkpoint50b_pre, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint48g_post, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint51b_pre, checkpoint47g_post, checkpoint51h_pre, checkpoint50g_post, checkpoint51g_post, checkpoint51f_post, checkpoint48b_post, checkpoint50b_post, checkpoint48a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47d_pre, checkpoint51d_post, checkpoint48c_pre, checkpoint50h_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre
Branch point for: branch-genmake2, branch-exfmods-curt
Two packages:  bulk_force (Bulk forcing)
and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
Earlier integration from Stephaine Dutkiewicz
and Patrick Heimbach.

Two ifdef statements for compile time,
ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE

Two switches in data.pkg to turn on at run-time:

cat data.pkg
# Packages
 &PACKAGES
 useBulkForce=.TRUE.,
 useThermSeaIce=.TRUE.,
 &

WARNING:  useSEAICE and useThermSEAICE are mutually exclusive.

The bulk package requires an additional parameter file
with two namelists, data.ice and data.blk.

c ADAPTED FROM:
c LANL CICE.v2.0.2
c-----------------------------------------------------------------------
c.. thermodynamics (vertical physics) based on M. Winton 3-layer model
c.. See Bitz, C. M. and W. H. Lipscomb, 1999:  "An energy-conserving
c..       thermodynamic sea ice model for climate study."  J. Geophys.
c..       Res., 104, 15669 - 15677.
c..     Winton, M., 1999:  "A reformulated three-layer sea ice model."
c..       Submitted to J. Atmos. Ocean. Technol.

c.. authors Elizabeth C. Hunke and William Lipscomb
c..         Fluid Dynamics Group, Los Alamos National Laboratory
c-----------------------------------------------------------------------

1 cheisey 1.1 c $Header: /u/u0/gcmpack/models/MITgcmUV/pkg/exf/exf_bulkcdn.F,v 1.1 2001/05/14 22:08:40 heimbach Exp $
2    
3    
4     #include "CPP_OPTIONS.h"
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_BULK_FORCEE, ALLOW_ATM_TEMP
28     c swd, Apr 2002: modified to be used in bulkf package
29     c
30     c ==================================================================
31     c FUNCTION exf_BulkCdn
32     c ==================================================================
33    
34     implicit none
35    
36     #ifdef ALLOW_BULK_FORCE
37    
38     c == global variables ==
39    
40     #include "SIZE.h"
41     #include "BULKF_ICE_CONSTANTS.h"
42    
43     c == function arguments ==
44     _RL umps
45    
46     c == end of interface ==
47    
48     exf_BulkCdn = cdrag_1/umps + cdrag_2 + cdrag_3*umps
49    
50     #else
51    
52     c == function arguments ==
53     _RL umps
54    
55     exf_BulkCdn = 0.0
56     #endif
57    
58     end

  ViewVC Help
Powered by ViewVC 1.1.22