| 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 |
Parent Directory
| 
Revision Log
| 
Revision Graph