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

Contents of /MITgcm/pkg/bulk_force/exf_bulkrhn.F

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


Revision 1.1 - (show annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 6 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint47b_post, checkpoint51f_post, checkpoint48i_post, checkpoint47a_post, checkpoint48d_post, checkpoint50b_post, checkpoint47i_post, checkpoint48g_post, branchpoint-genmake2, checkpoint50c_pre, checkpoint50, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint50f_post, checkpoint50a_post, checkpoint48e_post, checkpoint47c_post, checkpoint50f_pre, checkpoint48b_post, checkpoint47j_post, checkpoint47d_pre, checkpoint50d_pre, checkpoint47h_post, checkpoint51d_post, checkpoint51, checkpoint48c_pre, checkpoint48c_post, checkpoint50d_post, checkpoint47f_post, checkpoint51b_pre, checkpoint50e_post, checkpoint47g_post, checkpoint50h_post, checkpoint50c_post, checkpoint51a_post, checkpoint47d_post, checkpoint50e_pre, checkpoint50b_pre, checkpoint48d_pre, checkpoint50i_post, checkpoint51e_post, checkpoint51b_post, checkpoint48a_post, checkpoint51h_pre, checkpoint48f_post, checkpoint51i_pre, checkpoint50g_post, checkpoint51c_post, checkpoint51g_post, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint48h_post
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 c $Header: /u/u0/gcmpack/MITgcm/pkg/exf/exf_bulkrhn.F,v 1.1 2001/05/14 22:08:40 heimbach Exp $
2
3 #include "CPP_OPTIONS.h"
4
5 _RL function exf_BulkRhn(
6 I stab
7 & )
8
9 c ==================================================================
10 c FUNCTION exf_BulkRhn
11 c ==================================================================
12 c
13 c o Compute the Dalton number as a function of stability.
14 c
15 c started: Christian Eckert eckert@mit.edu 27-Aug-1999
16 c
17 c changed: Christian Eckert eckert@mit.edu 14-Jan-2000
18 c
19 c - Restructured the code in order to create a package
20 c for the MITgcmUV.
21 c
22 c Patrick Heimbach heimbach@mit.edu 05-May-2000
23 c
24 c - Included appropriate CPP options
25 c ALLOW_BULK_FORCEE, ALLOW_ATM_TEMP
26 c
27 c swd, April 2002: changed to be used in bulkf package
28 c
29 c ==================================================================
30 c FUNCTION exf_BulkRhn
31 c ==================================================================
32
33 implicit none
34
35 #ifdef ALLOW_BULK_FORCE
36
37 c == global variables ==
38
39 #include "SIZE.h"
40 #include "BULKF_ICE_CONSTANTS.h"
41
42
43 c == function arguments ==
44 _RL stab
45 _RL cstanton_1, cstanton_2
46 cstanton_1 = 0.0327000
47 cstanton_2 = 0.0180000
48
49 exf_BulkRhn = (1.d0 - stab)*cstanton_1 + stab*cstanton_2
50
51 #else
52
53 c == function arguments ==
54 _RL stab
55
56 exf_BulkRhn = 0.0
57
58 #endif
59
60 end

  ViewVC Help
Powered by ViewVC 1.1.22