<div dir="ltr"><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL-experts,</div><div><br></div><div>I am calculating potential temperature using the following ncl commands:</div><div> theta = pot_temp(p,temp,1,False) ; [time | 365] x [level | 17] x [lat | 73] x [lon | 144]<br></div><div><br></div><div>However, my inputs <b>p</b> and <b>temp </b>have levels that are in <b>hPa/mbar</b>.</div><div>According to the documentation, this should be in <b>Pa (i.e., levels*100)</b>.</div><div><br></div><div>I would like to ask how do I force the variables to have levels in Pa in the code below?</div><div><br></div><div>  <i>temp_f = addfile("<a href="http://air.day.ltm.1981-2010.nc">air.day.ltm.1981-2010.nc</a>", "r"<br>   p  = temp_f->level<br>   temp = temp_f->air(:,{p},:,:)<br></i></div><div><i><br></i></div><div><i>;;Print contents of file<br>   printVarSummary(p)<br>   printVarSummary(temp)<br><br>;;**************************************************<br>;;Calculates Potential Temperature<br>;;**************************************************<br>  theta = pot_temp(p,temp,1,False)<br>  print(theta&level)</i></div><div><i><br></i></div><div>I'll appreciate any help on this.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div><div><br></div></div></div></div></div></div></div></div></div></div></div>