[ARPSSUPPORT] Possible getarps bug

[Yunheng Wang]

Eddie,

You are right. That is obviously a bug (typo).

Thanks for detecting it.

Yunheng.


-------- Original Message --------
Subject: [ARPSSUPPORT] Possible getarps bug
From: Eddie Natenberg <enatenberg at meso.com>
To: arpssupport at caps.ou.edu
CC: Tim Wait <twait at meso.com>
Date: 9/14/2010 4:00 PM

> 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
>
>
>