<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hey David,<div class="">yes, the image that I attached in the first email was supposed to represents the final result. I know the “real” coordinate bounds just because I wrote a binary file directly from the postprocessing of the model with the real 2-D lat-lon coordinates in the native grid. Then, I printed the minimum and maximum of those data. Funny thing is that I wasn’t able to use these coordinates for my plot since they were somewhat randomly organized, still don’t know why (probably some fortran way of writing data…). However, the numbers that I attached in the first plot were not correct. These are the correct ones:</div><div class=""><br class=""></div><div class=""><div class="">lon_min<span class="Apple-tab-span" style="white-space:pre"> </span>-78.64264</div><div class="">lon_max<span class="Apple-tab-span" style="white-space:pre"> </span>25.00719</div><div class="">lat_min<span class="Apple-tab-span" style="white-space:pre"> </span>27.55382</div><div class="">lat_max<span class="Apple-tab-span" style="white-space:pre"> </span>67.51249</div><div class=""><br class=""></div><div class="">It puzzles me that they are not the same that you chose by replicating the plot “by eye”. And also, you changed the center of the map: I thought that the rotated pole coordinates that I had in the file were the only correct data :) </div><div class="">Using the script that Karin wrote, I should be able to create lat2d-lon2d array and unrotate the data, so that there is no need to play around with the map projection, right? </div><div class=""><br class=""></div><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class="Apple-interchange-newline">Guido Cioni</div><div class=""><a href="http://guidocioni.altervista.org" class="">http://guidocioni.altervista.org</a> </div><div class=""><br class=""></div></div><br class="Apple-interchange-newline" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><br class="Apple-interchange-newline"></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">Il giorno 30 nov 2015, alle ore 23:54, David Brown <<a href="mailto:dbrown@ucar.edu" class="">dbrown@ucar.edu</a>> ha scritto:</div><br class="Apple-interchange-newline"><div class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Hi Guido,</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Here's my take on this so far. You enclosed an image labelled "Bolam</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Model". Is that plot supposed to represent the actual grid of your</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">data? I am asking because it does not seem possible to reconcile the</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">coordinate information that you have provided with this plot. I am</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">attaching a couple of plots and a modified version of your script (it</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">just draws the map, since we do not have the data). I can</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">approximately duplicate the map projection of the Bolam plot by eye,</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">but I have to use different numbers than those you have supplied. The</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">script first draws a map based on the rotated lat/lon extent that you</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">supplied, after conversion into "real" lat/lon values. Then it draws</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">an approximation of the Bolam map.</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Mary supplied you with the NCL script from Karin that converts from</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">rotated coordinates to regular lat/lon. I used a C code that does the</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">same thing. You could use either method to create your own 2d lat/lon</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">coordinates for this grid. If you would like the C code let me know.</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">I am curious how you obtained the "real" corner values you used:</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@tfDoNDCOverlay=True</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@mpLimitMode = "Corners"</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@mpLeftCornerLatF = 30.35458</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@mpLeftCornerLonF = -75.14440</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@mpRightCornerLatF = 65.1125</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Res1@mpRightCornerLonF = 13.87071</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">These do not seem to me to be correct either.</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">It would help if you listed the complete set of rotated lat/lon</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">values. You could do it using ncl_filedump:</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">ncl_filedump -v rlat,rlon complete.nc</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">-dave</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">On Mon, Nov 30, 2015 at 2:18 PM, Mary Haley <</span><a href="mailto:haley@ucar.edu" style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">haley@ucar.edu</a><span style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">> wrote:</span><br style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><blockquote type="cite" style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Correct, NCL can't automatically "unrotate" the coordinates for you based on<br class="">the rotated coordinates.<br class=""><br class="">Karin Meier-Fleischer of DKRZ wrote a couple of functions that "unrotate"<br class="">the coordinates for you, and I've taken this example and her data file to<br class="">show how this might work for your data.<br class=""><br class="">See the attached script. The data file can be downloaded from:<br class=""><br class="">wget <a href="http://ftp.ucar.edu" class="">ftp.ucar.edu</a>:/pub/scd/haley/tas_mon_1961-1990_rectilinear_grid_2D.nc<br class=""><br class="">The rotated grid actually doesn't look that "rotated", because the center<br class="">lat/lon actually ends up being 0.0. But, hopefully you can use the<br class="">"unrot_lat" and "unrot_lon" functions to generate your own 2D lat/lon<br class="">arrays.<br class=""><br class="">--Mary<br class=""><br class=""><br class="">On Mon, Nov 30, 2015 at 1:31 PM, Guido Cioni <<a href="mailto:guidocioni@gmail.com" class="">guidocioni@gmail.com</a>> wrote:<br class=""><blockquote type="cite" class=""><br class="">I don't have the original 2D lat/lon values, that's the problem. However,<br class="">your script won't take into consideration the fact that coordinates are<br class="">rotated, right? My rlon spans from approximately 340 to 392 degrees, while<br class="">rlat from - 10 to +22.. Of course these are not the real coordinates to be<br class="">used for plotting the map, since they are expressed in the rotated<br class="">framework. Do you see what I mean?<br class="">I would like to provide you with the data file but is fairly big (gigs)<br class=""><br class="">Guido,<br class=""><br class="">I looked at your script further, and I think the issue is with the map<br class="">corners you've selected. You are using the 1D rotated coordinate arrays to<br class="">choose the corners, but this is likely not correct. What you really need are<br class="">the original 2D lat/lon values, as the two examples illustrated, but I don't<br class="">see that you have this information on the file.<br class=""><br class="">Did you try plotting the data using the coordinate arrays? I've created a<br class="">script for you to try, although it's untested because I don't have your data<br class="">file.<br class=""><br class="">Thanks,<br class=""><br class="">--Mary<br class=""><br class=""><br class="">On Mon, Nov 30, 2015 at 11:45 AM, Guido Cioni <<a href="mailto:guidocioni@gmail.com" class="">guidocioni@gmail.com</a>><br class="">wrote:<br class=""><blockquote type="cite" class=""><br class="">Sure, here it is. I think that the problem is that ncl is expecting 2-D<br class="">latitudes and longitudes while I have 1-D rlon and rlat..<br class=""><br class="">Yes I tried commenting the centers.<br class=""><br class="">m300382@mlogin103% ncl_filedump complete.nc<br class="">Copyright (C) 1995-2015 - All Rights Reserved<br class="">University Corporation for Atmospheric Research<br class="">NCAR Command Language Version 6.3.0<br class="">The use of this software is governed by a License Agreement.<br class="">See <a href="http://www.ncl.ucar.edu/" class="">http://www.ncl.ucar.edu/</a> for more details.<br class=""><br class="">Variable: f<br class="">Type: file<br class="">filename: complete<br class="">path: complete.nc<br class=""> file global attributes:<br class=""> CDI : Climate Data Interface version 1.4.7<br class="">(<a href="http://code.zmaw.de/projects/cdi" class="">http://code.zmaw.de/projects/cdi</a>)<br class=""> Conventions : CF-1.0<br class=""> history : Fri Nov 27 13:33:35 2015: cdo -f nc copy<br class="">2014011700_000.grib2 2014011700_001.grib2 2014011700_002.grib2<br class="">2014011700_003.grib2 2014011700_004.grib2 2014011700_005.grib2<br class="">2014011700_006.grib2 2014011700_007.grib2 2014011700_008.grib2<br class="">2014011700_009.grib2 2014011700_010.grib2 2014011700_011.grib2<br class="">2014011700_012.grib2 2014011700_013.grib2 2014011700_014.grib2<br class="">2014011700_015.grib2 2014011700_016.grib2 2014011700_017.grib2<br class="">2014011700_018.grib2 2014011700_019.grib2 2014011700_020.grib2<br class="">2014011700_021.grib2 2014011700_022.grib2 2014011700_023.grib2<br class="">2014011700_024.grib2 2014011700_025.grib2 2014011700_026.grib2<br class="">2014011700_027.grib2 2014011700_028.grib2 2014011700_029.grib2<br class="">2014011700_030.grib2 2014011700_031.grib2 2014011700_032.grib2<br class="">2014011700_033.grib2 2014011700_034.grib2 2014011700_035.grib2<br class="">2014011700_036.grib2 2014011700_037.grib2 2014011700_038.grib2<br class="">2014011700_039.grib2 2014011700_040.grib2 2014011700_041.grib2<br class="">2014011700_042.grib2 2014011700_043.grib2 2014011700_044.grib2<br class="">2014011700_045.grib2 2014011700_046.grib2 2014011700_047.grib2<br class=""> CDO : Climate Data Operators version 1.4.7<br class="">(<a href="http://code.zmaw.de/projects/cdo" class="">http://code.zmaw.de/projects/cdo</a>)<br class=""> dimensions:<br class=""> ncl_scalar = 1<br class=""> rlon = 706<br class=""> rlat = 434<br class=""> lev = 20<br class=""> time = 73 // unlimited<br class=""> variables:<br class=""> double rlon ( rlon )<br class=""> long_name : longitude in rotated pole grid<br class=""> units : degrees<br class=""> standard_name : grid_longitude<br class=""> axis : X<br class=""><br class=""> double rlat ( rlat )<br class=""> long_name : latitude in rotated pole grid<br class=""> units : degrees<br class=""> standard_name : grid_latitude<br class=""> axis : Y<br class=""><br class=""> character rotated_pole ( ncl_scalar )<br class=""> grid_mapping_name : rotated_latitude_longitude<br class=""> grid_north_pole_latitude : 44.999996<br class=""> grid_north_pole_longitude : 144<br class=""><br class=""> double lev ( lev )<br class=""> long_name : pressure<br class=""> units : Pa<br class=""> axis : Z<br class=""><br class=""> double time ( time )<br class=""> units : hours since 2014-01-17 00:00:00<br class=""> calendar : proleptic_gregorian<br class=""><br class=""> float gh ( time, lev, rlat, rlon )<br class=""> long_name : Geopotential Height<br class=""> units : gpm<br class=""> code : 5<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float t ( time, lev, rlat, rlon )<br class=""> long_name : Temperature<br class=""> units : K<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float u ( time, lev, rlat, rlon )<br class=""> long_name : U velocity<br class=""> units : m s**-1<br class=""> code : 2<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float v ( time, lev, rlat, rlon )<br class=""> long_name : V velocity<br class=""> units : m s**-1<br class=""> code : 3<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float w ( time, lev, rlat, rlon )<br class=""> long_name : Vertical velocity<br class=""> units : Pa s**-1<br class=""> code : 8<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float pv ( time, lev, rlat, rlon )<br class=""> long_name : Potential vorticity<br class=""> units : K m**2 kg**-1 s**-1<br class=""> code : 14<br class=""> grid_mapping : rotated_pole<br class=""><br class=""> float mterh ( time, rlat, rlon )<br class=""> long_name : Model terrain height<br class=""> units : m<br class=""> code : 7<br class=""> grid_mapping : rotated_pole<br class=""><br class="">Il 30 nov 2015 6:28 PM, "Mary Haley" <<a href="mailto:haley@ucar.edu" class="">haley@ucar.edu</a>> ha scritto:<br class=""><blockquote type="cite" class=""><br class="">Guido,<br class=""><br class="">I'm hoping Dave will weigh in here, but it would help if we could at<br class="">least see the output from:<br class=""><br class="">ncl_filedump pv.nc<br class=""><br class="">Sometimes the global attributes of the file will give some valuable<br class="">information about the correct map projection parameters.<br class=""><br class="">Did you actually try commenting out the mpCenterLatF and mpCenterLonF<br class="">settings? In the script you included, these two resources are still being<br class="">set.<br class=""><br class="">It's possible that the data has already been rotated for you.<br class=""><br class="">--Mary<br class=""><br class=""><br class=""><br class=""><br class="">On Mon, Nov 30, 2015 at 1:33 AM, Guido Cioni <<a href="mailto:guidocioni@gmail.com" class="">guidocioni@gmail.com</a>><br class="">wrote:<br class=""><blockquote type="cite" class=""><br class="">Hi all,<br class="">I’m trying to produce a contour plot with data on a rotated grid. The<br class="">original data look like that<br class=""><br class=""><br class="">Since I have the rotated north pole latitude (44.999996) and longitude<br class="">(144.0) in my data file I used the second example of the support page on<br class="">rotated grid.<br class=""><br class="">Res1@tfDoNDCOverlay=True<br class="">Res1@mpLimitMode = "Corners"<br class="">Res1@mpLeftCornerLatF = lat(0)<br class="">Res1@mpLeftCornerLonF = lon(0)<br class="">Res1@mpRightCornerLatF = lat(nlat-1)<br class="">Res1@mpRightCornerLonF = lon(nlon-1)<br class="">Res1@mpCenterLatF = 90 - north_pole_lat<br class="">Res1@mpCenterLonF = 180 + north_pole_lon<br class=""><br class="">However I get this result ( I believe this is Australia):<br class=""><br class=""><br class="">After several attempts I was able to plot the data using as bounds the<br class="">“real” coordinates of the rotated projection, without specifying the center,<br class="">e.g. :<br class=""><br class="">Res1@tfDoNDCOverlay=True<br class="">Res1@mpLimitMode = "Corners"<br class="">Res1@mpLeftCornerLatF = 30.35458<br class="">Res1@mpLeftCornerLonF = -75.14440<br class="">Res1@mpRightCornerLatF = 65.1125<br class="">Res1@mpRightCornerLonF = 13.87071<br class=""><br class="">But the contour plot is still slightly shifted to the west as you can<br class="">see in this image:<br class=""><br class=""><br class="">I know that in the data file my coordinates span from<br class="">rlon={340,392.875} and rlat={-9.962499,22.512501}.<br class=""><br class="">What am I missing?<br class="">Is there any additional documentation on rotated grid?<br class=""><br class="">Attached is my NCL script (data are too big).<br class=""><br class=""><br class=""><br class="">Guido Cioni<br class=""><a href="http://guidocioni.altervista.org" class="">http://guidocioni.altervista.org</a><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">ncl-talk mailing list<br class="">ncl-talk@ucar.edu<br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""><br class=""></blockquote><br class=""></blockquote></blockquote><br class=""></blockquote><br class=""><br class="">_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class=""><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br class=""><br class=""></blockquote><span id="cid:60183A5E-33D4-4193-9C9D-C234CF85947F@mpi.zmaw.de"><pv_test.000001.png></span><span id="cid:3930CFB2-496B-4DFB-A95B-A92099296A41@mpi.zmaw.de"><pv_test.000002.png></span><span id="cid:BCB155D3-D01A-4156-9EDF-C3E19266EB09@mpi.zmaw.de"><plot_pv-mod.ncl></span></div></blockquote></div><br class=""></div></body></html>