90 likes | 369 Views
RIP. A tool to create plots from WRF model output. RIP (Read, Interpolate, Plot). Developed over many years by Mark Stoelinga of the University of Washington
E N D
RIP A tool to create plots from WRF model output
RIP(Read, Interpolate, Plot) • Developed over many years by Mark Stoelinga of the University of Washington • Not a gui-driven visualization tool. To use rip, you edit a configuration file and run the rip executable, then view the static images created by rip • Allows user fine control over nearly all aspects of the plot; easy* to make publication-quality graphics (*once you get a feel for the arcane commands) • Not so great as a “quick-look” tool
RIP capabilities include… • Cross-sections on various types of quasi-horizontal surfaces • model sigma, pressure, height, isentropic, etc. • Vertical cross-sections along arbitrary lines • Vertical profiles and skew-T diagrams • Contours, vectors, barbs, color-fill, maps, labels, etc., and almost any combination thereof • Difference fields between two model jobs
For the colloquium case studies… • Input datasets for RIP have been created for you… • i.e., you can skip the “ripdp” step • Sub-sampled and windowed for easier plotting on the classroom machines • Selected forecast hours to keep disk usage under control • Sample configuration files have been created for you… • to give you a starting point for creating your own configuration files • Plot files of many useful fields have been created for you… • lots of images are available for you to peruse before you even run RIP • Images on the full model grid, not sub-sampled or windowed • Images at hourly time interval
To use RIP • 1) Prepare the data • Largely done for you • 2) Prepare the instruction file • This is the hard part. We offer some samples to get you started • 3) Run the program • 4) View the results • idt command
&userin title="RIP_example2" idotitle=1,titlecolor='def.foreground', ptimes=18,24,30, ptimeunits='h',tacc=120,timezone=-7,iusdaylightrule=1, iinittime=1,ifcsttime=1,ivalidtime=1,inearesth=0, flmin=.09, frmax=.92, fbmin=.10, ftmax=.85, ntextq=0,ntextcd=0,fcoffset=0.0,idotser=0, idescriptive=1,icgmsplit=0,maxfld=10,itrajcalc=0,imakev5d=0 &end =========================================================================== ---------------------- Plot Specification Table --------------------- =========================================================================== # # Vertical cross sections # feld=the; ptyp=vc; vcor=z; vwin=0,5; crsa=25,20; crsb=120,180; cint=1.0 feld=uuu,vvv,www; ptyp=vv; feld=tic; ptyp=vb; =========================================================================== feld=the; ptyp=vc; vcor=z; vwin=0,5; crsa=25,20; crsb=120,180; cint=2.0 feld=uuu,vvv; ptyp=vw; vcmx=-1; feld=tic; ptyp=vb; =========================================================================== # # Color Fill (smooth) # feld=wsp; ptyp=hc; vcor=s; levs=1fb; cmth=fill; cint=0.5; cosq=0,white,5,yellow,10,green,15,blue,20,red; feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== title and ptimes may be the only things you change in the namelist section feld=<something>; ptyp=<something> will begin almost every line “==============“ separates one frame from the next
… =========================================================================== # # Color Fill (discontinuous thresholds) # feld=wsp; ptyp=hc; cmth=fill; cosq=0,white,5,white,5,yellow,10,yellow,10,green, > 15,green,15,blue,20,blue,20,red; feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== # # Difference fields between sensitivity studies # feld=tmc; ptyp=hc; vcor=s; levs=1fb; diff=myj/myj; cmth=fill; cint=0.25; cbeg=-8; > cosq=-8,cyan,-4,blue,-0.5,white,0.5,white,4,red,8,green feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== # # Wind vectors # feld=uuu,vvv; ptyp=hv; vcor=p; levs=850, 700,500; intv=8; vcmx=10; feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== # # Wind barbs # feld=uuu,vvv; ptyp=hv; vcor=z; levs=1.5; intv=10; vcmx=-1; fulb=10kts feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== “>” at the end of a line is a continuation mark vcor specifies vertical coordinate to use levs specifies levels to plot
… =========================================================================== # # Skew-T diagrams # feld=tic; ptyp=sb; sloc=39.31lat,-94.72lon; hodo; sndg; feld=tmc; ptyp=sc; linw=2; colr=red; feld=tdp; ptyp=sc; linw=2; colr=blue; feld=uuu,vvv; ptyp=sv; hodo; sndg =========================================================================== # # Vertical profiles of other variables: # feld=www; ptyp=pc; sloc=39.31lat,-94.72lon; =========================================================================== # # Lines, boxes, bullets, annotations # feld=line; ptyp=hb; crsa=120,10; crsb=39.31lat,-94.72lon; colr=red; linw=3 feld=bullet; ptyp=hb; crsa=KMCI; colr=blue; linw=3 feld=box; ptyp=hb; crsa=120,40; crsb=160,80; colr=dark.green; linw=3 feld=bullet; ptyp=hb; crsa=100,180; titl=Put_your_label_here; feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== # # Windowing # feld=uuu,vvv; ptyp=hv; vcor=s; levs=1fb; intv=1; vcmx=-1; fulb=10kts; > xwin=120,160; ywin=40,80; feld=map; ptyp=hb; ouds=solid; feld=tic; ptyp=hb; =========================================================================== Coordinates may be specified in terms of lat/lon, grid x/y, and station ID Semicolons separate keyword settings Commas separate multiple values for a single keyword setting Whitespace is ignored, but blank lines are a no-no # anywhere in a line comments out the line
RIP suggestions • Bookmark the online documentation and refer to it often • Beware of commas and semicolons, and avoid blank lines • Don’t plot too much at once • It can take a long time, and you just end up waiting • Don’t spend too much time getting a plot to look “just right” • Some keyword settings are remembered from frame to frame. Some are forgotten.