meteo_symbol_ncl

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;
;  Draw station model data.
;
begin

 print("x = " + x)
 print("filename = " + filename)
;
;  Define a color map and open a workstation.
;
 cmap = (/                      \
            (/ 1., 1., 1. /),    \  ; color index 0 - white
            (/ 0., 0., 0. /)     \  ; color index 1 - black
         /)

 wks = gsn_open_wks("ps", filename)
 gsn_define_colormap(wks,cmap)
 res                   = True
 res@gsnMaximize       = True

 data_test   = x

;NP-39 2011-09-29-00
;UFTA 29001 99840 71499 41997 72002 11092 21104       49989 54000       77172 87050 333 21093=
;                       41997 72002 11092 21104 30000 49989 54000 60000 77172 87050
;                       41997 72002 11092 21104 30000 49989 54000 60000 72002 87050
; imdat                 11212 83320 10011 20000 30000 40147 52028 60111 77060 86792
 np39dat = "41997720021109221104300004998954000600007717287050"

;Resolution
 wmsetp("wbs",0.017)

;Correct visibility
 wmsetp("VVC",1)

;Metric system
 wmsetp("UNT",1)
 wmstnm(wks,0.5,0.5,data_test)

 frame(wks)
end