<div dir="ltr"><div dir="ltr"><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif" style="">Hi GrADS users, </font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">I was wondering if anyone knew if there was a way to plot a large enough symbol on a global map in GrADS where a difference field exceeds a specific threshold?</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">For example, in my GrADS script I am plotting sea surface temperatures and masking out land areas:</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default" style=""><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'set_parea 2 2 2 1 -margin 'margin</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">*'color -levs -15 -10 -5 5 10 15 -kind blue->white->red'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif" color="#0000ff"><b>'color -20 20 5 -kind blue->white-(0)->white->red'</b></font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif" color="#0000ff"><b>'d (TMPsfc.2-TMPsfc.1)'</b></font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'q gxinfo'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">line3=sublin(result,3)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">x1=subwrd(line3,4)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">x2=subwrd(line3,6)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">say '    (x1,x2) = ('x1','x2')'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">line4=sublin(result,4)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">y1=subwrd(line4,4)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">y2=subwrd(line4,6)</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">say '    (y1,y2) = ('y1','y2')'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'xcbar 'x1' 'x2' 'y1-0.3' 'y1-0.2' -fw .06 -fh .06 -fs 1'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">*'cbarn'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'draw title FV3-GFS'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">* 'set rgb 73 80 80 80'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'set rgb 73 255 255 255'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">** Plots ocean temperatures and masks land **</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">'basemap L 73 1 M'</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">** Plots land temperatures and masks ocean **</font></div><div class="gmail_default" style=""><font face="arial, helvetica, sans-serif">*'basemap O 0 1 H'</font></div><div><br></div><div>I want to say in GrADS, look at all points on the global map (not including land) and plot a large enough star that the user can see anywhere in the globe where the difference in surface temperature exceeds 10 degrees Celsius, for example.  Does anyone know how I can do this?</div><div><br></div><div>Thank you,</div><div>Tracey Dorian</div></div></div></div>