<div dir="ltr"><div>Hi,</div><div><br></div><div>Given the aspect ratio of the underlying data (range of -74.25...-70.17 in lon, 36...46.9 in lat), it seems like the plot produced via gsn_csm_contour_map is about what one would expect ? I see you manually assigned those geographic coordinates to the grid. Are those correct? Do you have other reasons to believe the storm shape should be more circular, rather than elongated?</div><div><br></div><div>Not sure what to suggest here...</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Jan 2, 2019 at 10:46 AM 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"><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" target="_blank">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-m_3941588858086951207gmail_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>
_______________________________________________<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>