<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>b=addfile("../../wrfout_d05_2016-03-31_00:00:00","r")</div><div>a=addfile("../../../wrfout_d05_2016-03-31_00:00:00","r")</div><div>t_Ind=a->TA_MEAN<br></div><div>ta=t_Ind(32:728,0,48:52,73:78)</div><div><br></div><div>t_No_Ind=b->TA_MEAN</div><div>tb=t_No_Ind(32:728,0,48:52,73:78)</div><div>ta_diff=ta-tb<br></div><div>copy_VarCoords(ta,ta_diff)</div><div><br></div><div>This gives out </div><div><br></div><div><div>Variable: ta_diff</div><div>Type: float</div><div>Total Size: 83640 bytes</div><div>            20910 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:   [Time | 697] x [south_north | 5] x [west_east | 6]</div></div><div>I want to first average over time so that I have (24x29) values. </div><div><br></div><div>Then I want to calculate percentiles P25, P50, P75, P90 and max values of ta_diff. </div><div><br></div><div>Can someone help please </div><div><br></div><div><br></div></div></div></div></div></div></div>