[Met_help] [rt.rap.ucar.edu #92553] History for MODE error

John Halley Gotway via RT met_help at ucar.edu
Fri Oct 11 09:47:53 MDT 2019


----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

I'm trying to run MODE on some hail data, and am getting an error I haven't seen before:

ERROR  :
ERROR  :
ERROR  :   ShapeData::shrink(const int) -> value too large ... new grid is empty
ERROR  :
HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:
  #000: H5T.c line 1763 in H5Tclose(): not a datatype
    major: Invalid arguments to routine
    minor: Inappropriate type
terminate called after throwing an instance of 'netCDF::exceptions::NcHdfErr'
  what():  NcHdfErr: An error was reported by the HDF5 layer.
file: ncFile.cpp  line:15
Abort

Any idea what this means?  It seems like it doesn't like something in the netCDF file, though I'm not sure what.  I have some files in /scratch/WEEKLY/Tina

Thanks,
Tina

----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: MODE error
From: John Halley Gotway
Time: Thu Oct 10 11:41:25 2019

Tina,

I moved this over to met-help.  Whenever debugging MET's ability to
read gridded data, I always test using plot_data_plane:

/usr/local/met-8.1.1/bin/plot_data_plane \
news-e_MET_12_20190518_0000_0100.nc \
news-e_MET_12_20190518_0000_0100.ps \
'name="HAILCAST"; level="(0,*,*)";'

That also fails, but with a different error message:

ERROR  : Box::set_lrbt() -> width and/or height negative!

Based on that, I suspect there's a problem with the definition of the
grid:

                :Projection = "Lambert Conformal" ;
                :scale_lat_1 = "60.00000" ;
                :scale_lat_2 = "30.00000" ;
                :lat_pin = 35.07053f ;
                :lon_pin = -105.6247f ;
                :x_pin = "0.000000" ;
                :y_pin = "0.000000" ;
                :lon_orient = -100.5997f ;
                :d_km = "3.000000" ;
                :r_km = "6371.200000" ;
                :nx = 300LL ;
                :ny = 300LL ;

Those lines for nx and ny = 300LL look a little suspicious to me.
They're usually listed as strings, like this:
                :nx = "185" ;
                :ny = "129 grid_points" ;
So I tried changing them to strings:
ncatted -a nx,global,o,c,"300" \
        -a ny,global,o,c,"300" \
        news-e_MET_12_20190518_0000_0100.nc \
        -o news-e_MET_12_20190518_0000_0100_JHG.nc

And now plot_data_plane is happy again:
/usr/local/met-8.1.1/bin/plot_data_plane \
news-e_MET_12_20190518_0000_0100_JHG.nc \
news-e_MET_12_20190518_0000_0100.ps \
'name="HAILCAST"; level="(0,*,*)";'

And produces the attached image.  Presumably, that's the underlying
problem in MODE as well.

Thanks,
John





------------------------------------------------
Subject: MODE error
From: Minna Win
Time: Thu Oct 10 18:39:14 2019

Hi Tina,

I'm going to assign this ticket to Randy, who handles our MODE and MDT
issues.

Thanks,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax:   303-497-8401



On Thu, Oct 10, 2019 at 11:32 AM John Halley Gotway via RT <
met_help at ucar.edu> wrote:

>
> Thu Oct 10 11:32:13 2019: Request 92553 was acted upon.
> Transaction: Ticket created by johnhg
>        Queue: met_help
>      Subject: MODE error
>        Owner: johnhg
>   Requestors: kalb at ucar.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=92553 >
>
>
> I'm trying to run MODE on some hail data, and am getting an error I
> haven't seen before:
>
> ERROR  :
> ERROR  :
> ERROR  :   ShapeData::shrink(const int) -> value too large ... new
grid is
> empty
> ERROR  :
> HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:
>   #000: H5T.c line 1763 in H5Tclose(): not a datatype
>     major: Invalid arguments to routine
>     minor: Inappropriate type
> terminate called after throwing an instance of
> 'netCDF::exceptions::NcHdfErr'
>   what():  NcHdfErr: An error was reported by the HDF5 layer.
> file: ncFile.cpp  line:15
> Abort
>
> Any idea what this means?  It seems like it doesn't like something
in the
> netCDF file, though I'm not sure what.  I have some files in
> /scratch/WEEKLY/Tina
>
> Thanks,
> Tina
>

------------------------------------------------
Subject: MODE error
From: Christina Kalb
Time: Thu Oct 10 20:06:18 2019

Thanks!  I made that change, and MODE ran fine now :-D

Tina

On Thu, Oct 10, 2019 at 11:41 AM John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Tina,
>
> I moved this over to met-help.  Whenever debugging MET's ability to
read
> gridded data, I always test using plot_data_plane:
>
> /usr/local/met-8.1.1/bin/plot_data_plane \
> news-e_MET_12_20190518_0000_0100.nc \
> news-e_MET_12_20190518_0000_0100.ps \
> 'name="HAILCAST"; level="(0,*,*)";'
>
> That also fails, but with a different error message:
>
> ERROR  : Box::set_lrbt() -> width and/or height negative!
>
> Based on that, I suspect there's a problem with the definition of
the grid:
>
>                 :Projection = "Lambert Conformal" ;
>                 :scale_lat_1 = "60.00000" ;
>                 :scale_lat_2 = "30.00000" ;
>                 :lat_pin = 35.07053f ;
>                 :lon_pin = -105.6247f ;
>                 :x_pin = "0.000000" ;
>                 :y_pin = "0.000000" ;
>                 :lon_orient = -100.5997f ;
>                 :d_km = "3.000000" ;
>                 :r_km = "6371.200000" ;
>                 :nx = 300LL ;
>                 :ny = 300LL ;
>
> Those lines for nx and ny = 300LL look a little suspicious to me.
They're
> usually listed as strings, like this:
>                 :nx = "185" ;
>                 :ny = "129 grid_points" ;
> So I tried changing them to strings:
> ncatted -a nx,global,o,c,"300" \
>         -a ny,global,o,c,"300" \
>         news-e_MET_12_20190518_0000_0100.nc \
>         -o news-e_MET_12_20190518_0000_0100_JHG.nc
>
> And now plot_data_plane is happy again:
> /usr/local/met-8.1.1/bin/plot_data_plane \
> news-e_MET_12_20190518_0000_0100_JHG.nc \
> news-e_MET_12_20190518_0000_0100.ps \
> 'name="HAILCAST"; level="(0,*,*)";'
>
> And produces the attached image.  Presumably, that's the underlying
> problem in MODE as well.
>
> Thanks,
> John
>
>
>
>
>
>

--
*************************************************
Christina Kalb
Associate Scientist
National Center for Atmospheric Research
3450 Mitchell Lane
Boulder, CO 80301
kalb at ucar.edu
303-497-2819
*************************************************

------------------------------------------------


More information about the Met_help mailing list