Changeset a7167e4 in flexpart.git for src/calcfluxes.f90

Timestamp:

Sep 9, 2020, 10:29:08 AM (4 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

dev

Children:

03adec6

Parents:

969c1fe

Message:

bugfix: netcdffluxfield init, flux through upper, southern and western boundary

File:

• src/calcfluxes.f90 (modified) (3 diffs)

src/calcfluxes.f90

@@ -61 +61 @@
        (jyave.le.numygrid-1)) then
     do kz=1,numzgrid                ! determine height of cell
-      if (outheighthalf(kz).gt.zold) goto 11
+!      if (outheighthalf(kz).gt.zold) goto 11
+      if (outheight(kz).gt.zold) goto 11 !sec, use upper layer instead of mid layer
     end do
 11   k1=min(numzgrid,kz)
     do kz=1,numzgrid                ! determine height of cell
-      if (outheighthalf(kz).gt.ztra1(jpart)) goto 21
+!      if (outheighthalf(kz).gt.ztra1(jpart)) goto 21
+      if (outheight(kz).gt.ztra1(jpart)) goto 21
     end do
 21   k2=min(numzgrid,kz)
@@ -93 +95 @@
 
     if (abs(xold-xtra1(jpart)).lt.real(nx)/2.) then
-      ix1=int((xold*dx+xoutshift)/dxout+0.5)
-      ix2=int((xtra1(jpart)*dx+xoutshift)/dxout+0.5)
+      ix1=int((xold*dx+xoutshift)/dxout)              ! flux throught the western boundary (sec)
+      ix2=int((xtra1(jpart)*dx+xoutshift)/dxout)      ! flux throught the western boundary (sec)
+!      ix1=int((xold*dx+xoutshift)/dxout+0.5) 
+!      ix2=int((xtra1(jpart)*dx+xoutshift)/dxout+0.5)
       do k=1,nspec
         do ix=ix1,ix2-1
@@ -142 +146 @@
   if ((kzave.le.numzgrid).and.(ixave.ge.0).and. &
        (ixave.le.numxgrid-1)) then
-    jy1=int((yold*dy+youtshift)/dyout+0.5)
-    jy2=int((ytra1(jpart)*dy+youtshift)/dyout+0.5)
+    jy1=int((yold*dy+youtshift)/dyout)         ! flux throught the southern boundary (sec)
+    jy2=int((ytra1(jpart)*dy+youtshift)/dyout) ! flux throught the southern boundary (sec)
+!    jy1=int((yold*dy+youtshift)/dyout+0.5)
+!    jy2=int((ytra1(jpart)*dy+youtshift)/dyout+0.5)
 
     do k=1,nspec