All,<br><br> See examples below for computing chi/psi with fish() on Windows. <br><br> Arlindo<br><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">See Hai Ooi</b> <<a href="mailto:axl419@yahoo.com">axl419@yahoo.com</a>><br>
Date: Feb 3, 2008 5:16 PM<br>Subject: Re: building psichig.exe using pgf90 or g77 FORTRAN compiler<br>To: Arlindo da Silva <<a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>><br><br><br>Dear Mr Arlindo,<br>
<br> It works. Please see the attached <a href="http://ps.gs" target="_blank">ps.gs</a><br>(modified short version of your utFish.gs) and the<br>last output of velocity potential in VP.doc.<br><br> Thank you very much indeed.<br>
<div class="Ih2E3d"><br> Best regards.<br><br>Ooi<br><br><br><br>--- Arlindo da Silva <<a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>> wrote:<br><br></div><div class="Ih2E3d">> Hi,<br>><br>
> Look under<br>><br>> <a href="http://opengrads.org/pre-rel/win32/" target="_blank">http://opengrads.org/pre-rel/win32/</a><br>><br>> file grads-1.9.0-rc1.win32_gex.zip. Unzip, and<br>> replace the win32/gex folder<br>
> in your original PCGrADS with then one here. Also,<br>> copy utFish.gs included<br>> under lib/ here to your lib folder.<br>><br>> Start grads and type utFish. It shoudl work. Also,<br>> read the script<br>
> carefully.<br>><br>> Let me know if this works.<br>><br>> arlindo<br><br></div><div><div></div><div class="Wj3C7c">Send instant messages to your online friends <a href="http://uk.messenger.yahoo.com" target="_blank">http://uk.messenger.yahoo.com</a> </div>
</div><br>*<br> 'reinit'<br>*<br> 'sdfopen c:\cdc2\uvw\uwnd.1997.nc'<br> 'sdfopen c:\cdc2\uvw\vwnd.1997.nc'<br>*<br> 'set lev 200'<br>*<br> vars = '(uwnd,vwnd.2)'<br>*<br> 'set gxout shaded'<br>
*<br>* Vorticity and divervenge: enlarge x-domain to avoid grid undefs<br>* ---------------------------------------------------------------<br> xyrange()<br> xmin = _xmin - 1; xmax = _xmax + 1<br> 'set x ' xmin ' ' xmax<br>
'vort = hcurl'vars<br> 'div = hdivg'vars<br><br>* Reset x domain so that it is not x-wrapped<br>* ------------------------------------------<br> xmin = _xmin; xmax = _xmax<br> 'set x ' xmin ' ' xmax<br>
<br>* Display vorticity/divergence<br>* ----------------------------<br> 'd vort'; 'draw title Vorticity'<br> clear()<br> 'd div'; 'draw title Divergence'<br> clear()<br>*<br>* Stream function<br>
* ---------------<br> 'psi = fish(vort)'<br> 'd psi'<br> 'draw title Stream Function'<br> clear()<br>*<br>* Velocity Potential<br>* ------------------<br> 'chi = fish(div)'<br> 'set gxout shaded'<br>
'd chi'<br> 'draw title Velocity Potential'<br>*<br>function xyrange()<br><br> 'q file'<br> tmp = sublin ( result, 5 )<br> _xmin = 1<br> _xmax = subwrd(tmp,3)<br> _ymin = 1<br>
_ymax = subwrd(tmp,6)<br><br>return<br><br>function clear()<br><br> say 'Hit <CR> to continue...'; pull ans; 'clear'<br><br>return<br><br></div><br><br clear="all"><br>-- <br>Arlindo da Silva<br>
<a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>