<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">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">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_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>