The skewT was written by me 10 years ago.
I did it to learn various graphical aspects of NCL.
It was designed to reproduce the "USAF Skew-t, log p diagram.
(form dod-wpc 9-16-1)". When I donated it to the NCL project,
I added a few options. However, doing things like
eliminating the lower right hand corner and being abe to
only plot 1000 to 500 mb were not included.
The following shows how to eliminate the height scale
and, also, the wind barbs.
Good Luck
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
;**************************************************
begin
; --- Read Data; no winds; set to _FillValue------------;
diri = "./" ; input directory
fili = "sounding.testdata"
nlvl = 30
ncol = 16
TestData = asciiread (diri+fili , (/nlvl,ncol/), "float")
p = TestData (:,1)
z = TestData (:,2)
q = TestData (:,4)
tc = TestData (:,5) + 2. ; for demo purposes
tdc = TestData (:,9)
; create BOGUS winds
wspd = new ( nlvl, "float", -999.) ; wspd@_FillValue = -999.
wdir = new ( nlvl, "float", -999.) ; wdir@_FillValue = -999.
;*************************
; create plot
;*************************
wks = gsn_open_wks ("x11", "skewt_2")
; --- Create background skew-T; plot sounding ---------------------
skewtOpts = True
skewtOpts_at_DrawColAreaFill = True ; default is False
skewtOpts_at_tiMainString = "Raob Data; No Winds"
skewtOpts_at_DrawHeightScale = False
skewtOpts_at_DrawWind = False
skewt_bkgd = skewT_BackGround (wks, skewtOpts)
dataOpts = False ; no options
skewt_data = skewT_PlotData (wks, skewt_bkgd, p,tc,tdc,z \
, wspd,wdir, dataOpts)
draw (skewt_bkgd)
draw (skewt_data)
frame(wks)
Wolfgang Langhans wrote:
> HI,
> how can I not plot the wind barbs at the side of a skew_t plot, but only
> the temperatures. Is it then possible to ged rid of the extension of the
> diagram in the lower right corner?
>
> Thanks for helping,
> Wolfgang
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 18 2008 - 06:21:02 MDT
This archive was generated by hypermail 2.2.0 : Tue Sep 23 2008 - 14:12:01 MDT