<div dir="ltr"><span style="font-size:12.8px">GrADS users,</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I have come across a strange problem while running a gradscript. I am using this script to output average monthly data of different regions to a text file. Each block of commands contains a different mask, but uses the same datafile. The first block of commands is shown below:</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> &#39;reinit&#39;<div>&#39;open /data2/control/Region00_WLD_<wbr>mask_GLDAS.nc&#39;</div><div>&#39;set t 1 756&#39;</div><div>&#39;set dfile 1&#39;</div><div>&#39;set lon -179.5 179.5&#39;</div><div>&#39;set lat -89.5 89.5&#39;</div><div>&#39;sdfopen /data2/control/GLDAS_NOAH10_M.<wbr>020.ctl_1Deg_RainfallRateSfc.<wbr>nc&#39;</div><div>&#39;set dfile 2&#39;</div><div>&#39;set lon -179.5 179.5&#39;</div><div>&#39;set lat -89.5 89.5&#39;</div><div>&#39;set time 00Z01JAN1948&#39;</div><div>&#39;set dfile 1&#39;</div><div> nmonths=756</div><div> i=1</div><div> while (i&lt;=nmonths)</div><div> &#39;set t &#39;i</div><div>&#39;q file 1&#39;</div><div>region = sublin(result,2)</div><div>region = subwrd(region,2)</div><div>&#39;q file 2&#39;</div><div>dataset = sublin(result,2)</div><div>dataset = subwrd(dataset,2)</div><div> &#39;q dims&#39;</div><div> rectime = sublin(result,5)</div><div> recdate = subwrd(rectime,6) </div><div>&#39;d aave(maskout(rainfsfc.2*86400,<wbr>lterp(mask.1(time=<wbr>00Z01JAN1948),rainfsfc.2)),<wbr>lon=-179.5,lon=179.5,lat=-89.<wbr>5,lat=89.5)&#39; </div><div>  say &#39;t=&#39;i  </div><div> say result</div><div> rainfsfc = sublin(result,1)</div><div> rainfsfc = subwrd(rainfsfc,4)</div><div> string = dataset&#39; &#39;region&#39; &#39;rectime&#39; &#39;rainfsfc</div><div> rc=write(&#39;GLDAS_MASK_DATA.<wbr>txt&#39; ,  string ,&#39;append&#39;)</div><div> i = i + 1</div><div> endwhile</div><div><br></div><div>I repeat these commands in the same script over and over, only changing the second line (ex. the next one would read &#39;open /data2/control/Region01_ART_<wbr>mask_GLDAS.nc&#39;). The data is output perfectly when I have six masks in the script, but when I include the seventh, the data output from the third mask / block of commands is altered, and I get the following error:</div><div><br></div><div><div>QUERY FILE Error:  file 2 not open</div><div>Data Request Error:  File number out of range </div><div>  Variable = rainfsfc.2 </div><div>Operation Error:  Error from maskout function</div><div>Operation Error:  Error from aave function</div><div>  Error ocurred at column 1</div><div>DISPLAY error:  Invalid expression </div><div>  Expression = aave(maskout(rainfsfc.2*86400,<wbr>lterp(mask.1(time=<wbr>00z01jan1948),rainfsfc.2)),<wbr>lon=-179.5,lon=179.5,lat=-89.<wbr>5,lat=89.5)</div><div>t=703</div><div>Data Request Error:  File number out of range </div><div>  Variable = rainfsfc.2 </div><div>Operation Error:  Error from maskout function</div><div>Operation Error:  Error from aave function</div><div>  Error ocurred at column 1</div><div>DISPLAY error:  Invalid expression </div><div>  Expression = aave(maskout(rainfsfc.2*86400,<wbr>lterp(mask.1(time=<wbr>00z01jan1948),rainfsfc.2)),<wbr>lon=-179.5,lon=179.5,lat=-89.<wbr>5,lat=89.5)</div><div><br></div><div><br></div><div>I don&#39;t believe the error lies within the third mask file itself because it plots fine when it is run independently from the other masks. I also get the same output issue when I run, for example, masks 50-56 (one or two of these will come up with the error, but running masks 50-55 will yield no problems). </div><div><br></div><div>Any help would be appreciated to me solve this issue, as I&#39;m baffled.</div></div><div><br></div><div>-Justin</div></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><span style="font-size:12.8px"><span style="font-size:10pt;font-family:&quot;times new roman&quot;,serif;background-image:initial;background-repeat:initial"><font color="#000000"><b>Justin Hicks</b></font></span></span><br></div><div><i style="font-size:12.8px"><span style="font-size:10pt;font-family:&quot;times new roman&quot;,serif;background-image:initial;background-repeat:initial"><font color="#000000">Undergraduate<br>Atmospheric Sciences Major</font></span></i></div><div><i style="font-size:12.8px"><span style="font-size:10pt;font-family:&quot;times new roman&quot;,serif;background-image:initial;background-repeat:initial"><font color="#000000">GIS Minor<br>University of Maryland-College Park</font></span></i><font color="#000000" face="arial, helvetica, sans-serif" size="2"><b><br></b></font></div></div></div></div></div>
</div>