[ncl-talk] Sim reflectivity not plotting

[Ryan Connelly]

Hi Bill,

I had to explicitly define the sfXarray and sfYarray as lon and lat.  Which
is odd because I've never had to do that before - even the example script
from the MMM tutorial site doesn't.  But regardless, it works again.

I also implemented your suggestion to convert to linear Z for
interpolation, then back to dbz after.

Thanks for your help,
Ryan

On Thu, Oct 20, 2016 at 12:27 PM, Bill Ladwig <ladwig at ucar.edu> wrote:

> Hi Ryan,
>
> Here are a couple thoughts, that may or may not help.
>
> 1) Interpolating dbz (logarithmic data) using linear interpolation is
> probably not the best idea.  I would convert dbz to linear Z, then perform
> the interpolation, then convert back to dbz.  Note: Z = 10**(dbz/10.0)
>
> 2) Do a printVarSummary on dbz_plane and compare it to dbz to see if the
> coordinates got dropped.  I can't remember off the top of my head if that's
> going to happen.  Adding them back may fix your problem.
>
> 3) Note that the dbz produced by NCL is not identical to the result for
> REFL_10CM, and I would choose REFL_10CM over the NCL dbz result if you have
> it available.  This assumes you're using a microphysics scheme, like
> Thompson, that can produce simulated radar reflectivity.  Otherwise, the
> field will be 0s.
>
> 4) If all else fails, you can upload a file to our ftp via:
>
> ftp ftp.cgd.ucar.edu
> anonymous
> <use your email address for the password>
> cd incoming
> put ...
> put ...
> .
> .
> .
> quit
>
> After you upload it, you can send me a message offline to let me know the name of the file.  We can try to figure out what's going wrong.
>
> Hope this helps,
>
> Bill
>
>
> On Wed, Oct 19, 2016 at 9:09 PM, Ryan Connelly <rconne01 at gmail.com> wrote:
>
>> Hey Kacie,
>>
>> I am purposely interpolating it to a plane so that I can make an apples
>> to apples comparison to a 2 km interpolation of radar data without having
>> to deal with beam height uncertainty.  Otherwise you're right, you can just
>> call dbz directly, and I think it defaults to the lowest level?
>>
>> Ryan
>>
>> On Wed, Oct 19, 2016 at 9:52 PM, Kacie Shourd <Kacie.Shourd at dri.edu>
>> wrote:
>>
>>> Hey Ryan,
>>>
>>> Could you do a printVarSummary on dbz? I don't think it has to be
>>> interpolated to a plane. I could be wrong, but I know mdbz doesn't have to
>>> be. Change
>>>
>>> contour = wrf_contour(a,wks,dbz_plane(:,:),opts)
>>>
>>> to
>>>
>>> contour = wrf_contour(a,wks,dbz,opts)
>>>
>>> and see if that works.
>>>
>>> :)
>>> Best
>>> Kacie
>>>
>>> ------------------------------
>>> *From:* ncl-talk-bounces at ucar.edu [ncl-talk-bounces at ucar.edu] on behalf
>>> of Ryan Connelly [rconne01 at gmail.com]
>>> *Sent:* Wednesday, October 19, 2016 3:07 PM
>>> *To:* ncl-talk
>>> *Subject:* [ncl-talk] Sim reflectivity not plotting
>>>
>>> Hi,
>>>
>>> Any idea why sim reflec isn't plotting?  It's the baseline script which
>>> I've been using and modifying for many months, but suddenly today,
>>> nothing.  No error, but no contours generated.
>>>
>>> The WRF run has weather there: a Low, and accumulating rain, so it's
>>> surprising there's no sim reflec.
>>>
>>> I have do_radar_reflec = 1 in the namelist, and I checked that the
>>> values weren't outside the range being plotted.  I tried switching from
>>> "dbz" to "REFL_10CM", and neither produces any contours.
>>>
>>>      dbz = wrf_user_getvar(a,"dbz",it)
>>>      z = wrf_user_getvar(a,"z",it)
>>>
>>>      height = 2000.     ; 2 km
>>>      dbz_plane  = wrf_user_intrp3d(dbz,z,"h",height,0.,False)
>>>
>>>      opts = res
>>>      opts at cnFillOn = True
>>>      opts at ContourParameters = (/ 10., 50., 2./)
>>>      contour = wrf_contour(a,wks,dbz_plane(:,:),opts)     ; plot only
>>> lowest level
>>>      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
>>>
>>> PUBLIC RECORDS NOTICE: In accordance with NRS Chapter 239, this email
>>> and responses, unless otherwise made confidential by law, may be subject to
>>> the Nevada Public Records laws and may be disclosed to the public upon
>>> request.
>>>
>>
>>
>>
>> --
>> 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
ryan.connelly at stonybrook.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161020/8e723015/attachment.html