Dear Mr. Ian,<div>&nbsp; &nbsp; &nbsp; Thanks for your careful reply, again! I think it will be very helpful for me to calculate the zonal/meridional advection flux. Recently, i have another technique trouble. when i decomposed advection term into the contributions by method of perturbation, i&nbsp;<font size="2"><font face="Cambria Math">have no efficient way to&nbsp;distinguish the contributions from salinity anomaly or current anomaly, because&nbsp;</font></font>the total advection flux is represented by "<font size="2"><span style="font-style: italic; font-family: &quot;Cambria Math&quot;;">h ∇∙(</span><span style="font-style: italic; font-family: &quot;Cambria Math&quot;;">S</span><span style="font-style: italic; font-family: &quot;Cambria Math&quot;; font-weight: bold;">u</span><span style="font-family: &quot;Cambria Math&quot;;"><i>) " </i>at last, and &nbsp;</span></font><font size="2"><span style="font-family: &quot;Cambria Math&quot;;">the model output&nbsp;</span></font><font face="Cambria Math" size="2">such as&nbsp;</font>ADVx_SLT can not represent the current multiplied by the gradient of salinity. So i want to know about the detailed calculations and even matlab code if it is convenient. Could you give me some advices?</div><div>&nbsp; &nbsp; Thank you for your kindness again!</div><div><br></div><div><div>Best regards,</div><div>Michael</div></div>