<div dir="ltr">Justin,<div><br></div><div>It appears there may be other problems with your sample script. For example,</div><div><br></div><div>'q time' should yield only one line, not five. Perhaps you intended to get the results of 'q dim' instead? That one does produce five lines, and your subwrd(rectime,6) would indeed give you a date.</div><div><br></div><div>Also, after your 'd aave(maskout...' in the while loop, you have "evap = sublin(result,2)" which I believe produces only a single line in the result, not two as you code suggests. For the subsequent line, would you not want "evap = subwrd(evap,4)"?</div><div><br></div><div>Stephen Mc</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 5, 2016 at 11:13 AM, Justin Hicks <span dir="ltr"><<a href="mailto:jhicks2014@gmail.com" target="_blank">jhicks2014@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">GrADS users,<div><br></div><div>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."</div><div><br></div><div>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):</div><div><br></div><div><div>'reinit'</div><div>'open /data2/control/Region15_AFR_mask_TEST.nc'</div><div>'set t 1 last'</div><div>'set dfile 1'</div><div>'set lon -179.5 179.5'</div><div>'set lat -89.5 89.5'</div><div>'open /data2/control/GLDAS_VIC10_M.ctl'</div><div>'set dfile 2'</div><div>'set lon -179.5 179.5'</div><div>'set lat -89.5 89.5'</div><div>'set time 00Z01JAN1979 00Z01JAN2016'</div><div>'set dfile 1'</div><div>'d aave(maskout(evap.2,mask.1(time=00Z01JAN1979)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)'</div><div>'set gxout print'</div><div>'set prnopts %6.2f 1 1'</div><div>nmonths = 428</div><div>i = 1</div><div>while (i <= nmonths)</div><div> 'set t 'i</div><div> 'q time'</div><div>rectime = sublin(result,5)</div><div>recdate = subwrd(rectime,6) </div><div> 'd aave(maskout(evap.2,mask.1(time=00Z01JAN1979)),lon=-179.5,lon=179.5,lat=-89.5,lat=89.5)'</div><div>evap = sublin(result,2)</div><div>evap = subwrd(evap,1)</div><div>string = rectime' 'evap</div><div> 'write('GLDAS_MASK_DATA.txt' , string ,append)'</div><div> i = i + 1</div><div>endwhile</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div>I'm not sure why GrADS is giving me this error.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>-Justin</div></font></span></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>