[gradsusr] Cannot Locate Endwhile Statement Error

Davide Sacchetti davide.sacchetti at arpal.gov.it
Tue Jul 5 11:34:26 EDT 2016


there at least an error on write:

replace:
'write('GLDAS_MASK_DATA.txt' ,  string ,append)'
with:
rc=write('GLDAS_MASK_DATA.txt' ,  string ,'append')

hoping it may help
Davide

On Tue, 2016-07-05 at 11:13 -0400, Justin Hicks 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





More information about the gradsusr mailing list