<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
This may help:<div><br></div><div><br></div><div>The script is for my own needs, so it does not check for any exceptions. If you prefer to use hgtsfc converted to pressure instead of pressfc use the commented lines at the beginning. The script displays the model surface height as a grey mask across a transect with pressure values in the Y-axis.</div><div><br></div><div>Best,</div><div><br></div><div>Javier</div><div><br></div><div><div>*plotsfc.gs <a href="mailto:Javier.Corripio@uibk.ac.at">Javier.Corripio@uibk.ac.at</a></div><div><br></div><div>*'define earthg=9.80665' ;*acceleration due to gravity (m s-2)</div><div>*'define md=28.966' ;*molecular weight of dry air</div><div>*'define atmosr=8.31432' ;*universal gas constant (j deg-1 kmol-1)</div><div>*'define p0=prmslmsl(lev=1013)'</div><div>*'define T0=tmpprs(lev=1013)'</div><div>*'define lrate=(tmpprs-tmpprs(lev=1013))/hgtsfc(lev=1013)'</div><div>*'define base=(T0/(T0+lrate*hgtsfc))'</div><div>*'define exponent=((earthg*md)/(atmosr*lrate*1000))'</div><div>*'define pz=p0*POW(base,exponent)'</div><div><br></div><div>*get gx info</div><div>'q gxinfo'</div><div>data=sublin(result,3)</div><div>px0 = subwrd(data,4)</div><div>px1 = subwrd(data,6)</div><div>data=sublin(result,4)</div><div>py0 = subwrd(data,4)</div><div>py1 = subwrd(data,6)</div><div><br></div><div>say px0 ' ' px1 ' ' py0 ' ' py1</div><div><br></div><div>*get dim info</div><div>'q dim'</div><div>data=sublin(result,4)</div><div>lev0 = subwrd(data,6)</div><div>lev1 = subwrd(data,8)</div><div><br></div><div>say lev0 ' ' lev1</div><div><br></div><div>'set lev ' lev0</div><div><br></div><div>*set plotting area</div><div>'set parea ' px0 ' 'px1 ' ' py0 ' ' py1</div><div><br></div><div>*set Yaxis range</div><div>'set vrange ' lev0 ' ' lev1</div><div><br></div><div>'set gxout linefill'</div><div>'set lfcols 15 15'</div><div>'set cmark 0'</div><div>'d pressfc/100;const(pressfc/100,' lev0 ')'</div><div><br></div><div>'set gxout line'</div><div>'set parea off'</div><div>'set lev ' lev0 ' ' lev1</div><div><br></div></div><div><br></div><div><br><div><div>On 28 Aug 2008, at 08:18, Michael J Murphy wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><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></blockquote></div><br><div> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">-- </div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">***********************************************</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Javier G. Corripio</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Faculty of Geo- and Atmospheric Sciences</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">University of Innsbruck</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Innrain, 52</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">A-6020 Innsbruck</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Austria</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Tel.: +43 512 507 5418</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Fax: <span class="Apple-converted-space"> </span>+43 512 507 2895</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">email: <a href="mailto:Javier.Corripio@uibk.ac.at">Javier.Corripio@uibk.ac.at</a></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">url:<span class="Apple-converted-space"> </span><a href="http://www.uibk.ac.at/geographie/personal/corripio/">http://www.uibk.ac.at/geographie/personal/corripio/</a></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="http://web.mac.com/jgcorripio/">http://web.mac.com/jgcorripio/</a></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="http://www.meteoexploration.com">http://www.meteoexploration.com</a></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><br class="khtml-block-placeholder"></div></span></span><br class="Apple-interchange-newline"></span></span></span></span></span></span></span></span></span><br class="Apple-interchange-newline"></span> </div><br></div></body></html>