Changeset 0ecc1fe in flexpart.git for src/outgrid_init_nest.f90

Timestamp:

Nov 30, 2017, 4:04:54 PM (6 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:

d2a5a83

Parents:

c2bd55e

Message:

Changed handling of nested input fields to be consistent with non-nested

File:

• src/outgrid_init_nest.f90 (modified) (1 diff)

src/outgrid_init_nest.f90

@@ -69 +69 @@
   endif
 
+#ifdef USE_MPIINPLACE
+#else
 ! Extra field for totals at MPI root process
   if (lroot.and.mpi_mode.gt.0) then
-    ! allocate(griduncn0(0:numxgridn-1,0:numygridn-1,numzgrid,maxspec, &
-    !      maxpointspec_act,nclassunc,maxageclass),stat=stat)
-    ! if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
-
-    if (ldirect.gt.0) then
-      allocate(wetgriduncn0(0:numxgridn-1,0:numygridn-1,maxspec, &
-           maxpointspec_act,nclassunc,maxageclass),stat=stat)
-      if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
-      allocate(drygriduncn0(0:numxgridn-1,0:numygridn-1,maxspec, &
-           maxpointspec_act,nclassunc,maxageclass),stat=stat)
-      if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
-    endif
+! If MPI_IN_PLACE option is not used in mpi_mod.f90::mpif_tm_reduce_grid_nest(),
+! then an aux array is needed for parallel grid reduction
+    allocate(griduncn0(0:numxgridn-1,0:numygridn-1,numzgrid,maxspec, &
+         maxpointspec_act,nclassunc,maxageclass),stat=stat)
+    if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
+! allocate a dummy to avoid compilator complaints
+  else if (.not.lroot.and.mpi_mode.gt.0) then
+    allocate(griduncn0(1,1,1,1,1,1,1),stat=stat)
+  end if
+#endif
+!    if (ldirect.gt.0) then
+  if (lroot.and.mpi_mode.gt.0) then
+    allocate(wetgriduncn0(0:numxgridn-1,0:numygridn-1,maxspec, &
+         maxpointspec_act,nclassunc,maxageclass),stat=stat)
+    if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
+    allocate(drygriduncn0(0:numxgridn-1,0:numygridn-1,maxspec, &
+         maxpointspec_act,nclassunc,maxageclass),stat=stat)
+    if (stat.ne.0) write(*,*)'ERROR:could not allocate nested gridunc'
+!  endif
 ! allocate a dummy to avoid compilator complaints
   else if (.not.lroot.and.mpi_mode.gt.0) then