[gradsusr] Error: Axis sizes not the same

Tue Jul 12 15:00:27 EDT 2016

Jennifer,

Thank you, that makes sense! Your example was also successful.

If I wanted to find the maximum value of this masked, linear interpolated
areal average, how would I go about doing this?

I tried sticking in 'amax' before the aave, but an error prompted that it
was not a variable or function name. Sticking in a 'max' before the aave
yielded an error that said I had too many or too few arguments (I put in
starting and ending dimensions as well, but to no avail).

I also tried to use the 'amax' function without the aave since it already
accounts for the area:

define maximum=amax((maskout(prectot,lterp(mask.1(time=00:30z01feb1980),
prectot)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)

'd maximum'

However, I received this error:
Data Request Warning:  Request is completely outside file limits
  Entire grid contents are set to missing data
  Grid limits of file:     X = 1 361  Y = 1 180  Z = 1 1  T = 1 143  E = 1
1
  Grid limits of request:  X = -179 181  Y = 1 181  Z = 1 1  T = 0 0  E = 1
1
  Warning issued for variable = prectot.2

I'm not too sure where I went wrong.

-Justin

On Fri, Jul 8, 2016 at 9:42 AM, Jennifer M Adams <jadams21 at gmu.edu> wrote:

> Hi, Justin —
> If both grids have a lat/lon increment of 1.0 degree, but one starts at
> -180 and the other starts at 0, GrADS sees those as different and you have
> to use lterp() to perform operations between those grids.
>
> Your expression
>
> aave(maskout(prectot.2(x=1,y=1),mask.1(x=1,y=1
>>>> ,time=00:30z01feb1980)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)
>>>>
>>> ...isn’t quite right. You need to fix the X and Y dimensions *before*
> you invoke your aave() command, not as a local dimension override within
> your expression. Here’s an untested example of how it might look:
>
> set dfile 2
> set x 1
> set y 1
> set t 1
> d
> aave(maskout(prectot,lterp(mask.1(time=00:30z01feb1980),prectot)),lon=-179.5,
> … , … , lat=89.5)
>
>
> —Jennifer
>
>
> On Jul 7, 2016, at 4:43 PM, Justin Hicks <jhicks2014 at gmail.com> wrote:
>
> Jeff,
>
> This seems to work. Both of the files are gridded 1 degree, so I'm
> surprised they didn't work without this function.
>
> -Justin
>
> On Thu, Jul 7, 2016 at 4:09 PM, Justin Hicks <jhicks2014 at gmail.com> wrote:
>
>> Jeff,
>>
>> Thanks, I'll try the lterp function. I tried changing the xdef and ydef
>> before but I still had the same issue.
>>
>> -Justin
>>
>> On Thu, Jul 7, 2016 at 4:03 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:
>>
>>> Justin,
>>> As in my previous response to another thread, use the lterp function -
>>> it performs linear interpolation to match two grids. It is a fairly simple
>>> function. If for whatever reason you don't want to do that or it doesn't
>>> work, you can change the XDEF and YDEF lines of your control files to match
>>> each other, as that information is what GrADS uses to determine if you can
>>> perform operations on two grids. The PDEF entries do not have to match.
>>>
>>> Jeff Duda
>>>
>>> On Thu, Jul 7, 2016 at 2:38 PM, Justin Hicks <jhicks2014 at gmail.com>
>>> wrote:
>>>
>>>> GrADS users,
>>>>
>>>> I previously used GLDAS files and mask files to calculate monthly area
>>>> averages. However, I am trying to get MERRA2 data to work with the masks
>>>> now. I keep getting an error saying:
>>>>
>>>> ga-> run EXAMPLE_MASK_SCRIPT
>>>> Error in gagchk: axis sizes are not the same
>>>> Operation error:  Incompatable grids
>>>>   Dimension = 0
>>>>   1st grid range = -179 181   2nd = 1 360
>>>> Operation Error:  Error from maskout function
>>>> Operation Error:  Error from aave function
>>>>   Error ocurred at column 1
>>>> DISPLAY error:  Invalid expression
>>>>   Expression =
>>>> aave(maskout(prectot.2,mask.1(time=00:30z01feb1980)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)
>>>>
>>>>
>>>> I tried to fix x and y grid points using the following command, but it
>>>> did not work:
>>>>
>>>> aave(maskout(prectot.2(x=1,y=1),mask.1(x=1,y=1
>>>> ,time=00:30z01feb1980)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> The result of 'q dims' of the mask file is:
>>>>
>>>> Default file number is: 1
>>>> X is fixed     Lon = 0  X = 180.5
>>>> Y is fixed     Lat = -90  Y = 0.5
>>>> Z is fixed     Lev = 0  Z = 1
>>>> T is fixed     Time = 00:30Z01DEC1991  T = 143
>>>> E is fixed     Ens = 1  E = 1
>>>>
>>>> The result of 'q dims' of the MERRA2 file is:
>>>>
>>>> Default file number is: 2
>>>> X is fixed     Lon = 0  X = 1
>>>> Y is fixed     Lat = -90  Y = 1
>>>> Z is fixed     Lev = 0  Z = 1
>>>> T is fixed     Time = 00:30Z01DEC1991  T = 143
>>>> E is fixed     Ens = 1  E = 1
>>>>
>>>>
>>>> How can I get the x and y grid points to match up? I'm not even sure if
>>>> this is the issue, as the GLDAS and mask files didn't have matching x and y
>>>> grid points, yet they still worked.
>>>>
>>>> -Justin
>>>>
>>>> _______________________________________________
>>>> gradsusr mailing list
>>>> gradsusr at gradsusr.org
>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>
>>>>
>>>
>>>
>>> --
>>> Jeff Duda
>>> Post-doctoral research associate
>>> University of Oklahoma School of Meteorology
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> --
> Jennifer Miletta Adams
> Center for Ocean-Land-Atmosphere Studies (COLA)
> George Mason University
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160712/d3dac008/attachment.html 