<div dir="ltr"><div class="gmail_default" style="font-size:small">Michele,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When you set:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><div> res@tfDoNDCOverlay = True</div></div><br class="gmail-Apple-interchange-newline"><br></div><div class="gmail_default" style="font-size:small">then you are effectively telling NCL that you know the *exact* map projection parameters that the data was measured on, and thus NCL will not not do any calculations to transform your data to the map parameters that you've set.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It looks like you have the lat/lon arrays for your data, so as a debugging test, I suggest plotting your data and force NCL to do the calculations. To do this, you need to remove this line:</div><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"> res@tfDoNDCOverlay = True</span><br class="gmail-Apple-interchange-newline"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">and set these two lines:</div></div><br> res@sfXArray = lon ; this will force NCL to do the map projection <br> res@sfYArray = lat<div><br></div><div><div class="gmail_default" style="font-size:small">Let ncl-talk know if the plot looks any better. If you continue to have problems with this, then can you provide the data file needed by your script?</div><div class="gmail_default" style="font-size:small"><br>Thanks,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><br><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 15, 2018 at 3:37 AM, Michele Petrini - CITG <span dir="ltr"><<a href="mailto:M.Petrini@tudelft.nl" target="_blank">M.Petrini@tudelft.nl</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="word-wrap:break-word;line-break:after-white-space">
<div>Dear help desk,</div>
<div><br>
</div>
<div>I am trying to plot Greenland ice thickness from a CISM simulation. Everything works fine - except that it looks that the ice sheet is a little bit shifted to the right (see attached figure: the small brown ice patch in the North-East should be
entirely on land, and that’s the same for other bits in the eastern margin). </div>
<div><br>
</div>
<div>The projection settings I am using coincides with those used in CISM (I am using lon/lat and projection settings as in the input file for my simulation), and the resources I set are </div>
<div><br>
</div>
<div>
<div> res@mpLimitMode = "Corners"</div>
<div> res@mpLeftCornerLatF = lat(0,0)</div>
<div> res@mpLeftCornerLonF = lon(0,0)</div>
<div> res@mpRightCornerLatF = lat(ny-1,nx-1)</div>
<div> res@mpRightCornerLonF = lon(ny-1,nx-1)</div>
<div> res@trGridType = "TriangularMesh" </div>
<div> res@mpProjection = "Stereographic" ; Change the map projection.</div>
<div> res@mpDataBaseVersion = "mediumres" </div>
<div> res@mpCenterLonF = -45.0</div>
<div> res@mpCenterLatF = 70.0</div>
<div> res@tfDoNDCOverlay = True</div>
</div>
<div><br>
</div>
<div>Do you know the reason for - and if it is possible to solve - this discrepancy? I attach also the contour_movie.ncl file if you want to have a look (the figure attached is just one snapshot). </div>
<div><br>
</div>
<div>thanks,</div>
<div>Michele</div>
<div><br>
</div>
<div></div>
</div>
<div style="word-wrap:break-word;line-break:after-white-space">
<div></div>
<div><br>
</div>
<div></div>
</div>
<div style="word-wrap:break-word;line-break:after-white-space">
<div></div>
<div><br>
</div>
<br>
<div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
******************************<wbr>******************************<wbr>******<br>
Michele Petrini<br>
<br>
PhD in Earth Science and Fluid Mechanics<br>
Department of Geoscience and Remote Sensing, TUDelft<br>
Delft, The Netherlands<br>
<br>
<a href="mailto:M.Petrini@tudelft.nl" target="_blank">M.Petrini@tudelft.nl</a> <br>
+39 3398367372<br>
michele.petrins (Skype alias)<br>
<br>
<br>
</div>
</div>
<br>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>