1 |
heimbach |
1.2 |
C $Header: |
2 |
|
|
|
3 |
|
|
#ifdef ALLOW_SEAICE |
4 |
|
|
|
5 |
|
|
C /==========================================================\ |
6 |
|
|
C | SEAICE_PARAMS.h | |
7 |
|
|
C | o Basic parameter header for sea ice model. | |
8 |
|
|
C \==========================================================/ |
9 |
|
|
|
10 |
|
|
C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model. |
11 |
|
|
C |
12 |
dimitri |
1.4 |
C SEAICEwriteState - If true, write sea ice state to file; |
13 |
|
|
C default is false. |
14 |
|
|
C |
15 |
|
|
C SEAICEuseDYNAMICS - If false, do not use dynamics; |
16 |
|
|
C default is to use dynamics. |
17 |
|
|
C |
18 |
|
|
C SEAICEuseLSR - If true, use LSR rather than ADI solver; |
19 |
|
|
C default is ADI solver. |
20 |
|
|
C |
21 |
|
|
C SEAICEwindOnCgrid - By default wind files, uwind and vwind, |
22 |
|
|
C are defined on North-East B-grid U and V points. |
23 |
|
|
C When this flag is set, wind files are defined on |
24 |
|
|
C South-West C-grid U and V points. |
25 |
|
|
C |
26 |
|
|
LOGICAL SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseLSR, |
27 |
|
|
& SEAICEwindOnCgrid |
28 |
heimbach |
1.2 |
COMMON /SEAICE_PARM_L/ |
29 |
dimitri |
1.4 |
& SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseLSR, |
30 |
|
|
& SEAICEwindOnCgrid |
31 |
heimbach |
1.2 |
|
32 |
|
|
C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model. |
33 |
dimitri |
1.3 |
C LAD - time stepping used for sea-ice advection: |
34 |
|
|
C 1 = LEAPFROG, 2 = BACKWARD EULER. |
35 |
heimbach |
1.2 |
C IMAX_TICE - number of iterations for ice heat budget 10 |
36 |
|
|
C NPSEUDO - number of pseudo-timesteps used in dynsolver |
37 |
|
|
C DO PSEUDO-TIMESTEPS TO OBTAIN AN ACCURATE VISCOUS-PLASTIC SOLUTION |
38 |
|
|
C 10 PSEUDO-TIMESTEPS OR MORE ARE SUGGESTED FOR HIGH-RESOLUTION (~10KM) |
39 |
|
|
C 1 PSEUDO-TIMESTEP CAN BE USED FOR LOW-RESOLUTION GLOBAL MODELING |
40 |
|
|
C |
41 |
|
|
INTEGER LAD, IMAX_TICE, NPSEUDO |
42 |
|
|
COMMON /SEAICE_PARM_I/ LAD, IMAX_TICE, NPSEUDO |
43 |
|
|
|
44 |
|
|
C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters. |
45 |
dimitri |
1.4 |
C uwindFile - File containing |
46 |
|
|
C vwindFile - File containing |
47 |
|
|
C atempFile - File containing |
48 |
|
|
C aqhFile - File containing |
49 |
|
|
C lwfluxFile - File containing |
50 |
|
|
C swfluxFile - File containing |
51 |
|
|
C precipFile - File containing |
52 |
heimbach |
1.2 |
C evapFile - File containing |
53 |
|
|
C runoffFile - File containing |
54 |
|
|
C |
55 |
dimitri |
1.4 |
CHARACTER*(MAX_LEN_FNAM) uwindFile |
56 |
|
|
CHARACTER*(MAX_LEN_FNAM) vwindFile |
57 |
|
|
CHARACTER*(MAX_LEN_FNAM) atempFile |
58 |
|
|
CHARACTER*(MAX_LEN_FNAM) aqhFile |
59 |
|
|
CHARACTER*(MAX_LEN_FNAM) lwfluxFile |
60 |
|
|
CHARACTER*(MAX_LEN_FNAM) swfluxFile |
61 |
|
|
CHARACTER*(MAX_LEN_FNAM) precipFile |
62 |
heimbach |
1.2 |
CHARACTER*(MAX_LEN_FNAM) evapFile |
63 |
|
|
CHARACTER*(MAX_LEN_FNAM) runoffFile |
64 |
dimitri |
1.4 |
COMMON /SEAICE_PARM_C/ uwindFile, vwindFile, atempFile, aqhFile, |
65 |
|
|
& lwfluxFile, swfluxFile, precipFile, evapFile, runoffFile |
66 |
heimbach |
1.2 |
|
67 |
|
|
C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model. |
68 |
|
|
C SEAICE_dumpFreq - SEAICE dump frequency. (s) |
69 |
|
|
C SEAICE_taveFreq - SEAICE time-averaging frequency. (s) |
70 |
|
|
C SEAICE_drag - air-ice drag coefficient |
71 |
|
|
C SEAICE_waterDrag - water-ice drag coefficient * water density |
72 |
|
|
C SEAICE_albedo - summer albedo |
73 |
|
|
C SEAICE_drySnowAlb - dry snow albedo |
74 |
|
|
C SEAICE_wetSnowAlb - wet snow albedo |
75 |
|
|
C SEAICE_waterAlbedo - water albedo |
76 |
|
|
C SEAICE_strength - sea-ice strength |
77 |
|
|
C SEAICE_sensHeat - sensible heat transfer coefficient |
78 |
|
|
C SEAICE_latentWater - latent heat transfer coefficient for water |
79 |
|
|
C SEAICE_latentIce - latent heat transfer coefficient for ice |
80 |
|
|
C SEAICE_iceConduct - sea-ice conductivity |
81 |
|
|
C SEAICE_snowConduct - snow conductivity |
82 |
|
|
C SEAICE_emissivity - Stefan-Boltzman constant * emissivity |
83 |
|
|
C SEAICE_snowThick - cutoff snow thickness |
84 |
|
|
C SEAICE_shortwave - penetration shortwave radiation factor |
85 |
|
|
C LSR_ERROR - sets accuracy of LSR solver |
86 |
|
|
C DIFF1 - parameter used in advect.F |
87 |
|
|
C A22 - parameter used in growth.F |
88 |
|
|
C HO - parameter used in growth.F |
89 |
|
|
C WindForcingStart - Time of first wind forcing record (s) |
90 |
|
|
C WindForcingEnd - Time of last wind forcing record (s) |
91 |
|
|
C WindForcingPeriod - Period between wind forcing records (s) |
92 |
|
|
C FluxForcingStart - Time of first flux forcing record (s) |
93 |
|
|
C FluxForcingEnd - Time of last flux forcing record (s) |
94 |
|
|
C FluxForcingPeriod - Period between flux forcing records (s) |
95 |
|
|
C SSTForcingStart - Time of first SST forcing record (s) |
96 |
|
|
C SSTForcingEnd - Time of last SST forcing record (s) |
97 |
|
|
C SSTForcingPeriod - Period between SST forcing records (s) |
98 |
|
|
C SSSForcingStart - Time of first SSS forcing record (s) |
99 |
|
|
C SSSForcingEnd - Time of last SSS forcing record (s) |
100 |
|
|
C SSSForcingPeriod - Period between SSS forcing records (s) |
101 |
|
|
C StartingYear - Starting year of integration |
102 |
|
|
C EndingYear - Ending year of integration |
103 |
|
|
C |
104 |
|
|
_RL SEAICE_dumpFreq, SEAICE_taveFreq |
105 |
|
|
_RL SEAICE_drag, SEAICE_waterDrag, SEAICE_albedo |
106 |
|
|
_RL SEAICE_drySnowAlb, SEAICE_wetSnowAlb |
107 |
|
|
_RL SEAICE_waterAlbedo, SEAICE_strength |
108 |
|
|
_RL SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce |
109 |
|
|
_RL SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity |
110 |
|
|
_RL SEAICE_snowThick, SEAICE_shortwave |
111 |
|
|
_RL LSR_ERROR,DIFF1,A22,HO |
112 |
|
|
_RL WindForcingStart, WindForcingEnd, WindForcingPeriod |
113 |
|
|
_RL FluxForcingStart, FluxForcingEnd, FluxForcingPeriod |
114 |
|
|
_RL SSTForcingStart, SSTForcingEnd, SSTForcingPeriod |
115 |
|
|
_RL SSSForcingStart, SSSForcingEnd, SSSForcingPeriod |
116 |
|
|
_RL StartingYear, EndingYear |
117 |
|
|
COMMON /SEAICE_PARM_RL/ |
118 |
|
|
& SEAICE_dumpFreq, SEAICE_taveFreq, |
119 |
|
|
& SEAICE_drag, SEAICE_waterDrag, SEAICE_albedo, |
120 |
|
|
& SEAICE_drySnowAlb, SEAICE_wetSnowAlb, |
121 |
|
|
& SEAICE_waterAlbedo, SEAICE_strength, |
122 |
|
|
& SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce, |
123 |
|
|
& SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity, |
124 |
|
|
& SEAICE_snowThick, SEAICE_shortwave, |
125 |
|
|
& LSR_ERROR, DIFF1, A22, HO, |
126 |
|
|
& WindForcingStart, WindForcingEnd, WindForcingPeriod, |
127 |
|
|
& FluxForcingStart, FluxForcingEnd, FluxForcingPeriod, |
128 |
|
|
& SSTForcingStart, SSTForcingEnd, SSTForcingPeriod, |
129 |
|
|
& SSSForcingStart, SSSForcingEnd, SSSForcingPeriod, |
130 |
|
|
& StartingYear, EndingYear |
131 |
|
|
|
132 |
|
|
C-- COMMON /SEAICE_BOUND_RL/ Various bounding values |
133 |
|
|
C MAX_HEFF - maximum ice thickness (m) |
134 |
dimitri |
1.4 |
C MIN_ATEMP - minimum air temperature (deg C) |
135 |
|
|
C MIN_LWFLUX - minimum incoming longwave (W/m^2) |
136 |
heimbach |
1.2 |
C MAX_TICE - maximum ice temperature (deg C) |
137 |
|
|
C MIN_TICE - minimum ice temperature (deg C) |
138 |
|
|
C |
139 |
dimitri |
1.4 |
_RL MAX_HEFF, MIN_ATEMP, MIN_LWFLUX, MAX_TICE, MIN_TICE |
140 |
heimbach |
1.2 |
COMMON /SEAICE_BOUND_RL/ |
141 |
dimitri |
1.4 |
& MAX_HEFF, MIN_ATEMP, MIN_LWFLUX, MAX_TICE, MIN_TICE |
142 |
dimitri |
1.3 |
|
143 |
|
|
C-- Seaice time step |
144 |
|
|
COMMON /SEAICE_R/ SEAICE_DT |
145 |
|
|
_RL SEAICE_DT |
146 |
|
|
|
147 |
|
|
C-- Constants used by sea-ice model |
148 |
|
|
_RL ZERO , ONE , TWO |
149 |
|
|
parameter ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 ) |
150 |
|
|
_RL QUART , HALF |
151 |
|
|
parameter ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 ) |
152 |
heimbach |
1.2 |
|
153 |
|
|
#endif ALLOW_SEAICE |