find the location of the max value over a lat-long domain

[Roger Shi]

Hi,

I am having problems to find the location of the max value over a
lat-long domain. GrADS only provides the "maxloc" function which only
loop through one dimension, e.g.

maxloc(expr, dim1, dim2 <,tinc>)

expr    - any valid GrADS expression
dim1    - the starting dimension expression
dim2    - the ending dimension expression
tinc    - optional time increment

Is there any way to loop through both x and y dimensions at the same time?

Best regards,
Roger Shi