Re: error in computation / new function

Hi James,

The wrf_rip_dbz function hasn't been implemented yet, because we need
to talk to the WRF developers about the correct input values and such.
I'll see where this is.

>From looking at the wrf_user_intrp3d code, it doesn't look like this
function specifically checks for missing values. I'm not sure if it's
supposed to or not. I've CC'ed wrfhelp to see if they can provide
some more input on the expected behavior of this function with
respect to missing values.

--Mar

On Mon, 2 Jun 2008, James Correia wrote:

> Is there anyway to obtain one of the future release functions (wrf_rip_dbz)?
>
> I have been trying to make a simple radar reflectivity cross section work
> with NCL:
> a = qrain + 0.2*qsnow + 2.*graup
> b = a*264083.11
> radar = 17.8*log10(b)
>
> Now the results, consistent with a previous log10 discussion, gives nan for
> values less than zero and -inf for zero values. Looks like there are LOTS of
> missing values from my computation, but I am not sure how to get around
> this.
>
> but for this plot I am making a vertical cross section between points A and
> B.
>
> rad_plane = wrf_user_intrp3d(radar,p,"v",plane,0.,True)
>
> the maximum value for radar is 67, but the max value after interpolation is
> -28.
>
> Perhaps the interpolation function (wrf_user_intrp3d) doesnt handle missing
> values?
>
> --
> James Correia Jr.
> Post Doc
> Climate Physics Group, PNNL
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 02 2008 - 16:14:16 MDT