<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=utf-8" http-equiv=Content-Type>
<STYLE type=text/css>BODY {
        PADDING-BOTTOM: 0px; MARGIN: 0px; PADDING-LEFT: 0px; PADDING-RIGHT: 0px; PADDING-TOP: 0px
}
BLOCKQUOTE {
        MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px; MARGIN-LEFT: 2em
}
OL {
        MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px
}
UL {
        MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px
}
</STYLE>

<META name=GENERATOR content="MSHTML 8.00.6001.19222"></HEAD>
<BODY style="MARGIN-TOP: 0px; MARGIN-LEFT: 0px; MARGIN-RIGHT: 0px"><!-- 52_302903_e11415f126781bf -->
<TABLE style="FONT-FAMILY: Verdana, Arial, Helvetica, sans-serif" id=Table_01 
border=0 cellSpacing=0 cellPadding=0 width="99.99%">
  <TBODY>
  <TR>
    <TD vAlign=top background=cid:__2@Foxmail.net><IMG 
      src="cid:__0@Foxmail.net" width=138 height=1></TD>
    <TD vAlign=top background=cid:__3@Foxmail.net>&nbsp;</TD>
    <TD vAlign=top background=cid:__4@Foxmail.net><IMG 
      src="cid:__1@Foxmail.net" width=80 height=60></TD></TR>
  <TR>
    <TD vAlign=top background=cid:__5@Foxmail.net>&nbsp; </TD>
    <TD 
    style="PADDING-BOTTOM: 70px; PADDING-LEFT: 30px; PADDING-RIGHT: 0px; HEIGHT: 300px; PADDING-TOP: 0px" 
    vAlign=top background=cid:__3@Foxmail.net width="100%">
      <TABLE 
      style="MARGIN-TOP: 0px; MARGIN-LEFT: 0px; cellpadding: 0px; cellspacing: 0px" 
      width="100%">
        <TBODY>
        <TR>
          <TD style="LINE-HEIGHT: 160%" vAlign=top>
            <DIV 
            style="WIDTH: 100%; WHITE-SPACE: normal; WORD-BREAK: break-all"><FONT 
            color=#265800 size=2 face=Verdana>
            <DIV><FONT color=#000000>I think streamfunciton is to see cyclonic 
            and anticyclonic circulations in horizontal plane, not in the 
            vertical. For&nbsp;vertical circulation&nbsp;purposes you may use 
            horizontal and vertical winds or mass-streamfunction. I don't know 
            how to calcuate the later in grads.</FONT></DIV>
            <DIV><FONT color=#000000>Good Luck!</FONT></DIV>
            <DIV><FONT color=#c0c0c0>2012-05-22 </FONT></DIV>
            <DIV align=left>
            <HR style="WIDTH: 122px; HEIGHT: 2px" align=left SIZE=2>
            </DIV>
            <DIV><FONT color=#c0c0c0><SPAN>
            <DIV>
            <DIV><FONT size=2 face=Verdana></FONT><FONT size=4 
            face="Franklin Gothic Book">Muhammad Afzaal Karori</FONT></DIV>
            <DIV><FONT face="Franklin Gothic Book">Institute of Atmospheric 
            Physics</FONT></DIV>
            <DIV><FONT face="Franklin Gothic Book">Chinese Academy of 
            Sciences</FONT></DIV>
            <DIV><FONT face="Franklin Gothic Book">Beijing - P.R. 
            China</FONT></DIV>
            <DIV>&nbsp;</DIV></DIV></SPAN></FONT></DIV>
            <HR>

            <DIV><STRONG>发件人:</STRONG> resmi ea </DIV>
            <DIV><STRONG>发送时间:</STRONG> 2012-05-22&nbsp; 18:18:13 </DIV>
            <DIV><STRONG>收件人:</STRONG> GrADS Users Forum </DIV>
            <DIV><STRONG>抄送:</STRONG> </DIV>
            <DIV><STRONG>主题:</STRONG> [gradsusr] stream function in pressure 
            level </DIV></FONT>
            <DIV>&nbsp;</DIV>
            <DIV>Dear users<BR>&nbsp;i have calculated stream function using 
            open grads from 1000-10 hPa level using NCEP/NCAR reanalysis data 
            ex.1dec2009<BR>The grads script give the grd file containing three 
            variable like psi,upsi and vpsi.<BR>Here i have attached the figure. 
            i find mistake in the figure...the programme doesnot undergone 
            integration part, at each level it calculate three variable and 
            finally dispalys.<BR>&nbsp;i want to know that whether it is correct 
            or not. how can i modify the programme<BR>'reinit'<BR>'sdfopen 
            G:\windows\ncep\uwnd\<A 
            href="http://uwnd.2009.nc">uwnd.2009.nc</A>'<BR>'sdfopen 
            G:\windows\ncep\uwnd\<A 
            href="http://vwnd.2009.nc">vwnd.2009.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&lt;=1)<BR>'set t 'tt<BR>zz=1<BR>while (zz&lt;=17)<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 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><BR clear=all><BR>-- <BR>Resmi E A<BR>Research 
            Scholar<BR>Cochin university of science and 
            technology<BR>Cochin-16<BR><BR></DIV></DIV></TD></TR></TBODY></TABLE></TD>
    <TD vAlign=top background=cid:__4@Foxmail.net 
width=47>&nbsp;</TD></TR></TBODY></TABLE></BODY></HTML>