<div dir="ltr"><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><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>Dear NCL experts,</div><div><br></div><div>I would like to compute the potential vorticity tendency (dPV/dt) from a wrf file. </div><div>However, there is no time dimension in the variables (although the "time" variable can be obtained using the wrf_user_getvar)</div><div><br></div><div>Here's my script so far:</div><div><br></div><div>filename = "./wrfout_d01_2011-12-12_12_00_00"<br>a     = addfile(filename,"r")<br>times = wrf_user_getvar(a,"times",-1)                   ; get all times in the file<br>it = 89                                                   ; first time step<br>theta = wrf_user_getvar(a,"theta",it)<br>pv    = wrf_user_getvar(a,"pvo",it)                     ;potential vorticity<br></div><div><br></div><div>I would like to compute dPV/dt using centered finite difference (center_finite_diff_n).</div><div><br></div><div>Any ideas/solutions/examples on how can I do this in NCL?</div><div><br></div><div>I'll appreciate any help on this matter.</div><div><br></div><div>-Lyndz</div><div></div></div></div></div></div></div></div></div></div></div></div>