[gradsusr] Cannot Locate Endwhile Statement Error

Justin Hicks jhicks2014 at gmail.com
Tue Jul 5 11:13:53 EDT 2016


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160705/e7ffc7e6/attachment.html 


More information about the gradsusr mailing list