<div dir="ltr"><div>I have not tried your script. However. I scanned it<br><br></div>[snip]<br><div> plot = gsn_csm_contour_map(xwks,LH(0,nLH,:,:),res_LH)<br> area_avg = wgt_areaave(LH(:,nLH,{19.5:23.5},{84.5:89.5}),1.0,1.0,1)<br> area_avg@_FillValue = -9999.9<br></div><div>[snip]<br><br></div><div>Note that the last argument of wgt_areaave has two settings. You specified one (1). Please read the following description of the opt argument.<br></div><div>===<br><em>opt</em> <br><br>If <i>opt</i> = 0, the area average is calculated using available non-missing data.
If <i>opt</i> = 1, then if <b>any</b> point in <em>q</em> is missing, the area average
is not computed. In this case, it will be set to the missing value, which is indicated by
<em>q</em>@<a href="http://www.ncl.ucar.edu/Document/Language/fillval.shtml">_FillValue</a>, or the default missing value if <em>q</em>@_FillValue is not set<br>===<br><br></div><div>I think you should set opt=0<br><br>===<br></div><div>Further ... given the (likely) different number of non-missing values for different time steps, I suggest you view the areal averages with caution.</div><div><br>===<br></div><div>Also ... explore the data you are plotting<br><br>x = LH(:,nLH,{19.5:23.5},{84.5:89.5}<br></div><div>printVarSummary(x)<br><br>opt = True<br>opt@PrintStat = True<br>statb = <strong>stat_dispersion</strong>(x, opt )
<br><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 13, 2016 at 4:50 AM, Geeta Geeta <span dir="ltr"><<a href="mailto:geetag54@yahoo.com" target="_blank">geetag54@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,Sans-Serif;font-size:16px"><div><span>Hi all. <br></span></div><div dir="ltr"><span>I have Global data of latent heat at 12 different levels. I have to compute the areal average of each layer , for one sector and only over the Ocean area. Accordingly, I masked the land area and plotted the data for each layer. <br></span></div><div dir="ltr"><span>It is plotted fine, but when I try to take the areal average of a particular layer and particular sector over the ocean, I run into problems. My output shows me <br></span></div><div dir="ltr"><span>1999 7, 0.5000,-----------------,-9999.90039062 <br> (DATE MONTH , 1.0000, ,<br> , 1.5000, ,<br> , 2.0000, ,<br> , 2.5000, ,<br> , 3.0000, ,<br> , 3.5000, ,<br> , 4.0000, ,<br> , 4.5000, ,<br> , 5.0000, ,<br> , 5.5000, ,<br> , 6.0000, </span></div><div><div dir="ltr"> (layers) (areal Average Value)</div><div dir="ltr"><br></div><div dir="ltr">I am unable to proceed further>>>>. Attaching my script and the figures and output file. </div><div dir="ltr">when I run this, I get areal_ave = -9999. for all the layers...</div><div dir="ltr">umber Of Attributes: 10<br> longitude : degrees_east<br> latitude : degrees_north<br> comments : Unknown1 variable comment<br> long_name : latent heat (C/hr)<br> units : <br> grid_name : grid-1<br> grid_type : linear<br> time_statistic : instantaneous<br> missing_value : -9999.9<br> _FillValue : -9999.9<br><br><br>Variable: date_final<br>Type: string<br>Total Size: 8 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates: <br>(0) 1999 7<br><br><br>Variable: area_avg<br>Type: float<br>Total Size: 4 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates: <br>Number Of Attributes: 1<br> _FillValue : -9999.9<br>(0) -9999.9<span class="HOEnZb"><font color="#888888"><br><br></font></span></div></div><span class="HOEnZb"><font color="#888888"><div>Geeta.</div></font></span></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>