<div>Grads users,</div><div><br></div>I am trying to plot topography in a cross section with pressure as the vertical coordinate.  I have read a lot about this in the user archives and even found a script (orog.gs) that was written to do it. <div> <div>This script uses the variable psfc (press. at the surface).  Unfortunately the variables that could be used in the script that are available from my WRF-ARW model output (terrain height, psfc, slp) are only on the lowest level.  Because of this when try:</div><div><br></div><div>set lat 21.5</div><div>set lev 1000 700</div><div>d psfc</div><div><br></div><div>The result is nothing being visibly plotted (because psfc is only on the lowest level).  Now if I set the vertical axis to a single pressure value like so:</div><div><br></div><div><div>set lat 21.5</div><div>set lev 1000</div><div>d psfc</div><div><br></div><div>The result is an x-y graph of the surface pressure with longitude on the x-axis and pressure on the y-axis.  The result is the similar when I try to plot terrain height and slp as well.  </div><div><br></div><div>Because psfc is only on the lowest level I cannot do the  var = lev-psfc/100  command.  What I really need is terrain height as a function of pressure!!</div><div>Maybe someone has run into this problem with WRF output data and has figure this out?  Maybe I would be better off using a different program such as RIP?  Any advice will be greatly appreciated.</div><div><br></div><div>Thanks,</div><div>MJ</div></div><div> </div><div><br></div></div>