Fwd: building psichig.exe using pgf90 or g77 FORTRAN compiler

[Arlindo da Silva]

All,

  See examples below for computing chi/psi with fish() on Windows.

          Arlindo

---------- Forwarded message ----------
From: See Hai Ooi <axl419 at yahoo.com>
Date: Feb 3, 2008 5:16 PM
Subject: Re: building psichig.exe using pgf90 or g77 FORTRAN compiler
To: Arlindo da Silva <dasilva at alum.mit.edu>


Dear Mr Arlindo,

    It works.  Please see the attached ps.gs
(modified short version of your utFish.gs) and the
last output of velocity potential in VP.doc.

    Thank you very much indeed.

    Best regards.

Ooi



--- Arlindo da Silva <dasilva at alum.mit.edu> wrote:

> Hi,
>
> Look under
>
> http://opengrads.org/pre-rel/win32/
>
> file grads-1.9.0-rc1.win32_gex.zip. Unzip, and
> replace the win32/gex folder
> in your original PCGrADS with then one here. Also,
> copy utFish.gs included
> under lib/ here to your lib folder.
>
> Start grads and type utFish. It shoudl work. Also,
> read the script
> carefully.
>
> Let me know if this works.
>
> arlindo

Send instant messages to your online friends http://uk.messenger.yahoo.com

*
 'reinit'
*
 'sdfopen c:\cdc2\uvw\uwnd.1997.nc'
 'sdfopen c:\cdc2\uvw\vwnd.1997.nc'
*
 'set lev 200'
*
  vars = '(uwnd,vwnd.2)'
*
  'set gxout shaded'
*
* Vorticity and divervenge: enlarge x-domain to avoid grid undefs
* ---------------------------------------------------------------
 xyrange()
 xmin = _xmin - 1;  xmax = _xmax + 1
 'set x ' xmin ' ' xmax
 'vort = hcurl'vars
 'div  = hdivg'vars

* Reset x domain so that it is not x-wrapped
* ------------------------------------------
 xmin = _xmin;  xmax = _xmax
 'set x ' xmin ' ' xmax

* Display vorticity/divergence
* ----------------------------
 'd vort'; 'draw title Vorticity'
 clear()
 'd div';  'draw title Divergence'
 clear()
*
* Stream function
* ---------------
 'psi = fish(vort)'
 'd psi'
 'draw title Stream Function'
 clear()
*
* Velocity Potential
* ------------------
 'chi = fish(div)'
 'set gxout shaded'
 'd chi'
 'draw title Velocity Potential'
*
function xyrange()

     'q file'
     tmp = sublin ( result, 5 )
     _xmin = 1
     _xmax = subwrd(tmp,3)
     _ymin = 1
     _ymax = subwrd(tmp,6)

return

function clear()

 say 'Hit <CR> to continue...'; pull ans; 'clear'

return




--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080203/acb88198/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: VP.doc
Type: application/msword
Size: 560640 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080203/acb88198/attachment.doc