Changeset 38b7917 in flexpart.git for src/boundcond_domainfill_mpi.f90

Timestamp:

Mar 8, 2016, 9:54:38 AM (8 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:

9b53903

Parents:

db712a8

Message:

Parallelization of domain fill option (save/restart not implemented yet)

File:

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

src/boundcond_domainfill_mpi.f90

@@ -192 +192 @@
         rhox=(rhohl(1)*dt2+rhohl(2)*dt1)*dtt
 
-  ! Calculate mass flux
-  !********************
-
-        fluxofmass=windx*rhox*boundarea*real(lsynctime)
+  ! Calculate mass flux, divided by number of processes
+  !****************************************************
+
+        fluxofmass=windx*rhox*boundarea*real(lsynctime)/mp_partgroup_np
 
 
@@ -425 +425 @@
         rhox=(rhohl(1)*dt2+rhohl(2)*dt1)*dtt
 
-  ! Calculate mass flux
-  !********************
-
-        fluxofmass=windx*rhox*boundarea*real(lsynctime)
+  ! Calculate mass flux, divided by number of processes
+  !****************************************************
+
+        fluxofmass=windx*rhox*boundarea*real(lsynctime)/mp_partgroup_np
 
   ! If the mass flux is directed into the domain, add it to previous mass fluxes;
@@ -587 +587 @@
   !*****************************************************************************
 
+! :TODO: eso parallelize
   if ((ipout.gt.0).and.(itime.eq.loutend)) then
     open(unitboundcond,file=path(2)(1:length(2))//'boundcond.bin', &