vertical cross section

Greeting,

I have some problem to draw vertical cross section.

plot=gsn_csm_contour(wks,ugrd(:,{19.4:19.4},:),res)

It shows following error ..

fatal:NclOneDValGetClosestIndex: finish coordinate index out of range, can't
continue
fatal:Could not obtain coordinate indexes, unable to perform subscript
fatal:Execute: Error occurred at or near line 52 in file sim_wind.ncl

In addition, my data is isobaric level,
when i draw axis level 200 to 1000,, i would like to draw 1000 to 200..

Would you give me comment..?

Here is scripts

;/opt/ncl/5.0.0_bin/bin/
;/opt/ncl/5.0.0_bin/bin/ncl_filedump
;/opt/ncl/5.0.0_bin/bin/ncl_convert2nc
;/opt/ncl/5.0.0_intel_9.1/bin
;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

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/csm/shea_util.ncl"

begin

dir = "../DENNIS04L.2005070700/PSTP/postprd/"

 ifim = sprinti("wrfprs_d02.%0.2i"+".grb",ispan(24,30,6))
 ofim = sprinti("windx_d02h_%0.2i",ispan(24,30,6))
 titl = sprinti("%0.2i",ispan(24,30,6))

 do i= 0, dimsizes(ifim)-1

 grb_file = addfile(dir+ifim(i),"r")
 names = getfilevarnames (grb_file)
 print (names)

 psfc = grb_file->PRMSL_GDS0_MSL(:, :)
 psfc = psfc*0.01
 psfc_at_units = "hPa" ; Pa

 u10m = grb_file->U_GRD_GDS0_HTGL(:, :)
 ugrd = grb_file->U_GRD_GDS0_ISBL(:, :,:)
 u10m_at_units = "m/s" ; m/s

 v10m = grb_file->V_GRD_GDS0_HTGL(:, :)
 vgrd = grb_file->V_GRD_GDS0_ISBL(:, :,:)
 v10m_at_units = "m/s" ; m/s

 wind=sqrt(ugrd*ugrd+vgrd*vgrd)
************************************************
 wks =gsn_open_wks("ps",ofim(i))
 gsn_define_colormap(wks,"gui_default")

  res = True
  res_at_cnFillOn = True
  res_at_tiMainString = "Wind Speed (m/s) " + titl(i) + "h Fcst.";title
  res_at_cnLevelSpacingF = 5. ; contour spacing

; ({700},{-30:50},{ 135:270})
 plot=gsn_csm_contour(wks,ugrd(:,{19.4:19.4},:),res)
; plot=SahdeGtContour(plot,-999,10)
; opt = True
; opt_at_gsnShadeFillType = "Pattern"
; opt_at_gsnShadeHigh = 100

; plot=gsn_csm_contour_shade(plot,-999,25,opt)

 end do

end

-- 
Ph. D. Student
Jin-Young Kim (jin_at_climate.pknu.ac.kr)
Integrated Climate System Modeling Lab.
Department of Environmental Atmospheric Sciences
Pukyong National University, Busan 608-737, South Korea
Web : http:climate.pknu.ac.kr/~jin
Visiting Scholar
Jin-Young Kim (Jin-Young.Kim_at_noaa.gov)
Global Systems Division
Earth System Research Laboratory
R/GSD7
325 Broadway
Boulder, CO80305-3328
+1-303- 497-4422 (Office)
+1-303- 497-7262 (FAX)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Nov 07 2008 - 13:56:46 MST