<div dir="ltr">Hi JP,<div><br></div><div>it seems that your script omits increment and endwhile line</div><div><br></div><div>t=t+1</div><div>endwhile </div><div><br></div><div>on the end of your process.</div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Raymond</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 1, 2016 at 6:27 AM, JP Kalb <span dir="ltr"><<a href="mailto:jpkalb1990@yahoo.com" target="_blank">jpkalb1990@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"><div><span>Hi there,</span></div><div dir="ltr"><span>I am new to learning GrADS. I am having trouble with my if loops for my code.</span></div><div dir="ltr">"open NAM12Ztest.ctl"</div><div dir="ltr">"set display color white"</div><div dir="ltr">"clear"</div><div dir="ltr">"set mpdset hires"</div><div dir="ltr">"set grads off"</div><div dir="ltr">"set lev 925"</div><div dir="ltr">t=1</div><div dir="ltr">while (t <= 2)</div><div dir="ltr">'set t 't</div><div dir="ltr">'query dim'</div><div dir="ltr">diminfo=result</div><div dir="ltr">'run <a href="http://rgbset.gs">rgbset.gs</a>'</div><div dir="ltr"> 'set lat 36 39'</div><div dir="ltr"> 'set lon -125 -120'</div><div dir="ltr">'set clevs 0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1.0'</div><div dir="ltr">'set ccols 0 4 13 3 10 7 12 8 2'</div><div dir="ltr">'set gxout shaded'</div><div dir="ltr">'define lcl=((TMP2m-DPT2m)*400)'</div><div dir="ltr">'define tc=(TMPprs-273.15)'</div><div dir="ltr">'define es=(6.11*exp((17.625*tc)/(tc+243.05)))'</div><div dir="ltr">'define e=(es*RHprs/100)'</div><div dir="ltr">'define wr=(622*e/(925-e))'</div><div dir="ltr">'define ht=(HGTprs*3.28)'</div><div dir="ltr">'define expo=(exp((2257*wr)/(1004*TMPprs)))'</div><div dir="ltr">'define numer=(-2.06327*TMPprs/(6516))'</div><div dir="ltr">'define denom=(pow((RHprs/100),0.46))'</div><div dir="ltr">'define dthetae=(expo*numer/denom)'</div><div dir="ltr">'define b=((-6.67*dthetae-0.667)*(-0.1*VVELprs))'</div><div dir="ltr">'define c=(0)'</div><div dir="ltr">'d b'</div><div dir="ltr"> if (ht <= lcl)</div><div dir="ltr"> ' d c'</div><div dir="ltr"> endif</div><div dir="ltr"> if (VVELprs >= 0)</div><div dir="ltr"> ' d c'</div><div dir="ltr"> endif</div><div dir="ltr"></div><div dir="ltr">'set grads off'</div><div dir="ltr"><br></div><div dir="ltr">Any help would be appreciated.</div><div dir="ltr">Thanks,</div><span class="HOEnZb"><font color="#888888"><div></div><div> </div><div>JP Kalb</div></font></span></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>-------------------------------------------------------------------------------------------------------------------------------</div>Raymond C. Ordinario<br>MS Atmospheric Science<br>Ateneo De Manila University<br>Quezon City, Philippines<br><br>Weather Specialist I<br><div>Technique Application and Meteorological Satellite Section<br>Weather DIvision<br>Philippine Atmospheric, Geophysical and Astronomical Services Admistration<br>Quezon City, Philippines</div></div></div>
</div>