Changeset a7167e4 in flexpart.git

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

Location:

src

Files:

• calcfluxes.f90 (modified) (3 diffs)
• netcdf_output_mod.f90 (modified) (1 diff, 1 prop)

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

src/netcdf_output_mod.f90

@@ -1655 +1655 @@
   call nf90_err(nf90_close(ncid))
 
+  ! Reinitialization of grid
+  !*************************
+
+  do ks=1,nspec
+    do kp=1,maxpointspec_act
+      do jy=0,numygrid-1
+        do ix=0,numxgrid-1
+          do kz=1,numzgrid
+            do nage=1,nageclass
+              do i=1,6
+                flux(i,ix,jy,kz,ks,kp,nage)=0.
+              end do
+            end do
+          end do
+        end do
+      end do
+    end do
+  end do
+
+
 end subroutine fluxoutput_netcdf