Created
June 24, 2013 07:37
-
-
Save WxAnalyst/5848366 to your computer and use it in GitHub Desktop.
GRaDS script to run a time/height diagram from WRF output
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
********************************************************** | |
* 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