Hi users,<br><br>I need help.<br><br>I did this script below to calculate potential temperature but I don&#39;t know how to call inside my script.<br><br>tempc = variable from my file and 700 level I want theta.<br><br>Thanks,<br>


<br>Guilherme.<br><br>&#39;open se_g1.ctl&#39;<br>&#39;set gxout shaded&#39;<br><br>gtheta(tempc,700)<br><br>function gtheta(temp,niv)<br>theta1 = (temp+273.15)*pow(1000/niv,0.286)<br>return(theta1)<br><br clear="all"><br>

-- <br>
-----------------------------------------------------------<br>José Guilherme Martins dos Santos<br>Grupo de Modelagem da Atmosfera e Interfaces - GMAI<br>Divisão de Modelagem e Desenvolvimento - DMD<br>Centro de Previsão de Tempo e Estudos Climáticos - INPE<br>


Rodovia Presidente Dutra, km 39 - CEP 12630-000<br>Cachoeira Paulista, SP, Brasil<br>Telefone(CP): +55 12 3186 8502<br>Celular(SJC): +55 12 91678770 / +55 12 81528580<br>E-mail: <a href="mailto:guilherme.martins@cptec.inpe.br" target="_blank">guilherme.martins@cptec.inpe.br</a><br>


E-mail alternativo: <a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a><br>Homepage: <a href="http://guilhermemartins.hd1.com.br/" target="_blank">http://guilhermemartins.hd1.com.br/</a><br>-----------------------------------------------------------<br>