Plots station model data on the given workstation.
procedure wmstnm ( wks : graphic, x : float, y : float, imdat : string ; of 50 characters )
An NCL workstation identifier specifying where you want to draw your station model data. The wks identifier is one returned either from calling gsn_open_wks or calling create to create a Workstation object.
World coordinate values specifying the center position of the base of the wind barb shaft. These may be multi-dimensioned arrays. They must have the same shape.imdat
A string (or array) of 50 characters encoded as per the WMO/NOAA guidelines. If an array, it must have the same shape as the x and y arrays. In more detail (where the characters are numbered from left to right, starting at character number 0):
|character 0||=||iR||-||the precipitation data indicator|
|character 1||=||iX||-||weather data and station type indicator|
|character 2||=||h||-||height above ground of base of lowest cloud|
|characters 3-4||=||VV||-||visibility in miles and fractions|
|character 5||=||N||-||total amount of cloud cover|
|characters 6-7||=||dd||-||direction from which wind is blowing|
|characters 8-9||=||ff||-||wind speed in knots|
|character 11||=||sn||- sign of temperature|
|characters 12-14||=||TTT||- current air temperature|
|character 16||=||sn||- sign of temperature|
|characters 17-19||=||Td||- dew point|
|characters 21-24||=||PO||- station pressure (not plotted)|
|characters 26-29||=||PPPP||- pressure reduced to sea level|
|character 31||=||a||- characteristic of barograph|
|characters 32-34||=||ppp||- pressure change, last 3 hrs.|
|characters 36-38||=||RRR||- precipitation|
|character 39||=||tR||- time duration of precipitation|
|characters 41-42||=||ww||- present weather|
|character 43||=||W1||- most significant past weather|
|character 44||=||W2||- 2nd most sig. past weather|
|character 46||=||Nh||- Fraction of sky cover|
|character 47||=||CL||- cloud type, low clouds|
|character 48||=||CM||- cloud type, medium clouds|
|character 49||=||CH||- cloud type, high clouds|
The appearance of the wind barbs that are part of the station model data is controlled by the parameters that control wind barbs. See the procedure wmbarb for details.
In addition to the wind barb parameters, the parameters VVC, WBC, WBL, and UNT apply to the station model display. WBC specifies the diameter of the sky cover circle at the base of the wind barb, WBL specifies the size of the text labels in the station model display, and UNT specifies whether you want to use Imperial units (the default) or metric units.
The procedure wmstnm does not call frame.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" ; ; Draw station model data over a map. ; begin ; ; Define a color map and create an X11 workstation. ; cmap = (/ \ (/ 1., 1., 1. /), \ ; color index 0 - white (/ 0., 0., 0. /) \ ; color index 1 - black /) wks = gsn_open_wks("x11","test_wmap") gsn_define_colormap(wks,cmap) ; ; Specify the station model data. ; imdat="11212800201001120000300004014752028601117706086792" ; ; Create a map over the central United States, but do not draw it. ; resources = True resources@gsnFrame = False ; Turn off drawing resources@mpProjection = "Stereographic" resources@mpCenterLonF = -95.0 resources@mpCenterLatF = 90.0 resources@mpLimitMode = "Corners" ; Limit the map view. resources@mpLeftCornerLonF = -103.722 resources@mpLeftCornerLatF = 33.968 resources@mpRightCornerLonF = -82.5161 resources@mpRightCornerLatF = 49.054 resources@mpOutlineBoundarySets = "AllBoundaries" resources@mpUSStateLineThicknessF = 1.25 resources@mpGridAndLimbOn = False resources@vpXF = 0.15 resources@vpYF = 0.85 resources@vpWidthF = 0.7 resources@vpHeightF = 0.7 map = gsn_map(wks,"Stereographic",resources) ; ; Specify the size of the wind barb (this scales the entire ; station model data plot. ; wmsetp("wbs",0.12) ; ; Since we want to plot the station model data over a map ; projection, we have to inform NCL that that is what we ; are doing. Setting a value for the internal parameter ; "ezf" to something other that its default of -1 does this. ; In addition to this, setting "ezf" will also attenuate the ; direction of the wind barb based on its latitude. For ; example, a wind direction of due north (as in this example) ; will show the wind barb pointing in the direction of the ; north pole, rather than straight up. More accurately, the ; barb is drawn on the map so that it is tangent to the projection ; of the local meridian and points northward along that meridian. ; wmsetp("ezf",1) ; ; Station model. ; wmstnm(wks,40.5,-92.5,imdat) ; ; Draw the plot. ; frame(wks) end