Thank you<span></span><br><br>On Friday, September 12, 2014, Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">a=addfile("~/Documents/Snow_<div>cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc" target="_blank">nhsce_v01r01_19661004_20140602.nc</a>","r")<br>sce=a->snow_cover_extent<br><br>;print(dimsizes(sce))<br><br>sce_area=a->area<br><br>B = sce_area<br><br> do i = 0,2486<br> B(i,:,:)=(/ sce(i,:,:)*sce_area(:,:) /)<br> end do<br></div><div>B@long_name = "sce*sce_area"<br></div><div><br>xMonthAvg = calculate_monthly_values(B, "avg", 0, False)<br></div>
</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 12, 2014 at 1:01 PM, Ipshita Majhi <span dir="ltr"><<a href="javascript:_e(%7B%7D,'cvml','ipmajhi@alaska.edu');" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi,<br></div>I am trying to create a new matrix where snow cover extent (sce) which are in bytes is multiplied by area to get the real km2 values. In the process I am loosing the time information. How to save that so that I can create monthly avg.<br></div>Best Regards<br>Ipshita<br><div><div><div><br></div><div><br>; This program is to calculate the area average for snow for all the data period<br>;*******************************************<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>;*******************************************<br><br>;To create a new matrix with area values multiplied to the binary values so that we can get <br><br>a=addfile("~/Documents/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc" target="_blank">nhsce_v01r01_19661004_20140602.nc</a>","r")<br>sce=a->snow_cover_extent<br><br>;print(dimsizes(sce))<br><br>sce_area=a->area<br><br>B = new((/2487,88,88/),float)<br><br> do i = 0,2486<br> B(i,:,:)=sce(i,:,:)*sce_area(:,:)<br> end do<br><br>xMonthAvg = calculate_monthly_values(B, "avg", 0, False)<br></div></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote>