[ncl-talk] Sim reflectivity not plotting

[Bill Ladwig]

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161020/a68f27e0/attachment.html