<div>I am trying to do an advection scheme for olr values and am not sure how to go about doing it.  I have made a nc file or average wind patterns for certain climatic patterns in the Pacific Ocean called <a href="http://uninan.nc">uninan.nc</a>.  this file contains zonal winds uninan.</div>

<div> </div>
<div>here is the script for modeling OLR advection  the start dates on both files match.</div>
<div> </div>
<div>&#39;sdfopen <a href="http://uninan.nc">uninan.nc</a>&#39;<br>&#39;sdfopen <a href="http://olr.day.mean.nc">olr.day.mean.nc</a>&#39;<br>&#39;set gxout shaded&#39;<br>&#39;set t 1&#39;<br>&#39;define temp=olr.2&#39;<br>
count=1<br>while (count &lt; 10)<br>    &#39;set clevs 140 160 180&#39;<br>    &#39;set ccols 14 4 11&#39;<br>    &#39;dtx = cdiff(temp,x)&#39;<br>    &#39;dty = cdiff(temp,y)&#39;<br>    &#39;dx=cdiff(temp,x)*3.14/180&#39;<br>
    &#39;dy=cdiff(lat,y)&#39;<br>    &#39;define temp=temp-((uninan*dtx/(cos(lat*3.14/180)*dx)))/2000&#39;<br>    &#39;d temp&#39;<br>count=count + 1<br>endwhile<br></div>
<div>based on the zonal wind values olr anomalies should advect across the Pacific Ocean in 4 days count=4 but they are not doing so or OLR values are off by quite a bit.</div>
<div>Do you have any suggestions?</div>
<div> </div>
<div>Thanks,</div>
<div> </div>
<div>Jason</div>
<div> </div>