<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px"><div id="yui_3_16_0_1_1439121182863_5446">Dear Mary,</div><div id="yui_3_16_0_1_1439121182863_6533"><br></div><div id="yui_3_16_0_1_1439121182863_6532">Thanks so much, and it really helps a lot!!! I will do check the map rotation paremeters.</div><div id="yui_3_16_0_1_1439121182863_7161"><br></div><div id="yui_3_16_0_1_1439121182863_7163">With my best regards</div><div id="yui_3_16_0_1_1439121182863_7162">Mark<br></div><div id="yui_3_16_0_1_1439121182863_5434"><span></span></div> <br><div class="qtdSeparateBR"><br><br></div><div style="display: block;" class="yahoo_quoted"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"> <font face="Arial" size="2"> On Friday, August 7, 2015 6:47 AM, Mary Haley <haley@ucar.edu> wrote:<br> </font> </div> <br><br> <div class="y_msg_container"><div id="yiv3868092052"><div><div dir="ltr"><div class="yiv3868092052gmail_default" style="font-size:small;">Dear Mark,</div><div class="yiv3868092052gmail_default" style="font-size:small;"><br clear="none"></div><div class="yiv3868092052gmail_default" style="font-size:small;">I don't think your "rot1_plot_good.ncl" is plotting the data in the correct location. I say this because it looks like your contours have a rough outline of what's supposed to be the Tibetan Plateau, but it doesn't look like that's where it is on the map.</div><div class="yiv3868092052gmail_default" style="font-size:small;"><br clear="none"></div><div class="yiv3868092052gmail_default" style="font-size:small;">Also, If I modify your code so that it actually uses the lat/lon arrays you read off the file to do the plotting, then the plot looks different. It does take longer to draw, because now you have to do the map transformation calculation behind the scenes, but the plot looks more correct to me. I just plotted it in a standard lat/lon projection, because when I use the parameters you provided, part of the contours go outside of the plot boundary.</div><div class="yiv3868092052gmail_default" style="font-size:small;"><br clear="none"></div><div class="yiv3868092052gmail_default" style="font-size:small;">Please compare the attached "rotated.png" image, which is from your original rot1_plot.ncl script, and the "non_rotated.png" image, which is the non-native image plotted using lat/lon arrays. I've called the script that generates the non-native plot "rot1_plot_good_mod.ncl" and attached it here as well.</div><div class="yiv3868092052gmail_default" style="font-size:small;"><br clear="none"></div><div class="yiv3868092052gmail_default" style="font-size:small;">Are you certain that you have the rotated parameters correct? This is usually the cause of why native data doesn't appear in the correct location.</div><div class="yiv3868092052gmail_default" style="font-size:small;"><br clear="none"></div><div class="yiv3868092052gmail_default" style="font-size:small;">--Mary</div><div class="yiv3868092052yqt5381827749" id="yiv3868092052yqt36469"><div class="yiv3868092052gmail_extra"><br clear="none"><div class="yiv3868092052gmail_quote">On Wed, Aug 5, 2015 at 4:16 AM, Mark Chan <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:cym263@yahoo.com" target="_blank" href="mailto:cym263@yahoo.com">cym263@yahoo.com</a>></span> wrote:<br clear="none"><blockquote class="yiv3868092052gmail_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>Dear Mary,</div><div><br clear="none"></div><div>Thanks so much for your response to my question!</div><div>The attached files are my ncl scripts and input data, and the script 'rot1_plot_good.ncl' works but 'rot2_plot_failed.ncl' always failed. Can you help to make the second script work?</div><div><br clear="none"></div><div dir="ltr">Actually I am trying to plot the rotated projection from WRF simulation ( see the attached namelist.wps), and wondering how to use the rotated message (such as: pole_lat = 77.61, pole_lon = -64.78) in my plot, so the plot could be exactly same as the model simulation (same as the ncview shows).</div><div dir="ltr"><br clear="none"></div><div dir="ltr">In addition, using the rotated message (77.61, -64.78), how is it possible to plot a global map, which shows the 'rotated map' looks like globally?</div><div dir="ltr"><br clear="none"></div><div dir="ltr">Thanks very much!</div><span class="yiv3868092052HOEnZb"><font color="#888888"></font></span><div dir="ltr">Mark<br clear="none"></div></div></div></blockquote></div><br clear="none"></div></div></div></div></div><br><br></div> </div> </div> </div></div></body></html>