* example GrADS plot script for model task 3
* ATM562
* http://www.atmos.albany.edu/facstaff/rfovell/ATM562/plot_vert_profile.gs
* version of 10/9/2018

* set background color white and clear screen
'set display color white'
'clear'

* set map projection off
'set mproj off'


* this formats the virtual page
'set vpage 0 8.5 0.5 8.5'

* smooth contours. enabled until switched off.
'set csmooth on'

* set location at domain center
'set lon 0'

* set range for horizontal plot
'set vrange -5 5'
'set grads off'

* plot temperature perturbation profile. color it blue
'set ccolor 4'
'd thpr'

* plot pressure perturbation profile (in mb). color ir red
'set ccolor 2'
'd pprmb'

* titles and labels
'draw title temperature and pressure pert profiles'
'draw xlab perturbation (K or mb)'
'draw ylab z (km)'

* ----------- make a PNG plot ---------------------------------
'gxprint init_tp_pp_profile.png'