Alter handling of big particles

metno · Apr 28, 2023 · 60f11e8 · 60f11e8
1 parent 74ff745
commit 60f11e8
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/common/drydep.f90 b/src/common/drydep.f90
@@ -147,7 +147,10 @@ subroutine drydep2(tstep, part)
     ! gas
       deprate = 1.0 - exp(-tstep*(0.008)/h)
     else if (def_comp(m)%radiusmym <= 10.0) then
-    ! particle 0.05<r
+    ! particle 0.05<r<10
+      deprate = 1.0 - exp(-tstep*(0.002+part%grv)/h)
+    else
+    ! particle r>=10.0
       deprate = 1.0 - exp(-tstep*(0.002+part%grv)/h)
     ! complete deposition when particle hits ground
       if (part%z == vlevel(1)) deprate = 1.