<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'>Thanks, Bob.<br><br>I read right over that statement from the link you sent and didn't even notice sh_fish. Appreciate your response.<br><br><div><hr id="stopSpelling">From: rhart@fsu.edu<br>To: gradsusr@gradsusr.org<br>Date: Thu, 24 Mar 2016 02:43:47 +0000<br>Subject: Re: [gradsusr] Streamfunction (fish_psi)<br><br>
<style style="display:none;">
.ExternalClass P {
}
</style>
<div id="ecxdivtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
Hi James,<BR>
<br>
<BR>
When the GFS went to 0.25 degree output last year, the fishpak routines became unreliable for the GFS:<BR>
<br>
<BR>
<a title="Ctrl+Click or tap to follow the link" id="ecxLPlnk49616" href="http://opengrads.org/doc/udxt/fish/fish.html" target="_blank">http://opengrads.org/doc/udxt/fish/fish.html</a><BR>
<br>
<BR>
"Notice that <code>fish()</code> becomes numerically unstable for horizontal resolutions finer that 1/2 degrees or so. In such cases use the spherical harmonic based
<strong>sh_fish</strong> given in extension <em>shfilt</em>."<BR>
<br>
<BR>
Indeed, I have found sh_fish is a very reliable replacement at these resolutions.<BR>
<br>
<BR>
Bob Hart<br>
<BR>
<br>
<br>
<div style="color:rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block;width:98%;">
<div id="ecxdivRplyFwdMsg" dir="ltr"><font style="font-size:11pt;" color="#000000" face="Calibri, sans-serif"><b>From:</b> gradsusr-bounces@gradsusr.org <gradsusr-bounces@gradsusr.org> on behalf of James Spinks <jspin881@hotmail.com><br>
<b>Sent:</b> Wednesday, March 23, 2016 4:22 PM<br>
<b>To:</b> gradsusr@gradsusr.org<br>
<b>Subject:</b> [gradsusr] Streamfunction (fish_psi)</font>
<div> </div>
</div>
<div>
<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>
</div>
</div>
</div>
<br>_______________________________________________
gradsusr mailing list
gradsusr@gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr</div> </div></body>
</html>