<div dir="ltr">GrADS users,<div><br></div><div>I'm having a few issues trying to find the maximum value of a variable for a shapefile. I found some documentation online of how to find the lat/lon location of the maximum value (with no shapefile), using: </div><div><br></div><div><pre><font color="#ff0000">*helps display lat</font></pre></div><div><pre><font color="#ff0000">'d maxloc(max('var',x=1,x=360),y=1,y=180)' </font></pre><pre><pre style="color:rgb(0,0,0)"><br></pre><pre><font color="#ff0000">*helps display lon</font></pre><font color="#ff0000">'d maxloc(max('var',y=1,y=180),x=1,x=360)' </font></pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)">This code works, but I'm having issues using the 'maskout' function with this code. </pre><pre style="color:rgb(0,0,0)">This is the display code I currently have that uses the maskout function with respect to the variables 'rainfsfc' and 'mask'...it is tested and works perfectly:</pre><pre style="color:rgb(0,0,0)"><br></pre><pre><font color="#ff0000">d maskout(rainfsfc.1*86400,lterp(mask.2(time=00Z01JAN1948),rainfsfc.1))</font><font color="#000000"><br></font></pre><pre><font color="#ff0000"><br></font></pre><pre><font color="#000000">I want to combine the two sets of display code to get the location of the maximum value within that shapefile, so I tried: </font></pre><pre style="color:rgb(0,0,0)"><br></pre><pre><pre style="color:rgb(34,34,34)"><font color="#ff0000">d maxloc(max(maskout(rainfsfc.1*86400,lterp(mask.2(time=00Z01JAN1948),rainfsfc.1)),x=1,x=360),y=1,y=180)</font><font color="#000000"><br></font></pre><pre><font color="#000000">and</font></pre><pre><pre><pre style="color:rgb(34,34,34)"><font color="#ff0000">d maxloc(max(maskout(rainfsfc.1*86400,lterp(mask.2(time=00Z01JAN1948),rainfsfc.1)),y=1,y=180),x=1,x=360)</font><font color="#000000"><br></font></pre><pre><font color="#000000"><br></font></pre><pre><font color="#000000">H</font>owever, these two did not work and I ended up with this error: </pre><pre><font color="#ff0000">Error from LTERP: The destination grid is a constant
Operation Error: Error from lterp function</font></pre><div style="color:rgb(0,0,0)"><font color="#ff0000"><br></font></div><div><font color="#000000">Is there an easier way to get this to output a lat/lon coordinate of the maximum value for this shapefile?</font></div><div><font color="#000000"><br></font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Thanks,</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">-Justin</font></div><div><font color="#000000"><br></font></div><div><font color="#000000"><br></font></div></pre></pre><div style="color:rgb(0,0,0)"><font color="#ff0000"><br></font></div></pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"><br></pre><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:"Times New Roman",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:"Times New Roman",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:"Times New Roman",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></div>