<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 &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<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(&quot;~/Documents/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>&quot;,&quot;r&quot;)<br>sce=a-&gt;snow_cover_extent<br><br>;print(dimsizes(sce))<br><br>sce_area=a-&gt;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, &quot;avg&quot;, 0, False)<br></div></div></div></div>