Commit ae096108 authored by Jorn Bruggeman's avatar Jorn Bruggeman
Browse files

do not rely on PRODUCT returning 1 for empty array (not true for PGI)

parent 3bb23f30
......@@ -228,6 +228,9 @@ k__=45
#if _FABM_DIMENSION_COUNT_>0
domain_extent = (/ _LOCATION_ /)
interior_count = product(domain_extent)
#else
interior_count = 1
#endif
! Set defaults
......@@ -235,11 +238,10 @@ if (ntest == -1) then
if (mode == 1) then
ntest = 1
else
ntest = 50000000/product(domain_extent)
ntest = 50000000/interior_count
end if
end if
interior_count = product(domain_extent)
#ifdef _FABM_DEPTH_DIMENSION_INDEX_
horizontal_count = interior_count / domain_extent(_FABM_DEPTH_DIMENSION_INDEX_)
#else
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment