[ncl-talk] Overlay plot of vertical cross section of WRF derived U and precipitation using NCL
Rashed Mahmood
rashidcomsis at gmail.com
Mon Mar 1 08:59:06 MST 2021
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)?
On Mon, Mar 1, 2021 at 11:24 AM Prasad Kunjeer <pskunjeer at gmail.com> wrote:
> Sir,
>
> I have tried the solution, now there is single figure but precipitation is
> not plotted over it.
>
> Regards
> प्रसाद कुंजीर / Prasad Kunjeer
> वैज्ञानिक - सी / Scientist - C
> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power Research
> Station
> पुणे 411 024/ Pune 411 024
>
>
> On Mon, Mar 1, 2021 at 12:38 AM Rashed Mahmood <rashidcomsis at gmail.com>
> wrote:
>
>> You need set the following in the respective resources list:
>>
>> res at gsnDraw = False ; do not draw the plot
>> res at gsnFrame = False ; do not advance the frame
>>
>>
>> res_z at gsnDraw = False ; do not draw the plot
>> res_z at gsnFrame = False ; do not advance the frame
>>
>>
>> On Sun, Feb 28, 2021 at 7:47 PM Prasad Kunjeer via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> I am writing a script to overlay U wind and precipitation from the wrf
>>> run. My script is as below.
>>>
>>> ;----------------------------------------------------------------------
>>> ; wrf_interp_3.ncl
>>> ;----------------------------------------------------------------------
>>> ; Concepts illustrated:
>>> ; - Interpolating a vertical cross-section from a 3D WRF-ARW field.
>>> ;----------------------------------------------------------------------
>>> ; wrf_user_vert_cross and wrf_user_interp_level replace the
>>> ; deprecated wrf_user_intrp3d function.
>>> ;
>>> ; NCL V6.6.0 or higher is required to run this example.
>>> ;----------------------------------------------------------------------
>>>
>>> begin
>>> dir = "/media/prasad/C_Drive/2016/res/regrid_exp/"
>>> a = addfile(dir+"wrfout_d03_2016-06-30_17_00_00.nc","r")
>>>
>>> z = wrf_user_getvar(a, "z",0) ; grid point height
>>> qv = wrf_user_getvar(a, "QVAPOR",0) ; cloud field
>>> lats = wrf_user_getvar(a, "lat",0)
>>> lons = wrf_user_getvar(a, "lon",0)
>>> pr = wrf_user_getvar(a, "pressure",0)
>>> ua = wrf_user_getvar(a, "ua",0)
>>> va = wrf_user_getvar(a, "va",0)
>>> wa = wrf_user_getvar(a, "wa",0)
>>>
>>> rain_exp = wrf_user_getvar(a,"RAINNC",0) / 25.4
>>> rain_con = wrf_user_getvar(a,"RAINC", 0) / 25.4
>>> rain_tot = (rain_exp + rain_con)
>>> copy_VarMeta(rain_exp,rain_tot)
>>>
>>> printVarSummary(rain_tot)
>>> printVarSummary(rain_exp)
>>> printMinMax(rain_tot,0)
>>>
>>> xlon = wrf_user_getvar(a, "XLONG",0)
>>>
>>> start_lat = 29.75
>>> end_lat = 29.75
>>> start_lon = 77
>>> end_lon = 83
>>> opt = True
>>> opt at latlon = True
>>> opt at linecoords = True
>>> opt at file_handle = a
>>> qv_latlon =
>>> wrf_user_vert_cross(ua,pr,(/start_lon,start_lat,end_lon,end_lat/),opt)
>>>
>>> printVarSummary(ua) ; [bottom_top | 31] x [south_north | 546]
>>> x [west_east | 480]
>>> printVarSummary(qv_latlon) ; [vertical | 100] x [cross_line_idx |
>>> 171]
>>> printMinMax(ua,0) ; [bottom_top | 31] x [south_north | 546] x
>>> [west_east | 480]
>>> printMinMax(qv_latlon,0) ; [vertical | 100] x [cross_line_idx | 171]
>>>
>>>
>>> ;---Interpolate the 2D Rain field to the lat/lon line
>>> t2_line =
>>> wrf_user_interp_line(rain_tot,(/start_lon,start_lat,end_lon,end_lat/),opt)
>>> printVarSummary(t2_line)
>>> long_line =
>>> wrf_user_interp_line(xlon,(/start_lon,start_lat,end_lon,end_lat/),opt)
>>>
>>> printMinMax(qv_latlon,0)
>>> printMinMax(t2_line,0)
>>> printMinMax(long_line,0)
>>>
>>> wks = gsn_open_wks("png","wrf_interp_u_wind")
>>>
>>> res = True
>>> res at gsnMaximize = True ; maximize plot in frame
>>> res at cnFillOn = True ; turn on contour fill
>>> res at cnFillOpacityF = 0.4
>>> res at cnLinesOn = False ; turn off contour lines
>>> res at cnLineLabelsOn = False ; turn off line labels
>>> res at lbOrientation = "Vertical"
>>> res at lbLabelFontHeightF = 0.01
>>> res at tiMainString = "Cross section from
>>> ("+start_lat+","+start_lon+ \
>>> ") to ("+end_lat +
>>> ","+end_lon+")"
>>> res at gsnStringFontHeightF = 0.015
>>> res at tmXBLabelFontHeightF = 0.01
>>> res at tmXBLabelAngleF = 45.
>>> res at trYReverse = True
>>> ;res at gsnDraw = False ; do not draw the plot
>>> ;res at gsnFrame = False ; do not advance the
>>> frame
>>>
>>> ;--Explicitly set lat/lon labels for X axis
>>> xvalues = ispan(0,dimsizes(qv_latlon(0,:))-1,1)
>>> ll_step = 15 ; step interval for
>>> tickmarks
>>> res at tmXBMode = "Explicit"
>>> res at tmXBValues = xvalues(::ll_step)
>>> res at tmXBLabels = sprintf("%6.2f",qv_latlon at lats(::ll_step)) +
>>> "~S~o~N~N~C~" + \
>>> sprintf("%6.2f",qv_latlon at lons(::ll_step)) +
>>> "~S~o~N~E"
>>>
>>> plot = gsn_csm_contour(wks,qv_latlon,res)
>>>
>>> res_z = res
>>> res_z at cnFillOn = False ; turn on contour fill
>>> res_z at cnLinesOn = True ; turn off contour lines
>>> res_z at cnLineLabelsOn = False ; turn off line labels
>>> res_z at trYReverse = False
>>> ;res_z at xyLineColors = "blue"
>>> res_z at xyLineThicknesses = (/3.0/) ; make line thicker
>>> res_z at xyLineColors = (/"green"/) ; change line color
>>>
>>> ;--Explicitly set lat/lon labels for X axis
>>> xvalues = ispan(0,dimsizes(long_line(:))-1,1)
>>> ll_step = 15 ; step interval for
>>> tickmarks
>>> res_z at tmXBMode = "Explicit"
>>> res_z at tmXBValues = xvalues(::ll_step)
>>> res_z at tmXBLabels = sprintf("%6.2f",qv_latlon at lats(::ll_step))
>>> + "~S~o~N~N~C~" + \
>>> sprintf("%6.2f",qv_latlon at lons(::ll_step)) +
>>> "~S~o~N~E"
>>>
>>> plot_ov = gsn_csm_xy(wks,long_line,t2_line,res_z)
>>>
>>>
>>> overlay(plot,plot_ov) ; overlay the U-wind plot
>>> on the temperature plot
>>> draw(plot) ; draw the temperature
>>> plot (with the rh plot overlaid)
>>> frame(wks) ; advance the frame
>>>
>>> end
>>>
>>>
>>> I am getting both precipitation and U wind plot but both on
>>> separate pages. I want to overlay them. The help is really appreciated
>>>
>>>
>>> Regards
>>> प्रसाद कुंजीर / Prasad Kunjeer
>>> वैज्ञानिक - सी / Scientist - C
>>> केन्द्रीय जल और विद्युत अनुसंधान शाला / Central Water and Power
>>> Research Station
>>> पुणे 411 024/ Pune 411 024
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210301/0cc4b573/attachment.html>
More information about the ncl-talk
mailing list