[gradsusr] Problem extracting HIRLAM gribfile

Tue Mar 15 11:41:32 EDT 2011

You've got another problem with this grib file. That error message  
"getdrstemplate: DRS..." is from the grib2 library, not GrADS. Data  
Representation Template 5.50002 is not in my GRIB2 doc, nor is it in a  
newer version of the doc from WMO.
--Jennifer

On Mar 15, 2011, at 11:24 AM, Jonas Kaufmann wrote:

> Hello Jennifer,
>
> thank you! This was the output that was generated by g2ctl.pl
> (including the negative sign), so I think the error is on the
> g2ctl-side.
>
> However, when I manually change this to a positive value, I get  
> another error:
> grib2map: scanning GRIB2 file: out.grib2
> getdrstemplate: DRS Template 5.50002 not defined.
> grib2map error: g2_getfld failed: ierr=12
>
> Is this a related problem or am I running into another problem with
> the grib file?
>
> Best regards,
> Jonas Kaufmann
>
> On Tue, Mar 15, 2011 at 15:14, Jennifer Adams <jma at cola.iges.org>  
> wrote:
>>
>> On Mar 15, 2011, at 10:05 AM, Jonas Kaufmann wrote:
>>
>> Hello Wesley,
>>
>> thank you for your help! I could convert the file to grib2 using  
>> ECMWF
>> grib_set, see the converted file as an attachment. I then used the
>> latest versions of g2ctl and wgrib2 to create a control file (also
>> attached). During this process I get the following error):
>> Argument "0.000000:winds(grid)" isn't numeric in multiplication (*)  
>> at
>> /usr/share/perl/5.10/Math/Trig.pm line 87.
>>
>> Now when I try to create an index file for the grib2-file, gribmap
>> complains with this error:
>> Open Error:  Missing or invalid dimension increment value
>>  --> The invalid description file record is:
>>  --> xdef 351 linear 245.437796901893 -0.459621342512909
>>  The data file was not opened.
>>
>> Not sure what the 1st problem mentioned above means, but this  
>> message is
>> telling you that XDEF cannot have a negative increment value. You  
>> must start
>> with the minimum longitude and have positive increment.
>> --Jennifer
>>
>> Is there a simple (manual?) way to fix this in the control file? Or  
>> is
>> there another option on how to display data from this file?
>>
>> Thank you for your help!
>>
>>
>> Best regards,
>> Jonas Kaufmann
>>
>> On Fri, Mar 11, 2011 at 20:08, Wesley Ebisuzaki
>> <Wesley.Ebisuzaki at noaa.gov> wrote:
>>
>> Jonas,
>>
>>   Your grib1 file is packed using "complex packing".  This kind of
>>
>> packing is
>>
>> not supported by grib2ctl and GrADS.
>>
>>   I think that you are getting the error message because the grid is
>>
>> (582 x 424) or 246768 points.
>>
>> However, the number of grid points (could be modified by undefineds  
>> and
>>
>> runs of identical values)
>>
>> is given by octets 19-20 in section 4. Since 246768 much too big to  
>> be
>>
>> stored in a 2 byte number,
>>
>> an error is flagged.  Try converting this file to grib2 using the  
>> ECMWF
>>
>> software.  Once in grib2,
>>
>> g2ctl and GrADS can read it.  I don't think that g2ctl will handle  
>> the
>>
>> projection right but at least you
>>
>> can display it as a x-y grid.
>>
>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/section4.html
>>
>>
>>    Wesley Ebisuzaki
>>
>> Jonas Kaufmann wrote:
>>
>> Hello,
>>
>> I ran into a problem trying to extract a HIRLAM data set on a rotated
>>
>> lat-lon grid. the file is in grib1-format. I could create a valid
>>
>> control file using the code from g2ctl.pl for grib2-files, merging
>>
>> this together with the code in grib2ctl.pl for grib1 files.  
>> Although a
>>
>> control file is created, I also get errors like:
>>
>> grib header at record 1: two values of nxny 246768 144639
>>
>>    LEN 162734 DataStart 14 UnusedBits 6 #Bits 9 nxny 246768
>>
>> When I then go ahead and try to plot a map using GrADS, I always get
>>
>> "Entire Grid undefined". I have attached the grib file, the index  
>> file
>>
>> and the control file.
>>
>> What am I doing wrong? How can I extract the data from this gribfile?
>>
>> Thank you for your help!
>>
>> Best regards,
>>
>> Jonas Kaufmann
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>>
>> gradsusr mailing list
>>
>> gradsusr at gradsusr.org
>>
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>> _______________________________________________
>>
>> gradsusr mailing list
>>
>> gradsusr at gradsusr.org
>>
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>> < 
>> out 
>> .grib2.ctl><out.grib2>_______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>> --
>> Jennifer M. Adams
>> IGES/COLA
>> 4041 Powder Mill Road, Suite 302
>> Calverton, MD 20705
>> jma at cola.iges.org
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110315/5dc5f69f/attachment-0003.html 