Skip to content

Instantly share code, notes, and snippets.

@WxAnalyst
Created June 24, 2013 07:37
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save WxAnalyst/5848366 to your computer and use it in GitHub Desktop.
Save WxAnalyst/5848366 to your computer and use it in GitHub Desktop.
GRaDS script to run a time/height diagram from WRF output
**********************************************************
* Displays a time-height diagram for stations
* Author: John Hinsberger
*********************************************************
*Bald Mountain*
'open d01_template.ctl'
'set background 5'
'clear'
'set background 5'
'set grads off'
'set zlog on'
'set lev 750 300'
'set lat 43.6612'
'set lon -114.4021'
'set t 0 36'
# RH
'set gxout shaded'
'run colorRH_Time.gs'
'd rhprs'
# Winds
'set gxout barb'
'set ccolor 0'
'd ugrdprs;vgrdprs'
# temps
'set gxout contour'
'set cmax 20'
'set cmin -40'
'set cint 2'
'set cthick 8'
'set ccolor 2'
'set clopts 1 6 0.11'
'set clskip 2'
'd tmpprs-273'
# Omega
'set cint 0.5'
'set clopts 14 8 0.11'
'set cmax -0.5'
'set cmin -30'
'set rgb 40 130 0 0'
'set ccolor 40'
'd vvelprs'
'run cbarnskip 1 0.75 0'
'draw title \\Bald Mountain RH(shaded)/winds/Omega/Temp (C)'
'printim d01_PIH_timeheight_bald.png x1024 y768'
'clear'
*Soldier Mountain*
'set grads off'
'set lev 750 300'
'set lat 43.48455'
'set lon -114.90712'
'set t 0 36'
# RH
'set gxout shaded'
'run colorRH_Time.gs'
'd rhprs'
# Winds
'set gxout barb'
'set ccolor 0'
'd ugrdprs;vgrdprs'
# temps
'set gxout contour'
'set cmax 20'
'set cmin -40'
'set cint 2'
'set cthick 8'
'set ccolor 2'
'set clopts 1 6 0.11'
'd tmpprs-273'
# Omega
'set cint 0.5'
'set clopts 14 5 0.11'
'set cmax -0.5'
'set cmin -30'
'set ccolor 40'
'd vvelprs'
'run cbarnskip 1 0.75 0'
'draw title \\Soldier Mountain RH(shaded)/winds/Omega/Temp (C)'
'printim d01_PIH_timeheight_soldier.png x1024 y768'
'clear'
*Sawtooth Range*
'set grads off'
'set lev 750 300'
'set lat 44.1'
'set lon -115.0'
'set t 0 36'
# RH
'set gxout shaded'
'run colorRH_Time.gs'
'd rhprs'
# Winds
'set gxout barb'
'set ccolor 0'
'd ugrdprs;vgrdprs'
# temps
'set gxout contour'
'set cmax 20'
'set cmin -40'
'set cint 2'
'set cthick 8'
'set ccolor 2'
'set clopts 1 6 0.11'
'd tmpprs-273'
# Omega
'set cint 0.5'
'set clopts 14 5 0.11'
'set cmax -0.5'
'set cmin -30'
'set ccolor 40'
'd vvelprs'
'run cbarnskip 1 0.75 0'
'draw title \\Sawtooth Range RH(shaded)/winds/Omega/Temp (C)'
'printim d01_PIH_timeheight_sawtooth.png x1024 y768'
'clear'
*Titus Peak*
'set grads off'
'set lev 750 300'
'set lat 43.8366'
'set lon -114.7104'
'set t 0 36'
# RH
'set gxout shaded'
'run colorRH_Time.gs'
'd rhprs'
# Winds
'set gxout barb'
'set ccolor 0'
'd ugrdprs;vgrdprs'
# Temps
'set gxout contour'
'set cmax 20'
'set cmin -40'
'set cint 2'
'set cthick 8'
'set ccolor 2'
'set clopts 1 6 0.11'
'd tmpprs-273'
# Omega
'set cint 0.5'
'set clopts 14 5 0.11'
'set cmax -0.5'
'set cmin -30'
'set ccolor 40'
'd vvelprs'
'run cbarnskip 1 0.75 0'
'draw title \\Titus Peak RH(shaded)/winds/Omega/Temp (C)'
'printim d01_PIH_timeheight_titus.png x1024 y768'
'clear'
'quit'
return
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment