[ARPSSUPPORT] Possible getarps bug

Eddie Natenberg enatenberg at meso.com
Tue Sep 14 16:00:54 CDT 2010


Hi,

I have been testing ext2arps with arps data and found a possible bug in
zpsoil allocated as:

zpsoil_ext(nx_ext,ny_ext,nz_ext)

in several places, when I believe it should be:

zpsoil_ext(nx_ext,ny_ext,nzsoil_ext)

This seems to only be a problem using MPI because of the bcast call:

CALL mpupdater(zpsoil_ext,nx_ext*ny_ext*nz_ext)

on line 738 of src/ext2arps/getextd3d.f90

which should be:

CALL mpupdater(zpsoil_ext,nx_ext*ny_ext*nzsoil_ext)

This does not seem to be a problem when running without MPI support or
when using ext2arps data other than arps.

Thanks,
Eddie



-- 
Eddie Natenberg
Scientific Programmer MESO Inc.
(518) 283-5169



More information about the ARPSSUPPORT mailing list