[ncl-talk] Overlay plot of vertical cross section of WRF derived U and precipitation using NCL
Prasad Kunjeer
pskunjeer at gmail.com
Mon Mar 1 10:26:02 MST 2021
My utimate aim is to create plot as attached in enclosed image
On Mon, 1 Mar, 2021, 9:29 PM Rashed Mahmood, <rashidcomsis at gmail.com> wrote:
> 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/b0b801ae/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: overlay_plot.jpg
Type: image/jpeg
Size: 910530 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210301/b0b801ae/attachment-0001.jpg>
More information about the ncl-talk
mailing list