<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hello<div>Dear Users, </div><div>I want to calculate anomaly precipitation with  chirps monthly data (It is one nc file with 477 t size since  jan. 1981 to oct. 2020 ). At first I calculated Climatic average for each month (Jan , Feb, ....) with this script:</div><div>'sdfopen F:\chrips\<a href="http://repre.nc">repre.nc</a>'</div><div>'define janclim=ave(repre,t=1,t=477,12)'</div><div>'set sdfwrite F:\chrips\<a href="http://janclim.nc">janclim.nc</a>'</div><div>'sdfwrite janclim'</div><div>So I calculated 12 Climatic average nc files (<a href="http://janclim.nc">janclim.nc</a> , Febclim.nc , .....). After that I wanted to calculate Anomaly precipitation each month.I mean, during this period 12 anomalies had to be calculated (Jananomaly.nc , Febanomaly.nc ,....) that t size in each file should be 40 . I tried the script below but the output files had 1 t size. </div><div>I would be very grateful if anyone could help me.</div><div><br></div><div><br></div><div>*****************************Script</div><div><div>'reinit'</div><div>'sdfopen F:\chrips\<a href="http://repre.nc">repre.nc</a>'</div><div>'sdfopen F:\chrips\meanclim\<a href="http://janclim.nc">janclim.nc</a>'</div><div>i = 1</div><div>while(i<=477)</div><div>'set t 'i</div><div>'define jananomaly=(repre.1-janclim.2)'</div><div>'set sdfwrite F:\chrips\<a href="http://jananomaly.nc">jananomaly.nc</a>'</div><div>'sdfwrite jananomaly'</div><div>i = i + 12</div><div>endwhile</div><div>'reinit'</div></div><div><br></div><div> </div><div dir="rtl">ه</div></div></div></div></div></div></div></div></div></div></div></div>