* * Display orography on vertical profile * * pressfc = surface pressure in Pa * lev is assumed to be in hPa * * display is assumed to be in "shaded" and vertical profile * 'define t1pm1=lev-pressfc/100' 'set gxout grfill' 'set clevs -1 1' 'set ccols 1 1' 'set ccols 0 0' 'd maskout(0,t1pm1)' 'set gxout shaded'