[ncl-talk] Attribute assignment type mismatch on wrf_map_overlays
Ryan Connelly
rconne01 at gmail.com
Mon Oct 31 20:43:03 MDT 2016
Hi Dennis,
That was it! Thanks!
Ryan
On Monday, October 31, 2016, Dennis Shea <shea at ucar.edu> wrote:
> I speculate the following:
>
> your code has
> mpres at mpRightCornerLonF = -69 ; <=== integer
> mpres at mpLeftCornerLonF = -75 ; "
>
> ===
>
> The error message states:
>
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 3351 in
> file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
>
> 3350 if ( opts at mpRightCornerLonF .lt. 0.0 ) then
> 3351 opts at mpRightCornerLonF = opts at mpRightCornerLonF + 360.0
> 3352 end if
>
> line number 3351 results in a *float* because of the 360.
>
> SO ...... You are trying to put a float into an integer variable.
>
> NCL is a *strongly typed language*. It will not allow placing a float into
> an integer variable .... Hence, the error message!!!
>
> ---
> I think a solution is:
>
>
> your code has
> mpres at mpRightCornerLonF = -69. ; <=== float
> mpres at mpLeftCornerLonF = -75. ; "
>
> On Mon, Oct 31, 2016 at 5:10 PM, Ryan Connelly <rconne01 at gmail.com
> <javascript:_e(%7B%7D,'cvml','rconne01 at gmail.com');>> wrote:
>
>> Hi Mary,
>>
>> Rick just asked the same thing. It's 6.3.0.
>>
>> Thanks,
>> Ryan
>>
>>
>> On Monday, October 31, 2016, Mary Haley <haley at ucar.edu
>> <javascript:_e(%7B%7D,'cvml','haley at ucar.edu');>> wrote:
>>
>>> What version of NCL are you running?
>>>
>>> ncl -V
>>>
>>> It doesn't really make sense that it's failing on the wrf_map_overlays
>>> line, because I don't think there's an attribute assignment happening
>>> there. But, if I know what version of NCL you have, then I can look at
>>> WRFUserARW.ncl for that version to see what's happening on the particular
>>> line numbers it's complaining about.
>>>
>>> --Mary
>>>
>>>
>>> On Mon, Oct 31, 2016 at 1:37 PM, Ryan Connelly <rconne01 at gmail.com>
>>> wrote:
>>>
>>>> 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
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>
>> --
>> 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 <javascript:_e(%7B%7D,'cvml','rconne01 at gmail.com');>
>> ryan.connelly at stonybrook.edu
>> <javascript:_e(%7B%7D,'cvml','ryan.connelly at stonybrook.edu');>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu <javascript:_e(%7B%7D,'cvml','ncl-talk at ucar.edu');>
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
--
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/d7674a2a/attachment.html
More information about the ncl-talk
mailing list