<div dir="ltr">Hi,<div><br></div><div>I'm trying to plot a zoomed in fixed region, rather than the entire WRF domain by default. I keep encountering this error:</div><div><br></div><div><div><font face="monospace, monospace">fatal:["NclAtt.c":262]:Attribute assignment type mismatch.</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 3351 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4550 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4619 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4702 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 111 in file wrf_dbz_discrete_area.ncl</font></div></div><div><br></div><div><br></div><div>Line 111 is the plot = wrf_map_overlays line. Here's the relevant snippets from my code:</div><div><br></div><div><div><font face="monospace, monospace"> pltres = True</font></div><div><font face="monospace, monospace"> mpres = True</font></div><div><font face="monospace, monospace"> </font></div><div><font face="monospace, monospace">; Plotting options for limited area consistent across domains</font></div><div><font face="monospace, monospace"> mpres@mpLimitMode = "Corners"</font></div><div><font face="monospace, monospace"> mpres@mpRightCornerLatF = 42.02</font></div><div><font face="monospace, monospace"> mpres@mpRightCornerLonF = -69</font></div><div><font face="monospace, monospace"> mpres@mpLeftCornerLatF = 45.88</font></div><div><font face="monospace, monospace"> mpres@mpLeftCornerLonF = -75</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">[...]</font></div><div><font face="monospace, monospace"><br></font></div><div><div><font face="monospace, monospace">; First get the variables we will need</font></div><div><font face="monospace, monospace">; Both dbz and mdbz will be calculated using constant intercept parameters</font></div><div><font face="monospace, monospace">; (as in early Reisner-2), with values of 8x10^6, 2x10^7, and 4x10^6 m^-4,</font></div><div><font face="monospace, monospace">; for rain, snow, and graupel, respectively.</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> mdbz = wrf_user_getvar(a,"mdbz",it)</font></div><div><font face="monospace, monospace"> dbz = wrf_user_getvar(a,"REFL_10CM",it)</font></div><div><font face="monospace, monospace"> z = wrf_user_getvar(a,"z",it)</font></div><div><font face="monospace, monospace"> lat2d = wrf_user_getvar(a,"XLAT",it)</font></div><div><font face="monospace, monospace"> lon2d = wrf_user_getvar(a,"XLONG",it)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> radar_Z = 10.^(dbz/10.)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> height = 2000. ; 2 km</font></div><div><font face="monospace, monospace"> radar_Z_plane = wrf_user_intrp3d(radar_Z,z,"h",height,0.,False)</font></div><div><font face="monospace, monospace"> dbz_plane = 10.*log10(radar_Z_plane)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"> opts = res</font></div><div><font face="monospace, monospace"> opts@sfXArray = lon2d</font></div><div><font face="monospace, monospace"> opts@sfYArray = lat2d</font></div><div><font face="monospace, monospace"> opts@cnFillOn = True</font></div><div><font face="monospace, monospace"> opts@ContourParameters = (/ 10., 50., 2./)</font></div><div><font face="monospace, monospace"> contour = wrf_contour(a,wks,dbz_plane(:,:),opts)</font></div><div><font face="monospace, monospace"> pltres@LatLonOverlay = True</font></div><div><font face="monospace, monospace"> plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)</font></div></div><div><br></div><div><br></div><div>Thanks,</div><div>Ryan</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font size="1">Ryan Connelly</font></div><div dir="ltr"><div dir="ltr" style="font-size:12.8px"><font color="#000000"><font size="1">M.S. Student in Atmospheric Sciences, </font><span style="font-size:x-small">Stony Brook University</span><br></font></div><div dir="ltr" style="font-size:12.8px"><div><font color="#000000"><font size="1">B.S. in Meteorology with Minors in Mathematics and GIS, </font><span style="font-size:x-small">Valparaiso University</span></font></div></div><div><a href="mailto:rconne01@gmail.com" style="font-size:x-small" target="_blank">rconne01@gmail.com</a><br></div><div><font size="1"><a href="mailto:ryan.connelly@stonybrook.edu" target="_blank">ryan.connelly@stonybrook.edu</a></font></div></div></div></div></div></div></div></div></div></div>
</div></div>