<div dir="ltr"><div>What does 'press1' look like?</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml"><b>printVarSummary</b></a>(press1)</div><div>-------------</div><div>Note: you can have the following WRF function calculate 'dbz' directly.</div><div><br></div><div> dbz = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_dbz.shtml"><b>wrf_user_getvar</b></a>(file1,"dbz",-1) <br> printVarSummary(dbz)<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 24, 2020 at 9:18 AM David W Reynolds via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;" lang="EN-US"><div class="gmail-m_3392284218114162836WordSection1"><p class="MsoNormal">When attempting to run my ncl script to execute wrf_dbz using this function I get the following error. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">fatal:wrf_dbz: The prs array must have at least 3 dimensions<u></u><u></u></p><p class="MsoNormal">fatal:["Execute.c":8575]:Execute: Error occurred at or near line 266 in file MAIN_plot_Reflectivity_Txt.ncl<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Line 266 is the line trying to run the function. Can anyone help with this error. I am very new to ncl so apologize for my ignorance. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">263 ; Calculate Reflectivity<u></u><u></u></p><p class="MsoNormal"> 264 ivarint = 0 ; or 1<u></u><u></u></p><p class="MsoNormal"> 265 iliqskin = 1 ; or 1<u></u><u></u></p><p class="MsoNormal"> 266 dbz = wrf_dbz( press1, tk1, qvapor1, qrain, qsnow, qgraup, ivarint, iliqskin)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Here is how I read the data. Trying to output a vertical profile of reflectivity over a given lat long. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">; Get fields from file<u></u><u></u></p><p class="MsoNormal"> pbtemp1 = file1->$PBTEMPNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> qvapor1 = file1->$QVAPORNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> bapress1 = file1->$BAPRESNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> pbpress1 = file1->$PBPRESNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> bageopot1 = file1->$BAGEOPNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> pbgeopot1 = file1->$PBGEOPNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> uwrf1 = file1->$UCOMPNAME1$(TIME1,:,:,:)<u></u><u></u></p><p class="MsoNormal"> vwrf1 = file1->$VCOMPNAME1$(TIME1,:,:,:)<u></u><u></u></p><p class="MsoNormal"> qrain = file1->$QRAINNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> qsnow = file1->$QSNOWNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"> qgraup = file1->$QGRAUPELNAME1$(TIME1,:,i1,j1)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Thanks and Happy Holidays!! <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Dave<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank">Mail</a> for Windows 10</p><p class="MsoNormal"><u></u> <u></u></p></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>