From 25b7b9457061c85b31f7dca3ac61d292f99fb0ef Mon Sep 17 00:00:00 2001
From: Bo Gao <83304414+gaobhub@users.noreply.github.com>
Date: Fri, 23 Feb 2024 13:57:25 -0500
Subject: [PATCH] Revise a bug for getting snow density (#239)

---
 .../land_cover/albedo_twocomponent_evaluator.cc                 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/pks/surface_balance/constitutive_relations/land_cover/albedo_twocomponent_evaluator.cc b/src/pks/surface_balance/constitutive_relations/land_cover/albedo_twocomponent_evaluator.cc
index 41002c4e2..eab547acf 100644
--- a/src/pks/surface_balance/constitutive_relations/land_cover/albedo_twocomponent_evaluator.cc
+++ b/src/pks/surface_balance/constitutive_relations/land_cover/albedo_twocomponent_evaluator.cc
@@ -92,7 +92,7 @@ AlbedoTwoComponentEvaluator::Evaluate_(const State& S, const std::vector<Composi
 
     for (auto c : lc_ids) {
       // albedo of the snow
-      albedo[1][c] = is_constant_snow_albedo_ ? a_snow_ : Relations::CalcAlbedoSnow(*(snow_dens)[0][c]); 
+      albedo[1][c] = is_constant_snow_albedo_ ? a_snow_ : Relations::CalcAlbedoSnow((*snow_dens)[0][c]); 
 
       double albedo_water =
         unfrozen_fraction[0][c] * a_water_ + (1 - unfrozen_fraction[0][c]) * a_ice_;