Re: vertical profile

Sima,

I believe you should contact wrfhelp@ucar.edu about this. I've CC'ed them.

The wrf_user_intrp3d routine doesn't return coordinates that would give you the values for sfXArray and sfYArray, but the WRF folks might be able to help you calculate them.

--Mary

On Feb 24, 2011, at 3:46 PM, sima sima wrote:

> Hi dear,
> I try to plot a vertical profile in an special point.
> I use gsn_contour function
> but its need too set sfXArrey and sfYArrey resources
> I don't know how set them
> i try to plot th on vertical levels in point(44,55)
>
> help me
>
> unfortunately I cant attached the script file because Im from Iran, then i copy it here:
>
> Thanks
>
> Sima.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
> a = addfile("/media/disk/wrfout/wrfout_d03_2009-05-05_10:00:00.nc","r")
>
> wks = gsn_open_wks("ps","vertical")
>
>
> ; Set some basic resources
> res = True
> res@MainTitle = "REAL-TIME WRF"
> res@Footer = False
>
> pltres = True
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> times = wrf_user_list_times(a) ; get times in the file
> ntimes = dimsizes(times) ; number of times in the file
> FirstTime = True
>
> mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file
> nd = dimsizes(mdims)
>
> ;---------------------------------------------------------------
>
> do it = 61, 61 ; TIME LOOP
>
> print("Working on time: " + times(it) )
> res@TimeLabel = times(it) ; Set Valid time to use on plots
>
> th = wrf_user_getvar(a,"th",it) ; T in C
> rh = wrf_user_getvar(a,"rh",it) ; relative humidity
> z = wrf_user_getvar(a, "z",it) ; grid point height
>
> ; max(z) = 3000
> if ( FirstTime ) then ; get height info for labels
> zmin = 0.
> zmax = max(z)
> nz = floattoint(zmax/2 + 1)
> FirstTime = False
> end if
>
> ;---------------------------------------------------------------
> angle = 0
> plane = new(2,float)
> plane = (/ 44,55/) ; pivot point is center of domain (x,y)
> opts = False ; start and end points not specified
>
> rh_plane = wrf_user_intrp3d(rh,z,"v",plane,angle,opts)
> th_plane = wrf_user_intrp3d(th,z,"v",plane,angle,opts)
> printVarSummary(th_plane)
> dim = dimsizes(th_plane) ; Find the data span - for use in labels
> zspan = dim(0)
>
> ; Options for XY Plots
> opts_xy = res
> opts_xy@tiYAxisString = "Height (m)"
> opts_xy@cnMissingValPerimOn = True
> opts_xy@cnMissingValFillColor = 0
> opts_xy@cnMissingValFillPattern = 11
> opts_xy@tmYLMode = "Explicit"
> opts_xy@tmYLValues = fspan(0,zspan,nz) ; Create tick marks
> opts_xy@tmYLLabels = sprintf("%.1f",fspan(zmin,zmax,nz)) ; Create labels
> opts_xy@tiXAxisFontHeightF = 0.020
> opts_xy@tiYAxisFontHeightF = 0.020
> opts_xy@tmXBMajorLengthF = 0.02
> opts_xy@tmYLMajorLengthF = 0.02
> opts_xy@tmYLLabelFontHeightF = 0.015
> opts_xy@PlotOrientation = th_plane@Orientation
>
>
> ; Plotting options for RH
> opts_rh = opts_xy
>
> ; Plotting options for Temperature
> opts_th = opts_xy
>
> ; Get the contour info for the rh and temp
> contour_th = wrf_contour(a,wks,th_plane,opts_th)
> contour_rh = wrf_contour(a,wks,rh_plane,opts_rh)
> printVarSummary(contour_th)
>
> ; MAKE PLOTS
> ; plot = wrf_overlays(a,wks,(/contour_rh,contour_tc/),pltres)
> plot = gsn_csm_y(wks,th_plane(:,44),opts_xy)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> end do ; END OF TIME LOOP
>
> end
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 25 08:37:37 2011