<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div id="yiv368885532"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: times new roman,new york,times,serif; font-size: 12pt;"><div>Thanks Kieren for your reply.</div><div>As you suggest, I tried only the following lines <br></div><div><br></div><div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white; color: rgb(0, 64, 127);"><span>'sdfopen
/home/precis/wind/uwnd.nc'<br>
'sdfopen /home/precis/precip/prate_Interp.nc'<br>
<br>
'set lat -90 90'<br>
'set lon -180 180'<br>
'set lev 1000 100'<br>
'set t 1 2'</span></div> 

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white; color: rgb(0, 64, 127);"> &nbsp;</div> 

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white; color: rgb(0, 64, 127);">‘set
dfile 2’</div> 

<div style="color: rgb(0, 64, 127);" class="yiv368885532MsoNormal"><span>'pre = prate.2*86400'</span></div><div><br><span></span></div><div><span>But I got the following error:</span></div><div><span><br></span></div><div><span><span style="color: rgb(0, 64, 127);">Data Request Warning:&nbsp; Request is completely outside file limits&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><br style="color: rgb(0, 64, 127);"><span style="color: rgb(0, 64, 127);">Define Error:&nbsp; Unable to allocate data memory&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><br style="color:
 rgb(0, 64, 127);"><span style="color: rgb(0, 64, 127);">&nbsp; Size of request was -19031830 grid elements&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><br style="color: rgb(0, 64, 127);"><span style="color: rgb(0, 64, 127);">*** glibc detected *** grads: double free or corruption (!prev): 0x08b07620 ***&nbsp;&nbsp;&nbsp;</span>&nbsp;&nbsp;&nbsp; </span></div><div>&nbsp;<br>What's wrong!!<br><br>After that, the GUI interface goes off. Wind data is 4D (t,lev,lat,lon) and precip 3D (t,lev,lat,lon). Dimensions of coordinates t,lat,lon are the same for both variables.<br><br>Thanks<br>Wilfied</div><div><br></div>  <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div
 style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1">  <b><span style="font-weight: bold;">De&nbsp;:</span></b> "Marsh, Kieren" &lt;k.marsh@geos-uae.com&gt;<br> <b><span style="font-weight: bold;">À&nbsp;:</span></b> wilfried pokam &lt;wpokam@yahoo.fr&gt;; GrADS Users Forum &lt;gradsusr@gradsusr.org&gt; <br> <b><span style="font-weight: bold;">Envoyé le :</span></b> Mardi 13 mars 2012 11h03<br> <b><span style="font-weight: bold;">Objet&nbsp;:</span></b> RE: [gradsusr] Performs array assignments based on a conditional array<br> </font> </div> <br><div id="yiv368885532">




 

<style>
<!--
#yiv368885532   
 filtered  {font-family:"Cambria Math";panose-1:2 4 5 3 5 4 6 3 2 4;}
#yiv368885532 filtered  {font-family:Calibri;panose-1:2 15 5 2 2 2 4 3 2 4;}
#yiv368885532 filtered  {font-family:Tahoma;panose-1:2 11 6 4 3 5 4 4 2 4;}
#yiv368885532   
 p.yiv368885532MsoNormal, #yiv368885532  li.yiv368885532MsoNormal, #yiv368885532  div.yiv368885532MsoNormal
        {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;font-family:"serif";}
#yiv368885532  a:link, #yiv368885532  span.yiv368885532MsoHyperlink
        {color:blue;text-decoration:underline;}
#yiv368885532  a:visited, #yiv368885532  span.yiv368885532MsoHyperlinkFollowed
        {color:purple;text-decoration:underline;}
#yiv368885532  span.yiv368885532EmailStyle17
        {font-family:"sans-serif";color:#1F497D;}
#yiv368885532  .yiv368885532MsoChpDefault
        {font-size:10.0pt;}
#yiv368885532 filtered  {margin:72.0pt 72.0pt 72.0pt 72.0pt;}
#yiv368885532  div.yiv368885532Section1
        {}
-->
</style>

<div>

<div class="yiv368885532Section1">

<div class="yiv368885532MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">You need to specify to GrADS which file to usesing set dfile:</span></div> 

<div class="yiv368885532MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> &nbsp;</span></div> 

<div class="yiv368885532MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">See my ‘set dfile’ additions below (this assumes you sdfopen the
files in the order you specify)</span></div> 

<div class="yiv368885532MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> &nbsp;</span></div> 

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">'sdfopen
/home/precis/wind/uwnd.nc'<br>
'sdfopen /home/precis/precip/prate_Interp.nc'<br>
<br>
'set lat -90 90'<br>
'set lon -180 180'<br>
'set lev 1000 100'<br>
'set t 1 2'</span></div> 

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;"> &nbsp;</span></div> 

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">‘set
dfile 2’</span></div> 

<div class="yiv368885532MsoNormal"><span style="color: black;">'pre = prate.3*86400'</span></div> 

<div class="yiv368885532MsoNormal"><span style="color: black;">‘set dfile 1’<br>
'p&nbsp;&nbsp; = const(maskout(uwnd.1,pre-1),0,-u)</span><span style="font-size: 11pt; color: rgb(31, 73, 125);"></span></div> 

<div class="yiv368885532MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> &nbsp;</span></div> 

<div>

<div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223); padding: 3pt 0cm 0cm;">

<div class="yiv368885532MsoNormal"><b><span style="font-size: 10pt;" lang="EN-US">From:</span></b><span style="font-size: 10pt;" lang="EN-US"> gradsusr-bounces@gradsusr.org
[mailto:gradsusr-bounces@gradsusr.org] <b>On Behalf Of </b>wilfried pokam<br>
<b>Sent:</b> 13 March 2012 13:06<br>
<b>To:</b> gradsusers; GrADS Users Forum<br>
<b>Subject:</b> [gradsusr] Performs array assignments based on a conditional
array</span></div> 

</div>

</div>

<div class="yiv368885532MsoNormal"> &nbsp;</div> 

<div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">Hello
Grads users,</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;"> &nbsp;</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">I want to
assign to the wind zonal wind (at all vertical levels) the value 0&nbsp; at
grid points where the precipitation is less than 10mm/day. I tried the
following, but I have an error</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">-------------------------------------------------------------<br>
'sdfopen /home/precis/wind/uwnd.nc'<br>
'sdfopen /home/precis/precip/prate_Interp.nc'<br>
<br>
'set lat -90 90'<br>
'set lon -180 180'<br>
'set lev 1000 100'<br>
'set t 1 2'</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;"> &nbsp;</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">'pre =
prate.3*86400'<br>
'p&nbsp;&nbsp; = const(maskout(uwnd.1,pre-1),0,-u)<br>
----------------------------------------------------</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;"> &nbsp;</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">Both
'uwnd.nc' and 'prate_Interp.nc' are on the same horizontal grids.</span></div> 

</div>

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;"> &nbsp;</span></div> 

</div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;" lang="EN-US">Any
help would be appreciated<br>
<br>
Regards<br>
Wilfried</span><span style="color: black;"></span></div> 

<div>

<div class="yiv368885532MsoNormal" style="background: none repeat scroll 0% 0% white;"><span style="color: black;">&nbsp; <br></span></div></div><div><div><div class="yiv368885532MsoNormal" style="margin-bottom: 12pt; background: none repeat scroll 0% 0% white;"><span style="color: black;"><br>
_______________________________________________<br>
gradsusr mailing list<br>
<a rel="nofollow" ymailto="mailto:gradsusr@gradsusr.org" target="_blank" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a rel="nofollow" target="_blank" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
</span></div> 

</div>

</div>

</div>

</div>

</div>

</div><br><br> </div> </div>  </div></div></div><br><br> </div> </div>  <img
 src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAABYAAAAUCAYAAACJfM0wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAK8AAACvABQqw0mAAAAB90RVh0U29mdHdhcmUATWFjcm9tZWRpYSBGaXJld29ya3MgOLVo0ngAAAAWdEVYdENyZWF0aW9uIFRpbWUAMDQvMDQvMDhrK9wWAAACLklEQVQ4jbXUP0wTcRQH8O/9ekdjkT8CUqpee00bRyNNmSRSV0PcJJoQg2i6ODTExEUHg04OaNSppqtCjQ4ukDSKSuLUwcm4NNZcQYsIGtD+u/f7MZSWXltqo/Ul7/JL7u7z3r3fLye53e5xj8ejoYWRSCSSstfr1YLBYHcr4XA4rMmMMciy3EoXjDHIjDEoivL/4fefrP1P3nYEvqzLajOIo8fQz5/cfH3cnVttCM8udQaODBxQFx44Ye9h4HxvdGWtgMlbSXV2SQoMHf0RNcGSJJlmvPLdos7fdyIWL+D5myx+ZwwUDAIRh2EU1wYRFItA6FwvIjdcGJr4qFYakiSBlTavlABwsJth7mUWmSyBOAfnAkQE4gKccxBx/MoYmHmcxuH+NgAwGQ03j3NeRjjnoGqcC/zcIgghAMBkNISJuKlbEy4EaKdoKerC5nNMxQdlgVx+t0siKhYQovwV1rbdtyoNxlhxxoqilBMA0uuES6Pt6NqP2hHsoDarhJuXD2F5NV/uuJR1T4XLzvTJ25/VyHUnzgzba0YkKq6pdB4T00m47EyvPhU1M54asy3ee5o55bvwQQWAr/PHMBfbQGhGrykCANqARZ8asy3+ccYjg/K3kcF9UQAYvrJ29dmrDUxHlnOxu72P+rpYrq5eFU39K649TCF0tnPB0WdtCt2z48rQHIp+8XTHu9ET7alm0aY6fnFHjda98a/w3wZjDJLP5xv3+/1aK+F
4PJ7cBm32CUNiyI2GAAAAAElFTkSuQmCC" style="position: absolute; visibility: visible; color: transparent; margin: 0px; border: medium none; z-index: 2147483647; left: 517px; top: 627px;" id="sm-target-image"></div></body></html>