> > real, dimension(size(ug_temp,1),size(ug_temp,2),size(ug_temp,3)) :: Z
> > real, dimension(size(ug_temp,1),size(ug_temp,2),size(ug_temp+1,3)) :: radius
> Assuming that ug_temp is an array, the expression ug_temp+1
> is also an array with the same shape.
Right. So it doesn't matter if it says +1, +2, +3,....
> It will have a lower bound
> in each dimension of 1 (regardless of the declared lower bound
> of ug_temp). But, that's a very subtle point. It won't have any
> consequences here.
Right---interesting point, but not an issue here.
Another fine point, possibly relevant. At least in the past, if ug_temp
is floating-point, then one can't use an expression involving it in a
declaration. IIRC. (I think floating-point operations were disallowed
to make it easier for cross-compilers or something.) I believe at least
some aspects of this restriction were lifted recently.
Since you (Catherine) thought it was an obvious typo, what do YOU think
the code should do?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|