Hi,
Please thanks for your response. For your convenience I indicated the error
line(119) in blue color(contour_ri =
wrf_contour(a,wks,ri_plane(0,:,:),opts_ri)). Your cooperation will be highly
appreciated.
Thanks again,
;*************************************
wks = gsn_open_wks("x11","ri")
gsn_define_colormap(wks,"BlWhRe")
res = True ; Set up some basic plot resources
res@gsnDraw = False ; don't draw
res@gsnFrame = False
pltres = True
;*******************************
FirstTime = res
rh = wrf_user_getvar(a,"rh",-1)
printVarSummary(rh)
;
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
;******************************
; Interpolate data vertically (in z)
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)
printVarSummary(ri_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@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@PlotOrientation = ri_plane@Orientation
; Plotting options for rh
opts_rh = opts_xy
opts_rh@ContourParameters = (/ 10., 90., 10. /)
opts_rh@cnFillOn = True
opts_rh@gsnSpreadColors = True
; Plotting options for ri
opts_ri = opts_xy
opts_ri@ContourParameters = (/.2 /)
; Get the contour info for the rh and ri
contour_rh = wrf_contour(a,wks,rh_plane(0,:,:),opts_rh)
contour_ri = wrf_contour(a,wks,ri_plane(0,:,:),opts_ri)
plot = wrf_overlays(a,wks,(/contour_rh,contour_ri/),pltres)
end
;************************************
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) min(U)=0.0100291 max(U)=62.1447
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: 1
_FillValue : -999
(0)
(0) min=0.0100291 max=62.1447
(0) min(dUdz)=-0.0857187 max(dUdz)=0.195239
(0) min(dUdz)=0.10001 max(dUdz)=0.195239
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.10001 max=0.195239
Variable: N
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) min(RI)=0.103217 max(RI)=1.95101
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.103217 max=1.95101
Variable: rh
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
description : Relative Humidity
units : %
FieldType : 104
MemoryOrder : XYZ
stagger :
coordinates : XLONG XLAT
Variable: rh_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: 4
_FillValue : -999
units : %
description : Relative Humidity
Orientation : Cross-Sesion: (0,84) to (135,84)
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)
fatal:The result of the conditional expression yields a missing value. NCL
can not determine branch, see ismissing function
fatal:Execute: Error occurred at or near line 1896 in file
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:Execute: Error occurred at or near line 2671 in file
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:Execute: Error occurred at or near line 119 in file
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Feb 13 18:49:02 2011
This archive was generated by hypermail 2.1.8 : Tue Feb 15 2011 - 09:43:19 MST