<div dir="ltr"><div><div><div><div>1) You're plotting meridional temperature advection rather than just the meridional gradient<br></div>2) units, units, units<br></div><div>3) The grid spacing of your data could impact the magnitudes of the final product. Models run on finer grids can resolve higher winds and temperature gradients. Taking the average might not fully cancel that out.<br><br></div><div>Jeff Duda<br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 11, 2017 at 1:19 AM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br clear="all"><div><div class="m_3069832306561089600gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear All,</div><div><br></div><div>I am trying to compute the meridional gradient of of potential temperature (925hPa). </div><div>Then create a latitude-time plot.</div><div>I am getting very small magnitudes of the gradient which is not consistent with the literature that I am following (<a href="http://wwwoa.ees.hokudai.ac.jp/people/yamazaki/D/Dthesis_yOkada.pdf" target="_blank">http://wwwoa.ees.hokudai.ac.<wbr>jp/people/yamazaki/D/Dthesis_<wbr>yOkada.pdf</a>).</div><div><br></div><div>Here's the script that I am using:</div><div><div>'sdfopen /mnt/dat/wind/vwind/10m/<a href="http://pentad_ensmean_vwnd_surf_1981-2010.nc" target="_blank">pentad<wbr>_ensmean_vwnd_surf_1981-2010.<wbr>nc</a>'</div><div><br></div><div>'set display color white'</div><div>'c'</div><div>'set gxout shaded'</div><div>'set mpdset hires'</div><div>'set grid off'</div><div>'set grads off'</div><div>'set parea 0.5 8 2 10.5'</div><div>'set font 1'</div><div>'set lat -90 90'</div><div>'set lon 0 360'</div><div>'set dfile 1'</div><div>'set t 1 73'</div><div>'define t=pottmp'</div><div>'set dfile 2'</div><div>'set t 1 73'</div><div>'define u=uwnd'</div><div>'set dfile 3'</div></div><div><div>'set t 1 73'</div><div>'define v=vwnd'</div><div><br></div><div><b>'define dty = cdiff(t,y)'</b></div><div><b>'define Re = 6.37e6'</b></div><div><b>'define dy = Re*cdiff(lat,y)*(3.1416/180)'</b></div><div><b>'define mergrad = -1*v*dty/dy'</b></div><div><br></div><div>'set t 1 73'</div><div>'set lat 15 45'</div><div>'set x 1'</div><div>'set xyrev on'</div><div>'set xaxis 1 73'</div><div>'define a = ave(mergrad,lon=110,lon=120)'</div><div>'d a'</div><div>'cbarn'</div><div>'printim test.png'</div><div>'quit'</div></div><div><br></div><div>Attached is the output image. In the literature,the unit is between -3 and 3, while in my case its x10-5.</div><div><br></div><div><br></div><div>Any suggestions on how I can do this correctly?</div><div><br></div><div><br></div><div>Lyndz</div><div><br></div></div></div></div></div></div></div></div></div>
</div>
<br>______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Jeff Duda<br>Post-doctoral research fellow<br>University of Oklahoma School of Meteorology<br></div></div></div></div></div>
</div>