<div dir="ltr"><div>A variation<br></div><div><br></div><div><a href="http://www.ncl.ucar.edu/Applications/2Dvertcoord.shtml"><b>http://www.ncl.ucar.edu/Applications/2Dvertcoord.shtml</b></a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 1, 2021 at 11:07 AM Rashed Mahmood via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</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>Ahh...Ok now I see one example on NCL page: <a href="https://www.ncl.ucar.edu/Applications/overlay.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/overlay.shtml</a></div><div>See Overlay_14.ncl.<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 1, 2021 at 7:03 PM Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com" target="_blank">rashidcomsis@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">Can you share your latest script and the data files? You can send directly to me, if you want.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 1, 2021 at 6:26 PM Prasad Kunjeer <<a href="mailto:pskunjeer@gmail.com" target="_blank">pskunjeer@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="auto">My utimate aim is to create plot as attached in enclosed image</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 1 Mar, 2021, 9:29 PM Rashed Mahmood, <<a href="mailto:rashidcomsis@gmail.com" rel="noreferrer" target="_blank">rashidcomsis@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">My mistake. I did not realize that you are trying to overlay an xy plot on top of a contour plot which would not work as it is. Can you be more specific what you are trying to do (what is the end goal here)?<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 1, 2021 at 11:24 AM Prasad Kunjeer <<a href="mailto:pskunjeer@gmail.com" rel="noreferrer noreferrer" target="_blank">pskunjeer@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>Sir,</div><div><br></div><div>I have tried the solution, now there is single figure but precipitation is not plotted over it.</div><div><br></div><div><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><span style="font-size:10pt;font-family:"Mangal","serif"" lang="HI">Regards<br>प्रसाद कुंजीर </span><span style="font-size:9pt;font-family:"Arial","sans-serif"">/ </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Prasad Kunjeer<br>
</span><span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">वैज्ञानिक -</span></span><span lang="hi"> सी</span><span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI"> / </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Scientist - C<br>
</span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">केन्द्रीय जल और विद्युत
अनुसंधान</span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"" lang="HI"> </span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">शाला / </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Central Water and Power Research Station<br>
</span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">पुणे </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">411 024/ Pune 411 024</span></span></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 1, 2021 at 12:38 AM Rashed Mahmood <<a href="mailto:rashidcomsis@gmail.com" rel="noreferrer noreferrer" target="_blank">rashidcomsis@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>You need set the following in the respective resources list:</div><div><br></div><div> res@gsnDraw              = False             ; do not draw the plot<br>res@gsnFrame             = False             ; do not advance the frame</div><div><br></div><div><br></div><div><div>res_z@gsnDraw              = False             ; do not draw the plot<br>res_z@gsnFrame             = False             ; do not advance the frame</div></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Feb 28, 2021 at 7:47 PM Prasad Kunjeer via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" rel="noreferrer noreferrer" target="_blank">ncl-talk@mailman.ucar.edu</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">I am writing a script to overlay U wind and precipitation from the wrf run. My script is as below. <div><br></div><div>;----------------------------------------------------------------------<br>; wrf_interp_3.ncl<br>;----------------------------------------------------------------------<br>; Concepts illustrated:<br>;   - Interpolating a vertical cross-section from a 3D WRF-ARW field.<br>;----------------------------------------------------------------------<br>; wrf_user_vert_cross and wrf_user_interp_level replace the <br>; deprecated wrf_user_intrp3d function. <br>;<br>; NCL V6.6.0 or higher is required to run this example.<br>;----------------------------------------------------------------------<br><br>begin<br>  dir = "/media/prasad/C_Drive/2016/res/regrid_exp/"<br>  a   = addfile(dir+"<a href="http://wrfout_d03_2016-06-30_17_00_00.nc" rel="noreferrer noreferrer" target="_blank">wrfout_d03_2016-06-30_17_00_00.nc</a>","r")<br><br>  z    = wrf_user_getvar(a, "z",0)          ; grid point height<br>  qv   = wrf_user_getvar(a, "QVAPOR",0)     ; cloud field<br>  lats = wrf_user_getvar(a, "lat",0)<br>  lons = wrf_user_getvar(a, "lon",0)<br>  pr   = wrf_user_getvar(a, "pressure",0)<br>  ua   = wrf_user_getvar(a, "ua",0)<br>  va   = wrf_user_getvar(a, "va",0)<br>  wa   = wrf_user_getvar(a, "wa",0)<br><br>  rain_exp = wrf_user_getvar(a,"RAINNC",0) / 25.4<br>  rain_con = wrf_user_getvar(a,"RAINC", 0) / 25.4<br>  rain_tot = (rain_exp + rain_con)<br>  copy_VarMeta(rain_exp,rain_tot)<br><br>  printVarSummary(rain_tot) <br>  printVarSummary(rain_exp)         <br>  printMinMax(rain_tot,0)          <br> <br>  xlon = wrf_user_getvar(a, "XLONG",0)<br><br>  start_lat       = 29.75<br>  end_lat         = 29.75<br>  start_lon       = 77<br>  end_lon         = 83<br>  opt             = True<br>  opt@latlon      = True<br>  opt@linecoords  = True<br>  opt@file_handle = a<br>  qv_latlon       = wrf_user_vert_cross(ua,pr,(/start_lon,start_lat,end_lon,end_lat/),opt)<br> <br>  printVarSummary(ua)          ; [bottom_top | 31] x [south_north | 546] x [west_east | 480]<br>  printVarSummary(qv_latlon)   ; [vertical | 100] x [cross_line_idx | 171]<br>  printMinMax(ua,0)          ; [bottom_top | 31] x [south_north | 546] x [west_east | 480]<br>  printMinMax(qv_latlon,0)   ; [vertical | 100] x [cross_line_idx | 171]<br><br><br>;---Interpolate the 2D Rain field to the lat/lon line<br>  t2_line  = wrf_user_interp_line(rain_tot,(/start_lon,start_lat,end_lon,end_lat/),opt)<br>  printVarSummary(t2_line)   <br>  long_line  = wrf_user_interp_line(xlon,(/start_lon,start_lat,end_lon,end_lat/),opt)<br><br>  printMinMax(qv_latlon,0)<br>  printMinMax(t2_line,0)<br>  printMinMax(long_line,0)<br><br>  wks = gsn_open_wks("png","wrf_interp_u_wind")<br><br>  res                      = True<br>  res@gsnMaximize          = True     ; maximize plot in frame<br>  res@cnFillOn             = True     ; turn on contour fill<br>  res@cnFillOpacityF       = 0.4<br>  res@cnLinesOn            = False   ; turn off contour lines<br>  res@cnLineLabelsOn       = False    ; turn off line labels<br>  res@lbOrientation        = "Vertical"<br>  res@lbLabelFontHeightF   = 0.01<br>  res@tiMainString         = "Cross section from ("+start_lat+","+start_lon+ \<br>                                           ") to ("+end_lat + ","+end_lon+")"<br>  res@gsnStringFontHeightF = 0.015<br>  res@tmXBLabelFontHeightF = 0.01<br>  res@tmXBLabelAngleF      = 45.<br>  res@trYReverse = True<br>  ;res@gsnDraw              = False             ; do not draw the plot<br>  ;res@gsnFrame             = False             ; do not advance the frame<br><br>;--Explicitly set lat/lon labels for X axis<br>  xvalues = ispan(0,dimsizes(qv_latlon(0,:))-1,1)<br>  ll_step = 15                                       ; step interval for tickmarks<br>  res@tmXBMode          = "Explicit"<br>  res@tmXBValues        = xvalues(::ll_step)<br>  res@tmXBLabels        = sprintf("%6.2f",qv_latlon@lats(::ll_step)) + "~S~o~N~N~C~" + \<br>                          sprintf("%6.2f",qv_latlon@lons(::ll_step)) + "~S~o~N~E"<br><br>  plot = gsn_csm_contour(wks,qv_latlon,res)<br><br>  res_z                      = res<br>  res_z@cnFillOn             = False    ; turn on contour fill<br>  res_z@cnLinesOn            = True   ; turn off contour lines<br>  res_z@cnLineLabelsOn       = False    ; turn off line labels<br>  res_z@trYReverse           = False<br>  ;res_z@xyLineColors         = "blue" <br>  res_z@xyLineThicknesses   = (/3.0/)              ; make line thicker<br>  res_z@xyLineColors        = (/"green"/)          ; change line color<br><br>;--Explicitly set lat/lon labels for X axis<br>  xvalues = ispan(0,dimsizes(long_line(:))-1,1)<br>  ll_step = 15                                       ; step interval for tickmarks<br>  res_z@tmXBMode          = "Explicit"<br>  res_z@tmXBValues        = xvalues(::ll_step)<br>  res_z@tmXBLabels        = sprintf("%6.2f",qv_latlon@lats(::ll_step)) + "~S~o~N~N~C~" + \<br>                          sprintf("%6.2f",qv_latlon@lons(::ll_step)) + "~S~o~N~E"<br><br>  plot_ov = gsn_csm_xy(wks,long_line,t2_line,res_z)<br><br><br>  overlay(plot,plot_ov)                       ; overlay the U-wind plot on the temperature plot<br>  draw(plot)                                  ; draw the temperature plot (with the rh plot overlaid)<br>  frame(wks)                                  ; advance the frame<br><br>end<br></div><div><br></div><div><br></div><div>I am getting both precipitation and U wind plot but both on separate pages. I want to overlay them. The help is really appreciated</div><div><br></div><div><br clear="all"><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><span style="font-size:10pt;font-family:"Mangal","serif"" lang="HI">Regards<br>प्रसाद कुंजीर </span><span style="font-size:9pt;font-family:"Arial","sans-serif"">/ </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Prasad Kunjeer<br>
</span><span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">वैज्ञानिक -</span></span><span lang="hi"> सी</span><span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI"> / </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Scientist - C<br>
</span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">केन्द्रीय जल और विद्युत
अनुसंधान</span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"" lang="HI"> </span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">शाला / </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">Central Water and Power Research Station<br>
</span><span style="font-size:10.5pt;font-family:"Mangal","serif"" lang="HI">पुणे </span><span style="font-size:10.5pt;font-family:"Arial","sans-serif"">411 024/ Pune 411 024</span></span></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" rel="noreferrer noreferrer" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer noreferrer noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>