<div>Dear all,</div><div>I'm dealing with the JRA four times per day grib data. And I want to average the four times data to get a daily data.</div><div>My script is listed below</div><div><div><b>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</b></div><div><b>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</b></div><div><b>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</b></div><div><b>fa=addfile("t.grib","r")</b></div><div><span style="line-height: 1.5;"><b>T=fa-&gt;TMP_GDSO_ISBL ;T(123,37,73,144)&nbsp;</b></span></div><div><b>Tpj=new((/31,37,73,144/),float)</b></div><div><b>do j=0,30&nbsp;</b></div><div><b>&nbsp; &nbsp; do i=0,123&nbsp;</b></div><div><b>&nbsp; &nbsp; &nbsp; &nbsp; Tpj(j,:,:,:)=(T(i,:,:,:)+T(i+1,:,:,:)+T(i+2,:,:,:)+T(i+3,:,:,:))/4</b></div><div><b>&nbsp; &nbsp; i=i+4</b></div><div><b>&nbsp; &nbsp;j=j+1</b></div><div><b>&nbsp; &nbsp;enddo;</b></div><div><b>enddo;</b></div></div><div>&nbsp;the error goes synatx error and error in statement. How should I do to get the daily average?</div><div>Any help would be appreciated!</div><div>Sincerely,</div><div>Yan</div>