[ncl-talk] Attribute assignment type mismatch on wrf_map_overlays

[Ryan Connelly]

Hi,

I'm trying to plot a zoomed in fixed region, rather than the entire WRF
domain by default.  I keep encountering this error:

fatal:["NclAtt.c":262]:Attribute assignment type mismatch.

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 3351 in
file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4550 in
file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4619 in
file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 4702 in
file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 111 in
file wrf_dbz_discrete_area.ncl


Line 111 is the plot = wrf_map_overlays line.  Here's the relevant snippets
from my code:

  pltres = True
  mpres = True

; Plotting options for limited area consistent across domains
  mpres at mpLimitMode = "Corners"
  mpres at mpRightCornerLatF = 42.02
  mpres at mpRightCornerLonF = -69
  mpres at mpLeftCornerLatF = 45.88
  mpres at mpLeftCornerLonF = -75

[...]

; First get the variables we will need
; Both dbz and mdbz will be calculated using constant intercept parameters
; (as in early Reisner-2), with values of 8x10^6, 2x10^7, and 4x10^6 m^-4,
; for rain, snow, and graupel, respectively.

     mdbz = wrf_user_getvar(a,"mdbz",it)
     dbz = wrf_user_getvar(a,"REFL_10CM",it)
     z = wrf_user_getvar(a,"z",it)
     lat2d = wrf_user_getvar(a,"XLAT",it)
     lon2d = wrf_user_getvar(a,"XLONG",it)

     radar_Z = 10.^(dbz/10.)

     height = 2000.     ; 2 km
     radar_Z_plane  = wrf_user_intrp3d(radar_Z,z,"h",height,0.,False)
     dbz_plane = 10.*log10(radar_Z_plane)

     opts = res
     opts at sfXArray = lon2d
     opts at sfYArray = lat2d
     opts at cnFillOn = True
     opts at ContourParameters = (/ 10., 50., 2./)
     contour = wrf_contour(a,wks,dbz_plane(:,:),opts)
     pltres at LatLonOverlay = True
     plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)


Thanks,
Ryan

-- 
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161031/86c59a40/attachment.html