[gradsusr] Cannot Locate Endwhile Statement Error

[Jennifer M Adams]

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<mailto: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<mailto: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<mailto: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<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr


_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr


_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer Miletta Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
George Mason University



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160705/86899520/attachment-0001.html