[ncl-talk] how can I avoid dividing by zero?
Sri Nandini
snandini at marum.de
Wed Aug 22 03:56:23 MDT 2018
Thank you
I used it and have seen that i have data points except some gridpoints which point to nan and the fillvalue both.
I think then there is no data over these points.
On Aug 22, 2018 10:46:56 AM, Barry Lynn wrote:
> Hello Sri:
> I would try it this way.
>
> Set > y> to a _FillValue value where it is equal to 0, and then can do the divide:
> > y> @_FillValue = 1.E36 ; or something like this.
> > print(where(y.ne.0,y,> @_FillValue)) ; I am assuming this works.
> > If nothing prints then maybe you have no values.
> > In any case, try the next two lines.
> > y = > where> (y.ne.0,y,y at _FillValue)
> yinv = 1. / y> However, just to check, I would
>
> On Wed, Aug 22, 2018 at 11:39 AM Sri Nandini <> snandini at marum.de> > wrote:
> > Hello
> >
> >
how can I remove the zero values to avoid dividing by zero the error I got is as the following:
> >
> >
fatal:divide: Division by 0, Can't continue
> >
fatal:Div: operator failed, can't continue
> >
> >
I have tried using yinv = 1. / where(y.ne.0, y, y at _FillValue) but now it gives me all nan values
> >
> >
Heres the bit of my code
> >
> >
T41 = f->Z3(:,{850},:,:)
> >
> >
printVarSummary(T41)
> >
T41 at _FillValue = -9.96921e+36
> >
> >
aveX = dim_avg_n_Wrap(T41,0)
> >
printVarSummary(aveX) ; (lat,lon)
> >
T41 at _FillValue = -9.96921e+36
> >
> >
aveX = 1. / where(aveX.ne.0, aveX, aveX at _FillValue)
> >
> >
aveX1 and aveX2 below are of same type as aveX.
> >
XLAT = conform ( aveX2 , xlat , 0 )
> >
printVarSummary(XLAT)
> >
> >
egr = eady_growth_rate(aveX2, aveX1, aveX, XLAT, 0, 1)
> >
> >
Output
> >
===========================
> >
Variable: egr
> >
Type: double
> >
Total Size: 110592 bytes
> >
13824 values
> >
Number of Dimensions: 2
> >
Dimensions and sizes: [lat | 96] x [lon | 144]
> >
Coordinates:
> >
lat: [ -90.. 90]
> >
lon: [ 0..357.5]
> >
Number Of Attributes: 3
> >
_FillValue : -9.969209968386869e+36
> >
long_name : maximum eady growth rate
> >
units :
> >
(0) maximum eady growth rate : min=nan max=nan
> >
> >
> >
Any advice is much appreciated
> >
_______________________________________________
> >
ncl-talk mailing list
> > ncl-talk at ucar.edu
> >
List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> --
> Barry H. Lynn, Ph.D> Senior Associate Scientist, Lecturer,
> > The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> > C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>
>
>
> >
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180822/a4f4ba8b/attachment.html>
More information about the ncl-talk
mailing list