1 | ! |
---|
2 | ! ice.F CLEANING REQUIRED |
---|
3 | ! |
---|
4 | |
---|
5 | ! tfsn Melting point temperature of the snow |
---|
6 | ! tfsg Melting point temperature of the ice |
---|
7 | ! xkn Conductivity of the snow |
---|
8 | ! xkg Conductivity of the ice |
---|
9 | ! rcpn Density times specific heat for the snow |
---|
10 | ! rcpg Density times specific heat for the ice |
---|
11 | ! rhog Density of the ice |
---|
12 | ! rhon Density of the snow |
---|
13 | ! emig Emissivity of the ice |
---|
14 | ! sglace Salinity of the ice |
---|
15 | ! hmelt Maximum melting at the bottom |
---|
16 | ! acrit(2) Minimum fraction for leads |
---|
17 | ! hgcrit(2) Ice thickness for lateral accretion |
---|
18 | ! hgmin Ice thickness corr. to max. energy stored in brine pocket |
---|
19 | ! hndif Computation of temp. in snow or not |
---|
20 | ! hgdif Computation of temp. in ice or not |
---|
21 | ! hglim Minimum ice thickness |
---|
22 | ! amax Maximum lead fraction |
---|
23 | ! uscomi =1.0/(1.0-amax) |
---|
24 | ! beta Numerical caracteritic of the scheme for diffusion in ice |
---|
25 | ! ddtb Time step for ice thermodynamics (s) |
---|
26 | ! swiqst Energy stored in brine pocket or not |
---|
27 | ! parlat Percentage of energy used for lateral ablation |
---|
28 | ! hakspl Slope of distr. for Hakkinen-Mellor's lateral melting |
---|
29 | ! hibspl Slope of distribution for Hibler's lateral melting |
---|
30 | ! exld Exponent for leads-closure rate |
---|
31 | ! hakdif Coefficient for diffusions of ice and snow |
---|
32 | ! hth Threshold thickness for comp. of eq. thermal conductivity |
---|
33 | ! hnzst Thickness of the surf. layer in temp. computation |
---|
34 | ! parsub Switch for snow sublimation or not |
---|
35 | ! cnscg ratio rcpn/rcpg |
---|
36 | ! nbits Number of time steps in Newton -Raphson procedure |
---|
37 | ! stefan Stefan-Boltzman constant |
---|
38 | ! vkarmn von Karman constant |
---|
39 | ! cevap Latent heat of evaporation of water |
---|
40 | ! zemise Emissivity of water |
---|
41 | ! rhoesn 1/rhon |
---|
42 | ! firg IR flux over the ice (only used for outputs) |
---|
43 | ! fcsg Sensible heat flux over the ice (only used for outputs) |
---|
44 | ! fleg Latent heat flux over the ice (only used for outputs) |
---|
45 | ! ts Surface temperature of the ice |
---|
46 | ! tfu Melting point temperature of sea water |
---|
47 | ! hnbq Snow thickness |
---|
48 | ! hgbq Ice thickness |
---|
49 | ! albq Leads fraction |
---|
50 | ! qstobq Energy stored in the brine pockets |
---|
51 | ! fbbq Heat flux at the ice base |
---|
52 | ! tbq Temperature inside the ice/snow layer |
---|
53 | ! dmnbq Variation of snow mass |
---|
54 | ! dmgbq Variation of ice mass |
---|
55 | ! qlbq heat balance of the lead (or of the open ocean) |
---|
56 | ! qcmbq Energy needed to bring the ocean surface layer until its freezing |
---|
57 | ! point (at a factor 2) |
---|
58 | ! thcm part of the solar energy used in the lead heat budget |
---|
59 | ! fstrbq Solar flux transmitted trough the ice |
---|
60 | ! ffltbq Array linked with the max heat contained in brine pockets (?) |
---|
61 | ! fscmbq Linked with the solar flux below the ice (?) |
---|
62 | ! fsbbq Also linked with the solar flux below the ice (?) |
---|
63 | ! qfvbq Array used to store energy in case of toral lateral ablation (?) |
---|
64 | ! xzo rugosity of the ice (no more used) |
---|
65 | ! dmgwi Variation of the mass of snow ice |
---|
66 | ! psbq Surface air pressure |
---|
67 | ! tabq Surface air temperature |
---|
68 | ! qabq Surface air humidity |
---|
69 | ! vabq Surface wind velocity |
---|
70 | ! hnplbq Snow precipitation |
---|
71 | ! fevabq Evaporation flux |
---|
72 | ! fsolcn Solar flux at the ocean surface |
---|
73 | ! fsolg Solar flux at the ice surface |
---|
74 | ! flecn Latent heat flux at the ocean surface |
---|
75 | ! fcscn Sensible heat flux at the ocean surface |
---|
76 | ! tenagx Wind stress at the ice surface (x) |
---|
77 | ! tenagy Wind stress at the ice surface (y) |
---|
78 | ! albg 03/08/2001 albedo obtenu du forcage simip2 |
---|
79 | ! albege Albedo of the snow or ice (only for outputs) |
---|
80 | ! tairox Wind stress at the ocean surface (x) |
---|
81 | ! tairoy Wind stress at the ocean surface (y) |
---|
82 | ! ratbqg Longwave downward radiation flux over the ice |
---|
83 | ! ratbqo Longwave downward radiation flux over the ocean |
---|
84 | ! cloud Cloud fraction |
---|
85 | ! tdew Air relative humidity |
---|
86 | ! albecn Albedo of the ocean (only for outputs) |
---|
87 | ! tauc Cloud optical depth |
---|
88 | ! runoff river runoff |
---|
89 | ! sdvt u*^2/(Stress/density) |
---|
90 | ! fcm1 Solar flux at the ocean surface |
---|
91 | ! fcm2 Non-solar flux at the ocean surface |
---|
92 | ! fwat Freshwater flux (change of definition between the routines) |
---|
93 | ! reslum Relative absorption of solar radiation in each ocean level |
---|
94 | ! |
---|
95 | !--COMMON blocs : |
---|
96 | !------------------------------------------------------------------------------ |
---|
97 | COMMON /ice_grid/ |
---|
98 | & deltaz_i_phy(maxnlay) , !: thicknesses of the physical ice layers |
---|
99 | & z_i_phy(maxnlay) , !: cotes of the physical ice layers |
---|
100 | & zb_i_phy(0:maxnlay) , !: cotes of the interfaces of the physical ice layers |
---|
101 | & deltaz_s_phy(maxnlay) , !: thicknesses of the physical snow layers |
---|
102 | & z_s_phy(maxnlay) !: cotes of the physical snow layers |
---|
103 | |
---|
104 | COMMON / ice_constants / |
---|
105 | & tfsn,tfsg,xkn,xkg,rcpn,rcpg,rhog,rhon, |
---|
106 | & emig,sglace,hmelt,acrit(2),hgcrit(2),hgmin,hndif, |
---|
107 | & hgdif,hglim,amax,uscomi,beta,ddtb,swiqst,parlat, |
---|
108 | & hakspl,hibspl,exld,hakdif,hth,hnzst,parsub,cnscg,nbits |
---|
109 | |
---|
110 | COMMON / fluxsf / |
---|
111 | & stefan,lvap,lsub,vkarmn,cevap,zemise,rhoesn |
---|
112 | |
---|
113 | COMMON / comdia / |
---|
114 | & firg(imax,jmax),fcsg(imax,jmax),fleg(imax,jmax) |
---|
115 | |
---|
116 | COMMON / comban / |
---|
117 | & ts(imax,jmax),tfu(imax,jmax),hnbq(imax,jmax), |
---|
118 | & hgbq(imax,jmax),albq(imax,jmax), |
---|
119 | & qstobq(imax,jmax),fbbq(imax,jmax), |
---|
120 | & dmnbq(imax,jmax),dmgbq(imax,jmax), |
---|
121 | & qlbq(imax,jmax),qcmbq(imax,jmax),thcm(imax,jmax), |
---|
122 | & fstrbq(imax,jmax),ffltbq(imax,jmax),fscmbq(imax,jmax), |
---|
123 | & fsbbq(imax,jmax),qfvbq(imax,jmax),xzo(imax,jmax), |
---|
124 | & dmgwi(imax,jmax),total(imax,jmax) |
---|
125 | c |
---|
126 | COMMON / comfor / |
---|
127 | & psbq(imax,jmax),tabq(imax,jmax), |
---|
128 | & qabq(imax,jmax),vabq(imax,jmax), |
---|
129 | & hnplbq(imax,jmax),fevabq(imax,jmax),fsolcn(imax,jmax), |
---|
130 | & hnpbq(imax,jmax), |
---|
131 | & fsolg(imax,jmax),flecn(imax,jmax),fcscn(imax,jmax), |
---|
132 | & tenagx(imax,jmax),tenagy(imax,jmax),albg(imax,jmax), |
---|
133 | & albege(imax,jmax),tairox(imax,jmax),tairoy(imax,jmax), |
---|
134 | & ratbqg(imax,jmax),ratbqo(imax,jmax),cloud(imax,jmax), |
---|
135 | & opt_dept, |
---|
136 | & tdew(imax,jmax), |
---|
137 | & albecn(imax,jmax),tauc(imax,jmax),runoff(imax,jmax), |
---|
138 | & sdvt(imax,jmax),fsolg2(imax,jmax) |
---|
139 | |
---|
140 | COMMON / comca / |
---|
141 | & fcm1(imax,jmax),fcm2(imax,jmax), |
---|
142 | & fwat(imax,jmax), |
---|
143 | & reslum(imax,jmax,0:kmax+1) |
---|
144 | |
---|
145 | c global characteristics of the ice pack |
---|
146 | COMMON / ice_global / |
---|
147 | & t_i(imax,jmax,maxnlay),t_s(imax,jmax,maxnlay),t_su(imax,jmax), |
---|
148 | & t_bo(imax,jmax),ht_s(imax,jmax),ht_i(imax,jmax),ht_w(imax,jmax), |
---|
149 | & s_i(imax,jmax,maxnlay) |
---|
150 | |
---|
151 | COMMON / heat_fluxes / |
---|
152 | & fc_int |
---|
153 | |
---|
154 | LOGICAL:: ln_evap |
---|
155 | |
---|
156 | LOGICAL:: ln_flo, ln_grd, ln_flu |
---|
157 | |
---|
158 | COMMON/snow_phy/ln_evap |
---|
159 | |
---|
160 | integer :: layer,layer_a,nconv,numofday |
---|
161 | integer :: modul0, modul1a, modul2, modul3, modul4 |
---|
162 | integer :: modul1b,modul5,nbot0,nbot1,ntop0,ntop1 |
---|
163 | |
---|
164 | c Number of layers in the ice and snow |
---|
165 | integer :: nlayi0, nlayi1, nlays0, nlays1 |
---|
166 | integer :: nlay_i, nlay_s |
---|
167 | |
---|
168 | COMMON/layers/ n_i, n_s |
---|
169 | |
---|
170 | COMMON/simip/zsim(0:maxnlay),tempsim(maxnlay),salsim(maxnlay), |
---|
171 | & tempint(maxnlay),salint(maxnlay), |
---|
172 | & thick0(maxnlay), |
---|
173 | & thick1(maxnlay), |
---|
174 | & hsold,hgold,hsnew,hgnew |
---|
175 | |
---|
176 | COMMON/vertres/zm0(0:maxnlay),zm1(0:maxnlay),qm0(maxnlay), |
---|
177 | & qm1(0:maxnlay+2), sal_new_layer |
---|
178 | |
---|
179 | COMMON/heateqcoe/lfus,cpg,cpw,gammac,betak,tmut,tpw, |
---|
180 | & cpoc,deltah, betak1, betak2, visc_br, beta_ocs |
---|
181 | |
---|
182 | COMMON/barrowconf/sal_read(11), hi_read(11), hgins, hnins, |
---|
183 | & tsuins, oce_sal, oce_flx, num_sal, nday1, |
---|
184 | & i_sal |
---|
185 | |
---|
186 | COMMON/snowprecip/numd_sn1, numd_sn2, numd_sn3, |
---|
187 | & sn_prec_1, sn_prec_2 |
---|
188 | |
---|
189 | COMMON/tuneforcing/sf_mult, tabq_ano |
---|
190 | |
---|
191 | COMMON/fluidtpt/flu_beta, rad_io, |
---|
192 | & frtr_si_phy, qsummer, d_br_mol, d_br_tur, |
---|
193 | & ra_c, ra_smooth, e_thr_flu, delta_cw, ini_swi, |
---|
194 | & s_ini, ln_flo, ln_flu, ln_grd |
---|
195 | |
---|
196 | CHARACTER*2 gravdr |
---|
197 | CHARACTER*4 alb_char, pre_char, sal_char |
---|
198 | CHARACTER*12 name_file_bar_par |
---|
199 | CHARACTER*17 name_file_bar_for |
---|
200 | CHARACTER*4 c_rad_scheme |
---|
201 | CHARACTER*3 c_rad_discr |
---|
202 | |
---|
203 | COMMON/chars/gravdr, alb_char, pre_char, sal_char, |
---|
204 | & name_file_bar_par, namefile_bar_for, |
---|
205 | & c_rad_scheme, c_rad_discr |
---|
206 | |
---|
207 | ! radiation transfer in sea ice |
---|
208 | REAL(8) :: |
---|
209 | & rad_inot_s_dry, rad_inot_s_wet, |
---|
210 | & rad_inot_i_dry, rad_inot_i_wet, |
---|
211 | & h_not_s, h_not_i, |
---|
212 | & fpar_fsw, qpar_fsw, |
---|
213 | & rad_kappa_s_su_d, rad_kappa_s_su_m, |
---|
214 | & rad_kappa_s_de_d, rad_kappa_s_de_m, |
---|
215 | & rad_kappa_i_su_d, rad_kappa_i_su_m, |
---|
216 | & rad_kappa_i_de_d, rad_kappa_i_de_m |
---|
217 | |
---|
218 | COMMON/icerad/ |
---|
219 | & rad_inot_s_dry, rad_inot_s_wet, |
---|
220 | & rad_inot_i_dry, rad_inot_i_wet, |
---|
221 | & h_not_s, h_not_i, |
---|
222 | & fpar_fsw, qpar_fsw, |
---|
223 | & rad_kappa_s_su_d, rad_kappa_s_su_m, |
---|
224 | & rad_kappa_s_de_d, rad_kappa_s_de_m, |
---|
225 | & rad_kappa_i_su_d, rad_kappa_i_su_m, |
---|
226 | & rad_kappa_i_de_d, rad_kappa_i_de_m |
---|
227 | |
---|
228 | CHARACTER*5 c_bio_model |
---|
229 | COMMON/limbio/ |
---|
230 | & c_bio_model |
---|
231 | ! |
---|
232 | !--fin du fichier "ice.com" |
---|
233 | !------------------------------------------------------------------------------- |
---|