<div dir="ltr"><div><div><div><div><div><div><div><div>Hi - I'm using the local_max function to find local maxima in a variable field. The variable covers all longitudes, but is calculated for a 30 deg latitude band only. All values outside this latitude band are set to the @_FillValue for the variable.<br><br></div>Where the local maximum for the variable is located well within the latitude band, there is no problem However, when the local maximum is located on the edge of the latitude band (i.e. at either the minimum or maximum latitude of the band), the function fails. This seems to be because the immediately neighbouring points are all "missing" values. However, the function documentation states that missing values are ignored. <br><br></div>This can be seen from a simple example:<br>*************************************************<br></div><div>begin<br><br></div>myvar = new((/3,3/),float)<br></div>myvar@_FillValue = -999.<br></div></div></div></div>myvar = (/ (/-999.0,-999.0,-999.0/),(/60,70,55/),(/43,45,23/)/)<br><div><div><div><div><div><div><div><div><div>print(myvar)<br></div><div>mymax = local_max(myvar, False, 0.0)<br></div><div>print(mymax)<br><br></div><div>mymaxval = mymax@maxval<br></div><div>print(mymaxval)<br><br></div><div>end<br></div><div><br></div><div>However, if the missing values are set to e.g. zero, the function works:<br>*********************************************************************<br>begin<br><br>myvar = new((/3,3/),float)<br>myvar@_FillValue = -999.<br>myvar = (/ (/0.0,0.0,0.0/),(/60,70,55/),(/43,45,23/)/)<br><div>print(myvar)<br></div><div>mymax = local_max(myvar, False, 0.0)<br></div>print(mymax)<br></div><div><div><div><div><br>mymaxval = mymax@maxvalmymaxval = mymax@maxval<br></div><div>print(mymaxval)<br><br></div><div>end<br><br></div><div>While it might seem to be an easy workaround to set the next highest and lowest latitude band values to some miniscule value, this shouldn't really be necessary for the function to work. Can anyone help?</div><div><br clear="all"><div><div><div dir="ltr"><div><div dir="ltr"><div>-- <br>Tess Parker<br>PhD Candidate<br>School of Earth, Atmosphere and Environment<br>Room 225, Building 28<br></div><div>9 Rainforest Walk<br></div><div>Monash University, Clayton VIC 3800<br>Phone: <a href="tel:%2B61%203%209905%204458" value="+61399054458" target="_blank">+61 3 9905 4458</a> <br></div><a href="http://users.monash.edu.au/~tess/index.html" target="_blank">http://users.monash.edu.au/~tess/index.html</a><br></div></div></div></div></div>
</div></div></div></div></div></div></div></div></div></div></div></div></div>