Re: Overlay for the aviation turbulence icons

From: Adam Phillips <asphilli_at_nyahnyahspammersnyahnyah>
Date: Mon, 02 Jun 2008 16:17:58 -0600

Hi Jung-Hoon,

The email client that I was using this morning mistakingly showed that
you didn't have an attachment. My mistake. I came up with an example NCL
script that shows you how to make use of NCL's font tables to create
PIREP symbols. You can see the font tables here:
http://www.ncl.ucar.edu/Document/Graphics/font_tables.shtml

There are a number of great examples on how to access the various
characters here:
http://www.ncl.ucar.edu/Applications/eqn.shtml

The script shown here creates the attached plot. For most symbols I had
to combine some characters, as most PIREP symbols are not available from
any font table.

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
   wks = gsn_open_wks("ps","pireps")
   res = True
   res_at_mpMinLatF = 24.0 ; map area
   res_at_mpMaxLatF = 50.0 ; latitudes
   res_at_mpMinLonF = -125.0 ; and
   res_at_mpMaxLonF = -70.0 ; longitudes
   res_at_gsnDraw = False
   res_at_gsnFrame = False
   res_at_mpDataBaseVersion = "Ncarg4_1"
   res_at_mpDataSetName = "Earth..1"
   res_at_mpOutlineBoundarySets = "AllBoundaries"
   res_at_mpFillOn = False
   map = gsn_csm_map(wks,res) ; Draw map.
        
   pres = True
   pres_at_txFontHeightF = 0.02
   pres_at_txFontColor = "green"
   pres_at_txFontThicknessF = 2.0
   gsn_text(wks,map,"~F8~f",260.0,40.,pres) ; null
   gsn_text(wks,map,"- - -",270.0,37.,pres) ; smooth-light
   gsn_text(wks,map,"~F34~Y",272.0,42.,pres) ; light
        
   lat = 42.4
   lon = 266.2
   pres_at_txFontColor = "blue"
   gsn_text(wks,map,"~F34~Y",lon,lat,pres)
   gsn_text(wks,map,"- - -",lon,lat-.4,pres) ; light-mod
        
   lat = 45.8
   lon = 266.9
   gsn_text(wks,map,"~F34~Y",lon,lat,pres)
   gsn_text(wks,map,"_",lon-.9,lat,pres)
   gsn_text(wks,map,"_",lon+.9,lat,pres) ; moderate
        
   lat = 44.8
   lon = 261.9
   pres_at_txFontColor = "purple"
   gsn_text(wks,map,"~F34~Y",lon,lat,pres)
   gsn_text(wks,map,"- - -",lon,lat-.4,pres)
   gsn_text(wks,map,"~F34~Y",lon,lat+.4,pres) ; mod-severe
        
   lat = 43.2
   lon = 263.2
   pres_at_txFontColor = "red"
   gsn_text(wks,map,"~F34~Y",lon,lat,pres)
   gsn_text(wks,map,"~F34~Y",lon,lat+.4,pres)
   gsn_text(wks,map,"_",lon-.9,lat,pres)
   gsn_text(wks,map,"_",lon+.9,lat,pres) ; severe
        
   lat = 41.2
   lon = 263.2
   pres_at_txFontColor = "brown"
   gsn_text(wks,map,"~F34~u",lon,lat,pres)
   gsn_text(wks,map,"~F34~Y",lon,lat+.4,pres)
   gsn_text(wks,map,"- - -",lon,lat-.25,pres) ; extreme
   draw(map)
   frame(wks)
end

Best regards,
Adam

김정훈 wrote:
> Dear developers for NCL
>
> Nice to see you. I'm Jung-Hoon Kim at Yonsei University in Seoul, South Korea. I'm Ph.D. Student for atmospheric dynamics.
> I usually use NCL for plotting all of my reseach results. While I plotted turbulence indices such as Ri (Richardson number), TI (Turbulence Index; Ellord and Knapp 1992) and so on, I would like to overlay the turbulence icons observed and reported by PIREPs (pilot reports) using NCL. I couldn't find how can I overlay turbulence icons (light, moderate, severe), even though I tried to find on the NCL web page. Is there something to solve this problem ? Please, let me know, if you have the solutions.
> following figure is the example of overlayed turbulence icons over the CONUS.
>
> Best regards,
>
> Jung-Hoon, Kim
>
> =======================================================
> Jung-Hoon, Kim
> Ph.D. Student
> Laboratory for Atmospheric Dynamics
> Department of Atmospheric Sciences, Yonsei University
> 262 Seongsanno, Seodaemun-gu, Seoul, 120-749, Korea
> Tel.: +82-2-2123-7643, Fax: +82-2-365-5163
> E-mail: kimjh99_at_yonsei.ac.kr
> =======================================================
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
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
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

pireps.gif
Received on Mon Jun 02 2008 - 16:17:58 MDT

This archive was generated by hypermail 2.2.0 : Thu Jun 05 2008 - 16:44:32 MDT