<div dir="ltr">Hi Mauro,<div>Please remember to always include ncl-talk on any response; that way others can assist. You stated that you did not have any error messages, but in the log file you sent me it ends like this:</div><div><br></div><div>(0) check_for_y_lat_coord: Warning: Data either does not contain<br>(0) a valid latitude coordinate array or doesn't contain one at all.<br>(0) A valid latitude coordinate array should have a 'units'<br>(0) attribute equal to one of the following values: <br>(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'<br>(0) check_for_lon_coord: Warning: Data either does not contain<br>(0) a valid longitude coordinate array or doesn't contain one at all.<br>(0) A valid longitude coordinate array should have a 'units'<br>(0) attribute equal to one of the following values: <br>(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'<br></div><div><br></div><div>The first 10 lines from the log file above all refer to latitude/longitude arrays not being correctly set, and are coming from one of these two lines:</div><div> plot_wr1 = gsn_csm_vector_scalar_map_polar(wks,u_avg(0,:,:),v_avg(0,:,:),geop_avg(0,:,:),res)<br> plot_er1 = gsn_csm_vector_scalar_map_polar(wks,u_regrid(0,:,:),v_regrid(0,:,:),z_regrid(0,:,:),res)<br></div><div>I would suggest that you use printVarSummary calls to examine the arrays that are being passed to the plotting functions in those two lines. As the vectors are being plotted correctly<br></div><div>in the 2nd plot, I would guess there's something amiss with your u_avg or v_avg arrays, and you should compare the output in printVarSummary between u_avg/v_avg and u_regrid/v_regrid to see if you see anything.</div><div>You are setting the 2-dimensional latitude and longitude arrays to the sfXArray/sfYArray resources (for the geopotential height scalar arrasy), but you are not setting vfXArray/vfYArray resources (for the wind component arrays).</div><div>This may or may not be the issue.</div><div><br></div><div>One can also set the 2-dimensional lat/lon arrays as attributes lat2d/lon2d attached to the arrays to be plotted. I do not know if you are doing that or not.</div><div><br></div><div>Note that setting gsnAddCyclic to True or False will not fix this particular issue.</div><div>Hope the above helps!</div><div>Adam</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Apr 7, 2021 at 7:43 AM Mauro Morichetti via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</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>
<div style="overflow-wrap: break-word;"><span style="font-size:14.6667px">Hi all,</span><br style="font-size:14.6667px">
<br style="font-size:14.6667px">
<span style="font-size:14.6667px">as the subject, I have a specific question about the function: gsn_csm_vector_scalar_map_polar.</span><br style="font-size:14.6667px">
<br style="font-size:14.6667px">
<span style="font-size:14.6667px">I am trying to plot a comparison between a WRF-Chem output and ERA5 data, for the geopotential height and the horizontal wind at 500 hPa, in arctic area.</span><br style="font-size:14.6667px">
<br style="font-size:14.6667px">
<span style="font-size:14.6667px">I have an issue plotting the wind vectors: in WRF-chem output (panel a - map attached) the wind doesn’t respect the domain as the values in the background. I had problem in the past with the gsnAddCyclic
function, but here I could not understand why the function doesn’t interfere with the vectors.</span><br style="font-size:14.6667px">
<br style="font-size:14.6667px">
<span style="font-size:14.6667px">Attached the script.</span><br style="font-size:14.6667px">
<br style="font-size:14.6667px">
<span style="font-size:14.6667px">Anyone can Help? Thanks</span><br>
<div><br>
</div>
<div>Regards</div>
<div><br>
---------<br>
<br>
Ph.D. Mauro Morichetti<br>
Department of Industrial Engineering and Mathematical Science<br>
Marche Polytechnic University, Ancona, Italy<br>
<a href="mailto:m.morichetti@pm.univpm.it" target="_blank">m.morichetti@pm.univpm.it</a>
</div>
<br>
<div></div>
</div>
<div><img src="cid:178af71643a6ea43201">
</div>
<div style="overflow-wrap: break-word;">
<div></div>
</div>
<div style="overflow-wrap: break-word;">
<div></div>
</div>
</div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div><br clear="all"><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><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>