@certik 
Yeah, I recall that it needed to be zero. Thank you for this PR! (and fortran90.org which I try to use as much as I can). I never got my head quite around assumed-shape arrays,
usually I'm interfacing older F77 code so passing dimensions is something I just tend to
do without thinking.