error in computation / new function

From: James Correia <jimmyc42_at_nyahnyahspammersnyahnyah>
Date: Mon, 2 Jun 2008 13:48:10 -0700

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 - 14:48:10 MDT

This archive was generated by hypermail 2.2.0 : Mon Jun 02 2008 - 16:43:00 MDT