arw graphics bruyere
TRANSCRIPT
-
7/26/2019 ARW Graphics Bruyere
1/78
WRF ARW Graphical ToolsWRF ARW Graphical ToolsWRF ARW Graphical Tools
Cindy Bruyre
-
7/26/2019 ARW Graphics Bruyere
2/78
outlineoutlineoutline
netCDF Data
NCL
RIP4
ARWpost ( o l d W RF- t o - Gr A DS & W RF- t o - V I S5 D )
-
7/26/2019 ARW Graphics Bruyere
3/78
netCDF datanetCDFnetCDF datadata
netCDF stands for netw o r k Co m m o n Da t aFo r m
netCDF is one of the current supporteddata formats chosen for WRF I/O API
WRF I/O supports netCDF (not CF compliant)/binary/GRIB/HDF
Most support graphical packages currentlyonly support netCDF file format
http://www.unidata.ucar.edu/
Documentation available at this site
http://www.unidata.ucar.edu/http://www.unidata.ucar.edu/ -
7/26/2019 ARW Graphics Bruyere
4/78
available graphicsavailable graphicsavailable graphics
NCL
RIP4
GrADS
Vis5D
IDL
Matlab
GMT
FX-Net Interactive
Visualization
Environment (I V E)
General MeteorologicalPackage /NationalAdvanced Weather
Interactive Processing
System (GEMPAK /NAWIPS, e.g.
garp andnmap)
Integrated Data
Viewer ( I D V )
Open-source DataExplorer ( O p e n D X )
-
7/26/2019 ARW Graphics Bruyere
5/78
NCL
NCL
NCL
Gr AD
S
Gr AD
S
Gr ADS
RIP
4
RIP
4
RIP
4
Vis5D
Vis5D
Vis5
D
?????
?
?????
?
?????
?
-
7/26/2019 ARW Graphics Bruyere
6/78
NCL RIP4 ARWpost(GrADS / Vis5D)
Directly ingest WRF
model output
c o n v e r t e r
/()c o n v e r t e r
Vertical Coordinate n / p / h n / p / h n / p / h
Model input & output i / o / w p s i / o / w p s i / o / w p s
Ideal 3 D / 2 D
3 D / 2 D
3 D / 2 D / 1 D
3 D / 2 D
WRF-Var (double precision)
Input format n e t CD F n e t CD F n e t CD F / GR I B 1
RIP4 easy, but need to know availablefields
GrADS - take some time to learn, but one
can get up to speed fast
Publish Quality Go o d Go o d Fa i r D i f f i c u l t
Easy of use NCL - need to do some scripting beyondstandard plots
ViS5D - need a large computer to run, asit is memory intensive
-
7/26/2019 ARW Graphics Bruyere
7/78
NCL RIP4 ARWpost(GrADS / Vis5D)
Software required
(All binaries are free)
H i g h
l e v e lN CL
NCAR
G r a p h i c s Gr AD S V i s 5 D
Diagnostics Som e > 1 0 0 Som e
Adding diagnostics F o r t r a n+ N CL
s c r i p t i n g
F o r t r a n Fo r t r a n + Gr A D S s c r i p t i n g
Intermediate files N O LOTS On e l a r g e f i l e
Interactive?
Can also be used with GRI B 1 + ASCI I
- - -GR I B 1 +
S t a t i o n - - -
Special Overlays S t a t i o n S t a t i o n 2 D + 3 D
Multiple Input Files
3D Capability
Cross-sections
-
7/26/2019 ARW Graphics Bruyere
8/78
NCL RIP4 ARWpost
(GrADS / Vis5D)
Skew-Ts
Plot sub-domains(zooming)
Animations?
Multiple images perpage?
Plotting format X - Y X - Y P r o j e c t e d 3 D
Create Trajectory data
Input to Vis5D
Output formats X1 1 , p s ,p d f , c gm
cg mX 1 1 , g i f ,
p n gX 1 1
-
7/26/2019 ARW Graphics Bruyere
9/78
examplesexamplesexamples
-
7/26/2019 ARW Graphics Bruyere
10/78
NCL
GrADS
-
7/26/2019 ARW Graphics Bruyere
11/78
NCLRIP4
-
7/26/2019 ARW Graphics Bruyere
12/78
NCL
GrADS
RIP4
-
7/26/2019 ARW Graphics Bruyere
13/78
RIP4
NCL
GrADS
-
7/26/2019 ARW Graphics Bruyere
14/78
GrADS
RIP4
NCL
-
7/26/2019 ARW Graphics Bruyere
15/78
GrADS
NCL
RIP4
-
7/26/2019 ARW Graphics Bruyere
16/78
Vis5D
NCL
-
7/26/2019 ARW Graphics Bruyere
17/78
NCL
-
7/26/2019 ARW Graphics Bruyere
18/78
NCLNCLNCL
NCL stands for NCAR Com m a n d L an g u a g e http://www.ncl.ucar.edu/
http://www.ncl.ucar.edu/http://www.ncl.ucar.edu/ -
7/26/2019 ARW Graphics Bruyere
19/78
downloading NCLdownloading NCLdownloading NCL
Go to:
http://www.ncl.ucar.edu/Download/
Read and agree to GPL licenseFill out short registration formDownload binaries
Set NCARG_ROOT environment variable: setenv NCARG_ROOT /usr/local/ncarg
Re com m e n d e d t o i n s t a l l N CAR Gr a p h i cs f i r s t a n d
t h e n NCL o n t o p o f N CAR Gr a p h i cs (most people dont- setenv NCARG_ROOT /usr/local/ncl)
W RF_ NCL ( Se p 2 0 0 5 v e r s i o n ) n e e d s NCL v e r s io n
4 . 2 . 0 . a 0 3 2 , o r h i g h e r
http://www.ncl.ucar.edu/Download/http://www.ncl.ucar.edu/Download/ -
7/26/2019 ARW Graphics Bruyere
20/78
WRF_NCLWRF_NCLWRF_NCL
Current release - Version 2.0
Release date September 2005
Re-write of functions and procedures
Cannot combine code from version 1 andversion 2
due to the re-write of WRF functions andprocedures
and upgrades made to newer versions of NCL
-
7/26/2019 ARW Graphics Bruyere
21/78
download NCL for WRFdownload NCL for WRFdownload NCL for WRF
From wrf-model.org web site(w r f _ n c l_ t a r .g z ):
http://www.mmm.ucar.edu/wrf/users/download/get_source.html
How to and examples:
http://www.mmm.ucar.edu/wrf/users/graphics/WRF_NCL/NCL.htm
http://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/graphics/WRF_NCL/NCL.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/WRF_NCL/NCL.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/WRF_NCL/NCL.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/WRF_NCL/NCL.htmhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.html -
7/26/2019 ARW Graphics Bruyere
22/78
README_FIRSTREADME_NCL
DOThluresfile
gsn_code.nclSkewTFunc.ncl
WRFOptions.nclWRFPlot.nclWRFUserARW.ncl
wrf_user_fortran_util_0.f
NCL for WRFNCL for WRFNCL for WRF
make_ncl_fortranmake_ncl_fortran.alphamake_ncl_fortran.ibm
make_ncl_fortran.ibm64make_ncl_fortran.linuxmake_ncl_fortran.sun
wrf_bwave.nclwrf_hill2d.nclwrf_grav2d.nclwrf_qss.ncl
wrf_squall_2d_x.nclwrf_squall_2d_y.nclwrf_real_input.nclwrf_real.nclwrf_cloud.ncl
-
7/26/2019 ARW Graphics Bruyere
23/78
NCL functions and proceduresNCL functions and proceduresNCL functions and procedures
gsn_code.nclHigh level NCL scripts
SkewTFunc.ncl
Skew-T plots
WRFOptions.nclBasic plotting options, like headers/footers
WRFPlot.nclmaps, contour plots, shaded plots, vectors
WRFUserARW.nclObtain native model variables
Basic diagnosticsInterface to FORTRAN code (wrf_fortran_user_util_0.f)
Must be loaded at top of NCL script, i.e.load WRFPlot.ncl
-
7/26/2019 ARW Graphics Bruyere
24/78
README_FIRSTREADME_NCL
DOThluresfile
gsn_code.nclSkewTFunc.ncl
WRFOptions.nclWRFPlot.nclWRFUserARW.ncl
wrf_user_fortran_util_0.f
NCL for WRFNCL for WRFNCL for WRF
make_ncl_fortranmake_ncl_fortran.alphamake_ncl_fortran.ibm
make_ncl_fortran.ibm64make_ncl_fortran.linuxmake_ncl_fortran.sun
wrf_bwave.nclwrf_hill2d.nclwrf_grav2d.nclwrf_qss.ncl
wrf_squall_2d_x.nclwrf_squall_2d_y.nclwrf_real_input.nclwrf_real.nclwrf_cloud.ncl
-
7/26/2019 ARW Graphics Bruyere
25/78
.hluresfile..hluresfilehluresfile
Required by NCL
Must be in your ~/ directory
Control color table
font
white/black background size of plot
control characters
h t t p : / / w w w .n c l .u c a r .e d u / D o cu m e n t / Gr a p h i c s / h l u r e s .s h t m l
cp DOThluresfile ~/.hluresfile NBNB
http://www.ncl.ucar.edu/Document/%0BGraphics/hlures.shtmlhttp://www.ncl.ucar.edu/Document/%0BGraphics/hlures.shtmlhttp://www.ncl.ucar.edu/Document/%0BGraphics/hlures.shtmlhttp://www.ncl.ucar.edu/Document/%0BGraphics/hlures.shtml -
7/26/2019 ARW Graphics Bruyere
26/78
README_FIRSTREADME_NCL
DOThluresfile
gsn_code.nclSkewTFunc.ncl
WRFOptions.nclWRFPlot.nclWRFUserARW.ncl
wrf_user_fortran_util_0.f
NCL for WRFNCL for WRFNCL for WRF
make_ncl_fortranmake_ncl_fortran.alphamake_ncl_fortran.ibm
make_ncl_fortran.ibm64make_ncl_fortran.linuxmake_ncl_fortran.sun
wrf_bwave.nclwrf_hill2d.nclwrf_grav2d.nclwrf_qss.ncl
wrf_squall_2d_x.nclwrf_squall_2d_y.nclwrf_real_input.nclwrf_real.ncl
wrf_cloud.ncl
-
7/26/2019 ARW Graphics Bruyere
27/78
wrf_user_fortran_util_0.fwrf_user_fortran_util_0.f
Contains FORTRAN routines for diagnostic
calculations
Code must be in F77
Code must contain special NCL START and ENDmarkers
Must create a shared object library from theseroutines in order for NCL to recognize the code.
C NCLFORTSTARTb ti t tk ( tk th t )
-
7/26/2019 ARW Graphics Bruyere
28/78
wrf_user_fortran_util_0.fwrf_user_fortran_util_0.fsubroutine compute_tk ( tk, pressure, theta, nx, ny, nz )implicit noneinteger nx,ny,nz
real pireal pressure(nx,ny,nz)real theta(nx,ny,nz)real tk(nx,ny,nz)
C NCLEND
integer i,j,kreal p1000mb, r_d, cp
parameter (p1000=100000., r_d=287., cp=7.*r_d/2.)
do k=1,nzdo j=1,ny
do i=1,nxpi=(pressure(i,j,k)/p1000)**(r_d/cp)tk(i,j,k) = pi*theta(i,j,k)enddoenddo
enddo
end
C NCLFORTSTARTb ti t tk ( tk th t )
-
7/26/2019 ARW Graphics Bruyere
29/78
wrf_user_fortran_util_0.fwrf_user_fortran_util_0.fsubroutine compute_tk ( tk, pressure, theta, nx, ny, nz )implicit noneinteger nx,ny,nz
real pireal pressure(nx,ny,nz)real theta(nx,ny,nz)real tk(nx,ny,nz)
C NCLEND
integer i,j,kreal p1000mb, r_d, cp
parameter (p1000=100000., r_d=287., cp=7.*r_d/2.)
do k=1,nzdo j=1,ny
do i=1,nxpi=(pressure(i,j,k)/p1000)**(r_d/cp)tk(i,j,k) = pi*theta(i,j,k)enddoenddo
enddo
end
-
7/26/2019 ARW Graphics Bruyere
30/78
wrf_user_fortran_util_0.sowrf_user_fortran_util_0.so
WRAPIT wrf_user_fortran_util_0.f
setenv NCARG_ROOT/usr/local/ncl
/usr/local/ncl/bin/WRAPIT
wrf_user_fortran_util_0.f
USE in NCL script
wrf_user_fortran_util_0 ::compute_tk(tk, pressure, theta, nx, ny, nz )
(WRFUserARW.ncl)
-
7/26/2019 ARW Graphics Bruyere
31/78
README_FIRSTREADME_NCL
DOThluresfile
gsn_code.nclSkewTFunc.ncl
WRFOptions.nclWRFPlot.nclWRFUserARW.ncl
wrf_user_fortran_util_0.f
NCL for WRFNCL for WRFNCL for WRF
make_ncl_fortranmake_ncl_fortran.alphamake_ncl_fortran.ibm
make_ncl_fortran.ibm64make_ncl_fortran.linuxmake_ncl_fortran.sun
wrf_bwave.nclwrf_hill2d.nclwrf_grav2d.nclwrf_qss.ncl
wrf_squall_2d_x.nclwrf_squall_2d_y.nclwrf_real_input.nclwrf_real.ncl
wrf_cloud.ncl
-
7/26/2019 ARW Graphics Bruyere
32/78
wrf_user_fortran_util_0.sowrf_user_fortran_util_0.so
IF WRAPIT DOES NOT WORK
Try one of the makefiles to create sharedobject library
make_ncl_fortran wrf_user_fortran_util_0
Potential Problem
wrapit77 is part of NCAR Graphics, and thepath to this function needs to be correctly
specified (.cshrc file)
-
7/26/2019 ARW Graphics Bruyere
33/78
a scripta scripta script
load "WRFOptions.ncl ; set up opt i ons resload "gsn_code.ncl"
load "WRFPlot.ncl"
load "WRFUserARW.ncl"
begin
; Open a f i l ea = addfile("wrfout_d01_2000-01-24_12:00:00.nc","r")
type = "pdf" ; Wi l l cr eat e pdf out put f i l es; open wor kst at i on + set out putwks = gsn_open_wks(type,"wrfout")
res@MainTitle = "ARW Real Data" ; Set pl ot t i t l e
times = wrf_user_list_times(a) ; get t i mes i n t he f i l entimes = dimsizes(times) ; number of t i mes i n t he f i l e
Not part of file name, but
needs to be added insidethe ncl script
-
7/26/2019 ARW Graphics Bruyere
34/78
a scripta scripta scriptdo it = 0,ntimes-1 ; l oop t hr ough f i l e f or al l t i mes
res@TimeLabel = times(it) ; Set Val i d t i me i nf omap = wrf_map(wks,a,res) ; Cr eat e a map background
ter = wrf_user_getvar(a,"HGT",time) ; get t er r ai nt2 = wrf_user_getvar(a,"T2",time) ; get T2
; Set some opt i ons f or t er r ai n pl ot
opts_ter = resopts_ter@cnFillOn = True
contour_ter = wrf_contour(a,wks,ter,opts_ter)
; Set some opt i ons f or t 2 pl otopts_t2 = res
contour_t2 = wrf_contour(a,wks,t2,opts_t2)
wrf_map_overlay(wks,map,(/contour_ter,contour_t2/),False)
end do
endOne image with map, shaded terrain and contour t2
-
7/26/2019 ARW Graphics Bruyere
35/78
a script - a l t e r n a t i v e a scripta script -- a l t e r n a t i v e a l t e r n a t i v e do it = 0,ntimes-1 ; l oop t hr ough f i l e f or al l t i mes
res@TimeLabel = times(it) ; Set Val i d t i me i nf omap = wrf_map(wks,a,res) ; Cr eat e a map background
ter = wrf_user_getvar(a,"HGT",time) ; get t er r ai nt2 = wrf_user_getvar(a,"T2",time) ; get T2
opts_ter = res
opts_ter@cnFillOn = Truecontour_ter = wrf_contour(a,wks,ter,opts_ter)
opts_t2 = res
contour_t2 = wrf_contour(a,wks,t2,opts_t2)
wrf_map_overlay(wks,map,(/contour_ter/),False)
wrf_map_overlay(wks,map,(/contour_t2/),False)
wrf_map_overlay(wks,map,(/contour_ter,contour_t2/),False)
end do
end
-
7/26/2019 ARW Graphics Bruyere
36/78
a script i d e a l ca s e a scripta script i d e a l c a s e i d e a l ca s e do it = 0,ntimes-1 ; l oop t hr ough f i l e f or al l t i mes
res@TimeLabel = times(it) ; Set Val i d t i me i nf omap = wrf_map(wks,a,res) ; Cr eat e a map background
ter = wrf_user_getvar(a,"HGT",time) ; get t er r ai nt2 = wrf_user_getvar(a,"T2",time) ; get T2
opts_ter = res
opts_ter@cnFillOn = Truecontour_ter = wrf_contour(a,wks,ter,opts_ter)
opts_t2 = res
contour_t2 = wrf_contour(a,wks,t2,opts_t2)
wrf_map_overlay(wks,map,(/contour_ter,contour_t2/),False)
wrf_overlay(wks,/contour_ter,contour_t2/),False)
end do
end
-
7/26/2019 ARW Graphics Bruyere
37/78
tc_plane = 1.8*tc2+32.
opts_tc = res
opts_tc@FieldTitle = "Surface Temperature"
opts_tc@UnitLabel = "F"
opts_tc@ContourParameters = (/-20.,90.,5./)
opts_tc@cnFillOn = Truecontour_tc=wrf_contour(a,wks,tc_plane,opts_tc)
opts_psl = res
opts_psl@ContourParameters = (/900.,1100.,4./)
opts_psl@cnLineColor = "NavyBlue"
opts_psl@cnHighLabelsOn = Trueopts_psl@cnLowLabelsOn = True
contour_psl=wrf_contour(a,wks,slvl,opts_psl)
u_plane = u10*1.94386
v_plane = v10*1.94386
opts_vct = res
opts_vct@FieldTitle = "Winds"
opts_vct@UnitLabel = "kts"
opts_vct@NumVectors = 47
vector=wrf_vector(a,wks,u_plane,v_plane,opts_vct)
wrf_map_overlay(wks,map,(/contour_tc,contour_psl,vector/),False)
exampleexampleexample
-
7/26/2019 ARW Graphics Bruyere
38/78
-
7/26/2019 ARW Graphics Bruyere
39/78
how to run NCL?how to run NCL?how to run NCL?
cp DOThluresfile ~/.hluresfile( c h a n g e i f n e e d e d )
Make external functionW RA PI T w r f _ u s e r _ f o r t r a n _ u t i l _ 0 . f
Edit WRFOptions.ncl( s e t h e a d e r a n d f o o t e r o p t i o n s )
Edit NCL script( s e t i n p u t a n d o u t p u t f i l e s ; s e t M a i n T i t l e )
Run NCL script
ncl w r f _ r e a l .n c l
-
7/26/2019 ARW Graphics Bruyere
40/78
News FlashNews FlashNews Flash
N Fl h
-
7/26/2019 ARW Graphics Bruyere
41/78
README_FIRSTREADME_NCL
DOThluresfile
gsn_code.nclSkewTFunc.ncl
WRFOptions.nclWRFPlot.nclWRFUserARW.ncl
wrf_user_fortran_util_0.f
News FlashNews FlashNews Flash
make_ncl_fortranmake_ncl_fortran.alphamake_ncl_fortran.ibm
make_ncl_fortran.ibm64make_ncl_fortran.linuxmake_ncl_fortran.sun
wrf_bwave.nclwrf_hill2d.nclwrf_grav2d.nclwrf_qss.ncl
wrf_squall_2d_x.nclwrf_squall_2d_y.nclwrf_real_input.nclwrf_real.ncl
wrf_cloud.ncl
RIP4RIP4
-
7/26/2019 ARW Graphics Bruyere
42/78
RIP4RIP4RIP4
Re a d / I n t e r p o l a t e / P l o t ( V e r s i o n 4 )
Originally written for MM5 input and output,and recently rewritten and generalized toinclude WRF
Written by Mark StoelingaUW / NCAR
lll
-
7/26/2019 ARW Graphics Bruyere
43/78
generalgeneralgeneral
Requires NCAR Graphics low-level routines
http://ngwww.ucar.edu
Documentation:
In program tar file under the Doc/ directory http://www.mmm.ucar.edu/wrf/users/
docs/ripug.htm
d l d RIP4 f WRFd l d RIP4 f WRFd l d RIP4 f WRF
http://ngwww.ucar.edu/http://www.mmm.ucar.edu/wrf/users/docs/ripug.htmhttp://www.mmm.ucar.edu/wrf/users/docs/ripug.htmhttp://www.mmm.ucar.edu/wrf/users/docs/ripug.htmhttp://www.mmm.ucar.edu/wrf/users/docs/ripug.htmhttp://ngwww.ucar.edu/ -
7/26/2019 ARW Graphics Bruyere
44/78
download RIP4 for WRFdownload RIP4 for WRFdownload RIP4 for WRF
From wrf-model.org web site(r i p 4 . t a r . g z ):
http://www.mmm.ucar.edu/wrf/users/download/get_source.html
How to and examples:
http://www.mmm.ucar.edu/wrf/users/graphics/RIP4/RIP4.htm
RIP4 tRIP4 tRIP4 t
http://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/graphics/RIP4/RIP4.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/RIP4/RIP4.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/RIP4/RIP4.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/RIP4/RIP4.htmhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.html -
7/26/2019 ARW Graphics Bruyere
45/78
RIP4 on your computerRIP4 on your computerRIP4 on your computer
set RIP_ROOT environment variable
setenv RIP_ROOT /usr/$USER/RIP4
Edit Ma k e f i l e to define paths to netCDFlibrary and include file on your computer:
NETCDFLIB and NETCDFINC
make (itll make suggestions)
make linux (example)
RIP4 has 2 parts (RIPDP and RIP)
ripdp_wrfripdp_mm5ripdp
i d fi d fripdp wrf
-
7/26/2019 ARW Graphics Bruyere
46/78
ripdp_wrfripdp_wrfripdp_wrf
ripdp_wrf is R I P D a t a P r e p a r a t i o n f o r W RF
RIP does not read WRF data directly
r i p d p _ w r f converts WRF netCDF data into RIPinput format (format described in thedocument)
RIP puts each variable at each time into aseparate file LOTS of files
mkdir storm_case
i i d frunningrunning ripdp wrfripdp wrf
-
7/26/2019 ARW Graphics Bruyere
47/78
running ripdp_wrfrunningrunning ripdp_wrfripdp_wrf
ripdp_wrf [n namelist-file] \ [basic/all] \
Optionalripdp_sample.in
Example
ripdp_wrf storm_case/test basic wrfout
use directory as part of themodel_data_name
ripdp wrf namelistripdp wrfripdp wrf namelistnamelist
-
7/26/2019 ARW Graphics Bruyere
48/78
ripdp_wrf namelistripdp_wrfripdp_wrfnamelistnamelist
&userin
ptimes=0,-72,1,ptimeunites=h,tacc=90.,
discard=LANDMASK,H2SO4,
iexpandedout=0&end
Use namelist to add control
ptimes times for ripdp_wrfto process
0,1,2
0,-72,1
0, 3,-24,3, 48
discard fields if all is selected on thecommand line
ripriprip
-
7/26/2019 ARW Graphics Bruyere
49/78
ripriprip
read the output generated by r i p d p _ w r f
read User Input File (UIF) ( r i p _ s a m p l e . i n )
First section is a list of general parameters
(namelist format) Second section is a series of plots in the Plot
Specification Table (PST)
generate meta file
running riprunning riprunning rip
-
7/26/2019 ARW Graphics Bruyere
50/78
running riprunning riprunning rip
Edit the User Input File (UIF)
setenv NCARG_ROOT /usr/local/ncarg
setenv RIP_ROOT y o u r - r i p - d i r e c t o r y
running riprunning riprunning rip
-
7/26/2019 ARW Graphics Bruyere
51/78
running riprunning riprunning rip
rip [-f] model-data-set-name \rip-execution-name
created byripdp_wrf
User Input File (UIF)
Example
rip -f storm_case/test rip_sample.in
use directory as part of themodel_data_set_name
output ; metacode
rip_sample.out
rip_sample.ncgm
rip UIFrip UIFrip UIF
-
7/26/2019 ARW Graphics Bruyere
52/78
rip UIFrip UIFrip UIF
userin
.
end
trajcalc
.
end
========================================
---------- Plot Specification Table -----------
========================================
feld=
feld=
========================================
feld=
feld=
========================================
Plot
Specification
Table (PST)
Frame specification
group (FSG)
Plot specification line (PSL)
Namelist controlling general parameters
Namelist for trajectory calculations
Only used if itrajcalc=1 in userin namelist
rip namelist - userinriprip namelistnamelist -- userinuserin
-
7/26/2019 ARW Graphics Bruyere
53/78
rip namelist - userinriprip namelistnamelist -- userinuserin
Use namelist to control
processing times, intervals
title information text quality on a plot
whether to do time series, trajectory, or to
write output for Vis5D
Full explanation for namelist variables isavailable in the user document
rip namelist - userinriprip namelistnamelist -- userinuserin
-
7/26/2019 ARW Graphics Bruyere
54/78
rip namelist - userinriprip namelistnamelist - userinuserin
idotitle first part of first title line
titlecolor color of title lines
ptimes, ptimeunits times to process
tacc tolerance for processing data
timezone display of local time
iusedaylightrule 1 applied, 0 not applied iinittime plotting of initial time
ivalidtime plotting of valid time
inearsth plot times as 2 / 4 digits flmin, frmax, fbmin, ftmax frame size
ntextq text quality
rip namelist - userinriprip namelistnamelist -- userinuserin
-
7/26/2019 ARW Graphics Bruyere
55/78
rip namelist - userinriprip namelistnamelist userinuserin
ntextcd text font
fcoffset 12 means hour 12 of the MM5
forecast is considered hour 0 by you idotser generate time series output
idescriptive more descriptive titles
icgmsplit split metacode into several files maxfld reserve memory for RIP (10-15)
itrajcalc 0, 1 ONLY when doing trajectory
calculations ( u s e a l s o n a m e l i s t t r a j c a l c ) imakev5d 0, 1 generate Vis5D data
exampleexampleexample
-
7/26/2019 ARW Graphics Bruyere
56/78
exampleexampleexample
=====================================
feld=tmc; ptyp=hc; vcor=s; levs=1fb;
cint=2; cmth=fill; >cosq=32,light.violet,-16,blue,>
0,yellow,16,orange,32,light.gray
feld=slp; ptyp=hc; cint=2; linw=2feld=uuu,vvv; ptyp=hv; vcmx=1;>
colr=white;intv=5
feld=map; ptyp=hbfeld=tic; ptyp=hb
=====================================
-
7/26/2019 ARW Graphics Bruyere
57/78
ARWpostARWpostARWpost
-
7/26/2019 ARW Graphics Bruyere
58/78
ARWpostARWpostARWpost
Converter, creates
GrADS ( . d a t ) and ( . c t l ) files, orVis5D input files
download ARWpostdownload ARWpostdownload ARWpost
-
7/26/2019 ARW Graphics Bruyere
59/78
download ARWpostdownload ARWpostdownload ARWpost
From wrf-model.org web site(ARW p o s t . t a r ):
http://www.mmm.ucar.edu/wrf/users/download/get_source.html
How to and examples:
http://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htm
generalgeneralgeneral
http://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htmhttp://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htmhttp://www.mmm.ucar.edu/wrf/users/download/get_source.htmlhttp://www.mmm.ucar.edu/wrf/users/download/get_source.html -
7/26/2019 ARW Graphics Bruyere
60/78
generalgeneralgeneral
MUST have WRFV2 compiled
IF Vis5D files are being created,
Vis5D libraries needed for compilation Vis5D is free and can be downloaded from:
http://www.ssec.wisc.edu/~billh/vis5d.html
codecodecode
http://www.ssec.wisc.edu/~billh/vis5d.htmlhttp://www.ssec.wisc.edu/~billh/vis5d.html -
7/26/2019 ARW Graphics Bruyere
61/78
codecode
READMEarch/clean*compile*configure*fields.plt
gribinfo.txtgribmap.txtmyLIST
namelist.ARWpostscripts/src/util/
compilationcompilationcompilation
-
7/26/2019 ARW Graphics Bruyere
62/78
compilationpp
./configure
Will use NETCDF in dir: /usr/local/netcdf-pgi- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -Pl ease sel ect f r om among t he f ol l owi ng suppor t edpl at f or ms.1. PC Li nux i 486 i 586 i 686, PGI compi l er (novis5d)2. PC Li nux i 486 i 586 i 686, PGI compi l er (vis5d)3. PC Li nux i 486 i 586 i 686, I nt el compi l er (novis5d)
4. PC Li nux i 486 i 586 i 686, I nt el compi l er(vis5d)
Ent er sel ect i on [ 1- 4] :
compilationcompilationcompilation
-
7/26/2019 ARW Graphics Bruyere
63/78
ppp
configure.arwp, will be created
If your WRFV2 code is not compiled under../WRFV2, edit configure.arwp, and setWRF_DIR to the correct location of yourWRFV2 code
./compile
This will create ARWpost.exe
namelist.ARWpostnamelist.ARWpostnamelist.ARWpost
-
7/26/2019 ARW Graphics Bruyere
64/78
ppp
set times to be processed (&datetime)
set input and output file names andvariables to be processed (&io)
set levels to interpolate too (&interp)
namelist.ARWpost ( & i o ) namelist.ARWpostnamelist.ARWpost ( &( &i oi o))
-
7/26/2019 ARW Graphics Bruyere
65/78
p ( )pp (( ))
io_form_input 2=netCDF, 5=GRIB1
input_root_name Path and root name of files to use as input.Do not use wild characters in the input_root_name.
output_root_name Output root name.
output_root_name.dat & output_root_name.ctl, OR
output_root_name.v5d
output_type Options are 'grads' (default) or 'v5d'
mercator_defs Set to true if mercator plots are distorted
namelist.ARWpost ( & i o ) namelist.ARWpostnamelist.ARWpost ( &( &i oi o))
-
7/26/2019 ARW Graphics Bruyere
66/78
p ( )p ( )
plot Which fields to process. Options are: all, basic, list, file, basic_file,
basic_list, list_file, all_file, all_list, basic_list_file, all_list_file
Order has no effect, i.e., all_list and list_all are similar.
"list" - a list of variables must be supplied under "fields"
"file" - a list of variables must be added to a file, and the
filename supplied under "fields_file
fields Fields to plot.
Only used is list was used in the plot variable.
fields_file File name that contains list of fields to plot.
Only used is file was used in the plot variable.
A v a i l a b l e d i a g n o s t i c s : height ,theta ,tc, tk, td, rh, umet, vmet,
pressure, dbz, max_dbz, u10m, v10m, slp, mcape, mcin, lcl, lfc,cape, cin
namelist.ARWpostnamelist.ARWpostnamelist.ARWpost
-
7/26/2019 ARW Graphics Bruyere
67/78
p
interp_method 0 = sigma levels,
-1 = code defined "nice" height levels,
1 = user defined height or pressure levels
interp_levels Only used if interp_method=1
Supply levels to interpolate to, in hPa (pressure) or km (height)
Supply levels bottom to top
runrunrun
-
7/26/2019 ARW Graphics Bruyere
68/78
./ARWpost.exe
Will create either,
o u t p u t _ r o o t _ n a m e .dat & o u t p u t _ r o o t _ n a m e .ctl
OR
o u t p u t _ r o o t _ n a m e .v5d
GrADS Specific NotesGrADSGrADS Specific NotesSpecific Notes
-
7/26/2019 ARW Graphics Bruyere
69/78
To display images requires GrADS softwarefreely available from
http://grads.iges.org/grads/grads.html
Documentation:http://grads.iges.org/grads/gadoc/index.html
GrADS - projectionsGrADSGrADS -- projectionsprojections
http://grads.iges.org/grads/grads.htmlhttp://grads.iges.org/grads/gadoc/index.htmlhttp://grads.iges.org/grads/gadoc/index.htmlhttp://grads.iges.org/grads/grads.html -
7/26/2019 ARW Graphics Bruyere
70/78
GrADS conversion - questionGrADSGrADS conversionconversion -- questionquestion
-
7/26/2019 ARW Graphics Bruyere
71/78
Why is a converter needed if GrADS candisplay netCDF files?
Can only display model surface coordinates
Cannot interpolate to height or pressurelevels
All diagnostics must be added via GrADSscript files
GRIB1 model output can also be read directly byGrADS, but above issues are still valid
For GRIB1, there is also a stagger problem
GrADS conversion - questionGrADSGrADS conversionconversion -- questionquestion
-
7/26/2019 ARW Graphics Bruyere
72/78
Why is a converter needed if GrADS can displaynetCDF files?
Staggering v v v v
m m
v
mv
m
vm
v
v
mv
m
vm
v
m
v
mv
m
vm
v
u m u u u u
v
u m u u u uv
u m u u u u
vu m u u u u
v
staggeringstaggeringstaggering
-
7/26/2019 ARW Graphics Bruyere
73/78
shaded=T ; black=U ; red=V
staggeringstaggeringstaggering
-
7/26/2019 ARW Graphics Bruyere
74/78
Since GrADS version 1.9
a new gradsnc interface is availablecreated by developers for WRF
To USE
must create 4 .ctl files (M ; U ; V ; W) must open the all at once
Utility ARWpost/util/WRFnc2ctl.f
GrADS script exampleGrADSGrADS script examplescript example
-
7/26/2019 ARW Graphics Bruyere
75/78
'open real.ctl'set mpdset hires''set display color white''run rgbset.gs'
'define tf=1.8*tc + 32'set gxout shaded''set clevs -20 -15 -10 -5 0 5 10 15 20 25 30 35 40
45 50 55 60 65 70 75 80'
'set ccols 0 51 53 55 57 58 59 49 47 45 39 37 3634 22 23 24 25 27 29 29'
'set z 1''d tf
'run cbar.gs
'set gxout contour''set ccolor 1''set cint 4''d slvl'draw title Surface T(F,color), SLP(mb)'
http://www.mmm.ucar.edu/mm5/WRF_post/WRF2GrADS/real_surface/images/surface6_gif.jpg -
7/26/2019 ARW Graphics Bruyere
76/78
Vis5D Specific NotesVis5D Specific NotesVis5D Specific Notes
http://www.mmm.ucar.edu/mm5/WRF_post/WRF2GrADS/real_surface/images/surface6_gif.jpg -
7/26/2019 ARW Graphics Bruyere
77/78
Vis5D is a three-dimensional visualizationsoftware
Vis5D is free and can be downloaded from:http://www.ssec.wisc.edu/~billh/vis5d.html
Run
vis5d output_root_name.v5d
Graphical Interface
http://www.ssec.wisc.edu/~billh/vis5d.htmlhttp://www.ssec.wisc.edu/~billh/vis5d.htmlhttp://www.ssec.wisc.edu/~billh/vis5d.htmlhttp://www.ssec.wisc.edu/~billh/vis5d.html -
7/26/2019 ARW Graphics Bruyere
78/78
NCL
NCL
NCL
GrAD
S
GrAD
S
GrAD
S
RIP
4
RIP
4
RIP
4
Vis5D
Vis5D
Vis5D