Apply precip scaling to computatational domain

to avoid discontinuities.
pism · Nov 22, 2017 · 6b96570 · 6b96570
1 parent 162a38b
commit 6b96570
Showing 1 changed file with 0 additions and 3 deletions.
diff --git a/src/coupler/util/PLapseRates.hh b/src/coupler/util/PLapseRates.hh
@@ -166,8 +166,6 @@ protected:
     for (Points p(*Mod::m_grid); p; p.next()) {
       const int i = p.i(), j = p.j();
 
-      if (thk(i,j) > 0) {
-
         //as in paleo_precip
         const double correction = exp(-lapse_rate * (surface(i,j) - m_reference_surface(i,j)));
 
@@ -178,7 +176,6 @@ protected:
         //const double correction = pow(2.0, -lapse_rate/0.05/10.0 * ((*surface)(i,j) - m_reference_surface(i,j)));
 
         result(i,j) *= correction;
-      }
     }
   }