Changeset 4c64400 in flexpart.git for src/drydepokernel.f90

Timestamp:

Nov 8, 2016, 4:42:27 PM (7 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

c8fc724

Parents:

16b61a5

Message:

Bugfix for double precision dry deposition calculation. Added (hardcoded) option to not use output kernel. Version/date string updated.

File:

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

src/drydepokernel.f90

@@ -40 +40 @@
   !                                                                            *
   !*****************************************************************************
+  ! Changes:
+  ! eso 10/2016: Added option to disregard kernel 
+  ! 
+  !*****************************************************************************
+
 
   use unc_mod
@@ -47 +52 @@
   implicit none
 
-  real :: x,y,deposit(maxspec),ddx,ddy,xl,yl,wx,wy,w
+  real(dep_prec), dimension(maxspec) :: deposit
+  real :: x,y,ddx,ddy,xl,yl,wx,wy,w
   integer :: ix,jy,ixp,jyp,ks,nunc,nage,kp
 
@@ -74 +80 @@
   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. &
@@ -112 +133 @@
   endif
 
-    end do
+  end do
+end if
 
 end subroutine drydepokernel