<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear Everyone <div><br></div><div>I guess the differences in the plots were because elevation from wrfout is extracted by the following equation</div><div><br></div><div><b style="color:rgb(32,33,34);font-family:sans-serif;font-size:14px">ELEVATION_THETA(i,j,k) =0.5*(PHB(i,j,k) + PH(i,j,k) + PHB(i,j,k+1) + PH(i,j,k+1))/9.81</b><br></div><div><br></div><div><br></div><div>Further information is provided in the following threads.</div><div><br></div><div><a href="https://forum.mmm.ucar.edu/phpBB3/viewtopic.php?f=39&t=9831">https://forum.mmm.ucar.edu/phpBB3/viewtopic.php?f=39&t=9831</a><br></div><div><a href="https://wiki.openwfm.org/wiki/How_to_interpret_WRF_variables">https://wiki.openwfm.org/wiki/How_to_interpret_WRF_variables</a><br></div><div><br></div><div><br></div><div>Let me know if this is the correct reason for the differences I observe in the plots.</div><div><br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 24, 2022 at 6:59 PM ali mughal <<a href="mailto:mughalali655@gmail.com">mughalali655@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear NCL users <div><br></div><div>I reported this issue a few days ago.</div><div><br></div><div>if I subtract the terrain height from the code by adding the following step</div><div><br></div><div><br></div><div><br></div><div>ter=  wrf_user_getvar(a, "ter",-1)<br></div><div><div>nheight = conform(z,ter,(/0,2,3/))</div></div><div>z_agl = z-nheight<br></div><div><br></div><div><div>do level=0, nlevels-1</div><div><br></div><div> height=levels(level)</div><div><br></div><div><br></div><div>u_plane = wrf_user_intrp3d( u,z_agl,"h",height,0.,False)</div><div>v_plane = wrf_user_intrp3d( v,z_agl,"h",height,0.,False)</div><div>ta_mean_plane = wrf_user_intrp3d(ta_mean_a ,z_agl,"h",height,0.,False)</div></div><div><br></div><div>I get the following plot</div><div><br></div><div><img src="cid:ii_kyski7bj0" alt="image.png" width="412" height="254"><br></div><div><br></div><div><br></div><div>As you can see there is change in the field compared to the plots I shared earlier but white spaces are removed.</div><div><br></div><div>Now I am not able to understand what is causing the white spaces earlier and why they are removed when I subtract the terrain height. </div><div><br></div><div>In the following example, a similar procedure is followed, however, the only difference between my code and the one in the example is the wrf file used in the example is for multiple time steps while I am using a file which is generated with cdo as an hourly average. </div><div><br></div><div><a href="https://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_Height.ncl" target="_blank">https://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_Height.ncl</a><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 21, 2022 at 6:58 PM ali mughal <<a href="mailto:mughalali655@gmail.com" target="_blank">mughalali655@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Dear NCL community <div><br></div><div>I attach the output from two ncl scripts. The first one (Seabreeze_hourly_Ind_avg.png)  interpolates temperature and wind field to the first pressure level from WRF out file. If I convert the pressure level to the height in m it comes out to be 10 m. The 2nd output I obtain by using the following method instead of using the first pressure level</div><div><br></div><div><div>levels = (/30.,35.,40.,45.,50.,55.,60.,65.,70.,75.,80.,85.,90./)</div><div>nlevels=dimsizes(levels)</div><div><br></div><div>do level=0, nlevels-1</div><div><br></div><div>  height=levels(level)</div><div>  u_plane = wrf_user_intrp3d( u,z,"h",height,0.,False)</div><div>  v_plane = wrf_user_intrp3d( v,z,"h",height,0.,False)</div><div>  ta_mean_plane = wrf_user_intrp3d(ta_mean ,z,"h",height,0.,False)</div><div>  U_1=u_plane(0,:,:)<br></div><div>  V_1=v_plane(0,:,:)<br></div><div> ta_mean_1=ta_mean_plane(0,:,:)</div><div> tc2=ta_mean_1-273<br></div></div><div><br></div><div>As you can see that the first level is 30m, considerably higher than the first pressure level but in the first plot, there are no white spaces which mean every thing is plotted above the terrain while in the second plot there are white spaces meaning these are points below the terrain. </div><div>Can someone kindly let me know how I can remove this difference ?</div><div><br></div><div>Regards </div><div><br></div><div>Ali</div><div><br></div></div></div>
</blockquote></div>
</blockquote></div>