Re: Ri contours are not appeared in the figure

In addition to the suggested print statements,
you were requested to make several graphical additions.
Specifically: "You can debug your code by looking at
the "contour_ri" plot on its own". Did you do what
was requested? What were the results?

On 2/16/11 2:07 AM, wei huang wrote:
> Hi,
>
> Thanks for your response, the data file is heavy sorry I could not send
> it please... The details of printVarSummary, printMinMax, and image are
> attached. Hope that will be helped to figure out the problem.
>
> thanks
>
> wei
> ****************
>
> Copyright (C) 1995-2010 - All Rights Reserved
>
> University Corporation for Atmospheric Research
>
> NCAR Command Language Version 5.2.1
>
> The use of this software is governed by a License Agreement.
>
> See http://www.ncl.ucar.edu/ for more details.
>
> Variable: T
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType :104
>
> MemoryOrder : XYZ
>
> description : perturbation potential temperature (theta-t0)
>
> units :K
>
> stagger :
>
> coordinates : XLONG XLAT
>
> Variable: u
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType :104
>
> MemoryOrder : XYZ
>
> description : x-wind component
>
> units :m s-1
>
> stagger :
>
> coordinates : XLONG XLAT
>
> Variable: v
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType :104
>
> MemoryOrder : XYZ
>
> description : y-wind component
>
> units :m s-1
>
> stagger :
>
> coordinates : XLONG XLAT
>
> Variable: z
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType :104
>
> MemoryOrder : XYZ
>
> description : Height
>
> units :m
>
> stagger :
>
> coordinates : XLONG XLAT
>
> (0)
>
> (0)perturbation potential temperature (theta-t0): min=275.254max=502.357
>
> (0)
>
> (0)x-wind component: min=-27.9777max=49.9323
>
> (0)
>
> (0)y-wind component: min=-18.5696max=60.0988
>
> (0)
>
> (0)Height: min=130.62max=20156.4
>
> Variable: U
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> (0)
>
> (0)min=0.00228019max=62.1447
>
> Variable: dTdz
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 1
>
> _FillValue :-999
>
> (0)
>
> (0)min=0.001max=0.0938801
>
> Variable: dUdz
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 1
>
> _FillValue :-999
>
> (0)
>
> (0)min=0.001max=0.195239
>
> Variable: nsqr
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 1
>
> _FillValue :-999
>
> (0)
>
> (0)min=2.07257e-05max=0.00314973
>
> Variable: RI
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[Time | 25] x [bottom_top | 39] x [south_north |
> 195] x [west_east | 135]
>
> Coordinates:
>
> Number Of Attributes: 1
>
> _FillValue :-999
>
> (0)
>
> (0) min=0.0011193max=1793.91
>
> Variable: times
>
> Type: string
>
> Total Size: 200 bytes
>
> 25 values
>
> Number of Dimensions: 1
>
> Dimensions and sizes:[25]
>
> Coordinates:
>
> Number Of Attributes: 2
>
> description : times in file
>
> _FillValue :missing
>
> Variable: ri_plane
>
> Type: float
>
> Total Size: 1360000 bytes
>
> 340000 values
>
> Number of Dimensions: 3
>
> Dimensions and sizes:[Time | 25] x [Vertical | 100] x [Horizontal | 136]
>
> Coordinates:
>
> Number Of Attributes: 3
>
> _FillValue :-999
>
> description :
>
> Orientation : Cross-Sesion: (0,84) to (135,84)
>
> (0)
>
> (0): min=-2263.84max=1622.94
>
>
>
> On Tue, Feb 15, 2011 at 11:20 PM, Mary Haley <haley@ucar.edu
> <mailto:haley@ucar.edu>> wrote:
>
> Wei,
>
> If a part of a plot is not showing up, you must examine your data by
> using "printVarSummary", "print(min(x))", "print(max(x))", etc.
> Since you didn't provide a data file, we are unable to run your
> script and determine if there are data values out-of-range. You will
> have to debug this problem yourself.
>
> You can debug your code by looking at the "contour_ri" plot on its own.
>
> After:
>
>> contour_ri = wrf_contour(a,wks,ri_plane(16,:,:),opts_ri)
>>
>
> add:
>
> draw(contour_ri)
> frame(wks)
>
> This will cause the "ri" contour plot to be drawn on its own, and
> you can see if it looks correct before it gets overlaid onto the
> other plot.
>
> If this plot doesn't look okay, then examine the ri_plane(16,:,:)
> values closely, to make sure they are the right units, etc. You can
> do this by printing the min and max of this data, as a start.
>
> Also, look at "opts_ri" closely. Are you using the correct contour
> levels?
>
> If the plot is okay on its own, then look at the part of your code
> that's doing the overlay. Perhaps one plot is covering up the other one.
>
> --Mary
>
>
>
>
> On Feb 14, 2011, at 8:48 PM, wei huang wrote:
>
>> Hi,
>> Thanks for all. I’m interested to Ri, so for this I made some
>> assumption to calculating it. The code is working but the Ri
>> contours are not appeared in the figure, may be my code is not
>> correct. Please could anyone help me to figure out the problem? My
>> code and the result are attached.
>> thanks
>> wei
>> = = = = = = =
>>
>> 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("wrfout_d03_2000-03-20.nc
>> <http://wrfout_d03_2000-03-20.nc/>","r")
>>
>> T = a->T
>>
>> T = T + 300. ;theta K
>>
>> u = wrf_user_getvar(a,"ua",-1)
>>
>> v = wrf_user_getvar(a,"va",-1)
>>
>> z = wrf_user_getvar(a, "z",-1)
>>
>> rh = wrf_user_getvar(a,"rh",-1)
>>
>>
>> U = (u*u+v*v)^(0.5) ; m/sec
>>
>> copy_VarCoords(u, U )
>>
>> printVarSummary(U)
>>
>> printMinMax(U, True)
>>
>> ;
>>
>> dTdz = center_finite_diff_n(T,z,False,0,1)
>>
>> eps = 1e-02 ; eliminate near zero values
>>
>> dTdz@_FillValue = -999.
>>
>> dTdz = where(dTdz.lt.eps, dTdz@_FillValue, dTdz)
>>
>> copy_VarCoords(T, dTdz )
>>
>> printVarSummary(dTdz)
>>
>> printMinMax(dTdz, True)
>>
>> ;
>>
>> dUdz = center_finite_diff_n(U,z,False,0,1)
>>
>> eps = 1e-02 ; eliminate near zero values
>>
>> dUdz@_FillValue = -999.
>>
>> dUdz = where(dUdz.lt.eps, dUdz@_FillValue, dUdz)
>>
>> copy_VarCoords(u, dUdz )
>>
>> printVarSummary(dUdz)
>>
>> printMinMax(dUdz, True)
>>
>> ;
>>
>> nsqr = ((9.8/T)*(dTdz)); Brunt vaisala frequency
>>
>> copy_VarCoords(T, nsqr)
>>
>> printVarSummary(nsqr)
>>
>> printMinMax(nsqr, True)
>>
>>
>> RI = nsqr/(dUdz)^2 ; to calculate Richardson Number(RI)
>>
>> copy_VarCoords(u, RI )
>>
>> printVarSummary(RI)
>>
>> printMinMax(RI, True)
>>
>> ;***********************************
>>
>> wks = gsn_open_wks("pdf","ri")
>>
>> gsn_define_colormap(wks,"gui_default")
>>
>> res = True ; Set up some basic plot resources
>>
>> res@Footer = False
>>
>> res@NoHeaderFooter = True
>>
>> pltres = True
>>
>> ;*******************************
>>
>> times = wrf_user_list_times (a); get times in the file
>>
>> ntimes = dimsizes(times) ; number of times in the file
>>
>> FirstTime = res
>>
>> mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for
>> the file
>>
>> nd = dimsizes(mdims)
>>
>>
>> ;do it = 10,10;ntimes-1
>>
>> ;***************************************
>>
>> if ( FirstTime ) then ; get height info for labels
>>
>> zmin = 0.
>>
>> zmax = max(z)/1000.
>>
>> nz = floattoint(zmax/2 + 1)
>>
>> FirstTime = False
>>
>> end if
>>
>> ;**************************************
>>
>> ip = 1 ; Just do the one (constant y coord) plot
>>
>> opts = True ; setting start and end times
>>
>> plane = new(4,float)
>>
>> if(ip .eq. 1) then
>>
>> plane = (/ 0,84, 200,84 /) ; start x;y & end x;y point
>>
>> end if
>>
>> ;
>>
>> rh_plane = wrf_user_intrp3d(rh,z,"v",plane,0.,opts)
>>
>> ri_plane = wrf_user_intrp3d(RI,z,"v",plane,0.,opts)
>>
>> printVarSummary(rh_plane)
>>
>>
>> dim = dimsizes(rh_plane) ; Find the data span - for use in labels
>>
>> zspan = dim(0)
>>
>> ; Options for XY Plots
>>
>> opts_xy = res
>>
>> opts_xy@tiYAxisString = "Height (km)"
>>
>> opts_xy@AspectRatio = 0.75
>>
>> 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 = ri_plane@Orientation
>>
>>
>> ; Plotting options for RH
>>
>> opts_rh = opts_xy
>>
>> opts_rh@pmLabelBarOrthogonalPosF = -0.07
>>
>> opts_rh@ContourParameters = (/ 10., 90., 10. /)
>>
>> opts_rh@cnFillOn = True
>>
>>
>> ; Plotting options for ri
>>
>> opts_ri = opts_xy
>>
>> opts_ri@cnInfoLabelOrthogonalPosF = 0.00
>>
>> opts_ri@ContourParameters = (/10. /)
>>
>> 　
>>
>> ; Get the contour info for the rh and ri
>>
>> contour_ri = wrf_contour(a,wks,ri_plane(16,:,:),opts_ri)
>>
>> contour_rh = wrf_contour(a,wks,rh_plane(16,:,:),opts_rh)
>>
>> ; MAKE PLOTS
>>
>> plot = wrf_overlays(a,wks,(/contour_rh,contour_ri/),pltres)
>>
>> ; Delete options and fields, so we don't have carry over
>>
>> delete(opts_ri)
>>
>> delete(opts_rh)
>>
>> delete(ri_plane)
>>
>> delete(rh_plane)
>>
>> end
>>
>> <ri.pdf>_______________________________________________
>> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Feb 16 07:35:59 2011