<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><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 need help with contour map plot for a Hurricane wind field. At first, I used gsn_csm_contour function for making the plot and the shape of the storm looked correct. But when I tried to make a contour map plot using gsn_csm_contour_map, the shape got distorted. I tried to solve this by setting the res@mpShapeMode = "FreeAspect". While it improved the contour plot but the map got stretched. Similarly, I set up the map projection to cylindrical equidistant but still got the same result. I have added the code and PrintVarSummary. Any suggestions will be very helpful. </div><div><br></div><div>Wishing you a Happy New Year!</div><div><br></div><div>--------------------------------------------------------------------------------------------------------</div><div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>begin</div><div>  ntim = 318</div><div>  nlat = 1202</div><div>  mlon = 1202</div><div><br></div><div>  f1 = addfile("<a href="http://bound_param_diag_var_rough.nc">bound_param_diag_var_rough.nc</a>","r")</div><div>  u = f1->um_bot   ; 97 x 1202 x 1202</div><div>  v = f1->vm_bot</div><div>  ws = sqrt(u^2+v^2)</div><div> </div><div> lat  = fspan(36.07092, 46.86241,nlat)</div><div>  lon  = fspan(-74.25677, -70.17657,mlon)</div><div><br></div><div>  lat!0= "lat"                   ; name dimension</div><div>  lat@units = "degrees_north"    ; associate units attribute</div><div>  lat&lat   = lat                ; coordinate variable</div><div><br></div><div>  lon!0= "lon"</div><div>  lon@units = "degrees_east"</div><div>  lon&lon   = lon</div><div><br></div><div>  u!0 = "time"</div><div>  u!1 = "lat"</div><div>  u!2 = "lon"</div></div><div><div>  u&lat = lat</div><div>  u&lon = lon</div><div>  v!0 = "time"</div><div>  v!1 = "lat"</div><div>  v!2 = "lon"</div><div><br></div><div>  v&lat = lat</div><div>  v&lon = lon</div><div><br></div><div>  ws!0 = "time"</div><div>  ws!1 = "lat"</div><div>  ws!2 = "lon"</div><div><br></div><div>  ws&lat = lat</div><div>  ws&lon = lon</div><div><br></div><div> printVarSummary(lat)</div><div> print("-----")</div><div> printVarSummary(lon)</div><div> print("-----")</div><div> printVarSummary(u)</div><div> print("-----")</div></div><div><br></div><div><div>  ws = mask(ws,ws.lt.4,False)</div><div>  u = mask(u,u.lt.4,False)</div><div>  v = mask(v,v.lt.4,False)</div><div><br></div><div>  wks = gsn_open_wks("png","gsn_csm_contour_map")</div><div>  res                          = True<br></div><div>  res@gsnDraw        = False</div><div>  res@gsnFrame      = False</div><div>  res@gsnAddCyclic   = False</div><div>  res@gsnMaximize    = True</div><div>  res@cnFillOn            = True</div><div>  res@cnLinesOn        = False</div><div>  res@cnFillMode        = "AreaFill"     </div><div>  res@cnFillPalette      = "matlab_jet"       </div></div><div><br></div><div>;  res@mpShapeMode = "FreeAspect"<br></div><div><div>;  res@vpWidthF = 0.6</div><div>;  res@vpHeightF = 0.5</div></div><div><br></div><div><div>  res@tmXBOn = True</div><div>  res@tmYLOn = True</div><div>  res@tiMainString   = "gsn_csm_contour_map"</div><div>  res@tmXBMode = "Explicit"</div><div>  res@tmXBValues = (/-66, -67, -68, -69, -70, -71,  -72,  -73,   -74,  -75, -76, -77/)</div><div>  res@tmXBLabels = (/"66W","67W","68W","69W","70W","71W","72W","73W","74W","75W","76W","77W"/)</div><div><br></div><div>  res@tmYLMode = "Explicit"</div><div>  res@tmYLValues = (/34, 35,  36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47/)</div><div>  res@tmYLLabels = (/"34N","35N","36N","37N","38N","39N","40N","41N","42N","43N","44N","45N","46N","47N"/)</div><div><br></div><div>  res@cnLevelSelectionMode = "ManualLevels"</div><div>  res@cnMinLevelValF       = 10</div><div>  res@cnMaxLevelValF       = 80</div><div>  res@cnLevelSpacingF      = 2.5</div><div><br></div><div>  res@mpFillOn = False</div><div>  res@mpGeophysicalLineThicknessF = 2.5</div><div>  res@mpOutlineBoundarySets     = "geophysicalandusstates"</div><div>  res@mpDataBaseVersion = "MediumRes"</div><div>  res@mpDataSetName = "Earth..4"</div><div>  res@mpGeophysicalLineColor = "Black"</div><div><br></div><div>  res@mpLimitMode = "LatLon"</div><div>  res@mpMinLatF = 36.0</div><div>  res@mpMaxLatF = 46.0</div><div>  res@mpMinLonF = -76.0</div><div>  res@mpMaxLonF = -70.0</div><div><br></div><div><div>  plot = gsn_csm_contour_map(wks,ws(22,:,:),res)</div><div><br></div><div>draw(plot)</div><div><br></div><div>frame(wks)</div><div><br></div><div>end</div></div><div><br></div><div><br></div><div><div>Variable: lat</div><div>Type: float</div><div>Total Size: 4808 bytes</div><div>            1202 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:   [lat | 1202]</div><div>Coordinates:</div><div>            lat: [36.07092..46.86241]</div><div>Number Of Attributes: 1</div><div>  units :       degrees_north</div><div>(0)     -----</div></div><div><br></div><div><div>Variable: lon</div><div>Type: float</div><div>Total Size: 4808 bytes</div><div>            1202 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:   [lon | 1202]</div><div>Coordinates:</div><div>            lon: [-74.25677..-70.17657]</div><div>Number Of Attributes: 1</div><div>  units :       degrees_east</div><div>(0)     -----</div></div><div><br></div><div><br></div><div><div>Variable: u</div><div>Type: float</div><div>Total Size: 1837790688 bytes</div><div>            459447672 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:   [time | 318] x [lat | 1202] x [lon | 1202]</div><div>Coordinates:</div><div>            lat: [36.07091903686523..46.86241149902344]</div><div>            lon: [-74.25676727294922..-70.17656707763672]</div><div>Number Of Attributes: 1</div><div>  units :</div><div>(0)     -----</div></div><div><br></div><div><br></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></div></div></div></div></div></div>