<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:14pt">Hi,<br><br>Dear Jennifer and Mahakur,<br><br>Thank you for your reply. I am sorry that instead of abs(), I wrote mod. I have tried your suggestions. However, the error <br><br>SET error: Invalid operand<br> Operand = <br>SET error: Invalid operand<br> Operand = <br>:<br>:<br>DEFINE error: Invalid expression.<br><div><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>persists. Let me send you my script.</span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new
york,times,serif; background-color: transparent; font-style: normal;"><span>'reinit'<br>'open SI-Decadal.ctl'<br><br>'set gxout fwrite'<br>'set fwrite data-sorted.dat'<br><br>t=1<br> while(t<=120)<br> 'set t 't<br><br> lv=1<br> while(lv<=15)<br> 'set lv 'lv<br><br> lt=-90<br> while(lt<=90)<br> 'set lt 'lt<br><br> ln=0<br> while(ln<=360)<br> 'set ln 'ln<br><br><br>'define si_avg = aave(si,lon=0,lon=360,lat=-90,lat=90)'<br>'define dv = pow(si_avg-si,2)'<br>'define sd = sqrt(aave(dv,lon=0,lon=360,lat=-90,lat=89))'<br>'define x = 2*sd'<br>'define y =
abs(si_avg-si)'<br><br>* if (y>=x);<br>* 'set gxout fwrite'<br>* 'set fwrite data-sorted.dat'<br> 'display maskout(si,si-x)'<br>* endif<br><br> ln=ln+1<br> endwhile<br> lt=lt+1<br> endwhile<br> lv=lv+1<br> endwhile<br> t=t+1<br>endwhile<br><br>'disable fwrite'<br></span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div
style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>Thank you again for looking into my problem. Now, I am planning to write a fortran program to read netcdf data and solve my problem.</span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>Regards</span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>Pankaj</span></div><div style="color: rgb(0, 0, 0); font-size: 18.6667px; font-family: times new roman,new york,times,serif;
background-color: transparent; font-style: normal;"><span>EES, Hokkaido University<br></span></div><div style="display: block;" class="yahoo_quoted"> <br> <br> <div style="font-family: times new roman, new york, times, serif; font-size: 14pt;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> On Monday, December 16, 2013 10:51 PM, Jennifer Adams <jma@cola.iges.org> wrote:<br> </font> </div> <div class="y_msg_container"><div id="yiv2968991608"><div><br clear="none"><div><div>On Dec 16, 2013, at 3:32 AM, pankaj kumar wrote:</div><br class="yiv2968991608Apple-interchange-newline" clear="none"><blockquote type="cite"><div><div style="color:rgb(0, 0, 0);background-color:rgb(255, 255, 255);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:14pt;"><div><span style="font-family:times new
roman, new york, times, serif;">Hi,</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"><span style="font-family:times new roman, new york, times, serif;"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">I am writing a grads script to sort out values which differnece from mean is more than 3 times standard deviation. The netcdf data is spread in x,y,z,t dimesions. Though while loops for x,y,z,t are running well, I get error as below:</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"><span style="font-family:times new
roman, new york, times, serif;"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"> Operand = lat1 -90<br clear="none">SET error: Invalid operand<br clear="none"> Operand = lon1 0<br clear="none">SET error: Invalid operand<br clear="none"> Operand = lon1 1<br clear="none">SET error: Invalid operand<br clear="none"> Operand = lon1 2<br clear="none">:</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">:</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span
style="font-family:times new roman, new york, times, serif;">:</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"> Operand = lon1 360<br clear="none"></span><span style="font-family:times new roman, new york, times, serif;"><span style="font-family:times new roman, new york, times, serif;">SET error: Invalid operand</span></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;">:</div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;">:<br clear="none"><span style="font-family:times new roman, new york, times, serif;"><span style="font-family:times new roman, new
york, times, serif;"></span></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"><span style="font-family:times new roman, new york, times, serif;"></span>DEFINE error: Invalid expression. <br clear="none"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;">My script for defining variables are as follows:</div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"><span
style="font-family:times new roman, new york, times, serif;"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">'define si_avg = aave(si,lon=0,lon=360,lat=0,lat=180)'<br clear="none">'define dv = pow(siavg-si,2)'<br clear="none">'define sd = sqrt(aave(dv,lon=0,lon=360,lat=0,lat=179))'<br clear="none">'define x = 3*sd'</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">'define y = mod(si_avg-si)'<br clear="none"></span></div></div></div></blockquote><div><br clear="none"></div>What is the mod() function?</div><div><br clear="none"><blockquote type="cite"><div><div style="color:rgb(0, 0,
0);background-color:rgb(255, 255, 255);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:14pt;"><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"> if (y>=x);<br clear="none"></span></div></div></div></blockquote><div><br clear="none"></div>The 'if' statment in the scripting language doesn't evaluate defined variables like this. You need to use the maskout() function to create a grid with missing values where x < y. Something like this:</div><div>' d maskout(si,y-x)'</div><div>--Jennifer</div><div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd69687"><div><br clear="none"></div><div><br clear="none"></div></div><div><div class="yiv2968991608yqt9437638797"
id="yiv2968991608yqtfd17905"><br clear="none"></div><blockquote type="cite"><div><div style="color:rgb(0, 0, 0);background-color:rgb(255, 255, 255);font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;font-size:14pt;"><div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd46113"><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"> 'set gxout
fwrite'<br clear="none"> 'set fwrite data-sorted.dat'<br clear="none"> 'd 'si<br clear="none"> endif<br clear="none"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"><span style="font-family:times new roman, new york, times, serif;"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">Please see what kind of error I am facing here in defining variables though it seems okay as variables are defined properly. If someone has script or fortran program to sort data in x,y,z,t
which has much higher standard deviation, please help me
urgently.</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><br clear="none"><span style="font-family:times new roman, new york, times, serif;"></span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">Thanks in advance!</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;">Regards</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york,
times, serif;">Pankaj</span></div></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:times new roman, new york, times, serif;background-color:transparent;font-style:normal;"><span style="font-family:times new roman, new york, times, serif;"><div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd87225">EES, Hokkaido University</div><br clear="none"></span></div></div></div>_______________________________________________<br clear="none">gradsusr mailing list<br clear="none"><a rel="nofollow" shape="rect" ymailto="mailto:gradsusr@gradsusr.org" target="_blank" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br clear="none">http://gradsusr.org/mailman/listinfo/gradsusr<br clear="none"></blockquote></div><br clear="none"><div>
<span class="yiv2968991608Apple-style-span" style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;orphans:2;text-indent:0px;text-transform:none;white-space:normal;widows:2;word-spacing:0px;font-size:medium;"><span class="yiv2968991608Apple-style-span" style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"></span></span><div style="word-wrap:break-word;"><span class="yiv2968991608Apple-style-span" style="border-collapse:separate;color:rgb(0, 0,
0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"><span class="yiv2968991608Apple-style-span" style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"></span></span><div>--</div><div>Jennifer M. Adams<br clear="none">Center for Ocean-Land-Atmosphere Studies (COLA)<br clear="none">111 Research Hall, Mail Stop 2B3<br clear="none">George Mason University<br clear="none">4400 University Drive<br clear="none">Fairfax, VA 22030 <div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd00546"><br clear="none"><br clear="none"></div></div><div
class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd81698"><div><br class="yiv2968991608khtml-block-placeholder" clear="none"></div><br class="yiv2968991608Apple-interchange-newline" clear="none"></div></div><div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd60377"><br class="yiv2968991608Apple-interchange-newline" clear="none">
</div></div><div class="yiv2968991608yqt9437638797" id="yiv2968991608yqtfd84012">
<br clear="none"></div></div></div><br><br></div> </div> </div> </div> </div></body></html>