<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><div>Dear NCL users,</div><div>I have hourly wrf arw output and plotted t<span style="font-size: 12pt;">he hovmoller diagramme &nbsp;along the cross section of two points &nbsp;for july as&nbsp;</span></div><div><ul><li>monthly average of each hour (average value of 0 hour, 1 hour ........23 hour of each day) precipitation for whole month ( 24 values ).</li><li>ploted hovmoller for &nbsp;hourly time &nbsp;vs cross section grid of two points.</li><li>i have used reshape function however i am confused whether its producing right output or not.(use of reshape function for my purpose). <b>(array as 24, number of days, ncross_section_points and averaging for days&nbsp;</b></li><li>furthermore, I would like to convert the UTC time to local time (5:45 hour ahead of UTC) in my plot but could not get idea to do so.</li></ul></div><div>I look forward for getting feedback from you.</div><div><b>my script</b></div><div><b><br></b></div><div>; This script is used to plot&nbsp;</div><div>; get each hour average value for 24 hours</div><div>; get data along two point section</div><div>; plot hovmoller for time vs grid along section&nbsp;</div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div>; The WRF ARW input file. &nbsp;</div><div>; This needs to have a ".nc" appended, so just do it.</div><div>&nbsp;DATADir = "/geodata/data/private/fgvy024/Build_WRF/WRFV6/run/"</div><div>&nbsp;FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d03_2015-07*") ; contain hourly time stamp</div><div>&nbsp;numFILES = dimsizes(FILES)</div><div>&nbsp;a = addfiles(FILES+".nc","r")</div><div>&nbsp;</div><div>&nbsp;type = "png" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; map type</div><div>&nbsp;wks = gsn_open_wks(type,"july_hovmoller_hourly_ppt") &nbsp; &nbsp; &nbsp; &nbsp; ; Create a plot workstation</div><div>&nbsp;rain_exp = wrf_user_getvar(a,"RAINNC",-1)</div><div>&nbsp;rain_con = wrf_user_getvar(a,"RAINC",-1)</div><div>&nbsp;rain = rain_exp + rain_con</div><div>&nbsp;Remp = dimsizes(rain)</div><div>&nbsp;nt = Remp(0)</div><div>&nbsp;nnlat = Remp(1)</div><div>&nbsp;nlon = Remp(2)</div><div>; get hourly rainfall</div><div>&nbsp; Ree = rain(1:nt-1,:,:) - rain(0:nt-2,:,:) ; to get hourly precipitation</div><div>&nbsp;; get new variable with full time step to make it multiplable by 24</div><div>&nbsp; Re = new((/nt, nnlat, nlon/),float)</div><div>&nbsp; Re(1:nt-1,:,:) = Ree(:,:, :)&nbsp;</div><div>&nbsp; Re(0, :, :) = 0 ; first array will have zero value</div><div>&nbsp; dims = dimsizes(Re)</div><div>&nbsp; pA_lon = 86.2233 ; longitude of point A</div><div>&nbsp; pA_lat = 27.88163333 ; lat of point A</div><div>&nbsp; pB_lon = 86.541725 ; similar for point B</div><div>&nbsp; pB_lat = 27.89950278</div><div>&nbsp; &nbsp; pointA = wrf_user_ll_to_ij(a, pA_lon, pA_lat, True) ; get array point</div><div>&nbsp; &nbsp; pointB = wrf_user_ll_to_ij(a, pB_lon, pB_lat, True)</div><div>&nbsp; &nbsp; xA = pointA(0) - 1</div><div>&nbsp; &nbsp; yA = pointA(1) - 1 &nbsp;</div><div>&nbsp; &nbsp; xB = pointB(0) - 1</div><div>&nbsp; &nbsp; yB = pointB(1) - 1&nbsp;</div><div>&nbsp; opts = True&nbsp;</div><div>&nbsp; plane = new(4,float) ; for section along two points</div><div>&nbsp; plane = (/ &nbsp;yA,xA, yB,xB &nbsp;/)</div><div><br></div><div>; interpolate value along the section of A and B</div><div>&nbsp;R_plane = wrf_user_intrp2d(Re,plane,0.,opts)</div><div>&nbsp;Re_dims = dimsizes(R_plane)</div><div>&nbsp;ntim &nbsp; &nbsp; = Re_dims(0)</div><div>&nbsp;ncross &nbsp; &nbsp; = Re_dims(1)</div><div>&nbsp;ndays &nbsp; &nbsp; &nbsp; = ntim/24</div><div>&nbsp;Re3d &nbsp; &nbsp; &nbsp; = reshape(R_plane,(/24,ndays, ncross/)) ; &nbsp;hours, days, crosssection</div><div>&nbsp;copy_VarAtts(R_plane,Re3d)</div><div>&nbsp;Re3d_hr_avg = dim_avg_n_Wrap(Re3d,1) ; hourly average (monthly, 0 hours average, 1 hour average, etc)</div><div>&nbsp; res &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; plot mods desired</div><div>&nbsp; res@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; turn on color fill</div><div>&nbsp; res@cnFillPalette &nbsp; &nbsp; &nbsp; &nbsp;= "BlWhRe" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; set color map</div><div>&nbsp; res@tiMainString &nbsp; &nbsp; &nbsp; &nbsp; = "Hourly Precipitation" &nbsp; ; title</div><div>&nbsp; res@cnLevelSelectionMode = "ManualLevels" &nbsp; &nbsp; ; manual contour levels</div><div>&nbsp; res@cnMinLevelValF &nbsp; &nbsp; &nbsp; = 0.<span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp; &nbsp; &nbsp; &nbsp;; min level</div><div>&nbsp; res@cnMaxLevelValF &nbsp; &nbsp; &nbsp; = &nbsp;2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; max level</div><div>&nbsp; res@cnLevelSpacingF &nbsp; &nbsp; &nbsp;= 0.1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; contour level spacing</div><div>&nbsp; plot = gsn_csm_hov(wks, Re3d_hr_avg(:,:), res)&nbsp;</div><div>end</div><div><br></div><div><i style="font-size: 12pt;">Regards</i></div><i>Ram</i><style><!--
.ExternalClass .ecxhmmessage P {
padding:0px;
}

.ExternalClass body.ecxhmmessage {
font-size:12pt;
font-family:Calibri;
}

--></style>                                               </div></body>
</html>