Changeset 38b7917 in flexpart.git for src/init_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/init_domainfill_mpi.f90 (modified) (2 diffs)

src/init_domainfill_mpi.f90

@@ -163 +163 @@
       pp(1)=rho(ix,jy,1,1)*r_air*tt(ix,jy,1,1)
       pp(nz)=rho(ix,jy,nz,1)*r_air*tt(ix,jy,nz,1)
-      colmass(ix,jy)=(pp(1)-pp(nz))/ga*gridarea(jy)
+      ! Each MPI process is assigned an equal share of particles
+      colmass(ix,jy)=(pp(1)-pp(nz))/ga*gridarea(jy)/mp_partgroup_np
       colmasstotal=colmasstotal+colmass(ix,jy)
+
     end do
   end do
 
-  write(*,*) 'Atm. mass: ',colmasstotal
+  if (lroot) write(*,*) 'Atm. mass: ',colmasstotal
 
 
@@ -415 +417 @@
 !***************************************************************************
 
+! :TODO: eso: parallelize
   if (ipin.eq.1) then
     open(unitboundcond,file=path(2)(1:length(2))//'boundcond.bin', &