The script in the previous reply was very useful, thank you Javier.<div><br></div><div>I am still having an issue of blank data around the mountains in my cross sections.  I tried searching the archives for a previous post about this but did not find anything.</div><div><br></div><div>Using Javier's script I have been able to get a nice plot of the terrain and my chosen meteorological variables in a cross section but the problem is the data cuts off and is blank near to the mountains.  This is obviously caused by the fact that the model output is interpolated into a finite number of pressure levels.  I have tried dramatically increasing the number of pressure levels the data is written onto and this helps to make the blank data  areas smaller but they are still quite visible.</div><div><br></div><div>Does anyone have any experience with this kind of problem?  I know it sounds weird but, maybe there is some way to plot the topography in a cross section with sigma levels instead of pressure levels?</div><div><br></div><div> Thanks,</div><div>MJ<br></div><div><br></div><br>----- Original Message -----<br>From: "Javier G. Corripio" &lt;Javier.Corripio@UIBK.AC.AT><br>Date: Thursday, August 28, 2008 12:14 am<br>Subject: Re: topography plot from WRF output<br>To: GRADSUSR@LIST.CINECA.IT<br><br><span><p></p></span><tr><td style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><p><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font> This may help:</p><div><br></div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Best,</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Javier</div><div><br></div><div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*plotsfc.gs  <a href="javascript:main.compose('new','t=Javier.Corripio@uibk.ac.at')">Javier.Corripio@uibk.ac.at</a></div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define earthg=9.80665'     ;*acceleration due to gravity (m s-2)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define md=28.966'           ;*molecular weight of dry air</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define atmosr=8.31432'     ;*universal gas constant (j deg-1 kmol-1)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define p0=prmslmsl(lev=1013)'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define T0=tmpprs(lev=1013)'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define lrate=(tmpprs-tmpprs(lev=1013))/hgtsfc(lev=1013)'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define base=(T0/(T0+lrate*hgtsfc))'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define exponent=((earthg*md)/(atmosr*lrate*1000))'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*'define pz=p0*POW(base,exponent)'</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*get gx info</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'q gxinfo'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>data=sublin(result,3)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>px0 = subwrd(data,4)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>px1 = subwrd(data,6)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>data=sublin(result,4)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>py0 = subwrd(data,4)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>py1 = subwrd(data,6)</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>say px0 ' ' px1 ' ' py0 ' ' py1</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*get dim info</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'q dim'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>data=sublin(result,4)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>lev0 = subwrd(data,6)</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>lev1 = subwrd(data,8)</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>say lev0 ' ' lev1</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set lev ' lev0</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*set plotting area</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set parea ' px0 ' 'px1 ' ' py0 ' ' py1</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>*set Yaxis range</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set vrange ' lev0 ' ' lev1</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set gxout linefill'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set lfcols 15 15'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set cmark 0'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'d pressfc/100;const(pressfc/100,' lev0 ')'</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set gxout line'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set parea off'</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>'set lev ' lev0 ' ' lev1</div><div><br></div></div><div><br></div><div><br><div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>On 28 Aug 2008, at 08:18, Michael J Murphy wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Grads users,</div><div><br></div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>set lat 21.5</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>set lev 1000 700</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>d psfc</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>set lat 21.5</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>set lev 1000</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>d psfc</div><div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Thanks,</div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>-- </div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>***********************************************</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Javier G. Corripio</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Faculty of Geo- and Atmospheric Sciences</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>University of Innsbruck</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Innrain, 52</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>A-6020 Innsbruck</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Austria</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>Tel.: +43 512 507 5418</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>email: <a href="javascript:main.compose('new','t=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; "><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>url:<span class="Apple-converted-space"> </span><a href="http://www.uibk.ac.at/geographie/personal/corripio/" target="1">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/" target="1"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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" target="1"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px; ">> </font>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></td></tr>