On 10/10/07, <b class="gmail_sendername">Ioannis A. Vamvakas</b> &lt;<a href="mailto:giavvns@ath.hcmr.gr">giavvns@ath.hcmr.gr</a>&gt; wrote:<div><span class="gmail_quote"></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br><br>I&#39;ve been using lats4d to convert a netCDF file into grib1.<br>Everything works just fine but the Center name/id.</blockquote><div><br>Do you mean the originating or the processing center?&nbsp;</div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
 It won&#39;t<br>see the name I provide in the lookup table and will use the<br>default name (=pcmdi) and code (=100).<br><br>The command I use to do this is:<br><br>lats4d -i &lt;file&gt;.nc -o &lt;file&gt; -format grib -table mytable
<br><br>where, &#39;mytable&#39;, is a copy of &#39;.grads.lats.table&#39; with a<br>line at the center section added as follows:<br><br>hcmr | 96 | 0 | 0</blockquote><div><br>The description of the center record in the lats table reads:
<br><br># The format of each record is:<br>#&nbsp;&nbsp; center | GRIB_id | GRIB_center | GRIB_subcenter<br>#<br># center = mnemonic for the center<br># GRIB_id = GRIB generating process id (PDS octet 6)<br># GRIB_center = the id of center managing the data (for AMIP II this is PCMDI) -
<br>&nbsp;see GRIB Table 0<br># GRIB_subcenter = the id of the subcenter<br><br><br>&quot;0&quot; is not a valid grib center; I believe this is being changed into 100.&nbsp; Try something like this:<br><br>hcmr | 96 | 96 | 2<br><br>
Here is what &quot;wgrib -V&quot; reports:<br><br>rec 190:1133616:date 1987010500 DIST kpds5=8 kpds6=100 kpds7=100 levels=(0,100) grid=255 100 mb anl: bitmap: 3312 undef<br>&nbsp; DIST=Geometric height [m]<br>&nbsp; timerange 0 P1 0 P2 0 TimeU 2&nbsp; nx 72 ny 46 GDS grid 0 num_in_ave 0 missing 0
<br>&nbsp; center 96 subcenter 2 process 96 Table 128 scan: WE:NS winds(N/S) <br>&nbsp; latlon: lat&nbsp; -90.000000 to 90.000000 by 4.000000&nbsp; nxny 3312<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; long 0.000000 to 355.000000 by 5.000000, (72 x 46) scan 2 mode 128 bdsgrid 1
<br>&nbsp; min/max data 3.40282e+38 -3.40282e+38&nbsp; num bits 16&nbsp; BDS_Ref 3.40282e+38&nbsp; DecScale 0 BinScale -32766<br><br>&nbsp; Good Luck,<br><br></div></div>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu
</a>