diff --git a/seims/preprocess/sp_terrain.py b/seims/preprocess/sp_terrain.py
index 361a1e814..1dea26554 100644
--- a/seims/preprocess/sp_terrain.py
+++ b/seims/preprocess/sp_terrain.py
@@ -490,8 +490,8 @@ def parameters_extraction(cfg):
 
         status_output('Generating depression storage capacity...', 20, f)
         TerrainUtilClass.depression_capacity(cfg.maindb, cfg.spatials.landuse,
-                                             cfg.spatials.soil_texture,
                                              cfg.spatials.slope,
+                                             cfg.spatials.soil_texture,
                                              cfg.spatials.depression,
                                              cfg.imper_perc_in_urban)
 
diff --git a/seims/src/seims_main/modules/hydrology_longterm/DEP_LINSLEY/DepressionLinsley.cpp b/seims/src/seims_main/modules/hydrology_longterm/DEP_LINSLEY/DepressionLinsley.cpp
index c6af1a937..38d4ed207 100644
--- a/seims/src/seims_main/modules/hydrology_longterm/DEP_LINSLEY/DepressionLinsley.cpp
+++ b/seims/src/seims_main/modules/hydrology_longterm/DEP_LINSLEY/DepressionLinsley.cpp
@@ -80,6 +80,10 @@ int DepressionFSDaily::Execute() {
                 m_ed[i] = m_pet[i] - m_ei[i];
             } else {
                 m_ed[i] = m_sd[i];
+            }
+            if (m_ed[i] < 0) {
+                printf("[DepressionLinsley] Warning! m_ed[%d](%f) < 0! m_pet[%d](%f) m_ei[%d](%f)\n", i, m_ed[i], i, m_pet[i], i, m_ei[i]);
+                m_ed[i] = 0;
             }
 			// �ݵ���ˮ��� - �ݵ��������
             m_sd[i] -= m_ed[i];