<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">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>