undef("skewT_BackGround_Original") function skewT_BackGround_Original (wks:graphic, Opts:logical) local copt begin ; statement of purpose. ; --------------------- ; this program generates a skew-t, log p thermodynamic diagram. this ; program was derived to reproduce the USAF skew-t, log p diagram. ; (form dod-wpc 9-16-1 current as of march 1978) ; source. ; -------- ; Tom Schlatter and Joe Wakefield [NOAA/PROFS] supplied fortran codes. ; ; history. ; -------- ; don baker 01 jul 85 original version. [NOAA/PROFS] ; don baker 01 dec 85 updated for product version. ; dennis shea oct 98 created the NCL version ; dennis shea 9 feb 99 fix: MoistAdiabat labels at top of curves localOpts = True localOpts@DrawIsotherm = True localOpts@DrawIsobar = True localOpts@DrawMixRatio = True localOpts@DrawDryAdiabat = True localOpts@DrawMoistAdiabat = True ; aka: saturation or pseudo adibat localOpts@DrawWind = True localOpts@DrawStandardAtm = True localOpts@DrawColLine = True localOpts@DrawColAreaFill = False localOpts@DrawFahrenheit = True ; Fahrenheit "x" axis localOpts@DrawHeightScale = False localOpts@DrawHeightScaleFt = True ;default is feet [otherwise km] localOpts@DrawStandardAtmThk= 2.0 localOpts@UseMyColorMap = False localOpts@Font = "helvetica" ; NCL resource names localOpts@tiMainString = " " localOpts@vpXF = 0.07 localOpts@vpYF = 0.925 localOpts@vpWidthF = 0.85 localOpts@vpHeightF = 0.85 ; Override localOpts ; with input Opts if (Opts .and. .not.any(ismissing(getvaratts(Opts)))) then localAtts = getvaratts(localOpts) OptsAtts = getvaratts(Opts) nOpts = dimsizes (OptsAtts) do i=0,nOpts-1 ; loop and add/override attributes localOpts@$OptsAtts(i)$ = Opts@$OptsAtts(i)$ ; new att end do end if if (isatt(localOpts,"PrintOpts").and.localOpts@PrintOpts) then print(localOpts) end if if(.not.localOpts@UseMyColorMap) then cmap_orig = fix_skewt_colormap_for_ncgm(wks,localOpts@UseMyColorMap) end if ; --- declare isotherm [C] values and pressures [hPa] where isotherms ; --- intersect the edge of the skew-t diagram. temp = (/-100.,-90.,-80.,-70.,-60.,-50.,-40.,-30. \ , -20.,-10., 0., 10., 20., 30., 40., 50. /) lendt= (/ 132., 181., 247., 337., 459., 625., 855. \ ,1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050./) rendt= (/ 100., 100., 100., 100., 100., 100., 100. \ , 135., 185., 251., 342., 430., 500., 580., 730., 993./) ntemp= dimsizes (temp) ; isotherms if (dimsizes(lendt).ne.dimsizes(rendt) .or. \ dimsizes(lendt).ne.ntemp) then print ("Error: dimsizes temp, lendt, rendt do not match") exit end if ; --- declare pressure values [hPa] and x coordinates of the endpoints of each ; --- isobar. these x,y values are computed from the equations in the ; --- transform functions listed at the beginning of this program. ; --- refer to a skew-t diagram for reference if necessary. pres = (/1050.,1000.,850.,700.,500.,400.,300.,250.,200.,150.,100./) xpl = (/-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0 \ ,-19.0,-19.0,-19.0/) xpr = (/27.1,27.1,27.1,27.1,22.83,18.6,18.6,18.6,18.6,18.6,18.6/) npres= dimsizes (pres) ; pressures if (dimsizes(xpl).ne.dimsizes(xpr) .or. dimsizes(xpl).ne.npres) then print ("Error: dimsizes pres, xpl, xpr do not match") exit end if ; --- declare adiabat values [C] and pressures where adiabats intersect the ; --- edge of the skew-t diagram. refer to a skew-t diagram if necessary. theta= (/-30.,-20.,-10.,0.,10.,20.,30.,40.,50.,60.,70.,80.,90. \ ,100.,110.,120.,130.,140.,150.,160.,170. /) lendth= (/880.,670.,512.,388.,292.,220.,163.,119.,100.,100.,100. \ ,100.,100.,100.,100.,100.,100.,100.,100.,100.,100. /) rendth= (/1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050. \ ,1003.,852.,728.,618.,395.,334.,286., 245., 210. \ , 180.,155.,133.,115. /) ntheta= dimsizes (theta) ; dry adiabats if (dimsizes(lendth).ne.dimsizes(rendth) .or. \ dimsizes(lendth).ne.ntheta) then print ("Error: dimsizes theta, lendth, rendth do not match") exit end if ; --- declare moist adiabat values and pressures of the tops of the ; --- moist adiabats. all moist adiabats to be plotted begin at 1050 mb. pseudo = (/ 32., 28., 24., 20., 16., 12., 8. /) lendps = (/250.,250.,250.,250.,250.,250.,250. /) npseudo= dimsizes (pseudo) ; moist adiabats ; --- declare mixing ratio lines. all mixing ratio lines will begin ; --- at 1050 mb and end at 400 mb. mixrat= (/ 20.,12.,8.,5.,3.,2.,1. /) nmix = dimsizes (mixrat) ; mixing ratios ; --- declare local stuff: arrays/variables for storing x,y positions ; --- during iterations to draw curved line, etc sx = new (200, float) sy = new (200, float) xx = new ( 2, float) yy = new ( 2, float) label = new ( 1, string) m2f = 3.2808 ; meter-to-feet f2m = 1./3.2808 ; feet-to-meter ; --- define absolute x,y max/min bounds corresponding to the outer ; --- edges of the diagram. these are computed by inserting the appropriate ; --- pressures and temperatures at the corners of the diagram. yminP = get_res_value_keep(Opts,"skewtYMinF",100.) ymaxP = get_res_value_keep(Opts,"skewtYMaxF",1050.) ; xmin = skewtx (-33.6, skewty(1050.)) [t=deg C] xmin =-19.0 ; xmin = skewtx (-109.1,skewty(100.)) [t=deg C] xmax = 27.1 ; xmax = skewtx (51.75, skewty(1050.)) [t=deg C] ;ymax =-0.935 ; ymax = skewty ( ymaxP ) ;ymin =44.061 ; i ymin = skewty ( yminP ) ; print(ymin+" "+ymax) ; --- specify arrays to hold corners of the diagram in x,y space. xc = (/ xmin, xmin, xmax, xmax, 18.60, 18.6, xmin /) yc = (/ ymin, ymax, ymax, 9.0, 17.53, ymin, ymin /) ; ---- depending on how options are set create Standard Atm Info if (localOpts@DrawStandardAtm .or. \ localOpts@DrawHeightScale .or. localOpts@DrawWind) then ; U.S. Standard ATmosphere zsa = fspan (0. , 16., 17) ; (km) source Hess/Riegel psa = (/ 1013.25, 898.71, 794.90, 700.99, 616.29, 540.07 \ , 471.65, 410.46, 355.82, 307.24, 264.19 \ , 226.31, 193.93, 165.33, 141.35, 120.86,103.30 /) tsa = (/ 15.0 , 8.5 , 2.0 , -4.5 , -11.0 , -17.5 \ , -24.0 , -30.5 , -37.0 , -43.5 , -50.0 \ , -56.5 , -56.5 , -56.5 , -56.5 , -56.5 , -56.5 /) nlvl= dimsizes (psa) end if ; PLOT if (localOpts@DrawColLine) then colGreen = "Green" colBrown = "Brown" colTan = "Tan" else colGreen = "Foreground" colBrown = "Foreground" colTan = "Foreground" end if ; --- draw outline of the skew-t, log p diagram. proceed in the upper left ; --- corner of the diagram and draw counter-clockwise. the array locations ; --- below that are hardcoded refer to points on the background where the ; --- skew-t diagram deviates from a rectangle, along the right edge. remember, ; --- the set call defines a rectangle, so as long as the boundary is along ; --- the edge of the set call, the points plotted are combinations of the min ; --- and max x,y values in the set call. if (localOpts@DrawFahrenheit) then tf = ispan (-20, 100, 20) ; deg F [nice round numbers] tc = 0.55555*(tf-32.) ; deg C corresponding to tf else tc = ispan(-30, 40, 10)*1.0 ; deg C (nice round numbers) end if xyOpts = True xyOpts@gsnDraw = False ; Don't draw the plot or advance the xyOpts@gsnFrame = False ; frame in the call to gsn_xy. ; specify location/size of skewT diagram xyOpts@vpXF = localOpts@vpXF xyOpts@vpYF = localOpts@vpYF xyOpts@vpWidthF = localOpts@vpWidthF xyOpts@vpHeightF = localOpts@vpHeightF xyOpts@tmYLMode = "Explicit" ; Define y tick mark labels. xyOpts@tmYLValues = skewty ( pres(1:)) ; skip 1050 xyOpts@tmYLLabels = (/"1000","850","700","500","400","300" \ ,"250" , "200","150","100"/) xyOpts@tmXBMode = "Explicit" ; Define x tick mark labels. xyOpts@tmXBValues = skewtx (tc, skewty(1050.)) ; transformed vals if (localOpts@DrawFahrenheit) then xyOpts@tmXBLabels = tf ; plot the nice deg F else xyOpts@tmXBLabels = tc ; plot the nice deg C end if xyOpts@trXMinF = xmin xyOpts@trXMaxF = xmax xyOpts@trYMinF = ymax ; Note: ymin > ymax xyOpts@trYMaxF = ymin xyOpts@xyComputeXMin= False xyOpts@xyComputeXMax= False xyOpts@xyComputeYMin= False xyOpts@xyComputeYMax= False xyOpts@tmXTOn = False xyOpts@tmYROn = False xyOpts@tmXTBorderOn = False xyOpts@tmXBBorderOn = False xyOpts@tmYRBorderOn = False xyOpts@tmYLBorderOn = False ; "tune" the plot xyOpts@tmXBMajorLengthF = 0.01 xyOpts@tmXBMajorThicknessF = 1.0 ; default is 2.0 xyOpts@tmXBMajorOutwardLengthF = 0.01 xyOpts@tmXTMajorLengthF = 0. xyOpts@tmYLMajorLengthF = 0. xyOpts@tmXBLabelFontHeightF = 0.014 xyOpts@tmYLLabelFontHeightF = xyOpts@tmXBLabelFontHeightF if (localOpts@DrawFahrenheit) then xyOpts@tiXAxisString = "Temperature (F)" else xyOpts@tiXAxisString = "Temperature (C)" end if xyOpts@tiXAxisFont = localOpts@Font xyOpts@tiXAxisFontColor = "Foreground" ; colTan xyOpts@tiXAxisFontHeightF = 0.0125 xyOpts@tiYAxisFont = localOpts@Font xyOpts@tiYAxisString = "P (hPa)" xyOpts@tiYAxisOffsetXF = 0.0200 xyOpts@tiYAxisOffsetYF = -0.0200 xyOpts@tiYAxisFontColor = "Foreground" ; colTan xyOpts@tiYAxisFontHeightF = xyOpts@tiXAxisFontHeightF xyOpts@tiMainString = localOpts@tiMainString xyOpts@tiMainFont = localOpts@Font xyOpts@tiMainFontColor = "Foreground" xyOpts@tiMainFontHeightF = 0.025 if (isatt(localOpts,"tiMainFontHeightF")) then xyOpts@tiMainFontHeightF = localOpts@tiMainFontHeightF end if xyOpts@tiMainOffsetXF = -0.1 if (isatt(localOpts,"tiMainOffsetXF")) then xyOpts@tiMainOffsetXF = localOpts@tiMainOffsetXF end if xyOpts@tiMainOffsetYF = 0.0 if (isatt(localOpts,"tiMainOffsetYF")) then xyOpts@tiMainOffsetYF = localOpts@tiMainOffsetYF end if if (localOpts@DrawHeightScale) then ; special for rhs xyOpts@trXMaxF = skewtx (55. , skewty(1013.)) ; extra wide xyOpts@tmYUseLeft = False ; Keep right axis independent of left. xyOpts@tmYRBorderOn = True xyOpts@tmYROn = True ; Turn on right axis tick marks. xyOpts@tmYRLabelsOn = True ; Turn on right axis labels. xyOpts@tmYRLabelFontHeightF = xyOpts@tmYLLabelFontHeightF xyOpts@tmYRMajorThicknessF = xyOpts@tmXBMajorThicknessF xyOpts@tmYRMajorLengthF = xyOpts@tmXBMajorLengthF xyOpts@tmYRMajorOutwardLengthF = xyOpts@tmXBMajorOutwardLengthF xyOpts@tmYRMinorOn = False ; No minor tick marks. xyOpts@tmYRMode = "Explicit" ; Define tick mark labels. zkm = fspan (0. , 16., 17) pkm = ftcurv(zsa,psa,zkm) zft = (/ 0., 2., 4., 6., 8.,10.,12.,14.,16. \ ,18.,20.,25.,30.,35.,40.,45.,50. /) pft = ftcurv(zsa,psa,zft*f2m) ; p corresponding to zkm if (localOpts@DrawHeightScaleFt) then znice = zft pnice = skewty(pft) zLabel= "Height (1000 Feet)" else znice = zkm pnice = skewty(pkm) zLabel= "Height (Km)" end if xyOpts@tmYRValues = pnice ; At each "nice" pressure value, xyOpts@tmYRLabels = znice ; put a "height" value label. end if ; DrawHeightScale xy = gsn_xy (wks,xc,yc,xyOpts) ; draw outline; x and y axis ; right *label* MUST be added AFTER xy created if (localOpts@DrawHeightScale) then txOpts = True txOpts@txAngleF = 270. txOpts@txFontColor = "Foreground" ; colTan txOpts@txFontHeightF = xyOpts@tmYLLabelFontHeightF xlab = skewtx (53., skewty(1013.)) ylab = skewty (350.) gsn_text (wks,xy,zLabel,xlab,ylab,txOpts) ; label rhs delete (txOpts) end if ; DrawHeightScale ; --- Color Fill Area of plot if (localOpts@DrawColAreaFill) then color1 = "PaleGreen" ; "LightGreen" if (isatt(localOpts,"DrawColAreaColor")) then color1 = localOpts@DrawColAreaColor end if color2 = "White" gsOpts = True do i=0,ntemp-2 if (i%2 .eq. 0) then ; alternate colors gsOpts@gsFillColor = color1 else gsOpts@gsFillColor = color2 end if nx = 3 ; this handles most cases sy(0) = skewty( lendt(i) ) sx(0) = skewtx( temp(i) , sy(0) ) sy(1) = skewty( lendt(i+1) ) sx(1) = skewtx( temp(i+1), sy(1) ) sy(2) = skewty( rendt(i+1) ) sx(2) = skewtx( temp(i+1), sy(2) ) sy(3) = skewty( rendt(i) ) sx(3) = skewtx( temp(i) , sy(3) ) ; special cases if (temp(i).eq.-40.) then nx = 5 sy(0:nx) = (/ 3.00, ymax, -0.935, 38.32, 44.06, 44.06 /) sx(0:nx) = (/ -18.88, xmin,-17.05 , 18.55, 18.55, 18.36 /) end if if (temp(i).eq.0.) then nx = 4 sy(0:nx) = (/ -0.935,-0.935, 16.15, 17.53, 20.53 /) sx(0:nx) = (/ -0.850, 4.55 , 20.05, 18.55, 18.55 /) end if if (temp(i).eq.30.) then nx = 4 sy(0:nx) = (/-0.935, -0.935, 6.02, 9.0 , 10.42 /) sx(0:nx) = (/15.35 , 20.75 , 27.06, 27.06, 25.65 /) end if gsn_polygon(wks, xy, sx(0:nx),sy(0:nx),gsOpts) end do ; special cases gsOpts@gsFillColor = color2 sy(0:2) = (/ 44.06, 44.06, 38.75 /) ; upper left triangle sx(0:2) = (/-14.04,-18.96,-18.86 /) gsn_polygon(wks, xy, sx(0:2),sy(0:2),gsOpts) gsOpts@gsFillColor = color2 sy(0:2) = (/ ymax, 0.13, ymax /) ; lower right triangle sx(0:2) = (/ xmax, xmax, 26.15 /) ; skewtx(51.6,ymax) gsn_polygon(wks, xy, sx(0:2),sy(0:2),gsOpts) delete (gsOpts) end if ; --- draw diagonal isotherms. ; [brown with labels interspersed at 45 degree angle] ; http://www.ncl.ucar.edu/Document/HLUs/Classes/GraphicStyle.shtml if (localOpts@DrawIsotherm) then gsOpts = True gsOpts@gsLineDashPattern = 0 ; solid gsOpts@gsLineColor = colTan gsOpts@gsLineThicknessF = 1.0 ;gsOpts@gsLineLabelFontColor = colTan ;gsOpts@gsLineLabelFontHeightF = 0.0125 txOpts = True txOpts@txAngleF = 45. txOpts@txFontColor = gsOpts@gsLineColor txOpts@txFontHeightF = 0.0140 txOpts@txFontThicknessF = 1.0 do i=0,dimsizes(temp)-3 yy(1) = skewty(rendt(i)) xx(1) = skewtx(temp(i),yy(1)) yy(0) = skewty(lendt(i)) xx(0) = skewtx(temp(i),yy(0)) ;gsOpts@gsLineLabelString = toint(temp(i)) ; floattointeger(temp(i)) gsn_polyline (wks,xy,xx,yy,gsOpts) xlab = xx(1)+0.625 ylab = yy(1)+0.55 label = toint(temp(i)) ; floattointeger(temp(i)) gsn_text(wks,xy,label,xlab,ylab,txOpts) end do delete (gsOpts) delete (txOpts) end if ; DrawIsotherm ; --- draw horizontal isobars. if (localOpts@DrawIsobar) then gsOpts = True gsOpts@gsLineDashPattern = 0 ; solid gsOpts@gsLineColor = colTan gsOpts@gsLineThicknessF = 1.0 ;gsOpts@gsLineLabelFontColor = colTan ;gsOpts@gsLineLabelFontHeightF = 0.0125 do i=0,npres-1 xx(0) = xpl(i) xx(1) = xpr(i) ypl = skewty(pres(i)) yy(0) = ypl yy(1) = ypl gsn_polyline (wks,xy,xx,yy,gsOpts) end do ; end "i=1,npres" delete (gsOpts) end if ; DrawIsobar ; --- draw saturation mixing ratio lines. these lines run between 1050 and 400 ; --- mb. the 20 line intersects the sounding below 400 mb, thus a special case ; --- is made for it. the lines are dashed green. the temperature where each line ; --- crosses 400 mb is computed in order to get x,y locations of the top of ; --- the lines. if (localOpts@DrawMixRatio) then gsOpts = True ; polyline [graphic style] opts gsOpts@gsLineThicknessF = 1.0 gsOpts@gsLineDashPattern = 2 ; sat mix ratio only gsOpts@gsLineColor = colGreen ; " txOpts = True txOpts@txAngleF = 65. ; " txOpts@txFontColor = colGreen ; " txOpts@txFontHeightF = 0.0100 ; " yy(1) = skewty(400.) ; y at top [right end of slanted line] yy(0) = skewty(1000.) ; y at bottom of line [was 1050.] do i=0,nmix-1 if (mixrat(i).eq.20.) then yy(1) = skewty(440.) ;tmix = THERMO::tmr(mixrat(i),440.) tmix = tmr_skewt(mixrat(i),440.) else yy(1) = skewty(400.) ;tmix = THERMO::tmr(mixrat(i),400.) tmix = tmr_skewt(mixrat(i),400.) end if xx(1) = skewtx(tmix,yy(1)) ;tmix = THERMO::tmr(mixrat(i),1000.) ; was 1050 tmix = tmr_skewt(mixrat(i),1000.) ; was 1050 xx(0) = skewtx(tmix,yy(0)) gsn_polyline (wks,xy,xx,yy,gsOpts) ; dashed green xlab = xx(0)-0.25 ylab = yy(0)-0.45 label = toint(mixrat(i)) ; floattointeger(mixrat(i)) gsn_text(wks,xy,label,xlab,ylab,txOpts) end do ; end "i=0,nmix-1" delete (gsOpts) delete (txOpts) end if ; DrawMixRatio ; --- draw dry adiabats. iterate in 10 mb increments to compute the x,y ; --- points on the curve. if (localOpts@DrawDryAdiabat) then gsOpts = True gsOpts@gsLineDashPattern = 0 gsOpts@gsLineColor = colTan gsOpts@gsLineThicknessF = 1.0 txOpts = True txOpts@txAngleF = 300. txOpts@txFontColor = colTan txOpts@txFontHeightF = 0.01 txOpts@txFontThicknessF = 1.0 pinc = 10. do i=0,ntheta-1 p = lendth(i)-pinc do j=0,dimsizes(sy)-1 p = p+pinc > yminP ;print(""+p) if (p.gt.rendth(i)) then sy(j) = skewty(rendth(i)) ;t = THERMO::tda(theta(i),p) ; get temp on dry adiabat at p t = tda_skewt(theta(i),p) ; get temp on dry adiabat at p sx(j) = skewtx(t,sy(j)) break end if sy(j) = skewty(p) ;t = THERMO::tda(theta(i),p) t = tda_skewt(theta(i),p) sx(j) = skewtx(t,sy(j)) end do ; end "j=0,dimsizes(sy)-1" ;gsn_polyline (wks,xy,sx(:j-1),sy(:j-1),gsOpts) ; whole line if(j.gt.2) if (theta(i).lt.170.) then gsn_polyline (wks,xy,sx(1:j-1),sy(1:j-1),gsOpts); label room ylab = skewty(lendth(i)+5. > yminP+5) ;t = THERMO::tda(theta(i),lendth(i)+5.>yminP+5) t = tda_skewt(theta(i),lendth(i)+5.>yminP+5) xlab = skewtx(t,ylab) label = toint(theta(i)) ; floattointeger(theta(i)) gsn_text(wks,xy,label,xlab,ylab,txOpts) else ; no laabel gsn_polyline (wks,xy,sx(:j-1),sy(:j-1),gsOpts) ; whole line end if end if end do ; end "i=0,ntheta-1 delete (gsOpts) delete (txOpts) end if ; DryAdiabat ; --- draw moist adiabats up to 230 [was 250] mb. ; --- draw the lines. iterate in 10 mb increments from 1060 mb. if (localOpts@DrawMoistAdiabat) then gsOpts = True gsOpts@gsLineColor = colGreen gsOpts@gsLineThicknessF = 0.5 gsOpts@gsLineDashPattern = 0 txOpts = True txOpts@txAngleF = 0. txOpts@txFontColor = colGreen txOpts@txFontHeightF = 0.0125 txOpts@txFontThicknessF = 1.0 pinc = 10. do i=0,npseudo-1 p = 1060. do j=0,dimsizes(sy)-1 p = p-pinc if (p.lt. 230) then ; was "250" break end if sy(j) = skewty(p) ;t = THERMO::satlft(pseudo(i),p) ; temp on moist adiabat at p t = satlft_skewt(pseudo(i),p) ; temp on moist adiabat at p sx(j) = skewtx(t,sy(j)) ;print ("j="+j+" p="+p+" t="+t+" sy(j)="+sy(j)+" sx(j)="+sx(j) ) end do ; end "j=0,dimsizes(sy)-1" gsn_polyline (wks,xy,sx(:j-1),sy(:j-1),gsOpts) ylab = skewty(p+0.5*pinc) ;t = THERMO::satlft(pseudo(i),p+0.75*pinc) t = satlft_skewt(pseudo(i),p+0.75*pinc) xlab = skewtx(t,ylab) ;label = floattointeger(pseudo(i)) ; 9 Feb 99 fix label = toint(pseudo(i)) ; Jan 2013 gsn_text(wks,xy,label,xlab,ylab,txOpts) end do ; end "i-0,dimsizes(pseudo)-1" delete (gsOpts) delete (txOpts) end if ; DrawMoistAdiabat if (localOpts@DrawStandardAtm) then gsOpts = True gsOpts@gsLineColor = colTan gsOpts@gsLineThicknessF = localOpts@DrawStandardAtmThk gsOpts@gsLineDashPattern = 0 do i=0,nlvl-1 sy(i) = skewty (psa(i)) sx(i) = skewtx (tsa(i), sy(i) ) end do gsn_polyline (wks,xy,sx(0:nlvl-1),sy(0:nlvl-1),gsOpts) delete (gsOpts) end if ; DrawStandardAtm ; --- draw vertical line upon which to plot wind barbs. if (localOpts@DrawWind) then gsOpts = True gsOpts@gsLineColor = "Foreground" gsOpts@gsLineThicknessF = 0.5 gsOpts@gsLineDashPattern = 0 gsOpts@gsMarkerIndex = 4 ; "hollow_circle"=> std pres gsOpts@gsMarkerColor = "Foreground" presWind = pres presWind(0) = 1013. ; override 1050 xWind = skewtx (45. , skewty(presWind(0))) sx(0:npres-1) = xWind ; "x" location of wind plot sy(0:npres-1) = skewty(presWind) gsn_polyline (wks,xy,sx(0:npres-1),sy(0:npres-1),gsOpts) gsn_polymarker (wks,xy,sx(1:npres-1),sy(1:npres-1),gsOpts) ; zwind => Pibal reports zftWind = (/ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10. \ ,12.,14.,16.,18.,20.,25.,30.,35.,40.,45.,50. /) zkmWind = zftWind*f2m pkmWind = ftcurv(zsa,psa,zkmWind) nzkmW = dimsizes (zkmWind) sx(0:nzkmW-1) = xWind ; "x" location of wind plot sy(0:nzkmW-1) = skewty(pkmWind) gsOpts@gsMarkerIndex = 16 ; "circle_filled" -> Pibal gsOpts@gsMarkerSizeF = 0.0035 ; 0.007 is default gsOpts@gsMarkerThicknessF= 0.5 ; 1.0 is default gsn_polymarker (wks,xy,sx(0:nzkmW-1),sy(0:nzkmW-1),gsOpts) delete (gsOpts) end if ; DrawWind ;---Restore color map, if necessary if(.not.localOpts@UseMyColorMap) then gsn_define_colormap(wks,cmap_orig) end if return (xy) ; return object end