/[MITgcm]/MITgcm/pkg/bulk_force/BULKF_ICE_CONSTANTS.h
ViewVC logotype

Annotation of /MITgcm/pkg/bulk_force/BULKF_ICE_CONSTANTS.h

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: checkpoint47b_post, checkpoint51f_post, checkpoint48i_post, checkpoint51l_post, checkpoint47a_post, checkpoint51k_post, checkpoint48d_post, checkpoint50b_post, checkpoint47i_post, checkpoint51o_post, checkpoint48g_post, branchpoint-genmake2, checkpoint50c_pre, checkpoint50, checkpoint51j_post, branch-exfmods-tag, checkpoint51q_post, checkpoint47e_post, checkpoint50f_post, checkpoint50a_post, checkpoint48e_post, checkpoint47c_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint48b_post, checkpoint47j_post, checkpoint47d_pre, checkpoint50d_pre, checkpoint47h_post, checkpoint51d_post, checkpoint51, checkpoint51r_post, checkpoint48c_pre, checkpoint52, checkpoint52b_pre, checkpoint48c_post, checkpoint50d_post, checkpoint51o_pre, checkpoint47f_post, checkpoint51t_post, checkpoint51b_pre, checkpoint52a_post, checkpoint51i_post, checkpoint50e_post, checkpoint47g_post, checkpoint50h_post, checkpoint50c_post, checkpoint51a_post, checkpoint51n_pre, checkpoint47d_post, checkpoint50e_pre, checkpoint50b_pre, checkpoint48d_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51e_post, checkpoint51b_post, checkpoint48a_post, checkpoint51h_pre, checkpoint48f_post, checkpoint51i_pre, checkpoint51l_pre, checkpoint50g_post, checkpoint51u_post, checkpoint51c_post, checkpoint51g_post, checkpoint51m_post, ecco_c52_e35, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint51s_post, checkpoint48h_post
Branch point for: checkpoint51n_branch, branch-nonh, tg2-branch, branch-genmake2, branch-exfmods-curt
File MIME type: text/plain 
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----------------------------------------------------------------------------
2     c.. Common blocks for almost everything that the sea ice model passes around.
3     c----------------------------------------------------------------------------
4     c.. number layers of ice
5     integer nlyr ! maximum number of ice layers
6     _RL rnlyr ! real value of nlyr
7     parameter (nlyr = 2)
8     parameter (rnlyr = 2.)
9    
10     c.. densities
11     _RL rhoa ! density of air (kg/m^3)
12     _RL rhos ! density of snow (kg/m^3)
13     _RL rhoi ! density of ice (kg/m^3)
14     _RL rhosw ! density of seawater (kg/m^3)
15     _RL rhoiw ! ocean-ice density difference (kg/m^3)
16    
17     parameter (rhoa =1.3)
18     parameter (rhos = 330.)
19     parameter (rhoi = 900.)
20     parameter (rhosw = 1026.)
21     parameter (rhoiw = rhosw - rhoi)
22    
23     c.. specific heats
24     _RL cpair ! specific heat of air (J/kg/K)
25     _RL cpice ! specific heat of fresh ice (J/kg/K)
26     _RL cpwater ! specific heat of water (J/kg/K)
27     _RL cpwv ! specific heat of water vapour (J/kg/KC)
28     _RL cpvir
29    
30     parameter(cpair = 1004.)
31     parameter(cpice = 2106.)
32     parameter (cpwater = 4218.)
33     parameter (cpwv = 1.81e3)
34     parameter (cpvir = cpwv/cpair -1.)
35    
36     c .. thermal conductivity. QQ check units
37     _RL kice ! thermal conductivity of pure ice (W/m/K)
38     _RL ksnow ! thermal conductivity of snow (W/m/K)
39    
40     parameter (kice = 2.03) !QQ originally 2.03)
41     parameter (ksnow = 0.30) !QQ originally 0.30)
42    
43     c .. heat transfer coefficient
44     _RL transcoef ! transfer coef between ice and water (unitless)
45    
46     parameter (transcoef=0.006) !QQ originally 0.006
47    
48     c .. latent heat
49     _RL Lvap ! latent heat of vaporizn at 0 C (J/kg)
50     _RL Lfresh ! latent heat of melting of pure ice (J/kg)
51    
52     parameter(Lvap = 2.5e+6)
53     parameter (Lfresh = 3.34e5)
54    
55     c.. Enthalpy
56     _RL qsnow ! snow enthalpy (J/kg)
57     parameter (qsnow = Lfresh)
58    
59     c .. Albedo
60     _RL albsnodry ! albedo of dry snow (Tsfc < 0)
61     _RL albsnowet ! albedo of melting snow (Tsfc = 0)
62     _RL albicemax ! max albedo of bare ice
63     _RL albicemin ! minimum ice albedo (thin melting ice)
64     _RL halb ! melt rate parameter for albedo QQ units
65    
66     parameter (albsnodry = 0.85)
67     parameter (albsnowet = 0.75)
68     parameter (albicemax = 0.65)
69     parameter (albicemin = 0.20)
70     parameter (halb = .5)
71    
72     c.. Solar parameters
73     _RL i0 ! fraction of penetrating solar rad
74     _RL ksolar ! bulk solar abs coeff of sea ice (m-1)
75    
76     parameter (i0 = 0.3) !QQQQQQ was 0.3
77     parameter (ksolar = 1.5) !QQQQQQ was1.5
78    
79    
80     c .. Salinity
81     _RL saltice ! salinity of ice (o/oo)
82     _RL S_winton ! winton salinity of ice (o/oo)
83     _RL mu_Tf ! Tf:brine salinity ratio (C/ppt)
84    
85     parameter (saltice = 4.)
86     parameter (S_winton = 1.)
87     parameter (mu_Tf = 0.054)
88    
89     c .. melting
90     _RL Tf0kel ! Freezing temp of fresh ice in Kelvin = 273.15
91     _RL Tmlt1 ! melting temp; depends on S (C)
92    
93     parameter (Tf0kel = 273.15)
94     parameter (Tmlt1=-mu_Tf*S_winton)
95    
96     c .. wind drag
97     c cdrag_[n] - n = 1,2,3 coefficients used to evaluate
98     c drag coefficient
99     _RL cdrag_1, cdrag_2, cdrag_3
100     parameter ( cdrag_1 = 0.0027000,
101     & cdrag_2 = 0.0001420,
102     & cdrag_3 = 0.0000764)
103    
104     c .. constants
105     _RL stefan ! Stefan-Boltzmann constant (W/m^2 K^4)
106     _RL xkar ! Von Karman constant QQ units?
107     _RL Rvap ! gas constant for H2O vapor (J/kg/K)
108    
109     parameter (stefan = 5.67e-8)
110     parameter (xkar = 0.4 )
111     parameter(Rvap = 461.)
112    
113     c.. Miscellaneous
114     _RL p0 ! surface pressure (mb)
115    
116     parameter(p0 = 1013.)
117    
118     c.. Combinations used for efficiency
119     _RL rhoi_Lfresh ! rhoi*Lfresh
120     _RL Qcoef ! another constant for latent heat flux
121     _RL lvrrv ! (Lvap + Lfresh) / Rvap
122    
123     parameter (rhoi_Lfresh = rhoi*Lfresh)
124     parameter (Qcoef = 6.11*0.622/p0)
125     parameter (lvrrv = (Lvap + Lfresh) / Rvap)
126    
127     c .. Min/Max
128     _RL himin ! minimum thickness for ice (m)
129     _RL Terrmax ! temperature convergence precision (C)
130    
131     parameter (himin = 0.01) !was 0.01?
132     parameter (Terrmax = 5.0e-1) !was 5.0e-3????
133    
134     c .. for bulk formula
135     _RL humid_fac ! const. for the evaluation of the virtual temp.
136     _RL saltsat ! reduction of sat. vapor pressure over salt water
137     _RL gamma_blk ! adiabatic lapse rate
138    
139    
140     parameter (humid_fac = 0.606)
141     parameter (saltsat = 0.980)
142     parameter (gamma_blk = 0.010)
143    
144    
145    
146     c .. for Vince bulk formula QQ check units
147     _RL Lvap_ice
148     _RL Rgas
149     _RL Sha
150    
151     parameter (Lvap_ice = 2.83e6) ! latent heat from sublimation
152     parameter (Rgas = 287.e0) ! gas constant for dry air
153     parameter (Sha = Rgas/.286) !
154    
155     c .. emissivities QQQ find
156     _RL atm_emissivity
157     _RL ocean_emissivity
158     _RL snow_emissivity
159     _RL ice_emissivity
160    
161     cQQQQQQ find real values
162     parameter(atm_emissivity=.90d0) !QQ
163     parameter(ocean_emissivity=.985d0) !QQ
164     parameter(snow_emissivity=.98d0) !QQ
165     parameter(ice_emissivity=.98d0) !QQ
  ViewVC Help
Powered by ViewVC 1.1.22