<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Hello,<br><br>I'm having some trouble using the fish_psi function to plot streamfunction for the GFS data.<br><br>'open gfs.ctl'<br>'d fish_psi(ugrdprs,vgrdprs)'<br><br>The field is returning constant value 0 or -inf.<br><br>Interestingly enough, I use to be able to plot streamfunction using fish_psi with no problem for the GFS. So this 'bug' has occurred over the last few months or so.<br><br>Anyone have any ideas what may be causing the problem with the GFS data? I have use the function fish_psi just fine using the ECMWF-interim data. <br><br>I'm using grads-2.1.a2.oga.1<br><br>Thanks,<br>James Spinks<br>                                               </div></body>
</html>