[gradsusr] Cannot Locate Endwhile Statement Error

Justin Hicks jhicks2014 at gmail.com
Thu Jul 7 12:33:35 EDT 2016 

GrADS users,

The GLDAS data is working with the masks. However, I am trying to get
MERRA2 data to work with the masks now. I keep getting an error saying:

Data Request Error: Invalid dimension expression
  Attempt to set or modify varying dimension
  Variable = prectot.2, Dimension = 0
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(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)

As you can see, I tried to fix x and y grid points because I kept getting
that error, but this did nothing.

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?

-Justin

On Tue, Jul 5, 2016 at 3:34 PM, Jennifer M Adams <jadams21 at gmu.edu> wrote:

>
> On Jul 5, 2016, at 3:21 PM, Justin Hicks <jhicks2014 at gmail.com> wrote:
>
> Jennifer,
>
> I am still getting values of 0 in my text file (and on the terminal).
> Typing in a 'q dims' command gives LEV=0  Z=1    for the mask and LEV=0
>  Z=0  for the GLDAS file. Could the inconsistency be the issue?
>
> Maybe. Try this:
>
> evap.2(z=1),mask.1(z=1,time=00Z01JAN1979)
>>>>
>>>
>
> There aren't any errors displayed on the terminal. And just to note,
> taking out the while statement gives me monthly values on the terminal when
> the maskout function is in use, so I assume the issue is within the while
> statement.
>
> Put some more ‘say’ statements in there to make sure you’re parsing and
> pasting all the strings correctly.
> —Jennifer
>
>
> -Justin
>
> On Tue, Jul 5, 2016 at 2:47 PM, Jennifer M Adams <jadams21 at gmu.edu> wrote:
>
>>
>>
>> On Jul 5, 2016, at 1:51 PM, Justin Hicks <jhicks2014 at gmail.com> wrote:
>>
>> Jennifer,
>>
>> How would I go about making the lat/lon values fixed? I know the aave
>> function needs 4 different lat/lon arguments. Would this just be 'd
>> ave(.......,lon=-179.5,lon=-179.5,lat=-89.5,lat=-89.5)’
>>
>> Correct. These args provide the range for your average.
>>
>>    or     'd ave(.......,x=1,x=1,y=1,y=1)’?
>>
>> No. This will give you an average over only one grid point.
>>
>> The ‘set x 1’ and ’set y 1’ commands are there to restrict the dimension
>> environment to match the result that comes out of aave(), which does not
>> vary in X or Y. Those dimensions need to be fixed, and it doesn’t matter at
>> all what they are, it could just as easily be ‘set x 22’ and ‘set y 89’.
>>
>> —Jennifer
>>
>>
>>
>> -Justin
>>
>> On Tue, Jul 5, 2016 at 12:25 PM, Jennifer M Adams <jadams21 at gmu.edu>
>> wrote:
>>
>>> Hi, Justin —
>>> Your dimension environment isn’t set properly. When using the aave()
>>> function, make sure the lat and lon dimensions are fixed — you’ve got them
>>> varying. Also, the LEV dimension value at z=1 may be different for your two
>>> data sets. Put in a 'say result’ statement after each display and it should
>>> print out more than just the warning message; it should show you what the
>>> dimensions are for the variables in your expression, and that should guide
>>> you as to where the mismatch is. I also think you can remove most of the
>>> code before the 'while’ statement:
>>>
>>> ‘open /data2/control/Region15_AFR_mask_TEST.nc'
>>> 'open /data2/control/GLDAS_VIC10_M.ctl'
>>> ‘set dfile 2’
>>> ‘set x 1’
>>> ‘set y 1’
>>> ‘set gxout print'
>>> ‘set prnopts %6.2f 1 1’
>>> nmonths=428
>>> i=1
>>> while (i<nmonths)
>>>   ‘set t ‘i
>>>   ...
>>>   ‘d aave ….'
>>>   say ‘t=’t’
>>>   say result
>>>   ….
>>>   i=i+1
>>> endwhile
>>>
>>> —Jennifer
>>>
>>> On Jul 5, 2016, at 12:00 PM, Justin Hicks <jhicks2014 at gmail.com> wrote:
>>>
>>> To those who responded,
>>>
>>> Thank you for the responses! Those seemed to do the trick. However,
>>> instead of the 'endwhile' error, I am now receiving this error:
>>>
>>> Data Request Warning:  Request is completely outside file limits
>>>
>>> This error is repeated 9 times. The text file that is created also
>>> contains only one column with each value being 0.0.
>>>
>>> -Justin
>>>
>>> On Tue, Jul 5, 2016 at 11:38 AM, sushant puranik <
>>> sushantpuranik at gmail.com> wrote:
>>>
>>>> Along with the suggested changes simply press enter button once at
>>>> ending line of your script.
>>>> On 5 Jul 2016 8:45 pm, "Justin Hicks" <jhicks2014 at gmail.com> wrote:
>>>>
>>>>> GrADS users,
>>>>>
>>>>> I am running a script that should output values into a 2 column table.
>>>>> The first column will be the time (monthly in the form of 00Z01JAN1979) and
>>>>> the second column will be the value of the variable (evap) at that time.
>>>>> After running the code below, I receive an error saying "Unable to locate
>>>>> ENDWHILE statement for the WHILE statement at line 18."
>>>>>
>>>>> This is the script I'm running; everything before the while statement
>>>>> is setting up a mask (in this case, the mask only outputs values for
>>>>> Africa):
>>>>>
>>>>> 'reinit'
>>>>> 'open /data2/control/Region15_AFR_mask_TEST.nc'
>>>>> 'set t 1 last'
>>>>> 'set dfile 1'
>>>>> 'set lon -179.5 179.5'
>>>>> 'set lat -89.5 89.5'
>>>>> 'open /data2/control/GLDAS_VIC10_M.ctl'
>>>>> 'set dfile 2'
>>>>> 'set lon -179.5 179.5'
>>>>> 'set lat -89.5 89.5'
>>>>> 'set time 00Z01JAN1979 00Z01JAN2016'
>>>>> 'set dfile 1'
>>>>> 'd
>>>>> aave(maskout(evap.2,mask.1(time=00Z01JAN1979)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)'
>>>>> 'set gxout print'
>>>>> 'set prnopts %6.2f 1 1'
>>>>> nmonths = 428
>>>>> i = 1
>>>>> while (i <= nmonths)
>>>>>  'set t 'i
>>>>>  'q time'
>>>>> rectime = sublin(result,5)
>>>>> recdate = subwrd(rectime,6)
>>>>>  'd
>>>>> aave(maskout(evap.2,mask.1(time=00Z01JAN1979)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)'
>>>>> evap = sublin(result,2)
>>>>> evap = subwrd(evap,1)
>>>>> string = rectime' 'evap
>>>>>  'write('GLDAS_MASK_DATA.txt' ,  string ,append)'
>>>>>  i = i + 1
>>>>> endwhile
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> I'm not sure why GrADS is giving me this error.
>>>>>
>>>>> -Justin
>>>>>
>>>>> _______________________________________________
>>>>> 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
>>>>
>>>>
>>> _______________________________________________
>>> 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
>>>
>>>
>> _______________________________________________
>> 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
>>
>>
> _______________________________________________
> 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/20160707/4f8ad7af/attachment.html

More information about the gradsusr mailing list