Changeset a1f4dd6 in flexpart.git

Timestamp:

Dec 1, 2015, 4:59:07 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:

7216707, 4c4261c

Parents:

6b22af9

Message:

Fixed a bug where receptor arrays were not properly zeroed

Location:

src

Files:

• mpi_mod.f90 (modified) (2 diffs)
• timemanager_mpi.f90 (modified) (3 diffs)

src/mpi_mod.f90

@@ -105 +105 @@
 ! If setting this to .false., numwfmem must be set to 3
 !===============================================================================
-  logical :: lmp_sync=.false. 
+  logical :: lmp_sync=.true. 
 !===============================================================================
 
@@ -202 +202 @@
       write(*,FMT='(80("#"))')
       write(*,*) '#### mpi_mod::mpif_init> WARNING: ', &
-           & 'numwfmem should be set to 2 for syncronous reading'
+           & 'numwfmem should be set to 2 for syncronous'
+      write(*,*) ' reading. Results will still be valid, but unneccesary '
+      write(*,*) 'amount of memory is being allocated.'
       write(*,FMT='(80("#"))')
 ! Force "syncronized" version if all processes will call getfields

src/timemanager_mpi.f90

@@ -459 +459 @@
               endif
             else
+! zero arrays on non-root processes
               gridunc(:,:,:,:,:,:,:)=0.
+              creceptor(:,:)=0.
             end if
           else 
@@ -471 +473 @@
               end if
             else
+! zero arrays on non-root processes
               gridunc(:,:,:,:,:,:,:)=0.
+              creceptor(:,:)=0.
             endif
           endif
@@ -788 +792 @@
 ! eso :TODO: this not implemented yet (transfer particles to PID 0 or rewrite)
 ! the tools to do this are already in mpi_mod.f90
+! :CUR: check this
   if (lroot) then 
     do j=1,numpart