Changeset 7999df47 in flexpart.git for src/init_domainfill_mpi.f90

Timestamp:

Mar 8, 2016, 3:28:18 PM (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:

c2162ce

Parents:

9b53903

Message:

For parallel domain-filling version, changed info written to screen more consistent with serial version

File:

• src/init_domainfill_mpi.f90 (modified) (2 diffs)

src/init_domainfill_mpi.f90

@@ -181 +181 @@
     ylat=ylat0+real(jy)*dy
     do ix=nx_we(1),nx_we(2)      ! loop about longitudes
-      ncolumn=nint(0.999*real(npart(1))*colmass(ix,jy)/ &
+      ncolumn=nint(0.999*real(npart(1)/mp_partgroup_np)*colmass(ix,jy)/ &
            colmasstotal)
       if (ncolumn.eq.0) goto 30
@@ -341 +341 @@
 
   fractus=real(numcolumn)/real(nz)
-  write(*,*) 'Total number of particles at model start: ',numpart
-  write(*,*) 'Maximum number of particles per column: ',numcolumn
-  write(*,*) 'If ',fractus,' <1, better use more particles'
+  write(*,*) 'Total number of particles at model start: ',numpart*mp_partgroup_np
+  write(*,*) 'Maximum number of particles per column: ',numcolumn*mp_partgroup_np
+  write(*,*) 'If ',fractus*mp_partgroup_np,' <1, better use more particles'
   fractus=sqrt(max(fractus,1.))/2.
 
   do jy=ny_sn(1),ny_sn(2)      ! loop about latitudes
     do ix=nx_we(1),nx_we(2)      ! loop about longitudes
-      ncolumn=nint(0.999/fractus*real(npart(1))*colmass(ix,jy) &
+      ncolumn=nint(0.999/fractus*real(npart(1)/mp_partgroup_np)*colmass(ix,jy) &
            /colmasstotal)
       if (ncolumn.gt.maxcolumn) stop 'maxcolumn too small'