fatal:The result of the conditional expression yields a missing value

From: wei huang <whua27_at_nyahnyahspammersnyahnyah>
Date: Sun Feb 13 2011 - 08:44:54 MST

Hi,

I’m interested in vertically cross sectional plot using wrf simulated data.
I have tried to fix the problem but could not succeed. When I running the
script I got the following errors:

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

Here is a part of script and print. I would greatly appreciate if somebody
could help me to fix this problem.

Thanks in advance,

Wei

= = = = = = = = = = =

  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)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Feb 13 08:45:01 2011

This archive was generated by hypermail 2.1.8 : Tue Feb 15 2011 - 09:43:19 MST