Changes in src/timemanager_mpi.f90 [0c8c7f2:20963b1] in flexpart.git

File:

• src/timemanager_mpi.f90 (modified) (6 diffs)

src/timemanager_mpi.f90

@@ -113 +113 @@
   integer :: j,ks,kp,l,n,itime=0,nstop,nstop1,memstat=0
 ! integer :: ksp
-  integer :: ip,irec
+  integer :: ip
   integer :: loutnext,loutstart,loutend
   integer :: ix,jy,ldeltat,itage,nage,idummy
@@ -129 +129 @@
 ! Measure time spent in timemanager
   if (mp_measure_time) call mpif_mtime('timemanager',0)
-
 
 ! First output for time 0
@@ -533 +532 @@
                 end if
 
-              else
+              else  ! :TODO: check for zeroing in the netcdf module
                 call concoutput_surf_nest(itime,outnum)
               end if
@@ -594 +593 @@
 46      format(' Simulated ',f7.1,' hours (',i13,' s), ',i13, ' particles')
         if (ipout.ge.1) then
-          if (mp_measure_time) call mpif_mtime('iotime',0)
-          irec=0
           do ip=0, mp_partgroup_np-1
-            if (ip.eq.mp_partid) then
-              if (mod(itime,ipoutfac*loutstep).eq.0) call partoutput(itime) ! dump particle positions
-              if (ipout.eq.3) call partoutput_average(itime,irec) ! dump particle positions
-            endif
-            if (ipout.eq.3) irec=irec+npart_per_process(ip)
+            if (ip.eq.mp_partid) call partoutput(itime) ! dump particle positions
             call mpif_mpi_barrier
           end do
-          if (mp_measure_time) call mpif_mtime('iotime',1)
         end if
 
@@ -765 +757 @@
         if (mp_measure_time) call mpif_mtime('advance',1)
 
-  ! Calculate average position for particle dump output
-  !****************************************************
-
-        if (ipout.eq.3) call partpos_average(itime,j)
-
 
 ! Calculate the gross fluxes across layer interfaces
@@ -908 +895 @@
     do ip=0, mp_partgroup_np-1
       if (ip.eq.mp_partid) then 
-        if (mp_dbg_mode) write(*,*) 'call partoutput(itime), proc, mp_partid',ip,mp_partid
+        !if (mp_dbg_mode) write(*,*) 'call partoutput(itime), proc, mp_partid',ip,mp_partid
         call partoutput(itime)    ! dump particle positions
       end if