Changes in src/drydepokernel.f90 [4c64400:e200b7a] in flexpart.git

File:

• src/drydepokernel.f90 (modified) (4 diffs)

src/drydepokernel.f90

@@ -40 +40 @@
   !                                                                            *
   !*****************************************************************************
-  ! Changes:
-  ! eso 10/2016: Added option to disregard kernel 
-  ! 
-  !*****************************************************************************
-
 
   use unc_mod
@@ -52 +47 @@
   implicit none
 
-  real(dep_prec), dimension(maxspec) :: deposit
-  real :: x,y,ddx,ddy,xl,yl,wx,wy,w
+  real :: x,y,deposit(maxspec),ddx,ddy,xl,yl,wx,wy,w
   integer :: ix,jy,ixp,jyp,ks,nunc,nage,kp
 
@@ -80 +74 @@
   endif
 
-  ! If no kernel is used, direct attribution to grid cell
-  !******************************************************
-
-  if (lnokernel) then
-    do ks=1,nspec
-      if ((abs(deposit(ks)).gt.0).and.DRYDEPSPEC(ks)) then
-        if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
-             (jy.le.numygrid-1)) then
-          drygridunc(ix,jy,ks,kp,nunc,nage)= &
-               drygridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)
-        end if
-      end if
-    end do
-  else ! use kernel 
-
 
   ! Determine mass fractions for four grid points
   !**********************************************
-  do ks=1,nspec
+    do ks=1,nspec
 
-   if ((abs(deposit(ks)).gt.0).and.DRYDEPSPEC(ks)) then
+    if ((abs(deposit(ks)).gt.0).and.DRYDEPSPEC(ks)) then
 
    if ((ix.ge.0).and.(jy.ge.0).and.(ix.le.numxgrid-1).and. &
         (jy.le.numygrid-1)) then
-     w=wx*wy
-     drygridunc(ix,jy,ks,kp,nunc,nage)= &
-          drygridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w
-     continue
-   endif
+    w=wx*wy
+      drygridunc(ix,jy,ks,kp,nunc,nage)= &
+           drygridunc(ix,jy,ks,kp,nunc,nage)+deposit(ks)*w
+      continue
+  endif
 
   if ((ixp.ge.0).and.(jyp.ge.0).and.(ixp.le.numxgrid-1).and. &
@@ -133 +112 @@
   endif
 
-  end do
-end if
+    end do
 
 end subroutine drydepokernel