Re: vertical cross section

If the coordinate you are specifying is not
between or is not *exactly* 19.4, NCL's
coordinate subscripting will not 'find' it.

ugrd(:,{19.3:19.5},:) or ugrd(:,{19.35:19.45},:)

Good luck

Jin-Young Kim wrote:
> 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 <mailto: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 <http://climate.pknu.ac.kr/~jin>
>
> Visiting Scholar
> Jin-Young Kim (Jin-Young.Kim_at_noaa.gov <mailto: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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Nov 07 2008 - 15:30:57 MST