FLEXPART-WRF crashes for array bound mismatch

[harish]

Hi
The flexpart-wrf version 3.3.2 crashes after about five hours of backward mode run when it is run using mpirun.

Using -fbacktrace and -fbounds-check option, I found that it crashes at line 132 in subroutine sendint_mpi.f90

if (tag.eq.1) npoint(jj2+1:numpart2)=mpi_npoint(1:chunksize2)

with error

Fortran runtime error: Array bound mismatch for dimension 1 of array 'npoint' (2500/2511)

My input file is similar to flexwrf.input.backward1 provided in the flexwrf_v31_testcases.tar.gz except for the release locations and dates.

Also, the par_mod.f90 file has been modified to accommodate larger grid size of the WRF data being used in following manner

integer, parameter :: naxmax=721, nymax=361,nuvzmax=64, nwzmax=64, nzmax=64

The code was run using LSF command bsub < mpijobfile where the mpijobfile is

mpirun -np 12 ./flexwrf33_gnu_mpi flexwrf.input.backward1

The code runs without crashing for serial mode, however grid_time file does not contain other than zero values in that case with all file size 372 bytes. I am not sure if this two issues are related or separate.

Any help/hint to resolve this issue will be highly appreciated. The input file being used is attached. I will be glad to provide more information.

Harish

[harish]

input file used for the run