[6] | 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 | !------------------------------------------------------------------------------- |
---|