[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 11:03:36 MST 2021
Can you share your latest script and the data files? You can send directly
to me, if you want.
On Mon, Mar 1, 2021 at 6:26 PM Prasad Kunjeer <pskunjeer at gmail.com> wrote:
> 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/50d8b784/attachment.html>
More information about the ncl-talk
mailing list