diff --git a/src/full_energy.c b/src/full_energy.c index c429f6e08..19d19d2bb 100644 --- a/src/full_energy.c +++ b/src/full_energy.c @@ -311,7 +311,12 @@ int full_energy(int gridcell, prepare_full_energy(iveg, Nveg, options.Nnode, all_vars, soil_con, moist0, ice0); /** Compute Bare (free of snow) Albedo **/ - bare_albedo = veg_var[iveg][0].albedo; + if (iveg!=Nveg){ + bare_albedo = veg_var[iveg][0].albedo; + } + else { + bare_albedo = BARE_SOIL_ALBEDO; + } /************************************* Compute the aerodynamic resistance