Hello Jeff,<div><br></div><div>Thank you so much for the help. </div><div><br></div><div>I added your suggestions to <a href="http://skew.gs" target="_blank">skew.gs</a>. The SfcPlev is showing 1000, which I dont think is wrong. But the result of Sfcdewp is something I am not understanding. This is the following result I am getting on running <a href="http://skew.gs" target="_blank">skew.gs</a> after additions made of "say result" and "say SfcPlev". </div>
<div><br></div><div><div>ga-> <a href="http://skew.gs" target="_blank">skew.gs</a> </div><div>No hardcopy metafile open</div><div>All files closed; all defined objects released;</div><div>All GrADS attributes have been reinitialized</div>
<div>Initial Time is 00Z08JAN2007</div><div>Create gif images as well (1=yes ; 0=no)</div><div>0</div><div>Please enter coordinated for sounding</div><div> </div><div>Latitude values between: 4.96654 and 50.0419</div><div>
Enter Latitude</div><div>33</div><div> </div><div>Longitude values between: 48.3967 and 105.603</div><div>Enter Longitude</div><div>78</div><div> </div><div>Overlay two time periods 1=yes ; 0=no</div><div>0</div><div> </div>
<div>Please enter time to plot</div><div>Available times: 1 to 17</div><div>5</div><div>Drawing temperature sounding.</div><div>Drawing dewpoint sounding.</div><div>Data Type = grid</div><div>Dimensions = -1 -1</div><div>
I Dimension = -999 to -999</div><div>J Dimension = -999 to -999</div><div>Sizes = 1 1 1</div><div>Undef value = -9.99e+08</div><div>Undef count = 1 Valid count = 0</div><div>Min, Max = 1e+30 1e+30</div><div><br></div><div>
1000</div><div>getrh</div><div>13.743</div><div>1e+30</div><div>1000</div><div>0.00992721580569856</div><div>286.893</div><div>43.385364499893</div><div>1573660000</div><div>Non-numeric args to numeric operation</div><div>
Error occurred on line 2437</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div><div> Error occurred on line 1108</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div>
<div> Error occurred on line 69</div>
<div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div></div><div><br></div><div>I have also attached the new <a href="http://skew.gs" target="_blank">skew.gs</a> file along with my ctl file. I rechecked my ctl file, the file does have variable td. How am I to proceed further.</div>
<div><br></div><div>Thank you so much for taking the time to helping me out.</div><div><br></div><div>With regards,</div><div>Amulya Chevuturi</div><div><br></div><div><br><br><div class="gmail_quote">On Wed, Dec 19, 2012 at 10:40 PM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">As you can see from your output, the script has the value 1e+30 stored for dewpoint. That is clearly not correct. Therefore, there is a problem with the variable dewp getting passed into this function, which is named Thetae. The other error output tells you where the call to that function came from (line 1106). That line is<br>
<br> SfcThee=Thetae(Sfctemp,Sfcdewp,SfcPlev)<br><br>So there must be a problem with the variable Sfcdewp outside of the Thetae function. The variable Sfcdewp is computed in the line above 1106. Clearly something is not right there. The 10 or so lines above this are where the determination is made of how to obtain the surface dewpoint. See the line "d "snddewp? I would insert a 'say result' command directly after that command and see what the output shows. Also insert the command say SfcPlev in that same area. The error could be from either source.<br>
<br>Lastly, make sure that you have the variable td in your control file. I'm pretty sure I saw it there when I looked at it a few days ago, but I no longer have that file, so I can't check. If you don't have that, that is the source of the problem.<span><font color="#888888"><br>
<br>Jeff</font></span><div><div><br><br><div class="gmail_quote">On Wed, Dec 19, 2012 at 4:47 AM, Amulya Chevuturi <span dir="ltr"><<a href="mailto:amulya.chevuturi@gmail.com" target="_blank">amulya.chevuturi@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thank you so much for the help that you are providing me MR. Jeff.<div><br></div><div>So as you suggested, I added the say commands for the line 2424. This is for the script variable mix. What is the term "getrh" in that? Is it a function or a variable. As when I input say getrh the output is getrh. So I added these lines in my gs file before the script for mix variable:</div>
<div><br></div><div><div>say getrh</div><div>say temp</div><div>say dewp</div><div>say pres</div><div>say ws</div></div><div><br></div><div><br></div><div>This is the output I am getting now. </div><div><br></div><div><div>
ga-> <a href="http://skew.gs" target="_blank">skew.gs</a> </div><div>No hardcopy metafile open</div><div>All files closed; all defined objects released;</div><div>All GrADS attributes have been reinitialized</div><div>
Initial Time is 00Z08JAN2007</div>
<div>Create gif images as well (1=yes ; 0=no)</div><div>0</div><div>Please enter coordinated for sounding</div><div> </div><div>Latitude values between: 4.96654 and 50.0419</div><div>Enter Latitude</div><div>33</div><div>
</div><div>Longitude values between: 48.3967 and 105.603</div><div>Enter Longitude</div><div>78</div><div> </div><div>Overlay two time periods 1=yes ; 0=no</div><div>0</div><div> </div><div>Please enter time to plot</div>
<div>Available times: 1 to 17</div><div>5</div><div>Drawing temperature sounding.</div><div>Drawing dewpoint sounding.</div><div>getrh</div><div>13.743</div><div>1e+30</div><div>1000</div><div>0.00992721580569856</div><div>
286.893</div><div>43.385364499893</div><div>1573660000</div><div>Non-numeric args to numeric operation</div><div> Error occurred on line 2435</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div>
<div> Error occurred on line 1106</div>
<div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div><div> Error occurred on line 69</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div><div><br></div><div>Any help would be most appreciated.</div>
<div><br>
</div><div>With regards,</div><div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Amulya Chevuturi</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
Junior Research Fellow (Ph.D. Student)</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Jawaharlal Nehru University</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
New Delhi, India</div></div><div><br></div><br><div class="gmail_quote">On Tue, Dec 18, 2012 at 11:45 PM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Amulya,<br>Thank you for adding the new information. I can tell that the value of the script variable mix, which is given as 1573660000, cannot be right. Check line 2424 of your script code. That is where the variable mix is computed. Insert say statements for all the component variables there. With this kind of investigation you will eventually be able to discover the source of the error. My gut tells me one of the input variables you supplied to the script in the main function call plotskew() is not what the script is expecting it to be. I'm guessing it is the dewpoint variable. But I saw your control file and it appears there should not be a problem with that. Nonetheless, insert say statements to find out what the values are, then go from there. If you need help knowing what to do next or where to look, please ask.<span><font color="#888888"><br>
<br>Jeff</font></span><div><div><br><br><div class="gmail_quote">On Tue, Dec 18, 2012 at 1:21 AM, Amulya Chevuturi <span dir="ltr"><<a href="mailto:amulya.chevuturi@gmail.com" target="_blank">amulya.chevuturi@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thank you very much for replying Mr. Jeff. <div><br></div><div>I have added the lines as stated. In the output all of them are showing numeric values only. So I cannot make out what the problem is. Could you or anyone in the grads user list help me out further. </div>
<div><br></div><div>This the output I am getting at the terminal when I am running the <a href="http://skew.gs" target="_blank">skew.gs</a> with the changes you suggested. </div><div><br></div><div><div><div>ga-> <a href="http://skew.gs" target="_blank">skew.gs</a> </div>
<div>No hardcopy metafile open</div><div>All files closed; all defined objects released;</div><div>All GrADS attributes have been reinitialized</div><div>Initial Time is 00Z08JAN2007</div><div>Create gif images as well (1=yes ; 0=no)</div>
<div>0</div><div>Please enter coordinated for sounding</div><div> </div><div>Latitude values between: 4.96654 and 50.0419</div><div>Enter Latitude</div><div>33</div><div> </div><div>Longitude values between: 48.3967 and 105.603</div>
<div>Enter Longitude</div><div>78</div><div> </div><div>Overlay two time periods 1=yes ; 0=no</div><div>0</div><div> </div><div>Please enter time to plot</div><div>Available times: 1 to 17</div><div>5</div><div>Drawing temperature sounding.</div>
<div>Drawing dewpoint sounding.</div></div><div>286.893</div><div>43.385364499893</div><div>1573660000</div><div><div>Non-numeric args to numeric operation</div></div><div> Error occurred on line 2430</div><div>
<div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div>
<div> Error occurred on line 1106</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div><div> Error occurred on line 69</div><div> In file <a href="http://skew.gs" target="_blank">skew.gs</a></div>
</div></div><div><br></div><div>I am attaching the <a href="http://skew.gs" target="_blank">skew.gs</a> file with the changes added. Please let me know how I can further correct the error. Thank you for the help in advance.</div>
<div>
<div><br></div><div>With regards,</div><div><div style="color:rgb(34,34,34);direction:ltr;font-size:13px;font-family:arial,sans-serif;margin:5px 15px 0px 0px;padding-bottom:5px">
<div><div>Amulya Chevuturi</div><div>Junior Research Fellow (Ph.D. Student)</div><div>Jawaharlal Nehru University</div><div>New Delhi, India</div><div></div></div></div><div style="clear:both;color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif;margin:15px 0px">
</div></div></div><div><div><div><br><div class="gmail_quote">On Tue, Dec 18, 2012 at 1:06 AM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Amulya,<br>I looked through your script and I can't find the cause of the error. I've used this script on WRF-ARW output myself before (also postprocessed using ARWpost). I remember having some issues getting the script to run properly the first time, but I got around that somehow. The errors indicate that in this line<br>
<br>thetae=theta*exp((3.376/TLcl-0.00254)*mix*1.0+0.00081*mix)<br><br>one of theta, TLcl, or mix is not being computed correctly. I'm not sure which of those it is, although you could put SAY statements just above that line for each of the variables to see which it is. What I mean is, insert these lines<br>
say theta<br>say TLcl<br>say mix<br>Above the line I copied above. Then see what values get spit out. One of those is either a NaN of some sort or has characters where there should be numbers.<br><br><br>Jeff Duda<br><br>
<div class="gmail_quote"><div><div>On Mon, Dec 17, 2012 at 12:56 AM, Amulya Chevuturi <span dir="ltr"><<a href="mailto:amulya.chevuturi@gmail.com" target="_blank">amulya.chevuturi@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
<span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Dear GrADS users,</span><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">I have recently tried to use the <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a> script to plot Skew-T/Log-P graphs for my WRF (ARW) output. I use ARWPost as my post processing software. This is the problem that is cropping up. </div>
<div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<div>ga-> <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a> </div><div>No hardcopy metafile open</div><div>All files closed; all defined objects released;</div><div>All GrADS attributes have been reinitialized</div>
<div>Initial Time is 00Z08JAN2007</div><div>Create gif images as well (1=yes ; 0=no)</div><div>0</div><div>Please enter coordinated for sounding</div><div> </div><div>Latitude values between: 4.96654 and 50.0419</div><div>
Enter Latitude</div><div>33</div><div> </div><div>Longitude values between: 48.3967 and 105.603</div><div>Enter Longitude</div><div>78</div><div> </div><div>Overlay two time periods 1=yes ; 0=no</div><div>0</div><div> </div>
<div>Please enter time to plot</div><div>Available times: 1 to 17</div><div>5</div><div>Drawing temperature sounding.</div><div>Drawing dewpoint sounding.</div><div>Non-numeric args to numeric operation</div><div> Error occurred on line 2427</div>
<div> In file <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a></div><div> Error occurred on line 1106</div><div> In file <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a></div>
<div> Error occurred on line 69</div><div> In file <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a></div></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Though I think it is plotting the temperature and dewpoint. The rest of the parameters in the line 69 is not being plotted due to some errors in parameters. I am attaching my ctl file and <a href="http://skew.gs/" style="color:rgb(17,85,204)" target="_blank">skew.gs</a>. Please let me know how to correct the problem.</div>
<div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
Thanking you in advance. </div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
With regards,</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Amulya Chevuturi</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
Junior Research Fellow (Ph.D. Student)</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Jawaharlal Nehru University</div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
New Delhi, India</div>
<br></div></div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
<br>
</font></span><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><span><font color="#888888"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
<br>
</font></span></div></div><span><font color="#888888"><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></font></span></blockquote></div><br></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br><br>
</div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>