<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Dear Mansur,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">There are a couple of issues with your script.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The first is that you are setting:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div> res@tfDoNDCOverlay = True ; do not transform data</div><br class="gmail-Apple-interchange-newline"></div><div class="gmail_default" style="font-size:small">This tells NCL to ignore any lat/lon information attached to your data, and to instead assume that the map resources (res@mpXXXX) have been set up exactly to match the projection of your data. Remove that line.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The second is that you are providing conflicting information to the plotting program by adding too much lat/lon information to your data variable. You have:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"> lat2d = fspan(13,47,1202)</div><div class="gmail_default"> lon2d = fspan(-95,-58,1202)</div><div class="gmail_default"> printVarSummary(u)</div><div class="gmail_default"> printVarSummary(lat2d)</div><div class="gmail_default"> printVarSummary(lon2d)</div><div class="gmail_default"> u!1="lat2d"</div><div class="gmail_default"> u!2="lon2d"</div><div class="gmail_default"><br></div><div class="gmail_default"> u@lat2d=lat2d</div><div class="gmail_default"> u@lon2d=lon2d</div><div class="gmail_default"><br></div><div class="gmail_default"> u&lat2d@units = "degrees_north"</div><div class="gmail_default"> u&lon2d@units = "degrees_east"'</div><div class="gmail_default"><br></div><div class="gmail_default">The first issue is that these two lines should NOT be set at all:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"> u@lat2d=lat2d</div><div class="gmail_default"> u@lon2d=lon2d</div><div class="gmail_default"><br></div><div class="gmail_default">They are only to be used in the case where you have curvilinear data, which is data represented by two-dimensional lat/lon arrays. The way you are setting up your lat/lon arrays is for rectilinear data, which means you have one-dimensional latitudes for all your Y locations, and one-dimensional longitudes for your X locations.</div><div class="gmail_default"><br></div><div class="gmail_default">While it's not wrong, I suggest not using "lat2d" and "lon2d" for your lat / lon variable names, because this is misleading. Instead, I suggest:</div><div class="gmail_default"><br></div></div></div><div class="gmail_default" style="font-size:small"><div class="gmail_default"> lat = fspan(13,47,1202)</div><div class="gmail_default"> lon = fspan(-95,-58,1202)</div><div class="gmail_default"> printVarSummary(u)</div><div class="gmail_default"> printVarSummary(lat)</div><div class="gmail_default"> printVarSummary(lon)</div><div class="gmail_default"> u!1="lat"</div><div class="gmail_default"> u!2="lon"</div><div class="gmail_default"><br></div><div class="gmail_default"> u&lat@units = "degrees_north"</div><div class="gmail_default"> u&lon@units = "degrees_east"'</div></div></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To learn more about rectilinear versus curvilinear, you might want to see the "plotting data on a map" examples page:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml">http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml</a><br></div><div><br></div><div><div class="gmail_default" style="font-size:small">Regards,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><br></div><br><div class="gmail_quote"><div dir="ltr">On Sun, Dec 2, 2018 at 9:28 PM Mansur Ali Jisan <<a href="mailto:jisan.mansur@gmail.com">jisan.mansur@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear NCL Community,<div><br></div><div>I'm facing some problems in projecting my data properly on a map. I'm working with a NetCDF file from a parametric wind model for the case of Hurricane Irma. The variable that I'm trying to plot is u10. The problem is that it doesn't have any coordinate values attached to it. But from the input data, I know the extent of latitude and longitude. So, I defined that in the NCL script but from the visualization, the location of the center position was shown in a wrong position.</div><div><br></div><div>Following is the script that I tried:</div><div><br></div><div><b>load </b>"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div><b>load</b> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>begin</div><div>;---Read data</div><div> f1 = addfile("<a href="http://boundary_parametric.nc" target="_blank">boundary_parametric.nc</a>","r")</div><div> u = f1->u10 ; 97 x 1202 x 1202</div><div> lat2d = fspan(13,47,1202)</div><div> lon2d = fspan(-95,-58,1202)</div><div> printVarSummary(u)<br></div><div> printVarSummary(lat2d)</div><div> printVarSummary(lon2d)</div><div> u!1="lat2d"<br></div><div> u!2="lon2d"</div><div><br></div><div> u@lat2d=lat2d</div><div> u@lon2d=lon2d</div><div><br></div><div> u&lat2d@units = "degrees_north"</div><div> u&lon2d@units = "degrees_east"</div><div><br></div><div><div> printVarSummary(u)</div><div><br></div><div><br></div><div>;---Start the graphics</div><div> wks = gsn_open_wks("png","irma-u10")</div><div> gsn_define_colormap(wks,"matlab_jet")</div><div>;---Set some resources</div><div> res = True</div><div> res@gsnMaximize = True ; maximize plot in frame</div><div> res@gsnAddCyclic = False</div><div> res@cnFillOn = True ; turn on contour fill</div><div> res@cnLinesOn = False ; turn off contour lines</div><div> res@cnLineLabelsOn = False ; turn off contour line labels</div><div> res@gsnScalarContour = True ; contours desired</div><div> res@tiMainString = "Hurricane Irma"</div><div> res@trGridType = "TriangularMesh"</div><div> res@gsnScalarContour = True</div><div> res@cnLevelSelectionMode = "AutomaticLevels"</div><div> res@mpProjection = "CylindricalEquidistant" ; projection<br></div><div> res@mpDataBaseVersion = "MediumRes" ; use more detailed maps</div></div><div> res@mpLimitMode = "LatLon"<br></div><div><div> res@mpMinLatF = 13.0</div><div> res@mpMaxLatF = 47.0</div><div> res@mpMinLonF = -95.0</div><div> res@mpMaxLonF = -58.0</div><div> res@tfDoNDCOverlay = True ; do not transform data</div><div><br></div><div> plot = gsn_csm_contour_map(wks,u(3,:,:),res)</div><div><br></div><div>end</div></div><div><br></div><div><b>NCL_Filedump</b> for the input file is the following:</div><div><br></div><div><div>Variable: f</div><div>Type: file</div><div>filename: boundary_parametric</div><div>path: <a href="http://boundary_parametric.nc" target="_blank">boundary_parametric.nc</a></div><div> file global attributes:</div><div> time-units : start time: 2017-09-10 00:06:00 +00:00</div><div> dimensions:</div><div> y_axis = 1202</div><div> x_axis = 1202</div><div> time = 97 // unlimited</div><div> variables:</div><div> double y_axis ( y_axis )</div><div> units : unit-less</div><div><br></div><div> double x_axis ( x_axis )</div><div> units : unit-less</div><div><br></div><div> float time ( time )</div><div> units : hours</div><div><br></div><div> float utx ( time )</div></div><div><div> units :</div><div><br></div><div> float uty ( time )</div><div> units :</div><div><br></div><div> float rmax ( time )</div><div> units :</div><div><br></div><div> float wsmax ( time )</div><div> units :</div><div><br></div><div> float u10 ( time, y_axis, x_axis )</div><div> units : m/s</div><div><br></div><div> float v10 ( time, y_axis, x_axis )</div><div> units : m/s</div><div><br></div><div> float utop ( time, y_axis, x_axis )</div><div> units : m/s</div><div><br></div><div> float vtop ( time, y_axis, x_axis )</div><div> units : m/s</div><div> <br></div></div><div>Please let me know if any additional information needed. </div><div><br></div><div>Thanks in advance!</div><div><br></div><div>Mansur</div><div><br></div><div>-- <br><div dir="ltr" class="gmail-m_-8249571406907405157gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><b style="font-family:Cambria;font-size:12pt"><span style="font-size:9pt;font-family:Arial">Mansur Ali Jisan</span></b></div><div><span style="font-family:Cambria;font-size:12pt"><span style="font-size:9pt;font-family:Arial">Ph.D. Student</span></span></div><div><span style="font-family:Cambria;font-size:12pt"><span style="font-size:9pt;font-family:Arial">URI Graduate School of Oceanography, RI 02882</span></span></div><div><br></div><div><div><div><a rel="nofollow" href="mailto:hasancee@iut-dhaka.edu" target="_blank"></a></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a><br>
</blockquote></div></div></div>