<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
<font style="font-size: 12pt;" face="Arial" size="3">Hello,<br>I'm having some problem with cropping and subsetting GRIB files.<br>I found lats4d script very useful for that, but there is some annoying error it makes - in the resulting GRIB file it writes the wrong variable names (while in the resulting CTL file they are correct).<br><br>For example, I'm taking the </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p GRIB file (attached, together with its CTL and IDX),<br>and it has 5 fields inside: fi, omega, relhum, t, u, v.<br>You can see it in the attached</font><font style="font-size: 12pt;" face="Arial" size="3"> lfff00120000p.ctl</font><font style="font-size: 12pt;" face="Arial" size="3"> file, as well as by running "wgrib </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p"</font><font style="font-size: 12pt;" face="Arial" size="3">.<br>Then I'm trying to extract only 2 fields, e.g. fi and t, to reduce the domain size, and to dump the result into a new GRIB file.<br>To do that, I'm using the following command:<br></font><font style="font-size: 12pt;" face="Arial" size="3"><br></font><div>
<meta http-equiv="Content-Type" content="text/html; charset=unicode">
<meta name="Generator" content="Microsoft SafeHTML"><font style="font-size: 12pt;" face="Arial" size="3">
lats4d.sh -i lfff00120000p.ctl -o lfff00120000p_crop.grb -format grads_grib -lat 29 33.5 -lon 34 36 -vars fi t</font><font style="font-size: 12pt;" face="Arial" size="3"><br><br>This command generates 3 new files: </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p_crop.grb, </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p_crop.ctl and </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p_crop.gmp (attached as well)</font><font style="" face="Arial">.</font><font style="" face="Arial"><br></font><font style="font-size: 12pt;" face="Arial" size="3">The resulting GRIB file seems to be valid, and the CTL </font><font style="font-size: 12pt;" face="Arial" size="3">file looks absolutely fine!<br>But when I try to run "wgrib </font><font style="font-size: 12pt;" face="Arial" size="3">lfff00120000p_crop.grb", instead of fi & t, I see </font><font style="font-size: 12pt;" face="Arial" size="3">PRES & PRMSL</font><font style="font-size: 12pt;" face="Arial" size="3"> variables!<br>My guess is, that the lats4d script is in some way unaware of the fact that the original GRIB file is encoded according to Offenbach (center code 78) GRIB table, but I don't know how to define it.<br><br>Any help would be highly appreciated!<br>Thanks!<br>Alon.<br></font>
<font style="font-size: 12pt;" face="Arial" size="3"><br></font><font style="font-size: 12pt;" face="Arial" size="3">P. S. some aditional technical info, don't know if it's relevant:<br></font><font style="font-size: 12pt;" face="Arial" size="3">
lats4d Version 2.1.5</font><font style="font-size: 12pt;" face="Arial" size="3"><br></font><font style="font-size: 12pt;" face="Arial" size="3">
Grads v2.0.a9.oga.1</font><font style="font-size: 12pt;" face="Arial" size="3"><br></font><font style="font-size: 12pt;" face="Arial" size="3">
Ubuntu 10.04 x86_64</font><font style="font-size: 12pt;" face="Arial" size="3"><br></font></div> </div></body>
</html>