Re: vertical cross section

Hi Jin-Young,

Can you add this line to your code right before the plot creation?

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

and tell us what the printVarSummary output is?
Thanks, Adam

Jin-Young Kim wrote:
> Greetings,
>
> Thanks for you comment..
>
> But when i tried to apply your comment at the dimension setting.
> All does not applied well.
>
> 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
>
> Do you know other way?
> As a matter of fact, it look the data could not recognize with latitude
> and longitude.
> That's why it's working ugrd(:,19,:) instead of ugrd (:,{19.4},:)
> Do you think so?
>
> Best regards,
> Jin-Young Kim
>
> On Fri, Nov 7, 2008 at 3:47 PM, Adam Phillips <asphilli_at_cgd.ucar.edu
> <mailto:asphilli_at_cgd.ucar.edu>> wrote:
>
> Hi Jin-Young,
>
> Beyond what Dennis said, a few other things I noticed:
>
> 1) Instead of this syntax
>
> ugrd(:,{19.4:19.4},:)
> use this
> ugrd(:,{19.4},:)
>
> Otherwise you are submitting an array to the plotting function that
> is number of levels by 1 by number of longitudes. (3D array)
>
> 2) (:, :) is not needed.. You can change this
>
> psfc = grb_file->PRMSL_GDS0_MSL(:, :)
> to this
>
> psfc = grb_file->PRMSL_GDS0_MSL
>
> 3) These are commented out, but you say this:
> gsn_csm_contour_shade
> when I think you mean this:
> gsn_contour_shade
>
> and you say this:
> SahdeGtContour
> when I think you mean this:
> ShadeGtContour
>
> Adam
>
> 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> <mailto: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/%7Ejin>
> <http://climate.pknu.ac.kr/~jin <http://climate.pknu.ac.kr/%7Ejin>>
>
> Visiting Scholar
> Jin-Young Kim (Jin-Young.Kim_at_noaa.gov
> <mailto:Jin-Young.Kim_at_noaa.gov> <mailto: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
>
>
> --
> --------------------------------------------------------------
> Adam Phillips asphilli_at_ucar.edu
> <mailto:asphilli_at_ucar.edu>
> National Center for Atmospheric Research tel: (303) 497-1726
> ESSL/CGD/CAS fax: (303) 497-1333
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>
>
>
>
> --
> 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)

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk