<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div>
<style type="text/css">
<!--
@page { margin: 2cm }
P { margin-bottom: 0.21cm }
-->
</style>
<div style="background: none repeat scroll 0% 0% transparent; font-style: normal;"><a href="" name="yui_3_7_2_17_1362125096836_137"></a>
<font face="Arial" size="2">Dear Jayakrishna, </font></div><div style="background: none repeat scroll 0% 0% transparent; font-style: normal; color: rgb(0, 0, 0); font-size: 13px; font-family: Arial;"><font color="#000000"><font face="Arial"><font size="2">I've tried the script but
the problem persists: </font></font></font>
</div>
<div style="background: transparent; font-style: normal"><a href="" name="yui_3_7_2_17_1362125096836_1053"></a><a href="" name="yui_3_7_2_17_1362125096836_130"></a><a href="" name="yui_3_7_2_17_1362125096836_142"></a>
<font color="#000000"><font face="Arial"><font size="2">ga-> define
psi= fish_psi(u,v.2)<br>Error from fish: rc = 1000<br>Operation
Error: Error from fish function<br>Operation Error: Error
from fish_psi function<br> Error ocurred at column 1<br>DEFINE
error: Invalid expression. </font></font></font></div>
<div style="background: transparent"><a href="" name="yui_3_7_2_17_1362125096836_146"></a>
</div><div style="margin-bottom: 0cm;">If push comes to shove I'll try the
idea proposed by Bernd Becker (Thanks for the email !!!)</div><div style="margin-bottom: 0cm;"><br></div><div style="margin-bottom: 0cm; color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">Arlindo <br></div>
</div><div> </div><div><br><br></div> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Jayakrishnan P.R <prjayakrishnan@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> GrADS Users Forum <gradsusr@gradsusr.org> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quinta-feira, 28 de Fevereiro de 2013 19:11<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [gradsusr] How to calculate the velocity potential for a local domain ?<br> </font> </div> <br><div id="yiv210830478">Dear Arlindo,<br> I herewith attach a GS file in which velocity potential and stream function is computed from
u and v components. I attempted this long back and I am sorry to say that I could not remember the whole steps. Still I hope this may be useful. Also the following link may be useful.<br>
<a rel="nofollow" target="_blank" href="http://opengrads.org/doc/udxt/fish/fish.html#computing_velocity_potential_from_wind_components">http://opengrads.org/doc/udxt/fish/fish.html#computing_velocity_potential_from_wind_components</a><br><br>'reinit'<br>
'sdfopen /media/JKPRMUSIC/JKPR/DATA/ncep/d/uvw/<a rel="nofollow" target="_blank" href="http://uwnd.2000.nc/">uwnd.2000.nc</a>'<br>'sdfopen /media/JKPRMUSIC/JKPR/DATA/ncep/d/uvw/<a rel="nofollow" target="_blank" href="http://vwnd.2000.nc/">vwnd.2000.nc</a>'<br>'set lon 0 357.5'<br>
'set lat -90 90'<br>prompt 'key-in file name for psi:'<br>pull fname<br>'set fwrite 'fname<br>tt=1<br>while (tt<=10)<br>'set t 'tt<br>zz=1<br>while (zz<=2)<br>'set z 'zz<br>'define psi= fish_psi(uwnd,vwnd.2)'<br>
'define one=1+0*lat'<br>'define upsi=mvadv(one,psi)'<br>'define vpsi=-muadv(one,psi)'<br>*'set gxout shaded'<br>*'display psi/1e7'<br>*'set gxout stream'<br>*'display upsi;vpsi'<br>
'set gxout fwrite'<br>'d psi/1e7'<br>'d upsi'<br>'d vpsi'<br>zz=zz+1<br>endwhile<br>tt=tt+1<br>endwhile<br>'disable fwrite'<br><br>If any suggestions I will again look into my old archives. Best of Luck.<br clear="all">
<br>-- <br>Sincerely<br>***********************************************<br>Jayakrishnan.P.R<br>Teaching Associate (Atmospheric Science)<br>Academy of Climate Change Education and Research<br>Kerala Agricultural University<br>
Vellanikkara-680 656<br><br>and also<br><br>Research Fellow (Part Time)<br>Department of Atmospheric Sciences<br>CUSAT, Cochin-682 016<br>Kerala, India.<br>Mob: 09496026547<br><br>***********************************************
</div><br>_______________________________________________<br>gradsusr mailing list<br><a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br><br><br> </div> </div> </div></body></html>