/[MITgcm]/MITgcm/pkg/thsice/THSICE_PARAMS.h
ViewVC logotype

Annotation of /MITgcm/pkg/thsice/THSICE_PARAMS.h

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


Revision 1.1 - (hide annotations) (download)
Sun Nov 23 01:20:13 2003 UTC (20 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: branch-netcdf, checkpoint52d_pre, checkpoint52d_post, checkpoint52b_post, checkpoint52c_post
Branch point for: netcdf-sm0
File MIME type: text/plain
new pkg "thSIce" (replace therm_seaice).

1 jmc 1.1 C $Header: $
2     C $Name: $
3    
4     #ifdef ALLOW_THSICE
5     C *==========================================================*
6     C | THSICE_PARAMS.h
7     C | o Header file for Therm_SeaIce package parameters:
8     C | - basic parameter ( I/O frequency, etc ...)
9     C | - physical constants (used in therm_SeaIce pkg)
10     C *==========================================================*
11    
12     C----------------------------------------------------------------------------
13     C.. Common blocks for almost everything that the sea ice model passes around.
14     C----------------------------------------------------------------------------
15    
16     C-- COMMON / THSICE_PHYSPAR_R / physical (real) parameter
17     C.. densities
18     C rhos :: density of snow (kg/m^3)
19     C rhoi :: density of ice (kg/m^3)
20     C rhosw :: density of seawater (kg/m^3)
21     C rhofw :: density of fresh water (kg/m^3)
22     C rhoiw :: ocean-ice density difference (kg/m^3)
23     C.. specific heats
24     C cpice :: specific heat of fresh ice (J/kg/K)
25     C cpwater :: specific heat of water (J/kg/K)
26     C .. thermal conductivity. QQ check units
27     C kice :: thermal conductivity of pure ice (W/m/K)
28     C ksnow :: thermal conductivity of snow (W/m/K)
29     C .. heat transfer coefficient
30     C transcoef :: transfer coef between ice and water (unitless)
31     C .. latent heat
32     C Lfresh :: latent heat of melting of pure ice (J/kg)
33     C.. Enthalpy
34     C qsnow :: snow enthalpy (J/kg)
35     C .. Albedo
36     C albsnodry :: albedo of dry snow (Tsfc < 0)
37     C albsnowet :: albedo of melting snow (Tsfc = 0)
38     C albicemax :: max albedo of bare ice
39     C albicemin :: minimum ice albedo (thin melting ice)
40     C halb :: melt rate parameter for albedo QQ units
41     C.. Solar parameters
42     C i0 :: fraction of penetrating solar rad
43     C ksolar :: bulk solar abs coeff of sea ice (m-1)
44     C .. Salinity
45     C saltice :: salinity of ice (o/oo)
46     C S_winton :: winton salinity of ice (o/oo)
47     C mu_Tf :: Tf:brine salinity ratio (C/ppt)
48     C .. melting
49     C Tf0kel :: Freezing temp of fresh ice in Kelvin = 273.15
50     C Tmlt1 :: melting temp; depends on S (C)
51     C .. Min/Max
52     C himin :: minimum thickness for ice (m)
53     C Terrmax :: temperature convergence precision (C)
54     C hiMax :: Maximum thickness for ice (m)
55     C hsMax :: Maximum thickness for snow (m)
56     C .. for fractional ice
57     C iceMaskmax :: maximum Ice fraction (=1 for no fractional ice)
58     C iceMaskmin :: mimimum Ice fraction (=1 for no fractional ice)
59     C himin0 :: minimum ice height (=himin for no fractional ice)
60     C frac_energy:: fraction of energy going to melting (=0 for no fract. ice)
61     C hihig :: ice height above which freezing only occurs over open ocean
62     C (=large for no fractional ice)
63     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
64    
65     COMMON / THSICE_PHYSPAR_R /
66     & rhos, rhoi, rhosw, rhofw, rhoiw,
67     & cpice, cpwater,
68     & kice, ksnow,
69     & transcoef, Lfresh, qsnow,
70     & albsnodry, albsnowet, albicemax, albicemin, halb,
71     & i0, ksolar,
72     & saltice, S_winton, mu_Tf,
73     & Tf0kel, Tmlt1,
74     & himin, Terrmax, hiMax, hsMax,
75     & iceMaskmax, iceMaskmin, himin0,
76     & frac_energy, hihig
77    
78     _RL rhos
79     _RL rhoi
80     _RL rhosw
81     _RL rhofw
82     _RL rhoiw
83     _RL cpice
84     _RL cpwater
85     _RL kice
86     _RL ksnow
87     _RL transcoef
88     _RL Lfresh
89     _RL qsnow
90     _RL albsnodry
91     _RL albsnowet
92     _RL albicemax
93     _RL albicemin
94     _RL halb
95     _RL i0
96     _RL ksolar
97     _RL saltice
98     _RL S_winton
99     _RL mu_Tf
100     _RL Tf0kel
101     _RL Tmlt1
102     _RL himin
103     _RL Terrmax
104     _RL hiMax
105     _RL hsMax
106     _RL iceMaskmax
107     _RL iceMaskmin
108     _RL himin0
109     _RL frac_energy
110     _RL hihig
111    
112     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
113    
114     C-- COMMON / THSICE_PAR_I / ice model (integer) parameters
115     c startIceModel :: =1 : start ice model at nIter0 ; =0 : use pickup files
116     c nitMaxTsf :: maximum Nb of iter to find Surface Temp (Trsf)
117     COMMON / THSICE_PAR_I /
118     & startIceModel, nitMaxTsf
119    
120     INTEGER startIceModel
121     INTEGER nitMaxTsf
122    
123     C-- COMMON / THSICE_PAR_R / ice model (real) parameters
124     C thSIce_deltaT :: ice model time-step
125     C stressReduction :: reduction factor for wind-stress under sea-ice [0-1]
126     C thSIce_taveFreq :: Frequency^-1 for time-Aver. output (s)
127     C thSIce_diagFreq :: Frequency^-1 for diagnostic output (s)
128     C thSIce_monFreq :: Frequency^-1 for monitor output (s)
129     COMMON / THSICE_PAR_R /
130     & thSIce_deltaT,
131     & stressReduction,
132     & thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
133    
134     _RL thSIce_deltaT
135     _RL stressReduction
136     _RL thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
137    
138     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
139    
140     #endif /* ALLOW_THSICE */

  ViewVC Help
Powered by ViewVC 1.1.22